[R] generate

[David Winsemius]

On Mar 27, 2014, at 4:31 AM, IZHAK shabsogh wrote:

> i try to generate 27*5  matrix of observation using the following code but is given error
> kindly assist and correct the problem 
> 
> 
> x1<-c(-1,0,1,-1,0,1,-1,0,1,-1,0,1,-1,0,1,-1,0,1,-1,0,1,-1,0,1,-1,0,1)
> x2<-c(-1,-1,-1,0,0,0,1,1,1,-1,-1,-1,0,0,0,1,1,1,-1,-1,-1,0,0,0,1,1,1)
> x3<-c(-1,-1,-1,-1,-1,-1,-1,-1,-1,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1)
> mu<-50+5*(x1^2+x2^2+x3^2)
> sig<-100+5*((x1-0.5)^2+x2^2+x3^2)
> x.matrix<-matrix(0,nrow=27,ncol=5)
> for(i in 1:27){

Perhaps, instead:

 x.matrix[i, ]<-rnorm(5,mu[i],sig[i])  # since 'j' doesn't exist

> x.matrix[i,j]<-rnorm(5,mu[i],sig[i])
> }
> 
> 
> Thanks
> Ishaq
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius
Alameda, CA, USA