实验报告二答案.docVIP

实验报告二 课 程 统计软件分析与应用 学生姓名 学 号 学 院 数学与统计学院 专 业 统计学 指导教师 二Ｏ一四 年三月十日 1、向量 (1) 写出元素为3, -1.5, 3E-10的向量赋值给a。 a-c(3,-1,3E-10) (2) 写出从3开始每次增加3，长度为100的向量。 b-seq(length=100,from=3,by=3) (3) 写出(0, 2)重复10次的向量赋值给x。 x-rep(c(0,2),10) (4) 对向量x，选出其元素大于等于0小于1的x的值赋值给y。 y-x[x=0=1] (5) 对向量x，选出其元素都等于0的值并赋值给z。 z-x[x==0] (6) 设u为一个长100的整数向量。比如，u - floor(10*runif(100))。 显示u第21到30号元素。 u - floor(10*runif(100)) u[21:30] (7) 把u第31，35，39号元素赋值为0。 u - floor(10*runif(100)) u[c(31,35,39)]-0 (8) 显示u中除了第1号和第50号的元素之外的子集。 u-u[-c(1,50)] (9) 列出u中个位数等于3的元素。 u[u%%10==3] (10) 列出u中个位数等于3的元素的下标位置。 which(u%%10==3) (11) 给u的每一个元素加上名字，为u1到u100。 paste(u,1:100,sep= ) 2、建立一个4*4的矩阵A并求出其特征值、特征向量，行列式，QR分解，并对矩阵中每行求均值，标准差，以及标准化。 A-matrix(1:16,nrow=4,ncol=4) eigen(A) det(A) [1] 0 Aplus-qr(A) for(i in 1:4 ){ + y[i]-mean(A[i,]) + z[i]-sd(A[i,]) + x-scale(A[i,]) + } 或者 apply(A,1,mean) apply(A,1,sd) apply(A,1,scale) 3、建立一个数据框有15个样本，5个变量，变量分别是name(字符型)，身高(数值型)，体重(数值型)，BMI(体重指数)，等级(有序变量，因子)，其中有三个等级，前5个为1，中间为2，后五个等级是3.（除了最后的等级变量其他数值自己随便定义） df-data.frame(0) edit(df) name height weight BMI level 1 a 160 100 0.2 1 2 b 161 98 0.3 1 3 c 162 102 0.2 1 4 d 162 100 0.5 1 5 e 163 105 0.6 1 6 f 164 100 0.6 2 7 g 165 104 0.3 2 8 h 167 110 0.8 2 9 i 170 115 0.9 2 10 l 165 102 0.2 2 11 j 172 120 0.3 3 12 k 168 105 0.4 3 13 m 180 130 0.6 3 14 n 169 109 0.8 3 15 o 166 110 0.4 3 4、试画出参数为(0,1,0,1,0)的二元正态分布密度函数三维图。 x-y-seq(-5,5,0.1) f-function(x,y) {(1/2*pi)*exp((-x^2-y^2)/2)} 或者f-function(x,y){dnorm(x,0,1)*dnorm(y,0,1)} z-outer(x,y,f) persp(x,y,z,theta=30) 5、把user.txt数据中的性别、年龄、身高分别输入到R中。计算不同性别、不同年龄的人数，并计算每一组的平均身高。把这些变量组合成一个列表。把user.txt数据输入为R的数据框。 a-read.table(E:/R/us