R语言总和性试验详解.docx

中北大学理学院 实验报告 实验课程名称：R语言与统计分析实验类别：验证型专业：应用统计学班级号：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+B2、D=A*B3、E=A%*%B4、F-A[-3,][,1][,2][,3][,4][1,]15913[2,]261014[3,]481216G-B[,-3][,1][,2][,3][1,]124[2,]568[3,]91012[4,]131416H=F%*%G3、函数solve()有两个作用；solve（A,b）可用于求解线性方程组Ax=b,solve(A)可用于求解矩阵A的逆，用两种方法编程求解方程组Ax=b的解。A-matrix(1:9,3,3,byrow=TRUE)A[3,3]=10b=matrix(1:1,3,1)solve(A,b)[,1][1,]-1.000000e+00[2,]1.000000e+00[3,]3.806634e-164、用三种方法求解它们的內积与外积。1、x=c(1,2,3,4,5)y=c(2,4,6,8,10)a=t(x)%*%ye=x%*%t(y)2、b=crossprod(x,y)f=outer(x,y)3、c=x%*%yd=x%o%y5、编写一个用二分法求解非线性方程的根的函数，并求方程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=findingrootisfail!)else{repeat{if(abs(b-a)eps)breakx-(a+b)/2if(f(a)*f(x)0)b-xelsea-x}list(root=(a+b)/2,fun=f(x))}}f-function(x)x^3-x-1fzero(f,1,2,1e-5)$root[1]1.324718$fun[1]-1.405875e-05第三章：1、从1到100个自然数中随机不放回的抽取5个数，并求它们的和。sum(sample(1:100,5))[1]2052、从一副扑克牌（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-053、从正态分布N（100,100）中随机产生1000个随机数，（1）rnorm(1000,mean=100,sd=10)结果随机执行//产生1000个随机数作出这1000个正态随机数的直方图；hist(rnorm(1000,mean=100,sd=10))结果随机执行//产生对应直方图从这1000个随机数中随机有放回的抽取500个作出直方图（2）A-rnorm(1000,mean=100,sd=10)sample(A,500,replace=TRUE)或者sample(rnorm(1000,mean=100,sd=10),500)结果随机执行//产生500个对应随机数hist(sample(A,500,replace=TRUE))结果随机执行//产生对应直方图比较它们的样本均值和样本方差（3）mean(rnorm(1000,mean=100,sd=10))[1]100.0266var(rnorm(1000,mean=100,sd=10))[1]98.5499mean(sample(rnorm(1000,mean=100,sd=10),500))[1]99.9142var(sample(rnorm(1000,mean=100,sd=10),500))[1