Zelig April 2012

zelig@lists.gking.harvard.edu

7 participants
7 discussions

Re: MatchIt Zelig and multiple imputation

by Kosuke Imai

Unfortunately, I don't think we have an automated procedure for everything. You would have to multiply impute the data, do matching on each imputed data set, and then combine it in zelig() using mi() function. But this does not require any programming. You can simply run the same matching procedure on each data set via matchit() and then feed the resulting multiple matched data sets into zelig(). Good luck, Kosuke Department of Politics Princeton University http://imai.princeton.edu On Sep 13, 2011, at 6:02 PM, Pingaul jb wrote: > Dear Professor, > I’m a post-doctoral student at Montreal University. I’m actually in Columbia, working and propensity scores with a colleague and using MatchIt and Zelig. First, congratulations for your packages that are very flexible. > > My question is about multiple imputation and propensity scores with these softwares. From what I understand, combining both approaches would include: > > 1/ Doing multiple imputation and testing which variables to include. > > 2/ Propensity score analysis on each imputed data set and pooling the overall balance to check if it is ok (or on each data set?). > > 3/ Calculation of the quantities of interest for each data set > > 4/ Pooling the quantities across data sets. > > I would like to know if there is a written syntax to perform the MatchIt analysis for all of the imputed data set without having to do it manually and check the overall balance. Also, in theory, the number of individuals retained after propensity score matching and the weights can be different for each imputed data set. So that we have to perform the final analysis on each one and then pool the data with a specific procedure to take into account the eventual varying Ns? I normally use Mice package for multiple imputation but it seems that Zelig handle Amelia. My colleague seems to do be able to do all that in stata, but I’m not sure how to make all the three R packages work together. > > I would be very happy if you could indicate to me a reference or a place where I can find the syntax to do that (I’ve been using R for some times so I can use packages easily but I have no programming skills). > > > Best Regards! > > > > Jean-Baptiste > - Zelig Mailing List, served by Harvard-MIT Data Center Send messages: zelig(a)lists.gking.harvard.edu [un]subscribe Options: http://lists.gking.harvard.edu/?info=zelig Zelig program information: http://gking.harvard.edu/zelig/

10 years, 8 months

ZeligMultinomial

by Josue Guzman

*Please, let me know how can I install ZeligMultinomial package. I want to use mlogit, which according to the manual (page 50), is found in said package. I tried with the command, install.packages("ZeligMultinomial", repos="http://r.iq.harvard.edu/", type="source") But received the following message: Warning: dependency 'MNP' is not available trying URL ' http://r.iq.harvard.edu/src/contrib/ZeligMultinomial_0.5-4.tar.gz' Content type 'application/x-gzip' length 9730 bytes opened URL ================================================== downloaded 9730 bytes During startup - Warning messages: 1: Setting LC_CTYPE failed, using "C" 2: Setting LC_TIME failed, using "C" 3: Setting LC_MESSAGES failed, using "C" 4: Setting LC_PAPER failed, using "C" ERROR: dependency 'MNP' is not available for package 'ZeligMultinomial' * removing '/Library/Frameworks/R.framework/Versions/2.14/Resources/library/ZeligMultinomial' The downloaded packages are in '/private/var/folders/UL/ULLu+bi5GR8IM1G-7XZBJU+++TI/-Tmp-/Rtmp6jw6M9/downloaded_packages' Warning message: In install.packages("ZeligMultinomial", repos = "http://r.iq.harvard.edu/", : installation of package 'ZeligMultinomial' had non-zero exit status ***

11 years, 5 months

Gamma regression

by Salvador Becerra-Munoz

Please I need your help! This code sometimes predicts negative values because it is based on the inverse link. What is the modification to calculate de parameters using the loglink option? z.out = zelig(PT ~ Q + S + T + TDG, model = "gamma", data = PT) summary(z.out) x.out = setx(z.out) s.out = sim(z.out, x = x.out) plot(s.out) Thanks a lot for your advice

12 years

Conditional Frailty Models in Zelig

by Brenton Peterson

Hi, I am trying to estimate a conditional frailty survival model, as in Box-Steffensmeier, et al (2007 Political Analysis). I have run them previously in R but would like to implement them in Zelig for post-estimation purposes. Is there a way to do this? While the Cox PH with time-varying covariates is closest to my needs, it does not completely fit the bill. Is this an option in the current Zelig version? Thanks, B Peterson

12 years

readRDS

by Matthias Kretschmer

I tried to install Zelig on R 2.15 (on a Windows 7, 64bit machine) and keep getting an error message: *Error : .onAttach failed in attachNamespace() for 'Zelig', details: call: NULL error: '.readRDS' is defunct. Use 'readRDS' instead. See help("Defunct") Error: package/namespace load failed for ‘Zelig’ * Any idea what that means? Matthias

12 years

Cox Proportional Hazards Regression

by Johannes Castner

To whom this may concern, I am trying to run a Cox Proportional Hazards Regression in the following way, using Zelig: z.out <- zelig(Surv(duration, Censored) ~ Ethnic1 + Type + Year1, model = "coxph", data = data) but I get the following error: Error in fitter(X, Y, strats, offset, init, control, weights = weights, : NA/NaN/Inf in foreign function call (arg 6) In addition: Warning message: In fitter(X, Y, strats, offset, init, control, weights = weights, : Ran out of iterations and did not converge Does anyone know just what that means and how I can get around it? Johannes

12 years, 1 month

AUTO: Tom Aldenberg is out of the office (returning Fri 04/13/2012)

by Tom Aldenberg

I am out of the office from Sat 03/31/2012 until Fri 04/13/2012. Note: This is an automated response to your message "Zelig Digest, Vol 20, Issue 1" sent on 5-4-2012 18:00:02. This is the only notification you will receive while this person is away.

12 years, 1 month

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

2003

Zelig April 2012