作者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/cn.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