---------------------------------------------------------------------------------------------------
       log:  C:\mydocs\MPH\Text\SecondEdition\WebDoFiles\2.18.Funding.log
  log type:  text
 opened on:  28 Dec 2007, 10:55:27

. *  2.18.Funding.log
. *
. *  Explore the relationship between NIH research funds and disability-
. *  adjusted life-years lost due to the 29 diseases discussed by 
. *  Gross et al. (1999).  Look for transformed values of these variables 
. *  that are linearly related.  Perform a linear regression on these 
. *  transformed variables.  Replot this regression line as a function of
. *  the untransformed variables.
. *
. use C:\WDDtext\2.18.Funding.dta, clear

. scatter dollars disabil,  symbol(Oh) ylabel(0(.2)1.4, angle(0))         ///
>     xlabel(0(1)9) xsize(2.7) ysize(1.964) scale(1.5)

. more

. scatter dollars disabil if dollars <1                                   ///
>     , symbol(Oh) ylabel(0(.1).4) ymtick(.05(.1).35) xlabel(0(1)9)       ///
>     xsize(2.7) ysize(1.964) scale(1.5)

. more

. scatter dollars disabil if dollars <1                                   ///
>     , symbol(Oh) ylabel(0(.1).4, angle(0)) ymtick(.05(.1).35)           ///
>       xscale(log) xlabel(0.01 0.1 1 10)                                 ///
>       xmtick(.02 (.01) .09 .2 (.1) .9 2 (1) 9)                          ///
>       xsize(2.7) ysize(1.964) scale(1.5)

. more

. generate logdis = log(disabil)

. generate logdol = log(dollars)

. regress logdol logdis

      Source |       SS       df       MS              Number of obs =      29
-------------+------------------------------           F(  1,    27) =   18.97
       Model |  14.8027627     1  14.8027627           Prob > F      =  0.0002
    Residual |  21.0671978    27  .780266584           R-squared     =  0.4127
-------------+------------------------------           Adj R-squared =  0.3909
       Total |  35.8699605    28  1.28107002           Root MSE      =  .88333

------------------------------------------------------------------------------
      logdol |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
      logdis |   .4767575    .109458     4.36   0.000     .2521682    .7013468
       _cons |  -2.352205   .1640383   -14.34   0.000    -2.688784   -2.015626
------------------------------------------------------------------------------

. predict yhat,xb

. predict stdp, stdp

. generate ci_u = yhat +invttail(_N-2,.025)*stdp

. generate ci_l = yhat -invttail(_N-2,.025)*stdp

. sort logdis

. twoway rarea ci_u ci_l logdis, color(yellow)                            /// 
>     || line yhat logdis                                                 ///
>     || scatter logdol logdis,  symbol(Oh)                               ///
>     , ylabel(-4.61 "0.01" -2.3 "0.1" 0 "1", angle(0))                   ///
>       ymtick(-4.61 -3.91 -3.51 -3.22 -3.00 -2.81 -2.66 -2.53            ///
>       -2.41 -2.3 -1.61 -1.2 -.92 -.69 -.51 -.36 -.22 -.11 0)            ///
>       xtitle(Disability-Adjusted Life-Years Lost ($ millions))          ///
>       xlabel(-4.61 "0.01" -2.3 "0.1" 0 2.3 "10")                        ///
>       xmtick(-2.3 -1.61 -1.2 -.92 -.69 -.51 -.36 -.22 -.11              ///
>       0 .69 1.1 1.39 1.61 1.79 1.95 2.08 2.2 2.3) legend(off)           ///
>       xsize(2.7) ysize(1.964) scale(1.5)

. more            

. generate yhat2=exp(yhat)

. generate ci_u2= exp(ci_u)

. generate ci_l2= exp(ci_l)

. twoway  rarea ci_u2 ci_l2 disabil, color(yellow)                        ///
>     || line yhat2 disabil                                               ///
>     || scatter dollars disabil,  symbol(Oh)                             ///
>     , ytitle(NIH Research Funds ($ Billions))                           ///
>       ylabel(0(.2)1.4, angle(0)) ymtick(.1(.2)1.3) xlabel(0(1)9)        ///
>       legend(off) xsize(2.7) ysize(1.964) scale(1.5)

. 
. more

. twoway rarea ci_l2 ci_u2 disabil, color(yellow)                         ///
>     || line yhat2 disabil, sort                                         ///
>     || scatter dollars disabil if dollars<1, symbol(Oh)                 ///
>     , ytitle(NIH Research Funds ($ Billions))                           ///
>       ylabel(0 (.1) .5, angle(0)) xlabel(0 (1) 9) legend(off)           ///
>       xsize(2.7) ysize(1.964) scale(1.5)

. more

. log close
       log:  C:\mydocs\MPH\Text\SecondEdition\WebDoFiles\2.18.Funding.log
  log type:  text
 closed on:  28 Dec 2007, 10:56:18
-------------------------------------------------------------------------------------------------