|
|
由5000个标记（Marker1-Marker5000)和两个性状T1和T2组成的矩阵，想对每个标记的不同基因型间的性状T1和T2进行方差分析，并将P值赋予该标记。想利用for循环完成，但是colnames提取的每个标记的名字是字符形式，如何将这个应用在aov（T1~Marker（1-5000））中是我一直搞不定的。
> fitA <- aov(Trait1 ~ markersname[1],data=mydata)
Error in model.frame.default(formula = Trait1 ~ markersname[1], data = mydata, :
variable lengths differ (found for 'markersname[1]')
有什么办法可以将colnames获得的名字利用在aov中？请指点。谢谢 |
|
发表于 2014-6-3 08:50:24
