similar to: mice library / survival analysis

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

I am working with some survival data with missing values. I am using the mice package to do multiple imputation. I have found code in this thread which handles pooling of the MI results: https://stat.ethz.ch/pipermail/r-help/2007-May/132180.html Now I would like to plot a survival curve using the pooled results. Here is a reproducible example: require(survival) require(mice) set.seed(2) dt

Hey there I am estimating a multilevel model using lmer. I have 5 imputed datasets so I am using mitools to pool the estimates from the 5 > > datasets. Everything seems to work until I try to use > MIcombine to produced pooled estimates. Does anyone have any suggestions? The betas and the standard errors were extracted with no problem so everything seems to work smoothly up until

Hi everyone, Does anyone here has experience using MICE to impute missing value? I am having problem to pool the imputed dataset for a MANOVA test, could you give me some advice please? Here is my code: > library(mice) >

Hi, I am a long time SPSS user but new to R, so please bear with me if my questions seem to be too basic for you guys. I am trying to figure out how to analyze survey data using logistic regression with multiple imputation. I have a survey data of about 200,000 cases and I am trying to predict the odds ratio of a dependent variable using 6 categorical independent variables (dummy-coded).

Hello I would like to pass a model formula as an argument to the with.mids function from the mice package. The with.mids functon fits models to multiply imputed data sets. Here's a simple example library(mice) #Create multiple imputations on the nhanes data contained in the mice package. imp <- mice(nahnes) #Fitting a linear model with each imputed data set the regular way works

R-users & helpers: I am using Amelia, mitools and cmprsk to fit cumulative incidence curves to multiply imputed datasets. The error message that I get "Error in eval(expr, envir, enclos) : invalid 'envir' argument" occurs when I try to fit models to the 50 imputed datasets using the "with.imputationList" function of mitools. The problem seems to occur

Hello: I am working with a stratified survey dataset with sampling weights and I want to use multiple imputation to help with missingness. 1. Is there a way to run an ordered logistic regression using both a multiply imputed dataset (i.e. from mice) and adjust for the survey characteristics using the weight variable? The Zelig package is able to do binary logistic regressions for survey

I tried multiple imputation with aregImpute() and fit.mult.impute() in Hmisc 3.8-3 (June 2010) and R-2.12.1. The warning message below suggests that summary(f) of fit.mult.impute() would only use the last imputed data set. Thus, the whole imputation process is ignored. "Not using a Design fitting function; summary(fit) will use standard errors, t, P from last imputation only. Use

Hello! If I understand this listserve correctly, I can email this address to get help when I am struggling with code. If this is inaccurate, please let me know, and I will unsubscribe. I have been struggling with the same error message for a while, and I can't seem to get past it. Here is the issue: I am using a data set that uses -1:-9 to indicate various kinds of missing data. I changed

I'm trying to use either mice or norm to perform multiple imputation to fill in some missing values in my data. The data has some missing values because of a chemical detection limit (so they are left censored). I'd like to use MI because I have several variables that are highly correlated. In SAS's proc MI, there is an option with which you can limit the imputed values that are

Hi all, How can I calculate the mean from several imputed data sets with the package mice? I know you can estimate regression parameters with, for example, lm and subsequently pool those parameters to get a point estimate using functions included in mice. But if I want to calculate the mean value of a variable over my multiple imputed data sets with fit <- with(data=imp, expr=mean(y)) and

Dear R-Users, The following question is more of general nature than a merely technical one. Nevertheless I hope someone get me some answers. I have been using the mice package to perform the multiple imputations. So far, everything works fine with the standard regressions analysis. However, I am wondering, if it is theoretically correct to perform nonparametric models (GAM, spline

Hi, I am attempting to impute my data for missing values using the mice function. However everytime I run the function it freezes or lags. I have tried running it over night, and it still does not finish. I am working with 17000 observations across 32 variables. here is my code: imputed.data = mice(data, + m = 1, + diagnostics = F) Thank you in advance, Richard [[alternative HTML version

Hi folks, I have five imputed data sets and would like to apply the same recoding routines to each. I could do this sort of thing pretty easily in Stata using MIM, but I've decided to go cold turkey on other stats packages as a incentive for learning more about R. Most of the recoding is for nominal variables, like race, religion, urbanicity, and the like. So, for example, to recode race

Hi I am trying to impute missing values for my data.frame. As I intend to use the complete data for prediction I am currently measuring the success of an imputation method by its resulting classification error in my training data. I have tried several approaches to replace missing values: - mean/median substitution - substitution by a value selected from the observed values of a variable - MLE

R 2.2.0 windows XP I am beginning to explore the mitools package contributed by Thomas Lumley (thank you Thomas) and I have a few questions: (1) In the examples given in the mitools documentation, the only family argument used is family=binomial. Does the package support family=gaussian and other link functions? I ran the with function with family=gaussian and I obtained results, but I am not

Dear R users I want to combine multiply imputed data-sets generated from mice to do classfication. However, I have various questions regarding the use of mice library. For example suppose I want to predict the class in this data.frame: data(nhanes) mydf=nhanes mydf$class="pos" mydf$class[sample(1:nrow(mydf), size=0.5*nrow(mydf))]="neg" mydf$class=factor(mydf$class) First I

Dear all, I request your help to solve a problem I've encountered in using 'mice' for multiple imputation. I want to apply a logistic regression model. I need to extract information on the fit of the model. Is there any way to calculate a likelihood ratio or the McFadden-pseudoR2 from the results of the logistic model? I mean, as it is possible to extract pooled averaging and odds

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