Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: ICE with Interval Error #459 multiple imputation
From
A Loumiotis <[email protected]>
To
[email protected]
Subject
Re: st: ICE with Interval Error #459 multiple imputation
Date
Mon, 13 Aug 2012 09:49:13 +0300
Perhaps you could also try to use ologit instead of intreg to impute
sclifeca2. Make sure you use the o. prefix, ie,
ice o.sclifeca2 ..., saving(imputed_dataset) m(5)
On Mon, Aug 13, 2012 at 9:20 AM, A Loumiotis
<[email protected]> wrote:
> Hi Laura,
>
> I'm not sure if you are setting up your lower and upper variables
> correctly because you said that when the observation is not missing
> both upper and lower bounds take the value of 7. This is not correct.
> They must be equal to the non-missing value. Here is an example of
> how they should look:
>
> input sclifeca2 ll1 ul7
> 4 4 4
> . 1 7
> 2.5 2.5 2.5
> 3 3 3
> . 1 7
> end
>
> Best,
> Antonis
>
>
>
>
> On Fri, Aug 10, 2012 at 9:32 PM, Laura Kudrna <[email protected]> wrote:
>> Hi,
>>
>> Thank you for your response. Apologies about ICE, I'm new to the
>> listserve: ICE comes from //www.homepages.ucl.ac.uk/~ucakjpr/stata/
>>
>> To answer your question about the bounds, yes: when I created the
>> upper limit variable ul7 that has no missing values and I include it
>> in the ice code, the upper and lower bounds remain 7. But for
>> example, when I impute my variable sclifeca2 "Life Satisfaction"
>> which has a scale of 1 to 7 I get imputed values ranging from -0.5 to
>> 9. I am using multiple imputation to impute all the predictor and
>> explanatory variables at once.
>>
>> I replicated the STATA code precisely from the ICE help for interval.
>>
>> ICE help for interval states: ...intlist has the syntax varname:llvar
>> ulvar [, varname:llvar ulvar ...], where each varname is an
>> interval-censored variable, each llvar contains the lower bound (a)
>> for varname and each ulvar contains the upper bound (b)
>>
>> My syntax is: ice ll1 ul7 sclifeca2 .....(other variables).....,
>> saving(imputed_dataset) m(5) interval(sclifeca2:ll1 ul7).
>>
>> I still can't understand why I am getting the error message and I
>> would be most grateful if you have any other suggestions on how to fix
>> it.
>>
>> Best wishes,
>> Laura
>>
>>
>> On Fri, Aug 10, 2012 at 7:32 AM, A Loumiotis
>> <[email protected]> wrote:
>>> If the variable that you are imputing with interval regression is not
>>> missing then the lower and upper bounds should equal that non missing
>>> value. Is this true in your case? See -help ice- especially the
>>> interval option.
>>>
>>> -ice- is a user written program available from SJ-9-3. Please state
>>> where user written programs come from as explained in Statalist FAQ.
>>>
>>>
>>>
>>>
>>> On Fri, Aug 10, 2012 at 12:29 AM, Laura Kudrna <[email protected]> wrote:
>>>> Hi,
>>>>
>>>> I am having the same error message / problem as the post below from
>>>> the STATA message boards - Error #459. It occurs when using the
>>>> interval option with ICE.
>>>>
>>>> I generated two variables, ll1 and ul7, containing values of 1 and 7
>>>> respectively to set the upper and lower bound for my variable
>>>> sclifeca2.
>>>>
>>>> My STATA code is:
>>>>
>>>> ice ll1 ul7 sclifeca2 sclifeca3 sclifeca4 sclddr2 sclddr3 sclddr4
>>>> iincome2 iincome3 iincome4 education2 occupation2 unemployed2
>>>> unemployed3 sex2 age2 age3 age4 partner2 partner3 partner4 vehicles2
>>>> vehicles3 vehicles4 holiday2 holiday3 holiday4 secprop2 secprop3
>>>> secprop4 weight2 weight3 weight4, saving(imputed_dataset) m(5)
>>>> interval(sclifeca2:ll1 ul7)
>>>>
>>>> My error message is:
>>>> Error #459 encountered while running -uvis-
>>>> I detected a problem with running uvis with command intreg on response sclifeca2
>>>> and covariates sclifeca3 sclifeca4 sclddr2 sclddr3 sclddr4 iincome2
>>>> iincome3 iincome4 education2 occupation2 unemployed2 unemployed3 sex2
>>>> ag
>>>>> e2 age3 age4 partner2 partner3 partner4 vehicles2 vehicles3 vehicles4 holiday2 holiday3 holiday4 secprop2 secprop3 secprop4 weight2 weight
>>>>> 3 weight4.
>>>>
>>>> The offending command resembled:
>>>> uvis intreg sclifeca2 sclifeca3 sclifeca4 sclddr2 sclddr3 sclddr4
>>>> iincome2 iincome3 iincome4 education2 occupation2 unemployed2
>>>> unemployed3
>>>>> sex2 age2 age3 age4 partner2 partner3 partner4 vehicles2 vehicles3 vehicles4 holiday2 holiday3 holiday4 secprop2 secprop3 secprop4 weight2
>>>>> weight3 weight4 , gen([imputed])
>>>>
>>>> you may wish to try the -persist- option to persist beyond this error.
>>>> dumping current data to ./_ice_dump.dta
>>>> something that should be true of your data is not
>>>>
>>>> When I use the "persist" option, the data is imputed, but doesn't set
>>>> the upper and lower bounds of my variable.
>>>>
>>>> Any ideas?!
>>>>
>>>> Thanks a lot.
>>>>
>>>> Best,
>>>> Laura
>>>>
>>>> --------------------------------------------------------------------------------
>>>> [Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
>>>> st: ICE with intervall
>>>>
>>>> --------------------------------------------------------------------------------
>>>> From "Renger van Nieuwkoop" <[email protected]>
>>>> To <[email protected]>
>>>> Subject st: ICE with intervall
>>>> Date Wed, 27 Oct 2010 11:02:00 +0200
>>>>
>>>> --------------------------------------------------------------------------------
>>>>
>>>> Hi
>>>> I am trying to use ICE with interval imputation (the values are censored
>>>> and are not allowed to be below 0).
>>>> I have generated two variables co2_lb (lower bound) = 0 and co2_ub
>>>> (upper bound) = . and then tried
>>>> ice c01 c02 c02_lb c02_ub inc_earned, m(5) saving(imputed)
>>>> interval(c02:c02_lb c02_ub) replace
>>>> Without interval I have no problems with ICE, with interval I get the
>>>> following output (see below).
>>>> Any idea what I am doing wrongly?
>>>> Thanks
>>>> Renger
>>>>
>>>>
>>>>
>>>> . ice c01 c02 c02_lb c02_ub inc_earned, m(5) saving(imputed)
>>>> interval(c02:c02_lb c02_ub) replace
>>>>
>>>> #missing |
>>>> values | Freq. Percent Cum.
>>>> ------------+-----------------------------------
>>>> 0 | 25,181 85.37 85.37
>>>> 1 | 3,579 12.13 97.50
>>>> 2 | 735 2.49 100.00
>>>> 3 | 1 0.00 100.00
>>>> ------------+-----------------------------------
>>>> Total | 29,496 100.00
>>>>
>>>> Variable | Command | Prediction equation
>>>> ------------+---------+-------------------------------------------------
>>>> ------
>>>> c02_lb | | [Lower bound for c02]
>>>> c02_ub | | [Upper bound for c02]
>>>> c01 | regress | c02 inc_earned
>>>> inc_earned | regress | c01 c02
>>>> c02 | intreg | c01 inc_earned
>>>> ------------------------------------------------------------------------
>>>> ------
>>>>
>>>> Imputing [Note: in regression for c01, permanently removing c02 due to
>>>> collinearity]
>>>> [Note: in regression for inc_earned, permanently removing c02 due to
>>>> collinearity]
>>>>
>>>> Error #459 encountered while running -uvis-
>>>> I detected a problem with running uvis with command intreg on response
>>>> c02
>>>> and covariates c01 inc_earned.
>>>>
>>>> The offending command resembled:
>>>> uvis intreg c02 c01 inc_earned , gen([imputed])
>>>>
>>>> you may wish to try the -persist- option to persist beyond this error.
>>>> dumping current data to ./_ice_dump.dta
>>>> something that should be true of your data is not
>>>> r(459);
>>>> *
>>>> * For searches and help try:
>>>> * http://www.stata.com/help.cgi?search
>>>> * http://www.stata.com/support/statalist/faq
>>>> * http://www.ats.ucla.edu/stat/stata/
>>> *
>>> * For searches and help try:
>>> * http://www.stata.com/help.cgi?search
>>> * http://www.stata.com/support/statalist/faq
>>> * http://www.ats.ucla.edu/stat/stata/
>> *
>> * For searches and help try:
>> * http://www.stata.com/help.cgi?search
>> * http://www.stata.com/support/statalist/faq
>> * http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/