统计软件课堂练习.docVIP

2017年10月11日 第一次作业 1、四个元素a,b,c,d,e,f都是相互独立的服从离散均匀分布的随机变量，其元素均在1，...，9中取值。设法求ade-bcf的分布。 要求：（1）产生一个外积矩阵； （2）用一个table（）函数来估计每一个值的出现概率； （3）绘制一个频数分布图，一个频率分布图。 1、源程序： x- outer(1:9,1:9) y- outer(1:9,x) z- outer(y, y,-) z f- table(z) f plot(as.numeric(names(f)), f, type=h, xlab =行列式值, ylab =频数) f- f / 9^6 plot(as.numeric(names(f)), f, type=h, xlab =行列式值, ylab =频率) 2、（1）写出元素为3，-1.5,3E-10的向量。 源程序： x-c (3, -1.5, 3*exp(-10)) x 输出结果： （2）写出从3开始每次增加3，长度为100的向量。 源程序： seq(from=3, by=3, length=100) 输出结果： （3）写出(0, 2)重复10次的向量。 源程序： rep(c(0,2),10) 输出结果： （4）对向量x，写出其元素大于等于0小于1的条件。 源程序： x[(x=0)(x1)] （5）对向量x，写出其元素都等于0的条件。 源程序： x[ ]=0 （6）写出包含12个月份名称的向量。 源程序： Q-paste(1:12,月份,sep=) Q 3、设x为一个长100的整数向量。比如，x - floor(100*runif(100))。 （1）显示x第21到30号元素。 源程序： x - floor(100*runif(100)) x[21:30] 输出结果： （2）把x第31，35，39号元素赋值为0。 源程序： x[31]-x[35]-x[39]-0 x 输出结果： （3）显示x中除了第1号和第50号的元素之外的子集。 源程序： x[-c(1,50)] 输出结果： （4）列出x中个位数等于3的元素。 源程序： x[x%%10 == 3] 输出结果： （5）列出x中个位数等于3的元素的下标位置。 源程序： which(x%%10 == 3) 输出结果： （6）给x的每一个元素加上名字，为x1到x100。 源程序： names(x)=paste(x,1:100,sep=) x 输出结果： （7）求x的平均值并求每一个元素减去平均值后的离差，计算x元素的平方和及离差平方和。 源程序： mean(x) x-mean(x) sum(x^2) sum((x-mean(x))^2) 输出结果： （8）把x从大到小排序。计算x的10%分位数到90%分位数之间的距离。 源程序： y-sort(x, decreasing = T) y y[10]-y[90] 输出结果： 4、定义一个维数为(3, 4, 2)的数组，其第一层（第三下标为1）取从1开始的奇数，第二层取从2开始的偶数。显示每一层的第2行元素。把第(1,1,1)，(2,2,2), (2, 2, 1)号元素赋值为零。把第一层加上100，把第二层加上200。分别计算第一层和第二层的平均值。 4、源程序： x- c((1:12)*2-1,(1:12)*2) A- array(x,dim=c(3,4,2)) A t(a[2,,]) A[1,1,1]-A[2,2,2]-A[2,2,1]-0 A A[, , 1] - A[, , 1] + 100 A[, , 2] - A[, , 2] + 200 A mean(A[,,1]) mean(A[,,2]) 输出结果： 2017年10月18日 第二次作业 1、查询tapply的函数帮助信息，并用帮助文件中的案例进一步学习（自己给一个新的例子）。 1、源程序： n - 17 fac - factor(rep(1:3, length = n), levels = 1:5) fac table(fac) tapply(1:n, fac, sum) 输出结果： （注：51=1+4+7+10+13+16，是对数据为1的下标进行求和，57，45同理。 2、把user.txt数据中的性别、年龄、身高分别输入到R中。计算不同性别、不同年龄的人数，并计算每一组的平均身高。把这些变量组合成一个列表。 2、源程序： data- read.table(file=user.txt, header=T) names(data) table(data$Sex) table(data$Age) tapply(data$Height,data$Sex,mean) tapply(data$Height,data$Age