R做方差因素分析

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)//结果展示