模拟来自对数正态分布的数据R

Question

假设我想模拟来自对数正态分布的10个观测值并重复这100次。我写了一些R代码，但由于某种原因，它不起作用。这里是代码：

for(i in 1:100) 
{ 

x = rlnorm(10, meanlog = 0, sdlog = 1) 

} 

有什么想法？

Answer 1

这可能工作：

lapply(1:100, function(i) rlnorm(10, meanlog = 0, sdlog = 1)) 

编辑
计算平均值和SD使用：

lapply(1:100, function(i) { 
    x <- rlnorm(10, meanlog = 0, sdlog = 1) 
    c(mean=mean(x), sd=sd(x)) 
}) 

或者将其以矩阵形式返回（使用do.call）：

do.call(rbind, lapply(1:100, function(i) { 
    x <- rlnorm(10, meanlog = 0, sdlog = 1) 
    c(mean=mean(x), sd=sd(x)) 
})) 

而且还使原来的代码工作（见迪文注）使用方法：

x <- list() 
for(i in 1:100) { 
    x[[i]] <- rlnorm(10, meanlog = 0, sdlog = 1) 
}