非参数统计王星版第一章课后答案.pdf

非参数统计王星版第一章课后答案非参数统计王星版第一章课后答案 非参数统计第一章课后答案 #1.1 AGE=c(18,23,22,21,20,19,20,20,20) first3=sort(AGE)(length(AGE)-2):length(AGE) delete.1=order(AGE)(length(AGE)-2):length(AGE) except=AGE-delete.1 c(except1:2,19,except3:length(except) except2=22 #1.2 a1=rep(1:3,rep(2,3) a2=c(1,8,10,11) a3=seq(1,30,length(a2) a4=seq(1,5,2) a4=c(a1,a4,rep(0,2) a5=2:10 a6=c(a2,a3-(1:3),a4) a7=c(rep(1,10),rep(0,8) #1.3 library(MASS) data(geyser) a=geyser a1=subset(a,waiting70) a1=geysergeyser$waiting70,#两个 a1 等价 a2=subset(a,waiting70waiting!=57) a2=geysergeyser$waiting70geyser$waiting!=57,#两个 a2 等价 a3=subset(a,waiting70,c(duration) a4=subset(a,duration70,c(waiting) #1.3 library(MASS) a=geyser attach(a) b1=waitingwaiting70 b2=waitingwaiting70waiting!=57 b3=durationwaiting70 b4=waitingduration70 #1.4 x=c(0,1,1,2,3,4) num=function(x) r=0 p=c(rep(0,length(x)-1) q=c(rep(0,length(x)-1) for(i in 1:(length(x)-1) for(j in (i+1):length(x) pi=pi+I(xixj) qi=qi+I(xixj) r=r+pi-qi r h=rep(0,1000) for(i in 1:length(h) x=runif(5,-5,5) hi=num(x) y=as.factor(h) yy=levels(y) a=rep(0,11)#记数 记录-10:10:2 各种结果出现的次数 p=rep(0,11) for(i in 1:11) ai=sum(h=(-12+2*i) pi=ai/length(h) a p hist(h) h=rep(0,10000) for(i in 1:length(h) x=sample(1:5,5,replace=T)#-5 到 5 中间的整数随机抽样 hi=num(x) y=as.factor(h) yy=levels(y) a=rep(0,length(yy)#记数 记录-10:10 各种结果出现的次数 p=rep(0,length(yy) for(i in 1:length(yy) ai=sum(h=(-11+i) pi=ai/length(h) a p hist(h) #当随机取 10000 次的一个结果a=71 321 774 1255 1637 1825 1684 1256 743 338 96 # p=0.0071 0.0321 0.0774 0.1255 0.1637 0.1825 0.1684 0.1256 0.0743 0.0338 0.0096 #当随机取十万次数据的一个结果 #a=79534217553 12521 16771 18180 16538 125537418 3400850 #p=0.007950.034210.075530.125210.167710.18180 0.16538 0.12553 0.07418 0.03400 0.00850 #1.5 uniroot(f=function(x) 2*x3-4*x2+3*x-6, interval=c(-10,10) f=function(x)2*x3-4*x2+3*x-6 f(0) a=-10 b=10 root=function(a,b) c=(a+b)/2; while(abs(f(c)0.00001) if(f(c)*f(a)0)b=c; c=(a+b)/2; else a=c; c=(a+b)/2; c #1.6 x=seq(0,2*pi,0.2) y=sin(x)/(cos(x)+x) #1.7 chartonum=function(x) a=c(abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRS TUVWXYZ); b=strsplit(a,); for(i in 1:52) if(b1i=x)t=i;i=i+1 elsei=i+1 t #将字符转化为数字，小写为前 26 位，大写为后 26 位，输入为单个字 符 f7=function(x) y=strsplit(x,);#将输入分为单个字符 for(i in 1:length(y1) t=chartonum(y1i) if(t14)t=t+13;y1i=LETTERSt; else if(t=14t=26)t=t-13;y1i=LETTERSt; else if(t=27t=39)t=t+13- 26;y1i=letterst; elset=t-13-26;y1i=letterst; mima=y; mima #1.8 f1.81=function(a) for(i in 1:length(a) ci=alength(a)-i+1 c f1.82=function(a,b) for(i in 1:(length(a)+length(b) if(i%2=1)ci=a(i+1)/2; elseci=bi/2; c #%为求余符号 #1.9 f1.9=function(n,m) a=rep(1,n); i=1; k=0; t=1; while(sum(a)1|at=0) t=i%n; if(t=0)t=n; elset=i%n if(at!=0) k=k+1;i=i+1; elsei=i+1 if(k%m=0)k=0;at=0; t #1.10 student-read.table(C:Documentsand SettingsAdministrator桌面非参数统计配套数据各章数据第 一章student.txt,header = T) student1=as.data.frame(student) means=apply(student1,2:6,1,mean) b1=data.frame(student1,means) b2=student1which(b=max(b), b3=I(student1,2:660)#判断每个学生每门课程是否及格 b4=apply(b3,1,sum)#找出每个学生有几门课程不及格 b5=b1which(b41), b6=b1-which(b41), t.test(b5$means,b6$means) #不说明情况下认为两总体方差不同，特殊说明 var.equal=TRUE 置信区 间为均值差的置信区间 wilcox.test(b5$means,b6$means) wilcox.test(b5$means,b6$means,paires=FALSE) #1.11 basket-read.table(C:Documentsand SettingsAdministrator桌面非参数统计配套数据各章数据第 一章basket.txt,header = T) A=I(basket,2:6=A) A=basket,2:6=A Asum=apply(A,1,sum) B=basket,2:6=B Bsum=apply(B,1,sum) a1=basketwhich(Asum0Bsum0),1 length(a1)#求个数 a2=basketwhich(Asum0Bsum=3), length(a2$ID)#求个数 #1.12 x=seq(-10,10,0.05) y1=sin(x) y2=cos(x) y3=y1+y2 plot(x,y1,col=1,lty=1,axes=T,main=习题 1.12,sub=sin（）,cos（）,sin()+cos()对比图) points(x,y2,col=2,lty=3) points(x,y3,col=3,lty=4) plot(x,y1,col=1,lty=1,main=习题 1.12,sub=sin（）,cos（）,sin()+cos()对比图) lines(x,y1,col=1,lty=1) lines(x,y2,col=2,lty=3) lines(x,y3,col=3,lty=4) curve(sin(x),-10,10,col=1,lty=1,main=习题 1.12,sub=sin（）,cos（）,sin()+cos()对比图) curve(cos(x),add=TRUE,col=2,lty=3) curve(sin(x)+cos(x),add=TRUE,col=3,lty=4) #1.13 x=seq(-5,5,length=50) a=runif(500,-5,5) y=0.1*a*sin(2*a) f1=function(x,y)1-exp(-1/x2+y2) z1=outer(x/2,x/2,f1) persp(z1) f2=function(x,y)0.1*x*sin(2*y) z2=outer(x,x,f2) persp(z2) f3=function(x,y)sin(x)+cos(x) z3=outer(x,x,f3) persp(z3) plot(sin(3*x),sin(6*x),type=l) #1.14 a=rnorm(100,3,sqrt(5) b=rnorm(20,5,sqrt(3) c=c(a,b) hist