中北大学理学院
实验报告
实验课程名称： 实验类别： 专班学姓
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、
由1.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,] -1.000000e+00 [2,] 1.000000e+00 [3,] 3.806634e-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=\else{repeat{if(abs(b-a) 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] 1.324718 $fun [1] -1.405875e-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] 1.539077e-06 抽到的是同花顺 （2）9*choose(4,1)/choose(52,5) [1] 1.385169e-05 3、 从正态分布N（100,100）中随机产生1000个随机数， （1）rnorm(1000,mean=100,sd=10) 结果随机执行//产生1000个随机数 作出这1000个正态随机数的直方图； hist( rnorm(1000,mean=100,sd=10)) 结果随机执行//产生对应直方图