非参数统计(R软件)参考答案

非参数统计(R软件)参考答案用时可以删除sergeyserwhichwaiting#提取满足条件(waiting<70)的数据,which()，读取下标Sub2geyser=geyser[which((waiting<70)&(waiting!=57)),1];#提取满足条件(waiting<70&(waiting!=57)的数据.Subgeysergeyserwhichwaiting<70),2];#提取满足条件(waiting<70)的数据,which()，读取下标Subgeysergeyserwhichwaiting>70),1];#提取满足条件(waiting<70)的数据,which()，读取下标差异。eanATthphysicschemliteratenglishKatty28479Leo6455merep30);foriinxstudentimeimeanx;}以如此:student[,7]=medenteantmeathphysicschemliteratenglishmeanLiggle76Datastudent#赋初值foriinxstudentiakstudentiathphysicschemliteratenglishmeanRicky57Simon71675257Jed83100794150Jack55Jetty56Corner52Osten52Amon747995595901data2=student[1,];k=0;ofxmeanofy在一张图上，用取值(-10，10)之间间隔均等的1000个点，采用不同的线型一颜色给制sincossincos有主标题和副标题，标示出从坐标xlabxylabyylimctypelcolZZ01理由：>1000是批发商的意愿，违背这个意愿，也就是拒绝原假设H，他就购这批灯泡0了。不能轻易否定的事情应置于被保护地位H。这个问题的检验统计量为0==01H00H,怪不得我也。面此假设违返旧过程，这样的假设毫无意义。001001情况，不要联系(2)中的具体数据例子)，如果没有给定水平，如何用p值来做出结论12n H。则有100100a/200100a00100a装2装12n12n12n12na，BaxswaxHydrocarboncdfhistHydrocarbonrbonqqlineHydrocarbon图图MeltingPointcdfhistMeltingPointPointqqlineMeltingPoint在两条直线的附近。所以两近似正态分布(高斯分布)。ntP据呈正态分布)。程序如下：,,,2,81677984429设1e101Xz2.672612123456789011x=c(91,46,108,99,110,105,191,57,34,81);y=c(81,51,63,51,46,45,66,64,90,28);sg=sum(z>0);sl=sum(z<0);n1=sg+sl;k=min(sg,sl)Exactbinomialtestdata:kandn1numberofsuccesses=3,numberoftrials=10,p-value=alternativehypothesis:trueprobabilityofsuccessisnotequalto95percentconfidenceinterval:0.sampleestimates:probabilityofsuccesscoxonsignedranktestV=45,p-value=vehypothesistruelocationisnotequalto中，由于数据的分布不存在显著不对称的迹象，是可靠的，因而理好。事211154113413306您能否同意所声称的12岁的这种鱼的长度的中位数总是在69~72cm之间xc0,70,70,70,71,71,71,71,71,72,72,72,73,733);shNULLnlengthxforjininwxixjwalshcwalshwhashNULLnlengthxforjininwxixjwalshcwalsh,w);}}sortwalshwalsh的置信区间，对称地砍掉左尾和右尾forkinseqnwFpbinomnwknw-pbinom(k,nw,;ilklkuknwlk(2)其它置信区间，基于Bootstrap方差的枢轴区间是最好的，它是(69，73)，还是没有家欲了解抑郁症的发病率是否在一年时间随季节的不而不同，他使用了来年一所大医季节春季夏季秋季冬季合计50349150349195病率是否与季节有关0123411234V=c(495,503,491,581);p=1/4;n=sum(V);df=4-1;chi2=sum((V-n*p)^2/(n*p))pvalue=1-pchisq(chi2,df);pvalue;#请思考：为什么用右尾概率结论：在a=时拒绝原假设，认为发病率与季节有关。具体地说，冬天的发病率高(p=)。3当然，为了要得到科学的结论，应该要规范抽样，使得样本有代表性，毕竟一个医院的数据其代表性是值得商榷的。P106:;;.进行攻击性测试，测量得分显示在如下表中(得分越高表示攻击性越强)设其位置参数(均值或中位数)是无显著差异，即检验假设为：.0AB1ABAc0,8,12,16,5,9,7,11,6);Bc5,20,18,13,14,9,16);minmincAB);max=max(c(A,B));plotAtypebpchAxlim=c(0,9),ylim=c(min,max));箱线图，程序如下，图如图"A","B"),c(9,8)))cABgrouptdata:AtD=,p-value=1data:AandBalternativehypothesis:truedifferenceinmeansislessthan0percentconfidenceinterval:-InftimatesmeanofxmeanofywosampleKolmogorovSmirnovtestedianAandBmedianBD=,p-value=ativehypothesistwosidedWilcoxonranksumtestwithcontinuitycorrectionW=,p-value=alternativehypothesis:truelocationshiftislessthan0两个不同学院教师一年的课时量分别为(单位：学时)本，两个学院教师讲课的课时是否存在显著差异估计这些差异。从两个学院教情况是否类似”，再选择检验方法，推断是否存在显著Ac21,266,256,386,330,329,303,334,299,221,365,250,258,342,243,298,238,317);Bc593,507,428,807,342,512,350,672,589,665,549,451,492,514,391,366,469);(1)检验“教师完成讲课任务的情况是否类似”1ABAB01orrectionWp-value=tisnotequaltoplettestalue,,,,,,,,,,,,,,,,,,,,12m112n201211212m012m12nim+n2ii=1data:AandBdataSSAandSSBpleynyTVbyvarybVBycVByVby;MVxmeanVBx;MVy=mean(VBy);VarxyvarVBxvarVBy;ppnormZppnormZ);pvalueminpppvalue此结果与(2)相同。C解：(1)三个样本A、B、C均为独立随机样本，非区组试验数据，样本量不同，只能用H:三样本的中位数相同VSH:三样本的中位数不全相同01A=c(83,64,67,62,70);B=c(85,81,80,78);C=c(88,89,79,90,95);n1=length(A);n2=length(B);n3=length(C);BCgroupfactorrepcnn2,n3)));Kruskal-WallisranksumtestoupKruskal-Wallischi-squared=,df=2,p-value=C(2)进一步分析差异出自何处，请看箱线盒须图：txgroup(3)两两比较的程序和结果如下Ac3,64,67,62,70);B=c(85,81,80,78);C=c(88,89,79,90,95);nlengthAnlengthBn3=length(C);k=3;ncnnnalphaalphas=alpha/(k*(k-1));Z=qnorm(alphas,0,1);Nsumn;MST=N*(N+1)/12;xcABCRrankxRbarrep,k);upfactorrepcnnnriinkRbarimedianRgroupiforiinkforjinik){SE=sqrt(MST*(1/n[i]+1/n[j]));bariRbarjSElengthddsigrepndforiinndifdiZdsigi]=1}}112345ABCAc,;cngthAnlengthBnlengthCMmatrixxbyrowT;manranksumtestFriedmanchisquareddf=2,p-value=(2)进一步分析差异出自何处，请看箱线盒须图：p(3)两两比较的程序与结果如下A=c,,,,;B=c,,,,;Cc,,;n1=length(A);n2=length(B);n3=length(C);k=3;n=c(n1,n2,n3);alpha=;alphas=alpha/(k*(k-1));Z=qnorm(alphas,0,1);N=sum(n);MST=N*(N+1)/12;x=c(A,B,C);R=rank(x);Rbar=rep(0,k);group=factor(rep(1:3,c(n1,n2,n3)));for(iin1:k){Rbar[i]=median(R[group==i])}for(iin1:(k-1)){for(jin(i+1):k){SE=sqrt(MST*(