Can I also ask, is it also the case that the bootstrapping option in the sim() function is
not compatible with Zelig? I get this error:
Error in simpar[, 1:k] : incorrect number of dimensions
Thanks
________________________________________
From: owner-zelig_at_lists_gking_harvard_edu(a)mail.hmdc.harvard.edu
[owner-zelig_at_lists_gking_harvard_edu(a)mail.hmdc.harvard.edu] On Behalf Of Anneley
McMillan [AXM900(a)bham.ac.uk]
Sent: 03 September 2010 10:42
To: Kosuke Imai
Cc: zelig(a)lists.gking.harvard.edu
Subject: RE: [zelig] RE: Validation
Hello Kosuke,
Thanks for the info, however I'm afraid that didn't work. I read that there was
support for cv.glm through this document:
http://www.uvg.edu.gt/antropologia/estadistica/Zelig%20Documentation.pdf
using the example on page 33, but this was not with the weighted method.
Would you instead be able to suggest the best way to calculate the MSE from the zelig
relogit output, if this makes sense? I can then compare it to the prior method at least.
Thank you for your time,
Anneley
________________________________________
From: Kosuke Imai [kimai(a)Princeton.Edu]
Sent: 03 September 2010 00:17
To: Anneley McMillan
Cc: zelig(a)lists.gking.harvard.edu
Subject: Re: [zelig] RE: Validation
Zelig does not support cv.glm() but you might be able to trick Zelig by
specifying the class. That is,
class(z.out) <- c("relogit", "glm.robust", "glm")
I'm guessing this, and so use this at own your risk though!
Good luck,
Kosuke
--
Department of Politics
Princeton University
http://imai.princeton.edu
On Thu, 2 Sep 2010, Anneley McMillan wrote:
A quick update on this... I have run relogit with 'prior' instead of
'weighting' and the cv.glm() function now works. Specifically it seems to have
an issue with the glm.robust object. Can anyone explain to me why this might be the case?
Thanks in advance!
________________________________________
From: owner-zelig_at_lists_gking_harvard_edu(a)mail.hmdc.harvard.edu
[owner-zelig_at_lists_gking_harvard_edu(a)mail.hmdc.harvard.edu] On Behalf Of Anneley
McMillan [AXM900(a)bham.ac.uk]
Sent: 01 September 2010 08:48
To: zelig(a)lists.gking.harvard.edu
Subject: [zelig] Validation
Hi there,
I have a question about validation. First off, I would like to be able to use the
validate.lrm() function, however I believe that it is not supported to put zelig functions
into it yet... Is there any other options to use this using Zelig outputs at the moment?
i.e.
fit <- zelig(Y_bin ~ TMIN_WI + TMAX_SU, model = "relogit", tau = 0.0076,
case.control = "weighting", robust = TRUE, data = table, x=TRUE,y=TRUE)
where Y_bin is a dichotomous variable. You need to add the x=TRUE and y=TRUE, following
the methodology on the function manual.
Then I try to put it into a cross validation or bootstrap function with some test values
and I get the following errors:
validate(fit, method=”boot”,B=10)
Error:
unexpected input in " validate(fit, method=”"
validate(fit, B=10)
Error in
UseMethod("validate") :
no applicable method for 'validate' applied to an object of class
"c('relogit', 'glm.robust')"
My second question is for the cv.glm function, which I thought was supported:
cv.out <-cv.glm(z.out, data = table)
it again gives a similar error:
Error in UseMethod("predict") :
no applicable method for 'predict' applied to an object of class
"c('relogit', 'glm.robust')"
Could you offer any suggestions as to what may be missing here?
Thanks
-
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/
-
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/
-
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/
-
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/