similar to: the mice pool function and aov()

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

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

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

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

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

Hello everyone!I have a data set with missing observations that I am trying to impute. I am using MICE and I would like the imputed values to all be positive. I have two types of variables: prices (P1 to P136) and quantities (Q1 to Q136) and I also want the range of these two types to be different. Besides these variables. I am using the squeeze function but I am unable to set it such that I

Dear list, Is there any one use MICE package deal with multilevel missing values here? I have a question about the 2lonly.pmm() and 2lonly.norm(), I get the following error quite often. Here is the code the error, could you give me some advice please? Am I using it in the right way? > ini=mice(bhrm,maxit=0) > pred=ini$pred > pred V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15

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

Hi, I'm using the library MICE to make multiple imputations. I'can pool the results to show how the predicted values fit, but how to combine the five imputed datasets? I take the mean? for exemple : x1<-complete(imp) x2<-complete(imp, 2) x3<-complete(imp, 3) x4<-complete(imp, 4) x5<-complete(imp, 5) impcomp<-(x1+x2+x3+x4+x5)/5 Or there is an other way to do this?

Hello all, I am trying to impute some missing data using the mice package. The data set I am working with contains 125 variables (190 observations), involving both categorical and continuous data. Some of these variables are missing up to 30% of their data. I am running into a peculiar problem which is illustrated by the following example showing both the original data (blue) and the imputed

Hi I am unable to impute using the MICE command in R when imputing a binary variable using linear discriminant analysis. To illustrate my problem I have created a dataset, which consists of 1 continuous and 1 binary variable. The continuous variable is complete and the binary variable is partially observed. I am able to impute using the MICE command where the imputation methods is logistic

Dear list, I have a question about imputing 2 level data in MICE, could you give me some suggestions please? Thank you very much. The data set contains 35634 cases and 1007 variables, 280 of them are categorical variables, and the rest of them are continuous variables. On the second level, there are 198 units. I am trying to impute missing values for 270 categorical variables by using the

R 2.7.2 PPC Mac OS X 10.4.11 library mice 1.13.1 I try to use mice for multivariate data imputation. My variables are numeric, factors, count data, ordered factors. First I created a vector for the methods to use with each variable ImpMethMice<-c(rep("logreg", 62), rep("polyreg",1), rep("norm",12), rep("polyreg",12)) next step was

Saludo gente, antes que nada gracias por la ayuda que puedan aportarme, soy iniciante en R, estoy usando el paquete Mice para realizar imputaciones múltiples sobre variables en su mayoría categóricas. El problema está que cuando expresó este comando imp <- mice(dataset,method="polr",maxit=1) donde el dataset es un data.frame me tirá este error : iter imp variable 1 1 pial1a

Hi all, I have a question about multiple imputation within the MICE package. I want to use passive imputation for my variable called X, because it is calculated out of multiple variables, namely Y, Z. Let's give an example with BMI. I know, that if I want to use passive imputation for BMI, I can use the following command: meth["BMI"] <- "~I(weight/(height/100)^2)"

Saludo gente, antes que nada gracias por la ayuda que puedan aportarme, soy iniciante en R, estoy usando el paquete Mice para realizar imputaciones múltiples sobre variables en su mayoría categóricas. El problema está que cuando expresó este comando imp <- mice(dataset,method="polr",maxit=1) donde el dataset es un data.frame me tirá este error : iter imp variable 1 1 pial1a

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

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