R语言与统计分析第五章习题答案.pdfVIP

# x-c(3,5,7,9,11,13,15,17,19,21) y-c(21,16,15,26,22,14,21,22,18,25) e=sum(x*y)/sum(y) #样本期望 d=(sum(x*x*y)/sum(y))-e^2 #样本方差 a=(8*e+sqrt(64*e^2-4*4*(4*e^2-12*d)))/8 #估计结果 b=(8*e-sqrt(64*e^2-4*4*(4*e^2-12*d)))/8 a # x-c(0,1,2,3,4,5,6) y-c(17,20,10,2,1,0,0) e= f-function( λ)(e^(-50* λ)* λ^50)/(2^10*6^2*24) #似然函数 optimize(f,c(0,2),maximum=TRUE) # x-c(482,493,457,471,510,446,435,418,394,469) #置信区间 (x)$ ){ options(digits=4) result-list() n-length(x) v-var(x) result$var-v chi2-(n-1)*v/var result$chi2-chi2 p-pchisq(chi2,n-1) result$-p if(alternative==less) result$-pchaisq(chi2,n-1,=F) else if (alternative==) result$-2*min(pchaisq(chi2,n-1),pchaisq(chi2,n-1,=F)) result$-c( (n-1)*v/qchisq(alpha/2,df=n-1,=F), (n-1)*v/qchisq(alpha/2,df=n-1,=T) ) result } )$ #置信区间 # x-c(25,28,23,26,29,22) y-c(28,23,30,35,21,27) )$ #卷烟 A 方差置信区间 )$ #卷烟 B 方差置信区间 #方法二 两方差比的置信区间 (x,y) #样本方差作为整体方差 options(digits=4) m=length(x) n=length(y) xbar=mean(x)-mean(y) alpha= zstar=qnorm(1-alpha/2)*(sigma1/m+sigma2/n)^(1/2) xbar+c(-zstar,+zstar) } sigma1-var(x) sigma2-var(y) options(digits=4) m=length(x) n=length(y) xbar=mean(x)-mean(y) alpha= zstar=qnorm(1-alpha/2)*(sigma1/m+sigma2/n)^(1/2) xbar+c(-zstar,+zstar) } sigma1=2140 sigma2=3250 a #置信水平为的置信上限 b[1] #置信水平为的置信下限 # x-c,,,,,,, y-c,,,,,,,, (x,y) #x方差与 y 方差的比值极大，说明 x 方差大于 y 方差 # (224,400,=,correct=TRUE) # -function(s,alpha,d,m){ t0-qt(alpha/2,m,=FALSE) n0-(t0*s/d)^2 t1-qt(alpha/2,n0,=FALSE) n1-(t1*s/d)^2 while(abs(n1-n0){ n0-(qt(alpha/2,n1,=FALSE)*s/d)^2 n1-(qt(alpha/2,n0,=FALSE)*s/d)^2 } n1 } (10,,2,100) # -function(d,p,={ alpha= ((qnorm(1-alpha/2))/d)^2*p*(1-p) } ,,