R软件论文R软件论文.docVIP

R软件课程作业 指导老师： 毛树华 学院2014 —— 2015年 第 2 学期 随堂作业 1、R软件的基本操作 1.2（1）提取已婚人士收入数据： 程序：zxf1=subset(zxf,婚姻状况==已婚) v=zxf1[,10] v （2）计算其均值、方差等： 程序及结果：mean(v) [1] 393.7317 min(v) [1] 120 max(v) [1] 800 sd(v) [1] 195.4197 var(v) [1] 38188.85 hist(v) 1.3（1）高中：zxf2=subset(zxf1,文化程度==高中) v2=zxf2[,10] mean(v2) [1] 516.6667 var(v2) [1] 40845.28 （2）初中：zxf3=subset(zxf1,文化程度==初中) v3=zxf3[,10] mean(v3) [1] 364.92 var(v3) [1] 33786.19 （3）小学：zxf4=subset(zxf1,文化程度==小学) v4=zxf4[,10] mean(v4) [1] 284 var(v4) [1] 5817.931 2、数据描述性分析 习题3.1 setwd(e:/) book1=read.table(book1.txt,header=TRUE) book1 v=book1[,1] m=mean(v);m [1] 73.66162 s=sd(v);s [1] 3.95845 R=max(v)-min(v);R [1] 20 sm=s/sqrt(v);sm [1] 0.4459253 0.4772335 0.4482062 0.4717793 0.4411916 0.4411916 0.4741424…… [99] 0.4717793 cv=100*s/m;cv [1] 5.37383 css=sum((v-m)^2);css [1] 1535.594 uss=sum(v^2);uss [1] 538712.9 g1=n/((n-1)*(n-2))*sum((v-m)^3)/s^3;g1 [1] 0 g2=(n*(n+1))/((n-1)*(n-2)*(n-3))*sum((v-m)^4)/s^4;g2 [1] 3.102422 习题3.2 （1）hist(v) （2）par(mfrow=c(2,2)) hist(v,freq=FALSE) lines(density(v),col=blue) v=64:85 lines(v,dnorm(v,mean(v),sd(v)),col=red) （3）plot(ecdf(v),verticals=TRUE,do.p=FALSE) v=64:85 lines(v,pnorm(v,mean(v),sd(v))) （4）qqnorm(v);qqlines(v) 3、参数估计 习题4.5 x=c(54,67,68,78,70,66,67,70,65,69) source(E:/R软件/R上机/interval_estimate.R) interval_estimate(x) mean df a b 1 67.4 9 63.1585 71.6415 习题4.6 x=c(140,137,136,140,145,148,140,135,144,141) y=c(135,118,115,140,128,131,130,115,131,125) source(E:/R软件/R上机/interval_estimate2.R) interval_estimate2(x,y) mean df a b 1 13.8 13.01367 7.359713 20.24029 习题4.10 x=c(1067,919,1196,786,1126,936,918,1156,920,948) source(E:/R软件/R上机/interval_estimate3.R) interval_estimate3(x) mean df a_oneside b_oneside 1 997.2 9 921.048 1073.352 4、假设检验 习题5.4 x=c(-0.7,-5.6,2,2.8,0.7,3.5,4,5.8,7.1,-0.5,2.51,-1.6,1.7,3,0.4,4.5,4.6,2.5,6,-1.4) y=c(3.7,6.5,5,5.2,0.8,0.2,0.6,3.4,6.6,-1.1,6,3.8,2.1,1.6,2,2.2,1.2,3.1,1.7,-2) shapiro.test(x) shapiro.test(y