/*************************************************
 *  Kleinbaum, Chapter 7, Problem #1, p. 111     *
 *************************************************/

/* Read in the data */
capture clear
input age weight
 6 0.029
 7 0.052
 8 0.079
 9 0.125
10 0.181
11 0.261
12 0.425
13 0.738
14 1.130
15 1.882
16 2.812
end


/* Create log base 10 variable */
generate logweight=log10(weight)
format weight %9.3f
format logweight %9.3f
list

reg weight age
local f=e(mss)/(e(rss)/e(df_r))
local p= Ftail(e(df_m),e(df_r),e(F))
di _n _col(2) as text "F=" as result %8.2fc `f;' _col(15) as text "p=" as result %5.4f `p'


reg logweight age
local f=e(mss)/(e(rss)/e(df_r))
local p= Ftail(e(df_m),e(df_r),e(F))
di _n _col(2) as text "F=" as result %8.2fc `f;' _col(15) as text "p=" as result %5.4f `p'