R软件入门幻灯片.pptVIP

? x=1:3 x [1] 1 2 3 ? crossprod(x) [,1] [1,] 14 ? cbind(x) %*% x [,1] [,2] [,3] [1,] 1 2 3 [2,] 2 4 6 [3,] 3 6 9 ? x %*% rbind(x) [,1] [,2] [,3] [1,] 1 2 3 [2,] 2 4 6 [3,] 3 6 9 ? y=diag(x) y [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 2 0 [3,] 0 0 3 ? diag(y) [1] 1 2 3 ? diag(3) [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 1 0 [3,] 0 0 1 31 2020/3/20 矩阵合并与拉直 ? 函数 cbind() 把其自变量横向拼成一个大矩阵， rbind() 把其自变量纵向拼成一个大矩阵。 ? cbind() 的自变量是矩阵或者看作列向量的向 量，自变量的高度应该相等（对于向量，高 度即长度，对于矩阵，高度即行数）。 ? rbind 的自变量是矩阵或看作行向量的向量， 自变量的宽度应该相等（对于向量，宽度即 长度，对于矩阵，宽度即列数）。 ? 如果参与合并的自变量比其它自变量短则循 环补足后合并。例如： 32 2020/3/20 ? x1= rbind(c(1,2), c(3,4)) ? x1 [,1] [,2] [1,] 1 2 [2,] 3 4 ? x2 - 10+x1 ? x3 - cbind(x1, x2) ? x3 [,1] [,2] [,3] [,4] [1,] 1 2 11 12 [2,] 3 4 13 14 ? x4 = rbind(x1, x2) ? x4 [,1] [,2] [1,] 1 2 [2,] 3 4 [3,] 11 12 [4,] 13 14 ? cbind(1, x1) [,1] [,2] [,3] [1,] 1 1 2 [2,] 1 3 4 33 2020/3/20 ? a= rbind(c(1,2), c(3,4)) ? a [,1] [,2] [1,] 1 2 [2,] 3 4 ? as.vector(a) [1] 1 3 2 4 ? b=t(a) ? b [,1] [,2] [1,] 1 3 [2,] 2 4 ? c(a,b) [1] 1 3 2 4 1 2 3 4 34 2020/3/20 数组的外积 ? 两个数组 a 和 b 的外积是由 a 的每一个元素与 b 的每一个元 素搭配在一起相乘得到一个新元素，这样得到一个维数 向量等于 a 的维数向量与 b 的维数向量连起来的数组， ? 即若 d 为 a 和 b 的外积，则 dim(d)=c(dim(a), dim(b)) 。 ? a 和 b 的外积记作 a %o% b 。如 d - a %o% b 也可以写 成一个函数调用的形式： d - outer(a, b, *) 注意 outer(a, b, f) 是一个一般性的外积函数，它可以把 a 的任一 个元素与 b 的任意一个元素搭配起来作为 f 的自变量计算 得到新的元素值，外积是两个元素相乘的情况。函数当 然也可以是加、减、除，或其它一般函数。当函数为乘 积时可以省略不写。 35 2020/3/20 数组的外积 ? 例如，我们希望计算函数 在一个 x 和 y 的网格上的值用 来绘制三维曲面图，可以用如下方法生成网格及函数值： ? x - seq(-2, 2, length=20) ? y - seq(-pi, pi, length=20) ? f - function(x, y) cos(y)/(1+x^2) ? z - outer(x, y, f) ? persp(x,y,z) ? 用这个一般函数可以很容易地把两个数组的所有元素都 两两组合一遍进行指定的运算。 36 2020/3/20 apply 函数 ? 对于向量，我们有 sum 、 mean 等函数对其进行计算。对于 数组，如果我们想对其中一维（或若干维）进行某种计算， 可以用 apply 函数。其一般形式为： ?