Question

我使用for loop循环遍历两个不同的矩阵。代码如下所示：

x <- matrix(rnorm(1806),7,258)
x2 <- matrix(rnorm(1032),4,258)

samp_size <- 3
iter <- 1000
subs <- matrix(sample(1:nrow(x), samp_size*iter, replace=T),
               ncol=samp_size, byrow=T)
subs2 <- matrix(sample(1:nrow(x2), samp_size*iter, replace=T),
                ncol=samp_size, byrow=T)

for(j in 1:nrow(subs)){
  ad <- x[subs[j,],]
  ad1 <- x2[subs2[j,],]
  rd <- rbind(ad,ad1)
  dis <- dist(rd, method="euclidian")
  #CV <- sd(unlist(rd), na.rm=TRUE)/mean(unlist(rd), na.rm=TRUE)*100
  dis2 <- dis[dis!=0]
  list[[j]] <- mean(dis2)
  l <- unlist(list)
  l <- na.omit(l)
  sdis <- dist(ad, method="euclidian")
  sdis2 <- sdis[sdis!=0]
  res <- vector("list", 1000)
  res[[j]] <- mean(sdis2)
  l2 <- unlist(res)
  l2 <- na.omit(l2)
}

list[[j]]效果很好。

每次迭代的结果都会根据j的值作为列表元素添加到其中。

但是，res[[j]]最终只有NULL，只有最后一个元素包含值。

我不确定这里发生了什么。我希望res[[j]]的行为方式与list[[j]]相同。

非常感谢。

Answer 1

也许初始化res：

# before the loop
res <- rep(list(NA),nrow(subs))
# or in the loop
res <- rep(list(NA),1000) 
# or 
res <- list(numeric(1000))

我想在for循环代码中使用一个列表。 list [[i]]工作正常，但为什么不res＆lt; - vector（＆＃34; list＆＃34;，n）工作？

1 个答案: