2011-01-06 97 views
0

我写一些帮助R.最小二乘方法缺乏数据

我使用下面的脚本来 做一个简单的分析RCBD比较数据集基因型(名称) 为特质的“X”。

library(stats) 
data_1=read.table(file="test.txt", head=TRUE) 
result_X =aov(X~Block+Name, data=data_1) 
sink("result_X.txt") 
summary(result_X) 
sink() 

我的资料缺失(“不适用”)。所以,在计算LSD后,我想比较 基因型的降序。我不认为平均数 是好的,因为某些'名字'的某些 '块'缺失。

所以,问题是,什么是脚本打印出'最小 平均手段',其中我认为是最好的与LSD比较简单的平均数。 感谢你的帮助,

奥斯瓦尔德

+1

首先,请参阅[关于R帮助的讨论] [1],您还可以尝试R-forge上提供的[`lsmeans` package] [2]。 [1]:http://finzi.psych.upenn.edu/R/Rhelp02/archive/95809.html [2]:https://r-forge.r-project.org/projects/ LSMEANS / – caracal 2011-01-07 09:38:48

回答

0

考虑函数na.omit,它可以过滤掉包含NA数据行。