19 Jan
2010
19 Jan
'10
6:12 a.m.
Hi Matt
Yes, that was it.
Thanks
Peter
-----Original Message-----
From: Matt Blackwell <blackwel(a)fas.harvard.edu>
Sent: Jan 18, 2010 7:02 PM
To: Peter Flom <peterflomconsulting(a)mindspring.com>
Cc: "amelia(a)lists.gking.harvard.edu" <amelia(a)lists.gking.harvard.edu>
Subject: Re: [amelia] tscsPlot error
Hi Peter,
Quick question. Do any of these cross sections contain completely
missing data on a variable? Or do they only have one observation?
These situations can lead to problems when you use the "intercs"
option.
Cheers,
matt.
On Fri, Jan 15, 2010 at 8:55 AM, Peter Flom
<peterflomconsulting(a)mindspring.com> wrote:
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
-
Amelia mailing list served by Harvard-MIT Data Center
[Un]Subscribe/View Archive: http://lists.gking.harvard.edu/?info=amelia
More info about Amelia: http://gking.harvard.edu/amelia
Good morning,
I am still working with the data set described in my previous e-mails. Definitely making
progress by adding more info. But I now get an error that I do not understand
If I run:
<<<
susanMI.out <- amelia(susan2, m = 5, noms = "married", ts =
'time', cs = 'id',
intercs = T,
sqrts = "unprot_vag_sex")
-
Amelia mailing list served by Harvard-MIT Data Center
[Un]Subscribe/View Archive: http://lists.gking.harvard.edu/?info=amelia
More info about Amelia: http://gking.harvard.edu/amelia >>
and then (e.g.)
<<
tscsPlot(susanMI.out, cs = 11, var = "unprot_vag_sex")
>>
I get an error
Error in identical(m, vector(mode = "logical", length = length(m))) :
subscript out of bounds
This occurs with some values of cs but not others. I checked that those id numbers do,
in fact, exist in the data.
Any help appreciated.
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
-
Amelia mailing list served by Harvard-MIT Data Center
[Un]Subscribe/View Archive: http://lists.gking.harvard.edu/?info=amelia
More info about Amelia: http://gking.harvard.edu/amelia