R语言实验报告习题详解

1、言实验习题详解学院:班级:学号:姓名:导师:成绩:目录一、实验目的二、实验内容1.1问题叙述1.2问题求解1.3结果展示2.1问题叙述2.2问题求解2.3结果展示3.1问题叙述3.2问题求解3.3结果展示4.1问题叙述4.2问题求解4.3结果展示5.1问题叙述5.2问题求解5.3结果展示6.1问题叙述6.2问题求解6.3结果展示三、实验总结一、实验目的R是用于统计分析、绘图的语言和操作环境。R是属于GNU系统的一个自由、 免费、源代码开放的软件，它是一个用于统计计算和统计制图的优秀工具；本次试验要求掌握了解R语言的各项功能和函数，能够通过完成试验内容对 R语言有一定的了解，会运用软件对数据进行

2、分析；通过本实验加深对课本知识的理解以及熟练地运用 R语言软件来解决一些 复杂的冋题。二、实验内容1.1问题叙述将1,2，20构成两个4X 5阶的矩阵，其中矩阵A是按列输入，矩阵B 是按行输入，并做如下运算.C=A+B;D=A*B;F是由A的前3行和前3列构成的矩阵；G是由矩阵B的各列构成的矩阵，但不含 B的第3列.1.2问题求解1.2.1创建按列、行输入的4X5矩阵;Av-matrix(1:20,c(4,5)AB-matrix(1:20, nrow=4,byrow=TRUE)B程序求解A-matrix(1:20,c(4,5)AB-matrix(1:20, nrow=4,byrow=TRUE)

3、BC=A+BCD=A*BDF-A1:3,1:3FHv-matrix(c(1,2,4,5), nrow=1)HG-B,HG1.3结果展示2.1问题叙述已知有5名学生的数据，如下表所示.用数据框的形式读入数据.学生数据序号姓名性别年龄身高（cm）体重（kg）1张三女1415642.02李四男1516549.03王五女1615741.54赵六男1416252.05丁一女1515945.52.2问题求解StudentData 数据框Stude ntDatav-data.frame( name=c(zha ngsa n,lisi,wa ngwu,zhaoliu, din gyi),sex=c(F,M,F

4、,M,F),age=c(14,15,16,14,15),h eight=c(156,165,157,162,159),weight=c(42,49,41.5,52 ,45.5)程序求解Stude ntDatav-data.frame( name=c(zha ngsa n,lisi,wa ngwu,zhaoliu, din gyi),sex=c(F,M,F,M,F),age=c(14,15,16,14,15),h eight=c(156,165,157,162,159),weight=c(42,49,41.5,52 ,45.5)Stude ntData2.3结果展示3.1问题叙述某单位对100名

5、女生测定血清总蛋白含量（g/L），数据如下:74.379.968.878.070.480.580.569.771.273.579.575.675.078.872.072.072.074.371.272.075.073.578.874.375.865.074.371.269.768.073.575.072.064.375.880.369.774.373.573.575.875.868.876.570.471.281.275.070.468.070.472.076.574.376.577.667.372.075.074.373.579.573.574.765.076.581.675.472.772

6、.767.276.572.770.477.268.867.367.367.372.775.873.575.073.573.573.572.781.670.374.373.579.570.476.572.777.284.375.076.570.4绘制上述数据的直方图、密度估计曲线、经验分布图和QQ图.3.2问题求解直方图;hist(serumdata,freq=FALSE,col=purpl e,border=red,de nsity=3,a ngle=60,ma in=pas te(the histogram of serumdata),xlab=age,ylab=freque ncy)3.2

7、.2运用lines函数绘制密度估计曲线;lin es(de nsity(serumdata),col=blue)3.2.3运用plot函数绘制经验分布图;x lin es(x,d no rm(x,mea n( serumdata),sd(serumdata),col=gree n) p lot(ecdf(serumdata),verticals=TRUE,do. p=FALSE) lin es(x ,pno rm(x,mea n( serumdata),sd(serumdata),col=blue)3.2.4运用qqnorm函数绘制QQ图qqno rm(serumdata,col=purpl

8、e) qqli ne(serumdata,col=red)3.3结果展示直方图密度估计曲线经验分布图QQ图4.1问题叙述甲、乙两种稻谷分别播种在10块试验田中，每块实验田甲乙稻谷各种一半假设两稻谷产量X，Y均服从正态分布，且方差相等.收获后10块试验田的产量如下所示（单位：千克）甲种140137136140145148140135144141乙种135118115140128131130115131125求出两稻种产量的期望差12的置信区间（ o.05）.7.3结果展示4.2问题求解 x x1 140 137 136 140 145 148 140 135 144 141 y y1 135 1

9、18 115 140 128 131 130 115 131 125t.test函数求解 t.test(x,y,var.equal=TRUE)4.3结果展示由以上程序运行得两稻种产量的期望差12的95%置信区间为7.53626,20.06374.5.1问题叙述甲乙两组生产同种导线，现从甲组生产的导线中随机抽取 4根，从乙组生产的导线中随机抽取5根，它们的电阻值（单位:）分别为甲组0.1430.1420.1430.137乙组0.1400.1420.1360.1380.1402 工22假设两组电阻值分别服从正态分布N（和N（ 2,），未知.试求12的置信区间系数为0.95的区间估计.5.2问题求解

10、 x y t.test(x,y,var.equal=TRUE)Two Samp le t-testdata: x and y t = 1.198, df = 7, p-value = 0.2699 alter native hypo thesis: true differe nee in mea ns is not equal to 0 95 p eree nt con fide nee in terval:-0.0019963510.006096351 samp le estimates: mea n of x mea n of y0.141250.13920由以上程序运行甲乙两电阻的期望差12的95%置信区间为-0.001996351，0.006096351.6.1问题叙述已知某种灯泡寿命服从正态分布,在某星期所生产的该灯泡中随机抽取10只，测得其寿命（单位：小时）为1067 919119678511269369181156920948求这个星期生产出的灯泡能使用1000小时以上的概率.6.2问题求解x数据框 x=c(1067,919,1196,785,1126,936,918,1156,920,948) xpnorm函数求解 pn orm(1000,m