29 Apr
2008
29 Apr
'08
3:54 p.m.
Dear Zelig List,
I am trying to do estimate a random effects logit model using multiply imputed
data, but I get the following error message:
out<-zelig(cwarcons~rival_lag+weak10+civwar_nb+lref_totalnbl+lpopl
+ +lrgdp96l+polity2l+politysq+ethfrac+peace1+eeurop+lamerica+ssafrica+asia+
+ nafrme+lmtnest+elevdiff+Oil+instab+tag(1|ccode),data=mi(m1,m2,m3,m4,m5),
+ model="logit.mixed")
summary(out)
Error in apply(coef1, 1, mean) :
dim(X) must have a positive length
In addition: Warning messages:
1: In x$coef : $ operator is invalid for atomic vectors, returning NULL
2: In x$coef : $ operator is invalid for atomic vectors, returning NULL
3: In x$coef : $ operator is invalid for atomic vectors, returning NULL
4: In x$coef : $ operator is invalid for atomic vectors, returning NULL
5: In x$coef : $ operator is invalid for atomic vectors, returning NULL
Any ideas about how I might fix this?
Thanks,
Didi Kuo
-
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/
29 Apr
29 Apr
10:04 p.m.
Unfortunately, the mixed effects models do not work with multiply imputed
data sets yet. This is certainly on our to-do list. Right now, you would
have to do this by hand - apply zelig(), setx(), and sim() to each of the
data sets and then combine sim results.
Kosuke
---------------------------------------------------------
Kosuke Imai Office: Corwin Hall 041
Assistant Professor Phone: 609-258-6601
Department of Politics eFax: 973-556-1929
Princeton University Email: kimai(a)Princeton.Edu
Princeton, NJ 08544-1012 http://imai.princeton.edu/
---------------------------------------------------------
On Tue, 29 Apr 2008, Didi Kuo wrote:
Dear Zelig List,
I am trying to do estimate a random effects logit model using multiply imputed
data, but I get the following error message:
out<-zelig(cwarcons~rival_lag+weak10+civwar_nb+lref_totalnbl+lpopl
+ +lrgdp96l+polity2l+politysq+ethfrac+peace1+eeurop+lamerica+ssafrica+asia+
+ nafrme+lmtnest+elevdiff+Oil+instab+tag(1|ccode),data=mi(m1,m2,m3,m4,m5),
+ model="logit.mixed")
-
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/
summary(out)
Error in apply(coef1, 1,
mean) : dim(X) must have a positive length
In addition: Warning messages:
1: In x$coef : $ operator is invalid for atomic vectors, returning NULL
2: In x$coef : $ operator is invalid for atomic vectors, returning NULL
3: In x$coef : $ operator is invalid for atomic vectors, returning NULL
4: In x$coef : $ operator is invalid for atomic vectors, returning NULL
5: In x$coef : $ operator is invalid for atomic vectors, returning NULL
Any ideas about how I might fix this?
Thanks,
Didi Kuo
-
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/
1 May
1 May
8:08 a.m.
hi--
I guess one thing to keep in mind is that if you want to interpret the
results of the random effects estimation (meaning the variance
components), you will have to impute under that model as well. Doing so
is kind of a nuisance (you could use the package pan). One alternative
is to impute a wide file (meaning y1,y2, x,z ....), but that's not very
useful if the problem is very difficult.
i'd be curious how you imputed under a random-effects specification?
regards,
michael
Kosuke Imai wrote:
Unfortunately, the mixed effects models do not work
with multiply
imputed data sets yet. This is certainly on our to-do list. Right
now, you would have to do this by hand - apply zelig(), setx(), and
sim() to each of the data sets and then combine sim results.
Kosuke
---------------------------------------------------------
Kosuke Imai Office: Corwin Hall 041
Assistant Professor Phone: 609-258-6601
Department of Politics eFax: 973-556-1929
Princeton University Email: kimai(a)Princeton.Edu
Princeton, NJ 08544-1012 http://imai.princeton.edu/
---------------------------------------------------------
On Tue, 29 Apr 2008, Didi Kuo wrote:
--
E. Michael Foster
getmike.info
download vcard: www dot unc dot edu/~emfoster/Foster.vcf
Professor of Maternal and Child Health
of Health Policy and Administration
of Biostatistics (adjunct)
Associate Chair for Faculty Advancement
(W) 919-966-3773
(F) 919-966-0458
Office: 422D Rosenau Hall << use for FEDEX
School of Public Health
University of North Carolina, Chapel Hill
Rosenau Hall, Campus Box# 7445
Chapel Hill, NC 27599-7445
[the http is throwing off some spam filters)
visit the new FAQ section of my blog:
fosterfaq dot blogspot dot com
27) [Mollie Ivins] "What do you mean cancer is a gift. You're as a dumb as
Shrub. Cancer just killed me. Some gift."
32) If I have a vacant look on my face, what should you do?
33) If you were going to Brokeback mountain, whom would you take? (redux)
Hint:
"And I got the one I love beside me
My troubles behind me
I'm alive and I'm free
Who wouldn't wanna be me"
"I'm not quite dead yet, sir!" -- Concorde, Monty Python and the Holy
Grail
"It is a fearful thing to love what death can touch" -- Chaim Stern
"Manure is an art and I consider myself an accomplished artist" ~John Adams.
-
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/
Dear Zelig List,
I am trying to do estimate a random effects logit model using
multiply imputed
data, but I get the following error message:
out<-zelig(cwarcons~rival_lag+weak10+civwar_nb+lref_totalnbl+lpopl
+
+lrgdp96l+polity2l+politysq+ethfrac+peace1+eeurop+lamerica+ssafrica+asia+
+
nafrme+lmtnest+elevdiff+Oil+instab+tag(1|ccode),data=mi(m1,m2,m3,m4,m5),
+ model="logit.mixed")
-
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/ summary(out)
Error in apply(coef1, 1,
mean) : dim(X) must have a positive length
In addition: Warning messages:
1: In x$coef : $ operator is invalid for atomic vectors, returning NULL
2: In x$coef : $ operator is invalid for atomic vectors, returning NULL
3: In x$coef : $ operator is invalid for atomic vectors, returning NULL
4: In x$coef : $ operator is invalid for atomic vectors, returning NULL
5: In x$coef : $ operator is invalid for atomic vectors, returning NULL
Any ideas about how I might fix this?
Thanks,
Didi Kuo
-
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/
5862
days inactive
5864
days old
2 comments
3 participants
participants (3)
-
Didi Kuo -
E. Michael Foster -
Kosuke Imai