文档介绍：中北大学理学院
实验报告
实验课程名称:
R语言与统计分析
实验类别:
验证型
专业
应用统计学
:
班级:
13080441
学号:
1308044142
姓名:
吴庚雷
中北大学理学院
R语言与统计分析综合实验
【实验类型】验证性
【实验目的】
(1)掌握利用R语言实现数据处理并进行严格的统计分析;
(2)学会运用R语言进行程序的编写;
(3)熟练掌握R语言绘图功能;
(4)掌握R语言统计分析中的“参数估计”,“假设检验”,“方差分析”,“回归分析”, 等基本分析函数。
【实验要求】
(1)实验过程要求用R软件完成;
(2)实验结果逐个导入Word文档,并按问题作出解释;
(3)实验报告按照既定格式书写。
【实验仪器与软件平台】
计算机 R软件
【实验前的预备知识】
1、实验室电脑要求安装有R软件;
2、上实验课程的学生要对涉及到的统计概念有所了解;
3、要求学生事先查阅并熟悉R的相关命令。
【实验内容】
第二章:
1、
用rep()构造一个向量x,它由3个3,4个2,5个1构成。
x<-rep(c(3,2,1),c(3,4,5))
2、
...16构成两个方阵,其中矩阵A按列输入,矩阵B按行输入,并计算以下:
A<-matrix(1:16,4,4)
B<-matrix(1:16,4,4,byrow=TRUE)
1、
C=A+B
2、
> D=A*B
3、
> E=A%*%B
4、
F<-A[-3,]
[,1] [,2] [,3] [,4]
[1,] 1 5 9 13
[2,] 2 6 10 14
[3,] 4 8 12 16
> G<-B[,-3]
[,1] [,2] [,3]
[1,] 1 2 4
[2,] 5 6 8
[3,] 9 10 12
[4,] 13 14 16
> H=F%*%G
3、
函数solve()有两个作用;solve(A,b)可用于求解线性方程组Ax=b,solve(A)可用于求解矩阵A的逆,用两种方法编程求解方程组Ax=b的解。
A<-matrix(1:9,3,3,byrow=TRUE)
> A[3,3]=10
> b=matrix(1:1,3,1)
> solve(A,b)
[,1]
[1,] -+00
[2,] +00
[3,] -16
4、
用三种方法求解它们的內积与外积。
1、
> x=c(1,2,3,4,5)
> y=c(2,4,6,8,10)
> a=t(x)%*%y
e=x%*%t(y)
2、
> b=crossprod(x,y)
> > f=outer(x,y)
3、
> c=x%*%y
> d=x%o%y
5、
编写一个用二分法求解非线性方程的根的函数,并求方程x^3-x-1=0在区间[1,2]内的根,精度要求e=10^-5.
程序:
fzero<-function(f,a,b,eps=1e-5)
{
if(f(a)*f(b)>0)list(fail="finding root is fail!")
else{repeat{if(abs(b-a)<eps)
break
x<-(a+b)/2
if(f(a)*f(x)<0)
b<-x
else a<-x}
list(root=(a+b)/2,fun=f(x))}}
> f<-function(x) x^3-x-1
> fzero(f,1,2,1e-5)
$root
[1]
$fun
[1] --05
第三章:
1、
从1到100个自然数中随机不放回的抽取5个数,并求它们的和。
> sum(sample(1:100,5))
[1] 205
2、
从一副***牌(52张)中随机抽取5张求以下概率
抽到的是10,J,Q,K,A;
(1)> 4/choose(52,5)
[1] -06
抽到的是同花顺
(2)9*choose(4,1)/choose(52,5)
[1] -05
3、
从正态分布N(100,100)中随机产生1000个随机数,
(1)rnorm(1000,mean=100,sd=10)
结果随机执行//产生1000个随机数
作出这1000个正态随机数的直方图;
hist( rnorm(1000,mean=100,sd=10))
结果随机执行//产生对应直方图
从这1000个随机数中随机有放回的抽取500个