similar to: WinBUGS14 and R

Hi, I am trying to covert a Winbugs code into R code. Here is the winbugs code model{# model’s likelihoodfor (i in 1:n){time[i] ~ dnorm( mu[i], tau ) # stochastic componenent# link and linear predictormu[i] <- beta0 + beta1 * cases[i] + beta2 * distance[i]}# prior distributionstau ~ dgamma( 0.01, 0.01 )beta0 ~ dnorm( 0.0, 1.0E-4)beta1 ~ dnorm( 0.0, 1.0E-4)beta2 ~ dnorm( 0.0, 1.0E-4)#

Hello R and BUGS users, I am writing a heirarchical model in R to send to BUGS via R2WinBUGS and I am finding it difficult to get the model to run. I seem to be having two problems. 1) I can't seem to send variables classed as factors (Month), is there a way do this? 2) Checking the Log in WinBUGS I can see that the model is Syntactically correct, but Bugs is not able to recognise the the

Hi, I am trying to do run the following model saved in "C:/bugs/sus.bug" model { for (i in 1:n){ y[i] ~ dpois(lamdba[i]) log(lambda[i]) <- mu+bmale[male[i]]+bschn[schn[i]]+epsilon[i] # epsilon[i] ~ dnorm(0,tau.epsilon) } mu ~ dnorm(0,.0001) bmale ~ dnorm(0,.0001) tau.epsilon <- pow(sigma.epsilon, -2) sigma.epsilon ~ dunif(0,100) for (j in

I'm having difficulty running R2WinBugs on a setup that previously worked for me (Dell Laptop, Windows XP service pack 3, R 2.11.1, WinBugs 1.43) . When I issue the following command smds.sim <- bugs (data, inits, parameters, "SMDSbrandLoc2.bug", debug = T, n.thin = n.thin, n.chains = n.chains, n.burnin = n.burnin, n.iter = n.iter) I get the following error in R

Hi, I am trying to use WinBUGS using R, but i am having some problems using bugs function: out <- bugs(data = win.data, inits = inits, parameters.to.save = params, model.file = "model.txt", n.thin = nt, n.chains = nc, n.burnin = nb, n.iter = ni, debug = TRUE, DIC = TRUE, working.directory = getwd()) ls() WinBUGS file is in working directory=D:/Line/documents but R looks for

I'm writing a paper aimed at motivating the use of GEE within the field of economics. However, after computing using the geeglm function, I noticed there's one intercept in the summary output. I assume this means the function is pooling the data. That means my code is not what I want. I want a "fixed effects" model, meaning I want the intercept to vary by cluster. Here's

Hi all, I'm trying to do some object-oriented programming in R using the R.oo package. Right from the start I encountered a strange (at least for me) problem. I define a new class/constructor based on the R.oo documentation. However the first attempt to create an object fails: === Code: === library(R.oo); setConstructorS3("MyClass",function(param) { print(param);

Dear All 1) I'm fitting spatial CAR models using R2Winbugs and although everything seems to go reasonably well (or I think so) the next message appears from WINBUGS 1.4 window: gen.inits() Command #Bugs: gen.inits cannot be executed (is greyed out) The question is if this message means that something is wrong and the results are consequently wrong, or Can I assume it as a simple warning

R2WinBUGS users, I'm getting the error message "monitor could not be set" when I try to monitor (and output to R) a response variable that includes missing data (NAs)...ie.e imputing for y. Setting monitors for other parameters works OK and running the R-generated WinBUGS script file in WinBUGS, also works OK. Below is my R code and the log.txt file output by R2WinBUGS. As I want to

Hi All Does anyone know if I can avoid to use the write.model() function below? I dont want to do this. Can't bugs() do that automatically for me just by specifying the 4th argument 'model'? Just I like I am also using the 'inits' object! If I use 'model' in the same way as I use 'inits' I am getting the error: > sim <- bugs(data, inits, parameters,

Dear All, I wish to calculate the Gini index (ineq from same package) and some other indices for the diameter distribution of each plot (df dgtot). dgtot: IDPlot Diameter(cm) 1 4 34.0 2 4 23.0 3 4 38.0 ... 51 5 16.0 52 5 8.0 53 5 9.0 ... 5301 140 25.0 5302 140 12.0 5303 140 7.0 I use: >

Hi listers, Does it necessary to install any package in order to use the GINI or INEQ functions. If I use the following command the R tells me that didn't find the GINI function. x<-c(541, 1463, 2445, 3438, 4437, 5401, 6392, 8304, 11904, 22261) G<-gini(x) Thanks in advance, Marcio -- View this message in context:

Dear Users, I am new to both of things, so do not blame me too much... I am busy with semiparametric regression and use WinBUGS to sample posteriors. The code to call Winbugs is as follows: data <- list("y","X","n","m") #My variables inits.beta <- rep(0,K) inits.beta0 <- 0 inits <-

I came across the package 'ineq' that computes a variety of inequality measures (e.g. gini, theil etc). I want to compute the Theil index (racial segregation) and decompose the total into sub-components (by geog levels). I think the package doesn't report the decomposition (correct me if I'm wrong). Just wonder is that available elsewhere? K.

Appologies for cross-posting Dear R users, I am using R2WinBUGS to call WinBUGS from R. After loading data, model, and initial values I call this command res <- bugs(data = dfile, inits = list(ifile), parameters.to.save = c("beta"), model.file = mfile, working.directory = tdir, n.thin = nthin, n.chains = 1, n.iter = niterations*nthin, n.burnin = 0, DIC = F, debug = T,

Hi All I can specify whatever inits, it has no effect on the estimation. I am replicating a textbook example. The result is completely trash, having estimates of -58.7 (sd=59.3), where it should be closer to an ml estimate of 0.585 (SE=0.063). The two chains within one run are different, but with different inits for different runs, I get exactly the same chains, and I mean exactly. If I set

If I use the Ineq library and the Gini function in this way: >Gini(c(100,0,0,0)) I obtain the result 0.75 instead of 1 (that is the perfect inequality). I think Gini's formula in Ineq is based on a formula as reported here: http://mathworld.wolfram.com/GiniCoefficient.html but in the case of perfect inequality: x_1=.......=x_n-1 =0 x_n>0 these formula are equal to 1 - 1/n, not to

I'm back! I've just learned that, on a fully updated Fedora Core Linux5 sytem, the working solution to access Winbugs under wine via the R package "rbugs" no longer works. Here was my last post on this topic (with the formerly working solution) from January. http://finzi.psych.upenn.edu/R/Rhelp02a/archive/68497.html Currently, what happens is that WinBUGS starts up, but just

*Dear R-user, * I try to run Winbugs from R using bugs function in R2WinBUGS.My model works well in Winbugs except that I can't get DIC. Since I don't need DIC, when I try to run Winbugs from R , I set "DIC=FALSE". My model is as following: model { for (i in 1:N) { for(j in 1 : T ) { x[i, j] ~ dbin(p[i, j],n[i]) #Hier.prior p[i, j] ~ dbeta(alpha[i, j], beta[i, j]) alpha[i, j]

Hi, I??m trying to use the Win2BUGS package from R and I have a similar problem that reurns with the message: Error in FUN(X[[1]], ...) : symbol print-name too long But, there is no stray ` character in the file ( Sugestions given by: Duncan Temple Lang <duncan> Date: Mon, 26 Sep 2005 07:31:08 -0700 ) The progam in R is: library(R2WinBUGS) library(rbugs) dat <-