Zelig February 2010

zelig@lists.gking.harvard.edu

4 participants
6 discussions

summary with imputations bug ... any workarounds?

by Peter Flom

Good morning If I run <<< susan.lsmixed.out <- zelig(formula = unprot_vag_sex ~ married + age + TREATMENT.ARM*time + highest_grade + income + tag(1|id), data = susanMI.out$imputations, model = "ls.mixed") summary(susan.lsmixed.out) >>>> I get an error Error in x$coef : $ operator is invalid for atomic vectors Searching the archives, I see that others have had similar problems. Is there a workaround? summary(susan.lsmixed.out[[1]]) works fine; should I then average across the five imputed data sets? thanks! Peter Peter L. Flom, PhD Statistical Consultant Website: http://www DOT statisticalanalysisconsulting DOT com/ Writing; http://www.associatedcontent.com/user/582880/peter_flom.html Twitter: @peterflom - 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/

12 years, 1 month

Re: [zelig] setx generates "incorrect number of dimensions" error

by Donald Braman

I'm still having trouble with setx. Factoring the leftside prior doesn't help. Here's a link to the data & script: http://www.culturalcognition.net/storage/datasets/zelig.problem.zip > On Thu, Feb 11, 2010 at 12:09 AM, Kosuke Imai <kimai(a)princeton.edu> wrote: > >> Can you try doing as.factor() outside of the formula? That is, >> >> miset1$ssimportant <- as.factor(miset1$ssimportant) >> >> etc. for each data set. If this does not solve the problem, please let us >> know. >> >> Thanks, >> Kosuke >> >> Department of Politics >> Princeton University >> http://imai.princeton.edu >> >> On Feb 10, 2010, at 11:16 PM, Donald Braman wrote: >> >> > I've started to get this error when I run setx: >> > >> > Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] : >> > incorrect number of dimensions >> > >> > Any ideas as to why? >> > >> > >> > Here's an example: >> > >> > > z.out <- zelig(formula="as.factor(ssimportant) ~ egalitarianism + >> individualism ", model="ologit", >> data=mi(c(miset1,miset2,miset3,miset4,miset5))) >> > How to cite this model in Zelig: >> > Kosuke Imai, Gary King, and Oliva Lau. 2007. "ologit: Ordinal Logistic >> Regression for Ordered Categorical Dependent Variables" in Kosuke Imai, Gary >> King, and Olivia Lau, "Zelig: Everyone's Statistical Software," >> http://gking.harvard.edu/zelig >> > > >> > > summary(z.out) >> > Call: >> > zelig(formula = "as.factor(ssimportant) ~ egalitarianism + individualism >> ", >> > model = "ologit", data = mi(c(miset1, miset2, miset3, miset4, >> > miset5))) >> > >> > Coefficients: >> > Value Std. Error t value >> > egalitarianism -0.1078 0.06161 -1.749 >> > individualism -0.1054 0.06724 -1.567 >> > >> > Intercepts: >> > Value Std. Error t value >> > 1|2 -2.682 0.473 -5.673 >> > 2|3 -2.037 0.470 -4.334 >> > 3|4 -1.406 0.468 -3.002 >> > 4|5 -0.416 0.468 -0.890 >> > 5|6 0.701 0.469 1.495 >> > >> > Residual Deviance: 5568.83 >> > AIC: 5582.83 >> > > >> > > x.out <- setx(z.out) >> > Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] : >> > incorrect number of dimensions >> > >> > >> > -- >> > Donald Braman >> > phone: 413-628-1221 >> > http://www.culturalcognition.net/braman/ >> > http://ssrn.com/author=286206 >> > http://www.law.gwu.edu/Faculty/profile.aspx?id=10123 >> > >> >> >> >> >> > > > -- > Donald Braman > phone: 413-628-1221 > http://www.culturalcognition.net/braman/ > http://ssrn.com/author=286206 > http://www.law.gwu.edu/Faculty/profile.aspx?id=10123 > > > -- Donald Braman phone: 413-628-1221 http://www.culturalcognition.net/braman/ http://ssrn.com/author=286206 http://www.law.gwu.edu/Faculty/profile.aspx?id=10123 -- Donald Braman phone: 413-628-1221 http://www.culturalcognition.net/braman/ http://ssrn.com/author=286206 http://www.law.gwu.edu/Faculty/profile.aspx?id=10123 -- Donald Braman phone: 413-628-1221 http://www.culturalcognition.net/braman/ http://ssrn.com/author=286206 http://www.law.gwu.edu/Faculty/profile.aspx?id=10123

