〖R语言〗常用统计检验方法代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

源代码+代码解释+代码应用案例+参考

   

      

      R已经成为当前国际学术界最流行的统计和绘图软件之一，该语言较为简单易学，统计分析功能强大，且具有很强的绘图功能，能够绘制学术出版要求的多种图表.R语言在生物信息学，进化生物学、生态学与环境、经济学、语言学等领域有着极为广泛的应用。

      R软件是跨平台的，可以在Linux， MacOs， Windows等多种系统上运行。针对每个研究方向，有大量的科研人员编写了相关的程序包，可以导入到基本的程序平台上运行。现有的程序包已经超过了1800个，并且还在增加中。

      不仅如此，R是完全免费的，而且全部代码是公开的。

      读者可以到 http://cran.cnr.berkeley.edu/bin/windows/base/R-2.9.0-win32.exe

      下载windows版的R软件，安装程序仅为30M。

      学习并掌握R语言，对于需要用到统计学的研究人员和学生都是非常必要的。

      这里选取了R语言中若干操作实例，所有的命令行均可以在R中运行，并得到结果。

---

1、正态总体均值的假设检验——两总体T检验

---

X为旧炼钢炉出炉率，Y为新炼钢炉出炉率，问新的操作能否提高出炉率。

命令：

1. X<-c(78.1,72.4,76.2,74.3,77.4,78.4,76.0,75.5,76.7,77.3)
   
2. Y<-c(79.1,81.0,77.3,79.1,80.0,79.1,79.1,77.3,80.2,82.1)
   
3. t.test(X, Y, var.equal=TRUE, alternative = "less")

复制代码

2、正态总体方差的假设检验

---

从小学5年级男生中抽取20名，测量其身高（厘米）如下：
问，在0.05显著性水平下，
平均值是否等于149
sigma^2 是否等于 75

命令：

1. X<-scan()
   
2. 136 144 143 157 137 159 135 158 147 165
   
3. 158 142 159 150 156 152 140 149 148 155
   
4. var.test(X,Y)

复制代码

3、二项分布的总体检验   

---

案例一
有一批蔬菜种子的平均发芽率为P=0.85,现在随机抽取500粒，用种衣剂进行浸种处理，结果有445粒发芽，问种衣剂有无效果。

命令：

1. binom.test(445,500,p=0.85)

复制代码

案例二
按照以往经验，新生儿染色体异常率一般为1%，某医院观察了当地400名新生儿，有一例染色体异常，问该地区新生儿染色体是否低于一般水平？

命令：

1. binom.test(1,400,p=0.01,alternative="less")

复制代码

4、非参数检验   —— 数据是否正态分布的Neyman-Pearson 拟合优度检验-chisq

---

   案例一
5种品牌啤酒爱好者的人数如下
A 210
B 312
C 170
D 85
E 223
问不同品牌啤酒爱好者人数之间有没有差异？

命令：

1. X<-c(210, 312, 170, 85, 223)
   
2. chisq.test(X)

复制代码

案例二
检验学生成绩是否符合正态分布

命令：

1. X<-scan()
   
2. 25 45 50 54 55 61 64 68 72 75 75
   
3. 78 79 81 83 84 84 84 85 86 86 86
   
4. 87 89 89 89 90 91 91 92 100
   
5. A<-table(cut(X, br=c(0,69,79,89,100)))
   
6. p<-pnorm(c(70,80,90,100), mean(X), sd(X))
   
7. p<-c(p[1], p[2]-p[1], p[3]-p[2], 1-p[3])
   
8. chisq.test(A,p=p)

复制代码

案例三

大麦的杂交后代芒性状的比例 无芒：长芒： 短芒=9：3：4,而实际观测值为335：125：160 ,检验观测值是否符合理论假设？

命令：

1. chisq.test(c(335, 125, 160), p=c(9,3,4)/16)

复制代码

4、参考

---

   

内容来自 薛毅 陈立萍 《统计建模与Ｒ软件》 清华大学出版社 2006

本文引用地址：http://blog.sciencenet.cn/blog-255662-240107.html  此文来自科学网张金龙博客

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：统计检验 检验方法 R语言 Berkeley Windows 统计

现在咋全是空白？？