Matt, what happened to summary.mi() from previous Zelig releases? We used
to combine the MI output using the Rubin rules automatically in
summary(z.out-of-class-mi).
On Fri, Dec 14, 2012 at 12:04 PM, Matt Owen <mowen(a)iq.harvard.edu> wrote:
Hi Florian -
Zelig 4 fits each of these models separately at the moment, and as a
result . If you would like to merge these datasets into a single data
table, then the usual practice is with this line:
z.out <- zelig(democ ~ loglrgdpch + loglodac + lurban + logpop + lprimsch
+ llife +
lfuelx + alesinae + alesinar + ldemoc +
loglrgdpchldemoc +
loglodacldemoc + lurbanldemoc + logpopldemoc +
lprimschldemoc +
llifeldemoc + lfuelxldemoc + alesinaeldemoc +
alesinarldemoc,
model="probit",
data = rbind(outdata1, outdata2,
outdata3, outdata4, outdata5, outdata6, outdata7,
outdata8, outdata9, outdata10,
outdata11, outdata12, outdata13, outdata14,
outdata15, outdata16, outdata17, outdata18,
outdata19, outdata20, outdata21, outdata22,
outdata23, outdata24, outdata25, outdata26,
outdata27, outdata28, outdata29, outdata30,
outdata31, outdata32, outdata33, outdata34,
outdata35, outdata36, outdata37, outdata38,
outdata39, outdata40, outdata41, outdata42,
outdata43, outdata44, outdata45, outdata46,
outdata47, outdata48, outdata49, outdata50))
Note that the call to "mi" has been replaced with "rbind". This will
workd
if outdata1-50 have similar column structure.
A feature request has been opened to support this. This can be found here:
https://github.com/IQSS/Zelig/issues/32
On Dec 14, 2012, at 2:13 PM, Florian Reiche wrote:
Hello,
I have run my probit model in Zelig after using AmeliaII to impute the
data sets. When I use the "summary(z.out)" command, I receive the
coefficients for 50 separate data analyses, but no combined output for all
data sets together. Am I doing anything wrong?
Here is the code:
z.out <- zelig(democ ~ loglrgdpch + loglodac + lurban + logpop + lprimsch
+ llife +
lfuelx + alesinae + alesinar + ldemoc +
loglrgdpchldemoc +
loglodacldemoc + lurbanldemoc + logpopldemoc +
lprimschldemoc +
llifeldemoc + lfuelxldemoc + alesinaeldemoc +
alesinarldemoc,
model="probit",
data = mi(outdata1, outdata2,
outdata3, outdata4, outdata5, outdata6, outdata7,
outdata8, outdata9, outdata10,
outdata11, outdata12, outdata13, outdata14,
outdata15, outdata16, outdata17, outdata18,
outdata19, outdata20, outdata21, outdata22,
outdata23, outdata24, outdata25, outdata26,
outdata27, outdata28, outdata29, outdata30,
outdata31, outdata32, outdata33, outdata34,
outdata35, outdata36, outdata37, outdata38,
outdata39, outdata40, outdata41, outdata42,
outdata43, outdata44, outdata45, outdata46,
outdata47, outdata48, outdata49, outdata50))
summary(z.out)
Any help would be greatly appreciated.
All best
Florian
-
--
Zelig Mailing List, served by HUIT
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/mailman/listinfo/zelig
Zelig program information:
http://gking.harvard.edu/zelig/
Zelig mailing list
Zelig(a)lists.gking.harvard.edu
To unsubscribe from this list or get other information:
https://lists.gking.harvard.edu/mailman/listinfo/zelig
-
--
Zelig Mailing List, served by HUIT
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/mailman/listinfo/zelig
Zelig program information:
http://gking.harvard.edu/zelig/
Zelig mailing list
Zelig(a)lists.gking.harvard.edu
To unsubscribe from this list or get other information:
https://lists.gking.harvard.edu/mailman/listinfo/zelig