I tried the model
xi: glm lengthof stay lapscopy i.appdgn2 i.alderk prepermalign
precardioscleros preperdiabetes cons, eform link(power-1) nocons
and get an error message "unrecognized command: power"
Roland
2008/11/5 Nick Cox <[email protected]>:
> I doubt that anything is wrong with Tony's model except that
> -eform("exp(b)")- should just be -eform-.
>
> Nick
>
> roland andersson
>
> Peter and Maarten
>
> I am sorry Peter. Your model is not accepted by Stata. I tried
> different alternativ without success.
>
> However I tried Maarten GLM model
>
> xi: glm studytime i.drug c_age cons, family(gaussian) link(log) nocons
> eform
>
> on my data and got a different result compared to the regress of the
> lnLOS. Now laparoscopy has shorter LOS. Which method is correct?
>
>
> 2008/11/5 Lachenbruch, Peter <[email protected]>:
>
>> The issue seems to be that hospitals have a closure date on stay when
>> you are doing a study after patients are certain (or almost certain)
> to
>> have been discharged (e.g., all records are from admissions at least a
>> year old).
>>
>> An alternative model might fit the reciprocal of the mean rather than
>> the log of the observations (thus obviating problems with 0 days of
> stay
>> - e.g. an outpatient visit to the ER) in this case you could use
>> generalized linear models to get
>> xi: glm LOS lapscopic i.appdgn age agesq cons, eform("exp(b)")
>> link(power -1) nocons
>>
>>
>> Tony
>
> *
> * 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/
