Bookmark and Share

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: RE: Can't store estimates after logistic regression


From   Steve Samuels <[email protected]>
To   [email protected]
Subject   Re: st: RE: Can't store estimates after logistic regression
Date   Fri, 28 Jun 2013 18:54:12 -0400

Eduoard:

If you do

. des _est_mod1

you will see the variable label: "esample() from estimates store". 
esample() just records whether the observation was used in the analysis

It  not the kind of variable needed by -roccomp-, which has syntax:

. rocommp refvar classvar

According to the Manual:

"The two variables refvar and classvar must be numeric. The reference
variable indicates the true state of the observation, such as diseased
and nondiseased or normal and abnormal, and must be coded as 0 and 1.
The rating or outcome of the diagnostic test or test modality is
recorded in classvar, which must be at least ordinal, with higher values
indicating higher risk."

In short -estimates- is irrelevant for running -roccomp-. I suggest that you
reread the manual entry for -roccomp-.

Steve

On Jun 28, 2013, at 5:54 PM, Edouard Betsem wrote:

It's still not working.
by is not allowed here and using separate returned this error code

[__00000E] not found
r(111);

I'm wondering if "roccomp" can be used with "estimates"

I've even renamed the _est_mod1 to mod1. It's on the editor table but I don't know why Stata can't find it

           logistic tma whoall
           est store mod1
           rename _est_mod1 mod1

           logistic tma study_2
           est store mod2
           rename _est_mod2 mod2
           roccomp tma _est_mod1 _est_mod2, graph summary

the estimates variables return the value "1" for all my observations.
Wish someone can help me solve this.
Thanks

Edouard BETSEM, MD, PhD


On 28/06/2013 00:48, Bley N'Dede wrote:
> I am not sure but may be you should try something like
> 
> roccomp _est_mod1 _est_mod2, graph summary by(dis)
> ________________________________________
> From: [email protected] [[email protected]] on behalf of Edouard Betsem [[email protected] should try m]
> Sent: Friday, June 28, 2013 7:29 AM
> To: [email protected]
> Subject: st: Can't store estimates
> 
> Hi,
> I'm trying to use roccomp commande to compare AUC for two logistic
> models. but stata won't store my estimates.
> 
> here are the commande lines.
> 
> * all counts with binary
>              logistic dis ref
>              est store mod1
> 
>              logistic dis stud
>              est store mod2
>              roccomp dis mod1 mod2, graph summary
> 
> and this is the message i have:
> 
> variable mod1 not found
> r(111);
> 
> does anyone know the reason for this?
> Best
> Edouard
> 
> 
> 
> 
> --
> Edouard BETSEM, MD, PhD
> 
> 
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/faqs/resources/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/faqs/resources/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/faqs/resources/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/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index