search for: gelmans

Hi, A question. When I run gelman.diag and gelman.plot with mcmc lists obtained from MCMCregress, the results are following. > post.R <- MCMCregress(Size~Age+Status, data = data, burnin = 5000, mcmc = 100000, + thin = 10, verbose = FALSE, beta.start = NA, sigma2.start = NA, + b0 = 0, B0 = 0, nu = 0.001, delta = 0.001) > post1.R <- MCMCregress(Size~Age+Status, data

I would like to fit a hierarchical regression model from Witte et al. (1994; see reference below). It's a logistic regression of a health outcome on quntities of food intake; the linear predictor has the form, X*beta + W*gamma, where X is a matrix of consumption of 82 foods (i.e., the rows of X represent people in the study, the columns represent different foods, and X_ij is the amount of

Hi, all. I'm trying to automate some regression operations in R but am confused about how to evaluate expressoins that are expressed as character strings. For example: y <- ifelse (rnorm(10)>0, 1, 0) sex <- rnorm(10) age <- rnorm(10) test <- as.data.frame (cbind (y, sex, age)) # this works fine: glm (y ~ sex + I(age^2), data=test, family=binomial(link="logit"),

Dear Friends, I am trying to use the gelman-rubin convergence test. I generated a matrix samp[10,000x86] with the gibbs sampler. the test requires the creation of "mcmc" objects. Since I don't know how to define samp as a "mcmc" object, I tried to create one mcmc object by means of the mcmc() function. With this function I tried to create a mcmc object dul from samp but I

I want to write a function to standardize regression predictors, which will require me to do some character-string manipulation to parse the variables in a call to lm() or glm(). For example, consider the call lm (y ~ female + I(age^2) + female:black + (age + education)*female). I want to be able to parse this to pick out the input variables ("female", "age",

I decided to update my packages and then had a problem with loading the Matrix package http://cran.at.r-project.org/bin/windows/contrib/2.4/Matrix_0.9975-9.zip This is what happened when I tried to load it in: > library("Matrix") Error in importIntoEnv(impenv, impnames, ns, impvars) : object 'Logic' is not exported by 'namespace:methods' Error:

Hi, I am trying to recreate the right graph on page 524 of Gelman's 2006 paper "Prior distributions for variance parameters in hierarchical models" in Bayesian Analysis, 3, 515-533. I am only interested, however, in recreating the portion of the graph for the overlain prior density for the half-Cauchy with scale 25 and not the posterior distribution. However, when I try:

There are some interesting comments re Paul Murrell's new book, R, & SAS on Andrew Gelman's blog: http://www.stat.columbia.edu/~cook/movabletype/archives/2005/04/a_new_book_on_r.html -- Tony Plate

Dear, I'm trying to run diagnostics on MCMC analysis (fitting a log-linear model to rates data). I'm getting an error message when trying Gelman-Rubin shrink factor plot: >gelman.plot(out) Error in chol.default(W) : the leading minor of order 2 is not positive definite I take it that somewhere, somehow a matrix is singular, but how can that be remedied? My code: library(rjags)

Deal all, as MCMClogit does not allow for the specification of several chains, I have run my model 3 times with different random number seeds and differently dispersed multivariate normal priors. For example: res1 = MCMClogit(y~x,b0=0,B0=0.001,data=mydat, burnin=500, mcmc=5500, seed=1234, thin=5) res2 = MCMClogit(y~x,b0=1,B0=0.01,data=mydat, burnin=500, mcmc=5500, seed=5678, thin=5) res3 =

Thanks to some help by Doug Bates (and the updated version of the Matrix package), I've refined my question about fitting nested and non-nested factors in lmer(). I can get it to work in linear regression but it crashes in logistic regression. Here's my example: # set up the predictors n.age <- 4 n.edu <- 4 n.rep <- 100 n.state <- 50 n <- n.age*n.edu*n.rep age.id

Dear friends, I run the Gelman-Rubin Convergence test for a MCMC object I have and I got the following result Multivariate psrf 1.07+0i, What does this mean? I guess (if I am not mistaken) that I should get a psrf close to 1.00 but what is 1.07+0i? Is that convergence or something else? Jorge [[alternative HTML version deleted]]

Is there a quick way to fit student-t regressions (that is, a regression with t-distributed error, ideally with the degrees-of-freedom parameter estimated from the data)? I can do it easily enough in Bugs, or I can program the log-likelihood in R and optimize using optim(), but an R version (if it's already been written by somebody) would be convenient, especially for teaching purposes.

I am using -mice- for multiple imputation and would like to use the gelman diagnostic in -coda- to assess the convergence of my imputations. However, gelman.diag requires an mcmc list as input. van Buuren and Groothuis-Oudshoorn (2011) recommend running mice step-by-step to assess convergence (e.g. imp2 <- mice.mids(imp1, maxit = 3, print = FALSE) ) but this creates mids objects. How can I

Greetings: I'm going to encourage some students to try Bayesian ideas for hierarchical models. I want to run the WinBUGS and R examples in Tony Lancaster's An Introduction to Modern Bayesian Econometrics. That features MS Windows and "bugs" from R2WinBUGS. Today, I want to ask how people are doing this in Linux? I have found a plethora of possibilities, some of which are not

Dear R users, I've found bugs.R : the functions for running WinBUGs from R that is writen by Dr. Andrew Gelman who is a professor from Columbia University. The bugs.R would be very useful for me, and I think many of you know it as well. I followed the instuctions on Dr. Gelman's web to install all of documents that bugs.R needs, but when I try to run the school example the web posted in

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 While I'm not an R expert, I have used R on Windows XP. Now I've moved to Ubuntu (Intrepid), and I'm trying to configure R to work with the Gelman and Hill _Data Analysis Using Regression and Multilevel/Hierarchical Models_. So far, it's not working. I start by following the instructions for installing arm and BRugs at

Dear R users, I followed the instuctions on Dr. Gelman's web to install all of documents that bugs.R needs, but when I try to run the school example that the web posted in R, I got an error: couldn't find function "bugs", what's wrong? Thanks, Jia

I'm having another problem with lmer(), this time something simpler (I think) involving the coef() function for a model with varying coefficients. Here's the R code. It's a simple model with 2 observations per group and 10 groups: # set up the predictors n.groups <- 10 n.reps <- 2 n <- n.groups*n.reps group.id <- rep (1:n.groups, each=n.reps) # simulate the varying

I'm working with Andrew Gelman on a book project and we're having some difficulties getting coef() to work in some lmer() calls. Some versions of the model work and some do not. For example, this works (in that we can run the model and do coef() from the output): R2 <- lmer(y2 ~ factor(z.inc) + z.st.inc.full + z.st.rel.full + (1 + factor( z.inc) | st.num),