||
library(asreml) library(reshape) library(pastecs) data(oats) head(oats)
stat.desc(oats$yield,norm = T)
func <- function(x)(c(mean=mean(x,na.rm = T),sd=sd(x,na.rm = T),cv=sd(x,na.rm = T)/mean(x,na.rm = T)*100))cast(oats,Variety~. ,func) #分组因素为一个Variety cast(oats,Variety+Wplots~.,func) #分组因素为两个 Variety+Wplots
aggregate(oats$yield,by=list(Variety=oats$Variety),mean)#只能算一个平均数 aggregate(oats$yield, by=list(Variety=oats$Variety),func)#用自动以函数
by(oats$yield,oats$Variety,func) #结果不友好
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 00:23
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社