library(ISwR)
attach(juul)
4.3分组数据的汇总统计量
attach(red.cell.folate)
tapply(folate,ventilation,mean)#提取folate变量,根据ventilation分组,然后对每一组计算均值
N2O+O2,24h N2O+O2,op O2,24h
316.6250 256.4444 278.0000
tapply(folate,ventilation,sd)#标准差
N2O+O2,24h N2O+O2,op O2,24h
58.71709 37.12180 33.75648
tapply(folate,ventilation,length)#长度
N2O+O2,24h N2O+O2,op O2,24h
8 9 5
#较好的显示代码
xbar<-tapply(folate,ventilation,mean)
s<-tapply(folate,ventilation,sd)
n<-tapply(folate,ventilation,length)
cbind(mean=xbar,std.dev=s,n=n)#按列将均值,标准差,长度向量合在一起