统计软件期中考试A(2013级)答案.docx

学号： 姓名： 成绩： 统计软件期中试卷一、填空与选择（每小题2分，共40分）1. R软件编制的程序文件后缀名为 *.R 2. R软件命令窗口的命令提示符为 3. 安装程序包ggplot2的命令为 install.packages(ggplot2) 4.下列变量名中的 A 是合法变量（A）char_1,data_2,x_j (B)x*y ,a.1 (C)xy,a1234 (D)end,2bcx5.已知x=2(0:3),则命令2*x返回的结果为 2 4 8 16 6.命令2*3:7 返回的结果为 6 8 10 12 14 7.命令rep(c(1,3),each=3)返回结果为 1 1 1 3 3 3 8.命令sort(c(2, 10,7,9,5) 返回结果为 2 5 7 9 10 9. 命令matrix(2:7,nrow=2,byrow=F)返回结果为 2 4 6 3 5 7 10.已知A=1 32 4矩阵，则命令A%*%A返回的结果为 7 15 10 22 11. 已知A=1 32 4矩阵，则命令dim(A)返回的结果为 2 2 12. 已知矩阵A=1 32 4，则命令apply(A,2,sum)返回的结果为 3 7 13.命令rbind(1:3,1:4)返回的结果为 1 2 3 1 1 2 3 4 14执行命令b-list(name=Fred,wife=Mary,no.child=3,child.ages=c(4,7,9)后，则命令 b$name 或b1 返回的结果为Fred15. 执行命令b-list(name=Fred,wife=Mary,no.child=3,child.ages=c(4,7,9)后，则命,令 length(b)返回的结果为 4 16.执行命令xx=array(1:24,c(3,4,2)后，则命令xx2,1:3,2返回的结果为 1 4 17 20 17. 执行命令a=data.frame(class=c(1,1,2,2,3,3,3),sex=c(女,男,男,女,女,女,男),score=c(80,85,92,76,61,95,83)后，则命令a$scorea2=男返回的结果为 85,92, 83 18.已知变量x=c(2,3),则执行命令x5-5后变量x的结果为 2 3 NA NA 5 19.已知变量x=1:10+2,则命令x1:5返回的结果为 3 4 5 6 7 20.命令paste(x,1:5,sep=)的返回结果为 x1 x2 x3 x4 x5 二、论述题（10分） 谈谈你学习统计软件的真实体会、感受、建议与意见？要求不少于300字。三、编程题（共52分, 只需交程序，所有程序(每个小题前加上注释说明)写在一个R程序文件中，文件名以你的学号+姓名命名，最后传到06中的“统计软件R”文件夹下的“期中考试”目录下）1. （6+10+6+10=32分）现有7位成人的姓名、性别、年龄、身高（cm）、体重(kg)数据如下：Name Sex Age Height WeightAlice F 23 136.5 42.0 （你的学号为奇数，去掉第一行）Sandy F 21 141.3 50.5 Sharon F 25 152.5 60.5Alfred M 24 169.0 62.5Duke M 24 170.5 55.0Guido M 25 167.0 60.0Robert M 12 164.8 50.0 （你的学号为偶数，去掉最后一行）(1) 用数据框的形式读入数据，数据框名称为你的姓名拼音的首字母（如周丽：ZL）(2) 将上述数据写成一个纯文本文件（文件名为你的姓名拼音的首字母+学号后4位数字），并用函数read.table()读取该文件，筛选出体重指数BMI过小或过大的学生信息。注：体重指数（BMI）=体重（kg）身高（m）的平方，如：70kg（1.751.75）=22.86，BMI值介于2024之间为正常值。(3) 在一个窗口中绘制两幅图：（a）身高为横轴、体重为纵轴的散点图，横、纵坐标轴设置为中文的，标题为“散点图”（b）将上题中男生（M）的数据点、女生（F）的数据点，分别用绿色和红色的空心三角标志出来，并在图上绘制线性回归拟合直线（如果你是女生请用实线，线宽度为2，否则用虚线，线宽度为1.5）。程序：Name-c(Alice,Sandy,Shron,Alfred,Duke,Guido)Sex-c(F,F,F,M,M,M)Age-c(23,21,25,24,24,25)Height-c(136.5,141.3,152.5,169.0,170.5,167.0)Weight-c(42.0,50.5,60.5,62.5,55.0,60)ZL-data.frame(Name,Sex,Age,Height,Weight)write.table(ZL,file=E:/Rprogram/data/ZL1234.txt)D-read.table(E:/Rprogram/data/ZL1234.txt,header=TRUE)D$BMI-D$Weight/(D$Height/100)2subset(D,BMI24)plot(x=Height,y=Weight,xlab=身高,ylab=体重,main=散点图)ZL$mycolSex=F-redZL$mycolSex=M subset(D,BMI24) Name Sex Age Height Weight BMI2 Sandy F 21 141.3 50.5 25.293393 Shron F 25 152.5 60.5 26.014515 Duke M 24 170.5 55.0 18.919692. （8分）生成一个5阶的Hilbert矩阵H=hijnn,hij=1i+j-1,i,j=1,2,n(1) 计算Hilber矩阵H的行列式(2) 求H的逆矩阵程序:n-5H-matrix(ncol=n,nrow=n)for (i in 1:n) for (j in 1:n) Hi,j H ,1 ,2 ,3 ,4 ,51, 1.0000000 0.5000000 0.3333333 0.2500000 0.20000002, 0.5000000 0.3333333 0.2500000 0.2000000 0.16666673, 0.3333333 0.2500000 0.2000000 0.1666667 0.14285714, 0.2500000 0.2000000 0.1666667 0.1428571 0.12500005, 0.2000000 0.1666667 0.1428571 0.1250000 0.1111111 det(H)1 3.749295e-12 solve(H) ,1 ,2 ,3 ,4 ,51, 25 -300 1050 -1400 6302, -300 4800 -18900 26880 -126003, 1050 -18900 79380 -117600 567004, -1400 26880 -117600 179200 -882005, 630 -12600 56700 -88200 441003. （10分）鸡兔同笼问题:鸡和兔关在一个笼子里，已知共有头36个，脚100个，编制自定义函数，求笼内关了多少只兔和多少只鸡。JT-function(T,J) for(x