r<-function(n, mu,sigma,alpha,w)
+ {m<-length(w)
+  w<-cumsum(w)
+  U<-runif(n)
+  W<-matrix(0,n,m)
+  W[,1]<-U<w[1]
+  for (i in 2:m)
+  {W[,i]<-(U<w)&(U>=w[i-1])
+  }
+  z<-matrix(0,n,m)
+  for(i in 1:m){
+  z[,i]<-z[,i]+W[,i]*(alpha+exp(rnorm(n,mu,sigma)))
+  }
+  }
>  n<-100
>  mu=c(0.5,1.1,2.6)
>  sigma=c(0.2,0.3,0.2)
>  alpha=c(0.4,1.2,2.4)
>  w=c(0.33,0.33,0.34)
>  r(n,mu,sigma,alpha,w)
> summary(r)
错误于object[] : 类别为'closure'的对象不可以取子集

是为了生成混合的lognormal 分布的随机数。我怎么都不知道错在哪里。而且也不知道这个错误反馈是什么意思。有人知道的话麻烦帮我看看吧，感激不尽！！！