world_6520的博客

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,