原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
12营养2班生统复习资料
制作者：冯波
R语句集合：
变量- 数据
x - 1:5 #生成整数数列
x - seq(1, 5, 0.1) #生成数列
x - rep(5:6, times=10) #重复数列
x - rep(5:6, times=c(3,2)) #重复数列
x - rep(4:6, times=1:3) #重复数列
x - rnorm(10, 0, 1) #生成随机数字
x - c(480,490,500,505,510,520,535,560,585,600) #录入数据
median(x) #中位数
mean(x) #算术平均数
table(x) #众数
which.max(table(x)) #众数
range(x) #范围
var(x) #方差
sd(x) #标准差
summary(x) #简单统计描述
max(x),min(x),sum(x)
二项分布概率函数：dbinom(k, n, p)
二项分布分布函数：pbinom(k, n, p)
单个总体均数的假设检验
t检验： days - c(116,115,113,112,114, 117, 115, 116, 114, 113) #录入数据
t.test(days, mu = 114) #假设检验
t检验(单边)： vc - c(255, 260, 262, 248,244, 245,250, 238, 246, 248, 258, 270)
t.test(vc, mu = 246,alternative=greater) #算小于就把greater换成less两样本平均数的假设检验
1.非配对样本
bf1 - c(1.2, 1.32,1.1,1.28,1.35,1.08,1.18,1.25,1.3,1.12,1.19,1.05)
bf2 - c(2,1.85,1.6,1.78,1.96,1.88,1.82,1.7,1.68,1.92,1.8) #录入数据
t.test(bf1,y=bf2,var.equal=T) #假设检验1
var.test(bf1,bf2) #方差同质性检验
t.test(bf1,y=bf2,var.equal=F) #假设检验2
此案例假定两样本方差相同? 使用假设检验1
若两样本方差未知，需进行方差同质性检验，
p0.05 ? 方差不等? 假设检验2
p0.05 ? 方差相等? 假设检验1
2.配对样本
tmp1 - c(37.8,38.2,38.0,37.6,37.9,38.1,38.2,37.5,38.5,37.9)
tmp2 - c(37.9,39,38.9,38.4,37.9,39,39.5,38.6,38.8,39)
t.test(tmp1,y=tmp2,paired=TRUE)
单因素方差分析
ls - c(8,13,12,9,9, 7,8,10,9,7, 13,14,10,11,12, 13,9,8,8,10, 12,11,15,14,13)
bred - as.factor(rep(1:5,each=5))
case1 - aov(ls~bred)
summary(case1)
TukeyHSD(case1) #多重比较
数据转换
rate-c(94.3,64.1,47.7,43.6,50.4,80.5,57.8,26.7,9.4,42.1,30.6,40.9,18.6,40.9,18,35,20.7,31.6, 26.8,11.4,19.7)/100
area - as.factor(rep(1:3,each=7))
rate2 - asin(sqrt(rate)) #反正弦转换
case1 - aov(rate~area)
case2 - aov(rate2~area)
summary(case1)
summary(case2)
plot(density(rate),col=blue) #绘图
lines(density(rate2),col=red) #添加曲线
两因素方差分析
dog-c(45,69,53,51,51,50,62,68,52,18,25,32,54,72,69,66,53,63,67,70,51,59,47,42,91,87,89,91,69,73,77,74,66,68,70,64)
train -as.factor(rep(rep(c(‘A’,‘B’,‘C’),each=4),3))
prize - as.factor(rep(c(A,B,C),each=12))
case1 - aov(dog~train+prize+train*prize)
summary(case1)
TukeyHSD(cas
文档评论（0）