I'm having the exact same problem. I'm running R 2.0.1 on Windows XP
(fully updated) and I've updated Zelig, MASS, and boot using the
install.packages command. My classmates discovered that the problem
"goes away" if you go back to using R 1.9.1. Unfortunately, much of our
work is done in a lab where we cannot changes versions. The lab
computers are also running XP.
Thanks,
Stephen Haptonstahl
Washington University in St. Louis
=================================================
Could you let us know what operating system you are using? It looks like
you need to rebuild a few packages: MASS and boot. Try
install.packages("MASS")
install.packages("boot")
And then reload Zelig.
Hope this helps,
Kosuke
---------------------------------------------------------
Kosuke Imai Office: Corwin Hall 041
Assistant Professor Phone: 609-258-6601 Department of Politics
eFax: 973-556-1929
Princeton University Email: kimai@xxxxxxxxxxxxx
Princeton, NJ 08544-1012 http://www.princeton.edu/~kimai
<http://www.princeton.edu/%7Ekimai>
---------------------------------------------------------
On Thu, 10 Mar 2005, Mitch Morehart wrote:
> Dr. King,
>
> Sorry to bother you directly, but I am having trouble accessing the
> zelig listserv. I am trying to estimate the marginal effects for a
> multinomial logit model and the Zelig library provides some very useful
> functions towards that end. I installed the latest version of the
> library along with the latest version of the VGAM library and tried to
> run the example code for multinomial logit. The details of the R
> session and error for the sim procedure are provided below:
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 2.0.0 (2004-10-04), ISBN 3-900051-07-0
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R or R packages in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for a HTML browser interface to help.
> Type 'q()' to quit R.
>
>
> Welcome to R!
>
> > library(Zelig)
> Loading required package: MASS Loading required package: boot
> Please refer to http://gking.harvard.edu/zelig for full documentation
> or help.zelig() for help with commands and models supported by Zelig.
>
> Warning messages: 1: package 'MASS' was built under R version 2.0.1 2:
> package 'boot' was built under R version 2.0.1 3: S3 method(s)
> 'plot.zelig.beta' were declared in NAMESPACE but not
> found > data(mexico)
> > z.out <- zelig(vote88 ~ pristr + othcok + othsocok, model =
> "mlogit",
> + data = mexico)
> Loading required package: VGAM Loading required package: splines
> Loading required package: stats4 Creating a new generic function for
> "print" in "VGAM" Creating a new generic function for "coefficients"
> in "VGAM" Creating a new generic function for "fitted.values" in "VGAM"
> Attaching package 'VGAM':
>
>
> The following object(s) are masked from package:splines :
>
> bs ns
>
> The following object(s) are masked from package:boot :
>
> logit simplex
>
> The following object(s) are masked from package:stats :
>
> binomial Gamma gaussian inverse.gaussian poisson poly
> quasibinomial quasipoisson
>
> The following object(s) are masked from package:base :
>
> is.R jitter scale.default
> > x.weak <- setx(z.out, pristr = 1)
> > x.strong <- setx(z.out, pristr = 3)
> > s.out <- sim(z.out, x = x.strong, x1 = x.weak)
> Error in array("", c(n, 3), list(names(object), c("Length", "Class",
> "Mode"))) : length of dimnames [1] not equal to array extent
> >
> Any thoughts on a solution would be appreciated. I looked at the
> source code and my only thought is that there is some incompatability
> betwen the latest version of vglm and the zelig sim procedure.
>
>
> Mitchell Morehart
> Room 4121 N
> 1800 M Street NW
> Washington, DC 20036
> morehart@xxxxxxxxxxxx
>
> Be sure to visit the ERS website
> www.ers.usda.gov
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
Dear All,
I am new to zelig but not new to R, and I use Clarify with Stata. I
think the idea behind zelig is great. The only improvement I can see
is to make zelig work with nlme (and lme4). Since now panel data and
multilevel design have become quite common, but good multilevel
software that integrate multiple imputation, simulation, and
estimation in one intergrated package is still to come (GLLAMM has
simulation capabiltiy, but it is too slow to process complicated
models or large data sets; mlwin has experimental multiple imputation
macro, but it crashes all the time).
Thanks!
Shige Song
Assistant Professor
Institute of Sociology,
Chinese Academy of Social Sciences
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/
Could you let us know what operating system you are using? It looks like
you need to rebuild a few packages: MASS and boot. Try
install.packages("MASS")
install.packages("boot")
And then reload Zelig.
Hope this helps,
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://www.princeton.edu/~kimai
---------------------------------------------------------
On Thu, 10 Mar 2005, Mitch Morehart wrote:
> Dr. King,
>
> Sorry to bother you directly, but I am having trouble accessing the
> zelig listserv. I am trying to estimate the marginal effects for a
> multinomial logit model and the Zelig library provides some very useful
> functions towards that end. I installed the latest version of the
> library along with the latest version of the VGAM library and tried to
> run the example code for multinomial logit. The details of the R
> session and error for the sim procedure are provided below:
>
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 2.0.0 (2004-10-04), ISBN 3-900051-07-0
>
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type 'license()' or 'licence()' for distribution details.
>
> R is a collaborative project with many contributors.
> Type 'contributors()' for more information and
> 'citation()' on how to cite R or R packages in publications.
>
> Type 'demo()' for some demos, 'help()' for on-line help, or
> 'help.start()' for a HTML browser interface to help.
> Type 'q()' to quit R.
>
>
> Welcome to R!
>
> > library(Zelig)
> Loading required package: MASS
> Loading required package: boot
>
> Please refer to http://gking.harvard.edu/zelig for full documentation
> or help.zelig() for help with commands and models supported by Zelig.
>
> Warning messages:
> 1: package 'MASS' was built under R version 2.0.1
> 2: package 'boot' was built under R version 2.0.1
> 3: S3 method(s) 'plot.zelig.beta' were declared in NAMESPACE but not
> found
> > data(mexico)
> > z.out <- zelig(vote88 ~ pristr + othcok + othsocok, model =
> "mlogit",
> + data = mexico)
> Loading required package: VGAM
> Loading required package: splines
> Loading required package: stats4
> Creating a new generic function for "print" in "VGAM"
> Creating a new generic function for "coefficients" in
> "VGAM"
> Creating a new generic function for "fitted.values" in
> "VGAM"
>
> Attaching package 'VGAM':
>
>
> The following object(s) are masked from package:splines :
>
> bs ns
>
>
> The following object(s) are masked from package:boot :
>
> logit simplex
>
>
> The following object(s) are masked from package:stats :
>
> binomial Gamma gaussian inverse.gaussian poisson poly
> quasibinomial quasipoisson
>
>
> The following object(s) are masked from package:base :
>
> is.R jitter scale.default
>
> > x.weak <- setx(z.out, pristr = 1)
> > x.strong <- setx(z.out, pristr = 3)
> > s.out <- sim(z.out, x = x.strong, x1 = x.weak)
> Error in array("", c(n, 3), list(names(object), c("Length", "Class",
> "Mode"))) :
> length of dimnames [1] not equal to array extent
> >
>
> Any thoughts on a solution would be appreciated. I looked at the
> source code and my only thought is that there is some incompatability
> betwen the latest version of vglm and the zelig sim procedure.
>
>
> Mitchell Morehart
> Room 4121 N
> 1800 M Street NW
> Washington, DC 20036
> morehart(a)ers.usda.gov
>
> Be sure to visit the ERS website
> www.ers.usda.gov
>
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)latte.harvard.edu
[un]subscribe Options: http://lists.hmdc.harvard.edu/?info=zelig
Zelig program information: http://gking.harvard.edu/zelig/