Hi all,
I do not have a licence for stata to combine the imputations for amelia. I have thought about using an aggregate function in SPSS to generate means and standard deviations for cases across all the five files (My gut instinct is that this is way too simple). I have also thought about using NORM to combine the five file parameters (Again seems to simple). Any advice on how to combine the amelia data files without stata would be much appreciated.
Kind regards Paul
Dear Prof. Scheve,
Thanks a lot for replying to my message.
I got rid of those variables for which there was perfect linear dependence (there were 2 of them), and run Amelia but it did not work, the same message kept coming up:
Iteration Number:
sweep: elements of m cannot be zero
Currently active call: SWEEP [79]
When I tried to change the default start values it gives me either the option of a default with covariances zero, or file with starting values (I wanted to write 2 as you suggested me), however it asks me to enter the name of a file which contains the starting values and hit return, am I supposed to create such a file? is it a matrix?
WHat I did was:
I set AMpri equal to 2 as it is advised for high missingness data
I described my fully variables
I left AMmtho under the default (should I change it to 2?)
I did not define any idvar (I have tried with and without and it does not make any difference)
I did not define any time of sectional variable in the AMts and AMcs (or the AMusets, AMusecs and AMtstep)
I did not ask for any nominal variable to be returned as such
I asked for 1 ordinal variable to be returned as such
I got the same message again...
Sorry for taking so much time out of you and again, thank you for helping me put here,
Yours
Jose Merino
----- Original Message -----
From: Kenneth Scheve <kenneth.scheve(a)yale.edu>
Date: Saturday, February 28, 2004 4:03 pm
Subject: Re: [amelia] run but
>
> Actually this message is not about which variables are fully
> observed. It
> usually indicates there is something wrong with the data or at
> least that
> portion of the dataset being used in the initial calculations.
>
> Did you try setting the AMstart global to "2" rather than the default?
>
> Also check and make sure that you have not included data for which
> there is
> perfect linear dependence between some group of variables.
>
> Ken
>
> At 03:53 PM 2/28/2004, Jose Antonio Pena Merino wrote:
> >OK, the program is running now, I think I was asking way too much
> by
> >including too many variables and asking it to keep all the
> ordinal and
> >nominal as such.
> >
> >However, even if it is running it keeps returning the following
> message:>
> >Iteration Number:
> >sweep: elements of m cannot be zero
> >Currently active call: SWEEP [79]
> >
> >I define those variables which contain full information already,
> but it
> >looks like Amelia keeps detecting one variable with full
> information. In
> >order to detect which one was it I kept adding each variable to
> the list
> >of fullys, but all of them were detected as non-fully. So I
> really do not
> >know why I keep getting the same message. In sum, there is no
> other
> >variable besides those already defined as fully for which m=0
> (that is,
> >with no missing values)
> >
> >Thanks for clarifying this for me
> >
> >Yours
> >
> >Jose Merino
>
>
OK, the program is running now, I think I was asking way too much by including too many variables and asking it to keep all the ordinal and nominal as such.
However, even if it is running it keeps returning the following message:
Iteration Number:
sweep: elements of m cannot be zero
Currently active call: SWEEP [79]
I define those variables which contain full information already, but it looks like Amelia keeps detecting one variable with full information. In order to detect which one was it I kept adding each variable to the list of fullys, but all of them were detected as non-fully. So I really do not know why I keep getting the same message. In sum, there is no other variable besides those already defined as fully for which m=0 (that is, with no missing values)
Thanks for clarifying this for me
Yours
Jose Merino
Thanks so much for your help! I was hoping to get around the dummy variable
solution but it looks we will just have to use stata's xi with miest
somehow save the output.
Patrick
At 08:32 27/02/2004 -0500, Gary King wrote:
>you could use clarify or miest with regression and a set of dummy
>variables to predict the categories. the coefficients will return the
>fractions in each cell and the se's will be the se's.
>Gary
>
> : Gary King, King(a)Harvard.Edu http://GKing.Harvard.Edu :
> : Center for Basic Research Direct (617) 495-2027 :
> : in the Social Sciences Assistant (617) 495-9271 :
> : 34 Kirkland Street, Rm. 2 HU-MIT DC (617) 495-4734 :
> : Harvard U, Cambridge, MA 02138 eFax (928) 832-7022 :
>
>On Wed, 25 Feb 2004, Patrick Bernhagen wrote:
>
> > Hi,
> >
> > I am trying to obtain a table displaying proportions of a categorical
> > variable after multiple imputation using Amelia. I guess what I am looking
> > for is something equivalent to Stata's svyprop command that gives us
> > proportions with standard errors. Expressed differently, I would like to
> > get tables instead of summary statistics (as in misum) within the MI
> > command family, or be able to request something as simple as 'table' or
> > 'tab' etc., instead of the various more sophisticated regression
> estimators
> > offered for 'miest'.
> >
> > Many thanks in advance for your help!
> >
> > Patrick
> >
************************************************
Patrick Bernhagen M.A.
University of Dublin
Department of Political Science
Trinity College
Dublin 2, IRELAND
Tel.: +353-1-608-1651
Fax: +353-1-677-0546
mailto:bernhagp@tcd.ie
http://www.tcd.ie/Political_Science/P.Bernhagen
Hi,
I am trying to obtain a table displaying proportions of a categorical
variable after multiple imputation using Amelia. I guess what I am looking
for is something equivalent to Stata's svyprop command that gives us
proportions with standard errors. Expressed differently, I would like to
get tables instead of summary statistics (as in misum) within the MI
command family, or be able to request something as simple as 'table' or
'tab' etc., instead of the various more sophisticated regression estimators
offered for 'miest'.
Many thanks in advance for your help!
Patrick
I am currently working on a paper (with the help of Professor Nathaniel Beck)I've been trying to run Amelia for my data base, but it keeps crashing and giving me the following expression:
C:\GAUSS\SRC\SORTMC.SRC(89) : error G0057 : Procedure stack overflow - expression too complex
Currently active call: SORTMC [89]
By taking a look at the Amelia page, I took note of what others have done but things do not seem to improve here.
My data is composed of 66 variables with a lot of missings as well as ordinal and nominal variables, which tend to be those for which there is full information. After trying that data base, I decided to reduce it to 29 variables (I do not think I can go any lower). The data is mostly from the World Bank, but it also includes some variables from Przeworski's data and Profesor Bueno de Mesquita's data base. I suspect the problem is the huge amount of missingness (total of obs is supposed to be 7500 and for some variables there are only 1500-2000 obs). I really do not know how to proceed from here. If you have any suggestion, please do not hesitate to let me know.
Beforehand, I want to say thank you for the help you can provide me.
Yours
Jose Merino
Ph.D. Student
I think this is a really dumb question, but here goes:
In the past when using Amelia I have always computed parameter estimates and
standard errors across the five datasets by hand. This time I have too
many, and I want to do it in Stata. I just can't figure out how to tell
Stata which datasets to use. Giving the stub name doesn't seem to be enough
to get it to call up the right datasets from the right directory on my
computer if they're not open in Stata. Do I have to tell Stata to open all
of the datasets at once and then give the miest command syntax? If so, how
do you keep multiple datasets in memory at the same time?
Julie
Julia Lynch
Robert Wood Johnson Health Policy Scholar
Harvard University
CBRSS
34 Kirkland St
Cambridge MA 02138
tel 617 495 5366
fax 617 496 1636
email jlynch(a)latte.harvard.edu
and
Assistant Professor
Dept. of Political Science
University of Pennsylvania
ON LEAVE 2003-2005
Hi,
I downloaded the self-extracting executable for Windows from the website,
unpacked it into an Amelia directory in the C: root and got no error messages
on the unpack. However, when I invoke AMELIA.BAT I get an empty dos-type
window, and a program icon on my task bar that says "GAUSS" - but when I click
on that nothing happens.
This is on a Windows XP system with a very fast processor and over a gig of
memory. Any suggestions on how to get this program launching?
Many thanks,
Lori
yup, that was easy... serves me right for not using Stata more often...
Julia Lynch
Robert Wood Johnson Health Policy Scholar
Harvard University
CBRSS
34 Kirkland St
Cambridge MA 02138
tel 617 495 5366
fax 617 496 1636
email jlynch(a)latte.harvard.edu
and
Assistant Professor
Dept. of Political Science
University of Pennsylvania
ON LEAVE 2003-2005