非参第八次作业wdr.docVIP

非参数统计 第八次作业 王杜若 座机电话号码42 2012年月日星期 第三章习题 3.9 （1）计算Zi，Yi的相关系数 Y -c 576,635,558,578,666,580,555,661,651,605,653,575,545,572,594 Z -c 3.39,3.30,2.81,3.03,3.44,3.07,3.00,3.43,3.36,3.13,3.12,2.74,2.76,2.88,3.96 y as.matrix Y z as.matrix Z cor y,z [,1] [1,] 0.座机电话号码 （2）使用Bootsrap方法估计相关系数的标准误差 ##取有放回随机样本时，取样本量为5，即下n 5；取B 1000 B 1000 n 5 Tcor NULL SD NULL for i in 1:B + + Ysample sample Y,n,T + Zsample sample Z,n,T + TCor cor Ysample,Zsample + Tcor c Tcor,TCor + SD c SD,sd Tcor + 由上，输出结果： SD[1000] 0.座机电话号码 （3）计算置信度为0.95的相关系数Bootstrap Pivotal置信区间 Lcl 2*cor Y,Z -quantile Tcor,0.975 Ucl 2*cor Y,Z -quantile Tcor,0.025 PIVOTAL.interval c Lcl,Ucl PIVOTAL.interval 97.5% 2.5% 0.222143 1.958839 即所求置信区间为：[0.222143，1.958839] 3.10 比较三种Bootstrap置信区间的方法 ##先下载程序包moments，对偏度函数进行计算 y rnorm 50,0,1 x exp y sk.x skewness x sk.x [1] 3.38054 sk NULL sdsk NULL sd.sk NULL for i in 1:1000 xsample sample x,20,T sk0 skewness xsample sk c sk,sk0 sdsk c sdsk,sd sk alpha 0.05 ##正态方法 nLcl sk.x-qnorm 0.975,0,1 *sd sk nUcl sk.x+qnorm 0.975,0,1 *sd sk norm.inverval c nLcl,nUcl 得到区间为[0.座机电话号码， 2.座机电话号码] pLcl 2*sk.x-quantile sk,0.975 ##枢轴量法 pUcl 2*sk.x-quantile sk,0.025 pivotal.interval c pLcl,pUcl 得到区间为：[0.座机电话号码， 3.0148721] qLcl quantile sk,0.025 ##分位数法 qUcl quantile sk,0.975 得到区间为：[0.座机电话号码 ，2.座机电话号码] 3.12 ##首先，对总体的对称性进行检验。 median.long median long median.long [1] 71 wilcox.test long-median.long Wilcoxon signed rank test with continuity correction data: long - median.long V 269.5, p-value 0.6793 alternative hypothesis: true location is not equal to 0 由上可知，不拒绝原假设，可以认为总体是对称分布的。 ##下面，用Walsh平均法对所给长度样本求中位数95%置信区间 for i in 1:length long -1 for j in i+1 :length long walsh.long -c walsh.long, long[i]+long[j] /2 walsh.long -c walsh.long,long n length walsh.long alpha 0.05 for k in seq 1,n/2,1 F pbinom n-k,n,0.5 -pbinom k,n,0.5 if F 1-alpha a k-1 break sort.walsh.long sort walsh.long Lower sort.walsh.long[a] Upper sort.walsh.long[n-a+1] 得到[71,71.5] ,