R语言中 两个数组元素怎样堆叠

R语言中 两个数组元素怎样堆叠

x<-array(1:8,c(2,2,2))

x1<-array(2:9,c(2,2,2))

xx=array(,c(nrow(x)+nrow(x1),ncol(x),dim(x)[3]))

for(i in 1:dim(x)[3]) xx[,,i]=rbind(x[,,i],x1[,,i])

所谓数组(或向量)a和b的外积,指的是a的每一个元素和b的每一个元素搭配在一起相乘得到的新元素.当然运算规则也可自定义.外积运算符为 %o%(注意:百分号中间的字母是小写的字母o).例如: [plain] view plain copy > a <- 1:2 > b <- 3:5 > d <- a %o% b > d [,1] [,2] [,3] [1,] 3 4 5 [2,] 6 8 10 注意维数公式为: [plain] view plain copy dim(d) = c( dim(a) , dim(b) ) 实际上r语言提供了一个更为一般化得外积函数outer(x,y,fun),比较简单,看个示例吧. 下面的例子是为了绘制 z = x**2 +50cos(y)的3维表面图: [plain] view plain copy > x <- seq(from = -5, to = 5, by = 0.25) > y <- seq(from = -2, to = 2, by = 0.1) > f <- function(x,y){x**2+50*cos(y)} > z <- outer(x,y,f) > persp(x,y,z)