育种数据分析之放飞自我分享 http://blog.sciencenet.cn/u/yijiaobai 关注:生物统计,数量遗传,混合线性模型,生物信息,R,Perl,Python,GWAS,GS相关方法,文章及代码

博文

R语言汇总统计,包括分组数据的汇总,结果友好(代码分享)

已有 10091 次阅读 2016-7-5 09:00 |个人分类:农学统计|系统分类:科研笔记

这个包不是必须,如果没有这个包,感兴趣的人可以下载数据来演示:oats.csv


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函数来实现对分组数据的描述性统计

aggregate(oats$yield,by=list(Variety=oats$Variety),mean)#只能算一个平均数
aggregate(oats$yield, by=list(Variety=oats$Variety),func)#用自动以函数

第四种,可以用by函数直接处理

by(oats$yield,oats$Variety,func) #结果不友好



http://wap.sciencenet.cn/blog-2577109-988728.html

上一篇:ASReml-R多性状分析,怎么设置初始值?
下一篇:表型相关与遗传相关及其显著性检验(数据及代码分享)

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-10-29 02:07

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部