R语言总和性试验.docxVIP

中北大学理学院实验报告实验课程名称：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,] 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%*%G3、函数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-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%*%y d=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=finding root is fail!)else{repeat{if(abs(b-a)eps)breakx-(a+b)/2if(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] 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))结果随机执行//产生对应直方图