作者ronald736 (待業中 好想有工作)
看板Statistics
標題[程式] 請問R程式做多重比較(fisher法 LSD)
時間Sat Jan 3 13:06:27 2015
[軟體程式類別]:R
請填入軟體程式類別,例如:SAS、SPSS、R、EVIEWS...等
[程式問題]:目前在跑多重比較 (平均數差異比較)
資料處理、迴歸、敘述統計、logistic、probit...等
[軟體熟悉度]:低
請把以下不需要的部份刪除
新手(不到1個月)
低(1~3個月)
中(3個月到1年)
高(1年以上)
根據不同熟悉度
板上高手可以給予你不同指導
新手或許需要的還是語法的熟悉
而其他部分則可能往往是盲點上的點破
[問題敘述]: 我載入套件 library(asbio) 後 跑 LSD.test(y,treatment)
回應是: 沒有這個函數 沒有這個函數 "LSD.test"
請問是我套件少了嗎?
R版本是3.1.2版
請詳盡敘述遭遇到的問題,可能的話,分點敘述你要處理的流程
能夠分點敘述處理流程,基本上在撰寫程式或者跑統計時是相當有幫助的
甚至有經驗的人可以藉此告知你步驟是否錯誤
[程式範例]:
x=read.csv("D:/R/實設.csv",header=T)
x
y=x[,2]#y,反應變數值
y
treatment=x[,1]#處理值(因子)
treatment
lm.data=lm(y~treatment,data=x)#線性模型
anova(lm.data)#anova表
coef(lm.data)#各處理的影響程度 μi-μ
aov.data=aov(y~treatment,data=x)#用aov的anova表
aov.data
coef(aov.data)#各處理的影響程度 μi-μ
mean(y)#y平均為41 為μ
model.tables(aov.data,type="means")#為各個μi
library(asbio)
LSD.test(y,treatment)#沒有這函數? <-----在這邊出問題
雖然張貼程式很可怕,但基本上有些程式還是要張貼才能解決
-----------------------------------------------------------------------------
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.192.173.104
※ 文章網址: http://webptt.com/m.aspx?n=bbs/Statistics/M.1420261590.A.5E0.html
2F:→ andrew43: 上面是內建或agricolae package的自訂function 01/04 00:51
3F:推 celestialgod: library( agricolae) 01/07 12:05