pheno=read.table("XXLLpheno.txt",sep="",head=F)
y<-pheno[,8]//y是表型值
ffarm<-pheno[,2]//场固定效应
fsex<-pheno[,3]//性别固定效应
fyear<-pheno[,4]//年固定效应
fmonth<-pheno[,5]//月固定效应
fseason<-pheno[,6]//季节固定效应
fit<-aov(y~ffarm+fsex+fyear+fmonth+fseason,data=pheno)//方差因素分析
summary(fit)//结果展示
R做方差因素分析
最新推荐文章于 2024-07-07 17:31:45 发布