n=200
A=array(0,dim=c(n,n,n))
for(i in 1:n)
for(j in 1:n)
for(k in 1:n)
A[i,j,k]=sin(i+j+k)+i*j*k-i
A
这是可以在R中运行的代码，由于n=200时数组太大导致运行时间太长。请问各位老师上述数组赋值不使用循环怎么实现。
n <- 200
x <- rep(1:n,each=n^2);
y <- rep(rep(1:n,each=n),times=n);
z <- rep(1:n,times=n^2);
data <- cbind(x,y,z) ##建立(1,1,1),(1,1,2)...(1,200,200), ...(200,200,200)的数组
f <- function(x) {sin(x[1]+x[2]+x[3])+x[1]*x[2]*x[3]-x[1]}
results <- apply(data,1,f)
扫码加好友，拉您进群
全部版块
我的主页
收藏
好厉害！
