R语言第三章课后习题答案.docxVIP

##附加题斐波那契数列##：(1/√5)*{[(1+√5)/2]^n - [(1-√5)/2]^n}（又叫“比内公式”，是用无理数表示有理数的一个范例。）（√5表示根号5）rm(list=ls(all=TRUE))##方法1fei1-function(a,b,lim=5){ #递归，lim为要求输出的第lim个数字，PS：初始两个不算fun1-function(a,b,lim,fun){c=a+bif (lim==0) return(c)else {a=bb=creturn(fun(a,b,lim-1,fun))}}fun1(a,b,lim,fun1)}##方法2fei2-function(a,b,s){ #输出的第s个数for(i in 1:s){ c=a+ba=bb=c}return(c)}##3.1data-scan(C:\\3.1.txt)##均值mean(data) ##方差var(data)##标准差sd(data)^2##极差max(data)-min(data)##标准误sm-function(x){n-length(x)sm-sd(x)/sqrt(n)return(sm)}##变异系数CV-function(x){cv=(sd(x)/mean(x))*100return(cv)}##偏度G1-function(x){n-length(x)m-mean(x)s-sd(x)g1-n/((n-1)*(n-2))*sum((x-m)^3)return(g1)}##峰度G2-function(x){n-length(x)m-mean(x)s-sd(x)g2-((n*(n+1))/((n-1)*(n-2)*(n-3))*sum((x-m)^4)/s^4-(3*(n-1)^2)/((n-2)*(n-3)))return(g2)}##3.2data-scan(C:\\3.1.txt)##直方图hist(data,freq=F)##密度估计曲线lines(density(data),col=blue)x-60:88lines(x,dnorm(x,mean(data),sd(data)),col=red)##经验分布图plot(ecdf(data),verticals=T,do.p=F)x-60:88lines(x,pnorm(x,mean(data),sd(data)))##QQ图qqnorm(data)qqline(data)##3.3data-scan(C:\\3.1.txt)##茎叶图stem(data)##箱线图boxplot(data,notch=T,col=red)##五数总括fivenum(data)##3.4##W检验方法data-scan(C:\\3.1.txt)shapiro.test(data)##Kolmogorov-Smirnov检验方法ks.test(data,norm,mean(data),var(data))##3.5a-c(2,4,3,2,4,7,7,2,2,5,4)b-c(5,6,8,5,10,7,12,12,6,6)c-c(7,11,6,6,7,9,5,5,10,6,3,10)##1boxplot(a,b,c,notch=T,names=c(A,B,C),col=c(2,3,4))##2y-c(2,4,3,2,4,7,7,2,2,5,4,5,6,8,5,10,7,12,12,6,6,7,11,6,6,7,9,5,5,10,6,3,10)f-factor(c(rep(1,11),rep(2,10),rep(3,12)))plot(f,y)##生成日历rm=(list=ls(all=TRUE))##M=;W=;month-function(m){if(m==1|m==3|m==5|m==7|m==8|m==10|m==12)D=31else if(m==4|m==6|m==9|m==11)D=30elseD=28D}bt-c(Sun,Mon,Tue,Wed,Thr,Fri,Sat)yueli-function(m,w){a-month(m)if(w==7) w-1if(w!=1w7) w-w+1 n-1s-5if (w=6a==31) s-6if (w==7a!=28) s-6x-array(NA,dim=c(s,7))x[n,w]-1for(i in 2:a){w=w+1if(w==8) w-1if(w==1) n=n+1x[n,w]-i}rbind(bt,x)}