R语言实验报告-习题详解

R语言实验报告

习题详解

学院：

班级：

学号：

姓名：

导师：

成绩：

一、实验目的 (1)

二、实验内容 (1)

1.1问题叙述 (1)

1.2问题求解 (1)

1.2.1创建按列、行输入的4×5矩阵； (1)

1.2.2编写程序求解 (1)

1.3结果展示 (2)

2.1问题叙述 (2)

2.2问题求解 (3)

2.2.1创建StudentData数据框 (3)

2.2.2运行程序求解 (3)

2.3结果展示 (3)

3.1问题叙述 (3)

3.2问题求解 (4)

3.2.1运用hist函数绘制直方图； (4)

3.2.2运用lines函数绘制密度估计曲线； (4)

3.2.3运用plot函数绘制经验分布图； (4)

3.2.4运用qqnorm函数绘制QQ图 (4)

3.3结果展示 (4)

4.1问题叙述 (6)

4.2问题求解 (6)

4.2.1创建x、y数据框 (6)

4.2.2运用t.test函数求解 (6)

4.3结果展示 (7)

5.1问题叙述 (7)

5.2问题求解 (7)

5.2.1创建x、y数据框 (7)

5.2.2运用t.test函数求解 (7)

5.3结果展示 (7)

6.1问题叙述 (8)

6.2问题求解 (8)

6.2.1创建x数据框 (8)

6.2.2运用pnorm函数求解 (8)

6.3结果展示 (8)

三、实验总结 (8)

一、实验目的

R是用于统计分析、绘图的语言和操作环境。R是属于GNU系统的一个自由、免费、源代码开放的软件，它是一个用于统计计算和统计制图的优秀工具；

本次试验要求掌握了解R语言的各项功能和函数，能够通过完成试验内容对R语言有一定的了解，会运用软件对数据进行分析；

通过本实验加深对课本知识的理解以及熟练地运用R语言软件来解决一些复杂的问题。

二、实验内容

1.1问题叙述

将1,2，…，20构成两个4×5阶的矩阵，其中矩阵A是按列输入，矩阵B 是按行输入，并做如下运算.

C=A+B;

D=A*B;

F是由A的前3行和前3列构成的矩阵;

G是由矩阵B的各列构成的矩阵，但不含B的第3列.

1.2问题求解

1.2.1创建按列、行输入的4×5矩阵；

1.2.2编写程序求解

1.3结果展示

2.1问题叙述

已知有5名学生的数据，如下表所示.用数据框的形式读入数据.

序号姓名性别年龄身高（cm）体重（kg）1张三女1415642.0

2李四男1516549.0

3王五女1615741.5

4赵六男1416252.0

5丁一女1515945.5

2.2问题求解

2.2.1创建StudentData数据框

StudentData<-data.frame(name=c("zhangsan","lisi","wangwu","zhaoliu"," dingyi"),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"))

2.2.2运行程序求解

StudentData

2.3结果展示

3.1问题叙述

74.379.968.878.070.480.580.569.771.273.5 79.575.675.078.872.072.072.074.371.272.0 75.073.578.874.375.865.074.371.269.768.0 73.575.072.064.375.880.369.774.373.573.5 75.875.868.876.570.471.281.275.070.468.0 70.472.076.574.376.577.667.372.075.074.3 73.579.573.574.765.076.581.675.472.772.7 67.276.572.770.477.268.867.367.367.372.7 75.873.575.073.573.573.572.781.670.374.3 73.579.570.476.572.777.284.375.076.570.4

3.2问题求解

3.2.1运用hist 函数绘制直方图；

3.2.2运用lines

函数绘制密度估计曲线；

3.2.3运用plot 函数绘制经验分布图；

3.2.4运用

qqnorm 函数绘制QQ 图

3.3结果展示

直方图

the histogram of serumdata

age

f r e q u e

n c y

6570758085

0.000.020.040.060.080.1

密度估计曲线

the histogram of serumdata

age

f r e q u e n c y

6570758085

0.000.020.040.060.080.10

经验分布图

6570758085

0.00.20.4

0.60.81.0

ecdf(serumdata)

F n (x )

QQ 图

-2-1012

6570758085

Normal Q-Q Plot

Theoretical Quantiles

S a m p l e Q u a n t i l e

4.1问题叙述

甲、乙两种稻谷分别播种在10块试验田中，每块实验田甲乙稻谷各种一半.假设两稻谷产量X ，Y 均服从正态分布，且方差相等.收获后10块试验田的产量

求出两稻种产量的期望差

12μμ- 的置信区间（0.05α= ）.

4.2问题求解

4.2.1创建x 、y 数据框

4.2.2运用t.test 函数求解

4.3结果展示

由以上程序运行得两稻种产量的期望差12μμ-的95%置信区间为 [ 7.53626, 20.06374].

5.1问题叙述

甲乙两组生产同种导线，现从甲组生产的导线中随机抽取4根，从乙组生产的导线中随机抽取5根，它们的电阻值（单位：Ω ）分别为

甲组 0.143 0.142 0.143 0.137 乙组

0.140

0.142

0.136

0.138

0.140

假设两组电阻值分别服从正态分布

222

12(,)(,)N N μσμσσ和，未知.试求

μμ-的置信区间系数为0.95的区间估计.

5.2问题求解

5.2.1创建x 、y 数据框

> x<-c(0.143,0.142,0.143,0.137)

> y<-c(0.140,0.142,0.136,0.138,0.140) 5.2.2运用t.test 函数求解

> t.test(x,y,var.equal=TRUE)

5.3结果展示

Two Sample t-test

data: x and y

t = 1.198, df = 7, p-value = 0.2699

alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -0.001996351 0.006096351

sample estimates: mean of x mean of y 0.14125 0.13920

由以上程序运行甲乙两电阻的期望差12μμ-的95%置信区间为

[-0.001996351， 0.006096351].

6.1问题叙述

已知某种灯泡寿命服从正态分布，在某星期所生产的该灯泡中随机抽取10只，测得其寿命（单位：小时）为

1067 919 1196 785 1126 936 918 1156 920 948 求这个星期生产出的灯泡能使用1000小时以上的概率.

6.2问题求解

6.2.1创建x 数据框

6.2.2运用pnorm 函数求解

6.3结果展示

由以上程序运行得，x<=1000的概率为0.509,故x 大于1000的概率为0.491.

三、实验总结

在R 语言实验学习中，通过实验操作可使我们加深对理论知识的理解，学习和掌握R 语言的基本方法，并能进一步熟悉和掌握R 软件的操作方法，培养我们分析和解决实际问题的基本技能，提高我们的综合素质. 如有侵权请联系告知删除，感谢你们的配合！