假设我想从对数正态分布模拟10个观测值并重复100次。我写了一些R代码,但由于某种原因它没有用。这是代码:
for(i in 1:100)
{
x = rlnorm(10, meanlog = 0, sdlog = 1)
}
有什么想法吗?
答案 0 :(得分:5)
尝试此代码:
> 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)