r - R中的数据模拟

Question

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

for(i in 1:100) 
 {

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

 }

有什么想法吗？

Answer 1

也试试这个代码：

> x=matrix(0,nrow=10,ncol=100)
> for(i in 1:100) 
+  {
+ 
+ x[,i] = rlnorm(10, meanlog = 0, sdlog = 1)
+ 
+  }
> 
> apply(x,2,mean)
> apply(x,2,sd)
> library(moments)
> apply(x,2,skewness)