R软件及统计分析向量多维数组和矩阵.ppt

a=cbind(c(4,9,1),c(3,7,2)) a [,1] [,2] [1,] 4 3 [2,] 9 7 [3,] 1 2 apply(a,2,sum) [1] 14 12 apply(a,1,sum) [1] 7 16 3 apply(a,2,sort) [,1] [,2] [1,] 1 2 [2,] 4 3 [3,] 9 7 apply(a,1,sort ) [,1] [,2] [,3] [1,] 3 7 1 [2,] 4 9 2 t(apply(a,1,sort)) [,1] [,2] [1,] 3 4 [2,] 7 9 [3,] 1 2 apply(a, 2, sort) 但是，如果要对行排序，则apply(a, 1, sort)把a的每一行3个元素排序后的结果用第一维来引用，结果的维向量为c(N, dim(a)[1])=c(3, 4)，把原来的列变成了行，所以t(apply(a,1,sort)) 才是对a的每一行排序的结果。如： apply(a, 1, sort) t(apply(a,1,sort)) 上面我们只用了矩阵（二维数组）作为例子讲解apply的用法。实际上，apply可以用于任意维数的数组，函数FUN也可以是任意可以接收一个向量或数组作为第一自变量的函数。比如，设x是一个维数向量为c(2,3,4,5)的数组，则apply(x, c(1,3), sum)可以产生一个2行4 列的矩阵，其每一元素是x中固定第1维和第3维下标取出子数组求和的结果。 实验作业 1、六个元素 a, b, c, d , e, f都是相互独立的服从离散均匀分布的随机变量，其元素均在1,...,9中取值。设法求ade-bcf 的分布。 要求 (1)产生一个外积矩阵; (2)用一个table()函数来估计每一个值的出现概率; (3)绘制一个频数分布图，一个频率分布图。 要求：需给出程序、结果，存成word文档 发送到 3 用户名：rsoft 密码： 123456 2. (1) 写出元素为3, -1.5, 3E-10的向量。 (2) 写出从3开始每次增加3，长度为100的向量。 (3) 写出(0, 2)重复10次的向量。 (4) 对向量x，写出其元素大于等于0小于1的条件。 (5) 对向量x，写出其元素都等于0的条件。 (6) 写出包含12个月份名称的向量。 (7) 写出包含方程 的根的向量，并写出其幅角的余弦和正弦值。 3. 设x为一个长100的整数向量。比如，x - floor(100*runif(100))。 (1) 显示x第21到30号元素。 (2) 把x第31，35，39号元素赋值为0。 (3) 显示x中除了第1号和第50号的元素之外的子集。 (4) 列出x中个位数等于3的元素。 (5) 列出x中个位数等于3的元素的下标位置。Hint：use which( ) (6) 给x的每一个元素加上名字，为x1到x100。 (7) 求x的平均值并求每一个元素减去平均值后的离差，计算x元素的平方和及离差平方和。 (8) 把x从大到小排序。计算x的10%分位数到90%分位数之间的距离。 4. 定义一个维数为(3, 4, 2)的数组其第一层（第三 下标为1）取从1开始的奇数，第二层取从2开始 的偶数。显示每一层的第2行元素。把第(1,1,1)， (2,2,2), (2, 2, 1)号元素赋值为零。把第一层加上 100，把第二层加上200。分别计算第一层和第二 层的平均值。 5. 试画出任一给定的二元正态分布的密度函数图形。如二元标准正态分布等 三、取逻辑值的下标向量 v为和x等长的逻辑向量，x[v]表示取出所有v为真值的元素，如： x 输出： [1] 1.00 125.00 6.25 x10 输出： [1] TRUE FALSE TRUE x[x10] 输出： [1] 1.00 6.25 x[x0] 输出： numeric(0) 可见x[x10]取出所有小于10的元素组成的子集。 这种逻辑值下标是一种强有力的检索工具，例如x[sin(x)0]可以取出x中所有正弦函数值为正的元素组成的向量。 如果下标都是假值则结果是一个零长度的向量， 显示为numeric(0)