Hi Mitch and others who had a problem with VGAM earlier,
We just released a new version, 2.1-1, which fixes this problem. It will
be submitted to CRAN later today, but you can also try
install.packages("Zelig", CRAN="http://gking.harvard.edu")
to install it from our Zelig website. Your VGAM needs to be also updated:
the new version is 0.6-2.
Please let us know if you encounter any problem.
Thanks,
Zelig team
On Mon, 21 Mar 2005, Mitch Morehart wrote:
> Here is the same procedure under R 2.0.1. It does not seem to make a
> difference.
>
>
> R : Copyright 2004, The R Foundation for Statistical Computing
> Version 2.0.1 (2004-11-15), 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.
>
> > 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.
>
> > packageDescription("VGAM")
> Package: VGAM
> Version: 0.6-1
> Date: 2005-01-31
> Title: Vector Generalized Linear and Additive Models
> Author: Thomas W. Yee <t.yee(a)auckland.ac.nz>
> Maintainer: Thomas Yee <t.yee(a)auckland.ac.nz>
> Depends: R (>= 2.0.0), splines, methods, stats, stats4
> Description: Vector Generalized Linear and Additive Models, plus
> Reduced-Rank VGLMs and Quadratic RR-VGLMs
> License: Copyright TW Yee (2005)
> URL: http://www.stat.auckland.ac.nz/~yee
> LazyLoad: yes
> LazyData: yes
> Packaged: Tue Feb 1 13:40:14 2005; yee
> Built: R 2.0.0; i386-pc-mingw32; 2005-02-01 13:45:43; windows
>
> -- File: D:/R/rw2001/library/VGAM/DESCRIPTION
> > packageDescription("Zelig")
> Package: Zelig
> Version: 2.0-13
> Date: 2005-3-11
> Title: Zelig: Everyone's Statistical Software
> Author: Kosuke Imai <kimai(a)Princeton.Edu>, Gary King
> <king(a)harvard.edu>, Olivia Lau <olau(a)fas.harvard.edu>
> Maintainer: Olivia Lau <olau(a)fas.harvard.edu>
> Depends: R (>= 1.9.0), MASS, boot
> Suggests: VGAM
> Description: Zelig is an easy-to-use program that can estimate, and
> help interpret the results of, an enormous range of
> statistical models. It literally is ``everyone's statistical
> software'' because Zelig's simple unified framework
> incorporates everyone else's (R) code. We also hope it will
> become ``everyone's statistical software'' for applications
> and teaching, and so have designed Zelig so that anyone can
> easily use it or add their programs to it. Zelig also comes
> with infrastructure that facilitates the use of any existing
> method, such as by allowing multiply imputed data for any
> model, and mimicking the program Clarify (for Stata) that
> takes the raw output of existing statistical procedures and
> translates them into quantities of direct interest.
> License: GPL version 2 or newer
> URL: http://gking.harvard.edu/zelig
> Packaged: Fri Mar 11 16:47:33 2005; olau
> Built: R 2.0.1; ; 2005-03-14 12:24:09; windows
>
> -- File: D:/R/rw2001/library/Zelig/DESCRIPTION
> > packageDescription("boot")
> Package: boot
> Priority: recommended
> Version: 1.2-20
> Date: 2004-11-03
> Author: S original <http://statwww.epfl.ch/davison/BMA/library.html>
> by Angelo Canty <cantya(a)mcmaster.ca>. R port by Brian Ripley
> <ripley(a)stats.ox.ac.uk>.
> Maintainer: Brian Ripley <ripley(a)stats.ox.ac.uk>
> Description: functions and datasets for bootstrapping from the book
> "Bootstrap Methods and Their Applications" by A. C. Davison
> and D. V. Hinkley (1997, CUP).
> Title: Bootstrap R (S-Plus) Functions (Canty)
> Depends: R (>= 2.0.0), graphics, stats
> Suggests: survival
> LazyData: yes
> License: Unlimited distribution.
> Packaged: Wed Nov 3 07:58:39 2004; ripley
> Built: R 2.0.1; ; 2004-11-15 13:38:21; windows
>
> -- File: D:/R/rw2001/library/boot/DESCRIPTION
> > packageDescription("MASS")
> Package: MASS
> Description: The main library and the datasets
> Title: Main Package of Venables and Ripley's MASS
> LazyLoad: yes
> LazyData: yes
> Bundle: VR
> Priority: recommended
> Version: 7.2-10
> Date: 2004-11-03
> Depends: R (>= 2.0.0), graphics, stats
> Suggests: lattice, nlme, survival
> Author: S original by Venables & Ripley. R port by Brian Ripley
> <ripley(a)stats.ox.ac.uk>, following earlier work by Kurt Hornik
> and Albrecht Gebhardt.
> Maintainer: Brian Ripley <ripley(a)stats.ox.ac.uk>
> BundleDescription: Functions and datasets to support Venables and
> Ripley, 'Modern Applied Statistics with S' (4th edition).
> License: GPL (version 2 or later) See file LICENCE.
> URL: http://www.stats.ox.ac.uk/pub/MASS4/
> Packaged: Wed Nov 3 09:02:23 2004; ripley
> Built: R 2.0.1; i386-pc-mingw32; 2004-11-15 13:37:23; windows
>
> -- File: D:/R/rw2001/library/MASS/DESCRIPTION
> > 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
> >
> >
>
>
> 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
>
> >>> Kosuke Imai <kimai(a)Princeton.Edu> 3/21/2005 11:17:06 AM >>>
> I'm looking at your output.
>
> > 1: package 'MASS' was built under R version 2.0.1
> > 2: package 'boot' was built under R version 2.0.1
>
> This means that MASS and boot need to be reinstalled. I'm not sure if
> this
> will fix the problem, but it's still a good idea to fix these errors.
>
> Kosuke
>
> On Mon, 21 Mar 2005, Mitch Morehart wrote:
>
> > Here is the description for MASS and boot
> >
> > packageDescription("MASS")
> > Package: MASS
> > Description: The main library and the datasets
> > Title: Main Package of Venables and Ripley's MASS
> > LazyLoad: yes
> > LazyData: yes
> > Bundle: VR
> > Priority: recommended
> > Version: 7.2-12
> > Date: 2005-01-31
> > Depends: R (>= 2.0.0), graphics, stats
> > Suggests: lattice, nlme, survival
> > Author: S original by Venables & Ripley. R port by Brian Ripley
> > <ripley(a)stats.ox.ac.uk>, following earlier work by Kurt
> Hornik
> > and Albrecht Gebhardt.
> > Maintainer: Brian Ripley <ripley(a)stats.ox.ac.uk>
> > BundleDescription: Functions and datasets to support Venables and
> > Ripley, 'Modern Applied Statistics with S' (4th edition).
> > License: GPL (version 2 or later) See file LICENCE.
> > URL: http://www.stats.ox.ac.uk/pub/MASS4/
> > Packaged: Tue Feb 1 09:45:34 2005; ripley
> > Built: R 2.0.1; i386-pc-mingw32; 2005-02-02 12:27:57; windows
> >
> > -- File: D:/R/rw2000/library/MASS/DESCRIPTION
> >
> > packageDescription("boot")
> > Package: boot
> > Priority: recommended
> > Version: 1.2-22
> > Date: 2005-02-01
> > Author: S original <http://statwww.epfl.ch/davison/BMA/library.html>
> > by Angelo Canty <cantya(a)mcmaster.ca>. R port by Brian Ripley
> > <ripley(a)stats.ox.ac.uk>.
> > Maintainer: Brian Ripley <ripley(a)stats.ox.ac.uk>
> > Description: functions and datasets for bootstrapping from the book
> > "Bootstrap Methods and Their Applications" by A. C. Davison
> > and D. V. Hinkley (1997, CUP).
> > Title: Bootstrap R (S-Plus) Functions (Canty)
> > Depends: R (>= 2.0.0), graphics, stats
> > Suggests: survival
> > LazyData: yes
> > License: Unlimited distribution.
> > Packaged: Tue Feb 1 10:04:32 2005; ripley
> > Built: R 2.0.1; ; 2005-02-02 12:26:26; windows
> >
> > -- File: D:/R/rw2000/library/boot/DESCRIPTION
> >
> >
> >
> > 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
> >
> > >>> Kosuke Imai <kimai(a)Princeton.Edu> 3/21/2005 11:11:37 AM >>>
> > please update MASS and boot first.
> >
> > Kosuke
> >
> > On Mon, 21 Mar 2005, Mitch Morehart wrote:
> >
> > > I have the latest versions of both VGAM and Zelig installed. The
> > > example seems to fail on the "sim" command. Here is my output:
> > >
> > >
> > > 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.
> > >
> > > > 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
> > > > packageDescription("VGAM")
> > > Package: VGAM
> > > Version: 0.6-1
> > > Date: 2005-01-31
> > > Title: Vector Generalized Linear and Additive Models
> > > Author: Thomas W. Yee <t.yee(a)auckland.ac.nz>
> > > Maintainer: Thomas Yee <t.yee(a)auckland.ac.nz>
> > > Depends: R (>= 2.0.0), splines, methods, stats, stats4
> > > Description: Vector Generalized Linear and Additive Models, plus
> > > Reduced-Rank VGLMs and Quadratic RR-VGLMs
> > > License: Copyright TW Yee (2005)
> > > URL: http://www.stat.auckland.ac.nz/~yee
> > > LazyLoad: yes
> > > LazyData: yes
> > > Packaged: Tue Feb 1 13:40:14 2005; yee
> > > Built: R 2.0.0; i386-pc-mingw32; 2005-02-01 13:45:43; windows
> > >
> > > -- File: D:/R/rw2000/library/VGAM/DESCRIPTION
> > > > packageDescription("Zelig")
> > > Package: Zelig
> > > Version: 2.0-13
> > > Date: 2005-03-11
> > > Title: Zelig: Everyone's Statistical Software
> > > Author: Kosuke Imai <kimai(a)Princeton.Edu>, Gary King
> > > <king(a)harvard.edu>, Olivia Lau <olau(a)fas.harvard.edu>
> > > Maintainer: Olivia Lau <olau(a)fas.harvard.edu>
> > > Depends: R (>= 1.9.0), MASS, boot
> > > Description: Zelig is an easy-to-use program that can estimate,
> and
> > > help interpret the results of, an enormous range of
> > > statistical models. It literally is ``everyone's
> statistical
> > > software'' because Zelig's simple unified framework
> > > incorporates everyone else's (R) code. We also hope it
> will
> > > become ``everyone's statistical software'' for
> applications
> > > and teaching, and so have designed Zelig so that anyone
> can
> > > easily use it or add their programs to it. Zelig also
> comes
> > > with infrastructure that facilitates the use of any
> existing
> > > method, such as by allowing multiply imputed data for any
> > > model, and mimicking the program Clarify (for Stata) that
> > > takes the raw output of existing statistical procedures
> and
> > > translates them into quantities of direct interest.
> > > License: GPL version 2 or newer
> > > URL: http://gking.harvard.edu/zelig
> > > Packaged: Fri Mar 11 11:04:06 2005; king
> > > Built: R 2.0.0; ; 2005-03-11 11:04:08; windows
> > >
> > > -- File: D:/R/rw2000/library/Zelig/DESCRIPTION
> > > > 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
> > > >
> > > >
> > >
> > >
> > > 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
> > >
> > > >>> Kosuke Imai <kimai(a)Princeton.Edu> 3/21/2005 9:34:41 AM >>>
> > > Hi,
> > > These commands work for me. Please make sure you have the most
> > recent
> > >
> > > version of Zelig 2.0-13 and VGAM 0.6-1. For example, you can use
> > > packageDescription("VGAM") to get this info.
> > > Kosuke
> > >
> > > On Mon, 14 Mar 2005, Mitch Morehart wrote:
> > >
> > > > Traceback is included at the bottom:
> > > >
> > > >
> > > > 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(MASS)
> > > > Warning message:
> > > > package 'MASS' was built under R version 2.0.1
> > > > > library(boot)
> > > > Warning message:
> > > > package 'boot' was built under R version 2.0.1
> > > > >
> > > > > library(Zelig)
> > > >
> > > > Please refer to http://gking.harvard.edu/zelig for full
> > documentation
> > >
> > > > or help.zelig() for help with commands and models supported by
> > > Zelig.
> > > >
> > > > Warning message:
> > > > 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
> > > > > traceback()
> > > > 10: array("", c(n, 3), list(names(object), c("Length", "Class",
> > > > "Mode")))
> > > > 9: summary.default(object, corr = FALSE, dispersion =
> dispersion)
> > > > 8: summary(object, corr = FALSE, dispersion = dispersion)
> > > > 7: vcovvlm(object, ...)
> > > > 6: vcov(object)
> > > > 5: vcov(object)
> > > > 4: param.vglm(object, num = num, bootstrap = bootstrap)
> > > > 3: param(object, num = num, bootstrap = bootstrap)
> > > > 2: sim.default(z.out, x = x.strong, x1 = x.weak)
> > > > 1: sim(z.out, x = x.strong, x1 = x.weak)
> > > >
> > > >
> > >
> >
> ------------------------------------------------------------------------------------------------
> > > >
> > > >
> > > >
> > > > 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
> > > >
> > > > >>> Kosuke Imai <kimai(a)Princeton.Edu> 3/12/2005 10:51:23 AM >>>
> > > > could you send us the output of traceback() after you get the
> > error?
> > > > Thanks,
> > > > Kosuke
> > > >
> > > > On Fri, 11 Mar 2005, Mitch Morehart wrote:
> > > >
> > > > > I am running R 2.0 in Windows XP. I think the issue is with
> > the
> > > > new
> > > > > treatment of coefficients in VGAM.
> > > > >
> > > > > 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"
> > > > >
> > > > >
> > > >
> > >
> >
> --------------------------------------------------------------------------------------------
> > > > >
> > > > >
> > > > > MASS is part of the VR package which I updated from CRAN
> > > > >
> > > > > trying URL
> > > > >
> > `http://cran.r-project.org/bin/windows/contrib/2.0/VR_7.2-12.zip'
> > > > > Content type `application/zip' length 1551825 bytes
> > > > > opened URL
> > > > > downloaded 1515Kb
> > > > >
> > > > > bundle 'VR' successfully unpacked and MD5 sums checked
> > > > >
> > > > >
> > > > > install.packages("boot")
> > > > > trying URL
> > > > > `http://cran.r-project.org/bin/windows/contrib/2.0/PACKAGES'
> > > > > Content type `text/plain; charset=iso-8859-1' length 26875
> bytes
> > > > > opened URL
> > > > > downloaded 26Kb
> > > > >
> > > > > trying URL
> > > > >
> > >
> `http://cran.r-project.org/bin/windows/contrib/2.0/boot_1.2-22.zip'
> > > > > Content type `application/zip' length 767540 bytes
> > > > > opened URL
> > > > > downloaded 749Kb
> > > > >
> > > > > package 'boot' successfully unpacked and MD5 sums checked
> > > > >
> > > > > Here is the output of another attempt:
> > > > >
> > > > > > library(MASS)
> > > > > > library(boot)
> > > > > >
> > > > > > library(Zelig)
> > > > > > data(mexico)
> > > > > > z.out <- zelig(vote88 ~ pristr + othcok + othsocok, model =
> > > > > "mlogit",
> > > > > + data = mexico)
> > > > > > 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
> > > > >
> > > > >
> > > > >
> > > > > 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
> > > > >
> > > > > >>> Kosuke Imai <kimai(a)Princeton.Edu> 3/10/2005 9:21:19 PM >>>
> > > > > 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/