similar to: calculating gelman diagnostic for mice object

Dear R-users, Version V2.0 of the package mice is now available on CRAN for Windows, Linux and Apple users. Multivariate Imputation by Chained Equations (MICE) is the name of software for imputing incomplete multivariate data by Fully Conditional Specifcation (FCS). MICE V1.0 appeared in the year 2000 as an S-PLUS library, and in 2001 as an R package. MICE V1.0 introduced predictor selection,

Dear R-users, Version V2.0 of the package mice is now available on CRAN for Windows, Linux and Apple users. Multivariate Imputation by Chained Equations (MICE) is the name of software for imputing incomplete multivariate data by Fully Conditional Specifcation (FCS). MICE V1.0 appeared in the year 2000 as an S-PLUS library, and in 2001 as an R package. MICE V1.0 introduced predictor selection,

Hello, I have a general question about combining imputations as well as a question specific to the rms and Hmisc packages. The situation is multiple regression on a data set where multiple imputation has been used to give M imputed data sets. I know how to get the combined estimate of the covariance matrix of the estimated coefficients (average the M covariance matrices from the individual

Dear list, I'd like to use multiple imputations to try and save a somewhat badly mangled dataset (lousy data collection, worse than lousy monitoring, you know that drill... especially when I am consulted for the first time about one year *after* data collection). My dataset has 231 observations of 53 variables, of which only a very few has no missing data. Most variables have 5-10% of

Hi Folks, I don't often grumble, but this time I've found myself inconvenienced by a posting stored on R-help archives: https://stat.ethz.ch/pipermail/r-help/2009-September/211095.html This was Karin Groothuis-Oudshoorn & Stef van Buuren's message on 10 September about the new version of MICE. This has been sent by software which inserted no line-breaks. As a result, each

Hi, I have the following problem- I want to access a list whose elements are imp1, imp2, imp3 etc I tried theusing the paste comand in a for loop see the last for loop below. But I keep calling it df but df = imp1 (for the first run). Any ideas on how I can access the elements of the list? Isaac require(Amelia) library(Amelia) data.use <- read.csv("multiplecarol.CSV", header=T)

Hello R-Users, I would like to check the convergence of my imputations. However, when I use the function the plot.mids(), I always obtain the following error message Error in plot.new() : figure margins too large I read the same question in thread from November 2005 (see below). I actually have the same problem. Is it now possible to plot subsets of mids.objects. If yes, how? My

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)

Hello, Can anyone help me with the following: if one enters a function name in the R console then usually one sees the code of that function. But there are functions that one cannot see. For example I want to see the code of print.htest or t.test.default. These functions are non-visible. Is it possible to see the code anyway? Thanks in advance, Karin Groothuis-Oudshoorn [[alternative HTML

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

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

Hi everyone, I have carried out a multiple imputation in R using Amelia II and have created 5 multiply imputed datasets. The purpose of my research is to fit a Poisson Model to the data to estimate numbers of hospital admissions. Now that I have 5 completed datasets and I have to pool all the 5 datasets to get one combined output for a poisson model. I have checked previous queries about

R-helpers: I have a dataset that has 168 subjects and 12 variables. Some of the variables have missing data and I want to use the multiple imputation capabilities of the "mice" package to address the missing data. Given that mice only supports linear models and generalized linear models (via the lm.mids and glm.mids functions) and that I need to fit Cox models, I followed the previous

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]]

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:

Hi everyone! I have data consisting of several response variables and several explanatory variables. I wish to do a permanova on this using the vegan library and the adonis() function. However, my data had several missing values in it. In order to 'fix' this I used the mice() function from the mice library to make 5 imputations for all the missing values. To do analysis on the 5 datasets

Dear all, I am trying to impute data for a range of variables in my data set, of which unfortunately most variables have missing values, and some have quite a few. So I set up the predictor matrix to exclude certain variables (setting the relevant elements to zero) and then I run the imputation. This works fine if I use predictive mean matching for the continous variables in the data set. When I

Hello, I have installed R on my Apple Laptop. Next I wanted to install the package Rcmdr which requires the package tcltk. But then I get errors like: The downloaded packages are in /var/folders/0p/0pD8fDrwHouNDsQ+k8dGmU+++TI/-Tmp-//RtmpSp4q7p/ downloaded_packages Loading required package: tcltk Loading Tcl/Tk interface ... Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to

Hi, temp3<- read.table(text=" ID CTIME WEIGHT HM001 1223 24.0 HM001 1224 25.2 HM001 1225 23.1 HM001 1226 NA HM001 1227 32.1 HM001 1228 32.4 HM001 1229 1323.2 HM001 1230 27.4 HM001 1231 22.4236 #changed here to test the previous solution ",sep="",header=TRUE,stringsAsFactors=FALSE) ?tempnew<- na.omit(temp3) ?grep("\\d{4}",temp3$WEIGHT) #[1] 7 9 #not correct

Package: mcgibbsit Title: Warnes and Raftery's MCGibbsit MCMC diagnostic Version: 1.0 Author: Gregory R. Warnes <gregory_r_warnes at groton.pfizer.com> Description: mcgibbsit provides an implementation of Warnes & Raftery's MCGibbsit run-length diagnostic for a set of (not-necessarily independent) MCMC sampers. It combines the estimate error-bounding approach of Raftery