R软件如何不使用循环而实现数组赋值

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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时数组太大导致运行时间太长。请问各位老师上述数组赋值不使用循环怎么实现。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：r软件 不使用 array 运行时间 For 如何 软件

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)  

n=100时候是19s，n=200时候我的R会自动退出，你看看能用不
results是向量的格式，具体怎么对应你自己来吧

万人往LVR 发表于 2015-7-19 12:24
n

你的想法很好，请问你的计算机是台式机，还是笔记本呢，因为我用你的方法感觉时间不止19秒呢。

fun<-function(n){
  A=array(0,dim=c(n,n,n))
  N=1:n
  dat=expand.grid(N,N,N)
  dat=as.matrix(dat)
  A[dat]=sin(rowSums(dat))+dat[,1]*dat[,2]*dat[,3]-dat[,1]
  A
}
system.time(fun(200))
##用户 系统 流逝
##3.64 0.40 4.68

tjxxukai 发表于 2015-7-19 13:27
你的想法很好，请问你的计算机是台式机，还是笔记本呢，因为我用你的方法感觉时间不止19秒呢。

笔记本，n=100时候是19s,n=200我的R直接崩溃了

万人往LVR 发表于 2015-7-20 10:00
笔记本，n=100时候是19s,n=200我的R直接崩溃了

谢谢，感谢

ryoeng 发表于 2015-7-19 13:29
如果数据太大，还可以使用 ff 包嚄！

多谢，thanks

十里春风 发表于 2015-7-19 22:34
fun

非常感谢，有用

1. f <- function(n) {
   
2.   n <- n - 1
   
3.   k <- floor(n / 200^2) + 1
   
4.   n <- n %% 200^2
   
5.   j <- floor(n / 200) + 1
   
6.   i <- n %% 200 + 1
   
7.   sin(i + j + k) + i*j*k - i
   
8. }
   
9. 
   
10. system.time(a <- array(f(1:200^3), c(200, 200, 200)))
    
11. 
    
12.    user  system elapsed
    
13.   0.561   0.026   0.588

复制代码