A possible third problem is that 2:m+1 is the same as (2:m) + 1 and you probably want 2:(m+1)
Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Sarah Goslee > Sent: Monday, November 07, 2011 12:26 PM > To: Gyanendra Pokharel > Cc: R-help@r-project.org > Subject: Re: [R] Correction in error > > Hi, > > I see two problems right off: > > On Mon, Nov 7, 2011 at 3:10 PM, Gyanendra Pokharel > <gyanendra.pokha...@gmail.com> wrote: > > Hello R community, following is my code and it shows error, can some one > > fix this error and explain why this occurs? > > > > gibbs <-function(m,n, theta = 0, lambda = 1){ > > alpha <- 1.5 > > beta <- 1.5 > > gamma <- 1.5 > > x<- array(0,c(m+1, 3)) > > x[1,1] <- theta > > x[1,2] <- lambda > > x[1,3]<- n > > for(t in 2:m+1){ > > x[t,1] <- rbinom(x[t-1,3], 1, x[t-1,1]) > > The rbinom() command here returns a vector of values, but you're > trying to assign it to a single matrix element. You might want to > double-check the help for rbinom() again. > > > > x[t,2]<-rbeta(m, x[t-1,1] + alpha, tx[t-1,3] - x[t-1,1] + beta) > > tx isn't defined, but is probably a typo for x > > > x[t,3] <- rpois(x[t-1,3] - x[t-1,1],(1 - x[t-1,2])*gamma) > > } > > x > > } > > gibbs(100, 10) > > With a nice short function such as this, it's easy to set your > arguments to their default values and actually run the function line > by line at the command prompt so you can see whether what is happening > is what you expect. > > Sarah > > -- > Sarah Goslee > http://www.functionaldiversity.org > > ______________________________________________ > R-help@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. ______________________________________________ R-help@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.