14 years, 2 months

setx generates "incorrect number of dimensions" error

by Donald Braman

I've started to get this error when I run setx: Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] : incorrect number of dimensions Any ideas as to why? Here's an example: > z.out <- zelig(formula="as.factor(ssimportant) ~ egalitarianism + individualism ", model="ologit", data=mi(c(miset1,miset2,miset3,miset4,miset5))) How to cite this model in Zelig: Kosuke Imai, Gary King, and Oliva Lau. 2007. "ologit: Ordinal Logistic Regression for Ordered Categorical Dependent Variables" in Kosuke Imai, Gary King, and Olivia Lau, "Zelig: Everyone's Statistical Software," http://gking.harvard.edu/zelig > > summary(z.out) Call: zelig(formula = "as.factor(ssimportant) ~ egalitarianism + individualism ", model = "ologit", data = mi(c(miset1, miset2, miset3, miset4, miset5))) Coefficients: Value Std. Error t value egalitarianism -0.1078 0.06161 -1.749 individualism -0.1054 0.06724 -1.567 Intercepts: Value Std. Error t value 1|2 -2.682 0.473 -5.673 2|3 -2.037 0.470 -4.334 3|4 -1.406 0.468 -3.002 4|5 -0.416 0.468 -0.890 5|6 0.701 0.469 1.495 Residual Deviance: 5568.83 AIC: 5582.83 > > x.out <- setx(z.out) Error in dta[complete.cases(mf), names(dta) %in% vars, drop = FALSE] : incorrect number of dimensions -- Donald Braman phone: 413-628-1221 http://www.culturalcognition.net/braman/ http://ssrn.com/author=286206 http://www.law.gwu.edu/Faculty/profile.aspx?id=10123

14 years, 2 months

plot.ci() function options

by nela filimon

With apologies for cross posting Dear Zelig users, I'm working with the "sim()" function and I would need a hand with the following issues related to the graphical representation of the output obtained: 1. would it be possible to plot in a same graph, various outputs of the sim function? I do a graph for each case, but a single one would help a lot with comparisons among categorical variables and also save space in articles. 2. plot.ci() is graphing the confidence intervals; is it possible to have the representation of the mean values generated by the sim() instead of the confidence intervals? or both of them (in the same graph delivered by plot.ci) 3. Is it possible to selest for each variable a different type of line or symbol for the graphical representation? (like with colours; e.g. at page 52 in Ch. 5 . "Graphing commands" (Zelig documentation), college and high school education variables are represented with different colours but with the "same" type of line). I use text() option to put labels on the curves but I was thinking that maybe who knows, someone has already faced similar issues. Thanks a lot before hand for your help. Nela Filimon Unversitat de Girona, Spain - 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/

14 years, 2 months

plot function options

by nela filimon

With apologies for cross posting Dear Zelig users, I'm working with the "sim()" function and I would need a hand with the following issues related to the graphical representation of the output obtained: 1. would it be possible to plot in a same graph, various outputs of the sim function? I do a graph for each case, but a single one would help a lot with comparisons among variables and also save space in articles. 2. plot.ci() is graphing the confidence intervals; is it possible to have instead the representation of the "mean values" generated by the sim()? 3. selecting a different colour for each variable is very useful, especially in electronic format but when comes about printing the distinction becomes problematic (I use the text() options to put labels, and it works very well). Is it possible to selest for each variable a different type of line or symbol for the graphical representation? I've tried but I couldn't manage. Thanks a lot before hand for your help. Nela Filimon Unversitat de Girona, Spain - 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/

14 years, 3 months

error with poisson.survey

by Casey Klofstad

Hey all- I'm trying to run a weighted Poisson model (I did full matching and need to weight the model). The code looks like this: z.out <- zelig(polact ~ treat + pid_str + q44 + q52 + contact, weights="weights", data = mi(ds1, ds2, ds3, ds4, ds5), model = "poisson.survey") I get these errors: Error in 1/as.matrix(weights) : non-numeric argument to binary operator In addition: Warning message: In poisson(link = link) : use of poisson(link=link) is deprecated I am running version 3.4-5 in Windows 7. Any thoughts? Thanks much, -c -- Casey A. Klofstad University of Miami Department of Political Science Coral Gables, FL klofstad(a)gmail.com http://www.as.miami.edu/personal/cklofstad/ - 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/

14 years, 3 months

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

2003

Zelig February 2010