原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
非参数统计(第二版)习题R程序
P37.例2.1build.price-c(36,32,31,25,28,36,40,32,41,26,35,35,32,87,33,35);build.pricehist(build.price,freq=FALSE)#直方图lines(density(build.price),col=red)#连线#方法一：m-mean(build.price);m#均值D-var(build.price)#方差SD-sd(build.price)#标准差St=(m-37)/(SD/sqrt(length(build.price)));t#t统计量计算检验统计量t=[1] -0.1412332#方法二：t.test(build.price-37)#课本第38页例2.2binom.test(sum(build.price37),length(build.price),0.5)#课本40页例2.3P-2*(1-pnorm(1.96,0,1));P[1] 01-2*(1-pnorm(0.7906,0,1));P1[1] 0.4291774 例2.4 p-2*(pnorm(-1.96,0,1));p[1] 0 p1-2*(pnorm(-0.9487,0,1));p1[1] 0.3427732例2.5（P45）scores-c(95,89,68,90,88,60,81,67,60,60,60,63,60,92,60,88,88,87,60,73,60,97,91,60,83,87,81,90);length(scores)#输入向量求长度ss-c(scores-80);sst-0t1-0for(i in 1:length(ss)){if (ss[i]0) t-t+1#求小于80的个数else t1-t1+1求大于80的个数}t;t1 t;t1[1] 13[1] 15binom.test(sum(scores80),length(scores),0.75)p-value = 0.0014360.01Cox-Staut趋势存在性检验P47例2.6year-1971:2002;yearlength(year)rain-c(206,223,235,264,229,217,188,204,182,230,223,227,242,238,207,208,216,233,233,274,234,227,221,214,226,228,235,237,243,240,231,210)length(rain)#(1)该地区前10年降雨量是否变化？t1=0for (i in 1:5){if (rain[i]rain[i+5]) t1-t1+1}t1k-0:t1-1sum(dbinom(k,5,0.5))#=0.1875y-6/(2^5);y#=0.1875#(2)该地区前32年降雨量是否变化？t=0for (i in 1:16){if (rain[i]rain[i+16]) t-t+1}tk1-0:min(t,16-t)-1sum(dbinom(k1,16,0.5))#=0.0002593994pbinom(max(k1),16,0.5)#= 0.0002593994y1-(1+16)/(2^16);y1#=0.0002593994plot(year,rain)abline(v=(1971+2002)/2,col=2)lines(year,rain)anova(lm(rain~(year)))随机游程检验（P50）例2.8client-c(F,M,M,M,M,M,F,M,M,F,M,M,M,M,F,M,F,M,M,M,F,F,F,M,M,M);clientn-length(client);nn1-sum(client==M);n1n0-n-n1;n0t1-0for (i in 1:(length(client)-1)){if (client[i]==client[i+1]) t1-t1else t1-t1+1}R-t1+1;R#=12#find rejection region（不写）rl-1+2*n1*n0/(n1+n0)*(1-1.96/sqrt(n1+n0));rlru-2*n1*n0/(n1+n0)*(1+1.96/sqrt(n1+n0));ru#=15.33476（课本为ru=17）例2.9shuju39-data.frame(read.table(SHUJU39.txt,header=TRUE));shuju39attach(shuju39)sum.a=0sum.b=0sum.c=0for (i in 1:length(id)){if (pinzhong[i]==A) sum.a-sum.a+chanliang[i]
文档评论（0）