Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: Re: power analysis for panel data


From   "Christopher W. Ryan" <[email protected]>
To   [email protected]
Subject   st: Re: power analysis for panel data
Date   Tue, 07 Nov 2006 13:57:48 -0500

I was not very explicit in my original question.

To be more clear:

I'm conducting an observational study on the effect of simple office
advice by the physician on a patient's obesity (measured by BMI). I have
an unbalanced panel: roughly 100 patients observed for about two years.
 But they made their visits based on clinical needs, so the number of
visits varies from patient to patient.  I've "condensed" their visits
into one representative visit per quarter per patient--but not all
patients have a visit in any particular quarter.  qdate is the quarterly
date variable.

The patients are a 10% random sample of all the patients with BMI >= 30,
who made visits to our practice during a specified interval.  The
selected patients' charts were then abstracted, using every visit
subsequent to the index visit, for the ensuing 2.5 years.

So then I use -xtreg- with fe to examine the effect of a previous
visit's advice on the current visit's BMI.  So the independent variables
are lagged.  "Advice" can be to follow a certain diet, to exercise, to
see a dietician, to join a self-help group, etc.  Diet advice and
exercise advice are ordinal:  none, non-specific, or specific.  I also
look at the effects of charting the obesity in certain ways, and whether
the diagnosis of obsesity was "officially" entered into the chart.


Here is some of the output.

. xi:xtreg bmi qdate Lobsub Lobassess Lobicd9 i.Ldietadvice
i.Lexeradvice Ldietician Lbecouns
> Lmeds Lsurgery Lshgroup, fe

i.Ldietadvice     _ILdietadvi_0-2     (naturally coded; _ILdietadvi_0
omitted)
i.Lexeradvice     _ILexeradvi_0-2     (naturally coded; _ILexeradvi_0
omitted)

Fixed-effects (within) regression     Number of obs      =       272
Group variable (i): mrnumber          Number of groups   =        98

R-sq:  within  = 0.1510               Obs per group: min =         1
       between = 0.1323                              avg =       2.8
       overall = 0.0767                              max =         6

                                      F(13,161)          =      2.20
corr(u_i, Xb)  = 0.2001               Prob > F           =    0.0115

---------------------------------------------------------
         bmi |      Coef.   Std. Err.      t    P>|t|
-------------+-------------------------------------------
       qdate |   -.077571   .0607539    -1.28   0.204
      Lobsub |   .2856105   .3644859     0.78   0.434
   Lobassess |    -1.3555   .6527376    -2.08   0.039
     Lobicd9 |   1.408771   .6667927     2.11   0.036
_ILdietadv~1 |   .2555625   .6827685     0.37   0.709
_ILdietadv~2 |   1.912082   .7939111     2.41   0.017
_ILexeradv~1 |  -.1791865   .5411985    -0.33   0.741
_ILexeradv~2 |   .3961547   1.153241     0.34   0.732
  Ldietician |  -.5891818   .8675808    -0.68   0.498
    Lbecouns |  -1.180993   1.698886    -0.70   0.488
       Lmeds |   .3907819   1.052359     0.37   0.711
    Lsurgery |    4.65996   1.598309     2.92   0.004
    Lshgroup |   3.550525   1.706046     2.08   0.039
       _cons |   51.30236   11.05319     4.64   0.000
-------------+-------------------------------------------
     sigma_u |  5.4841536
     sigma_e |  1.3633508
         rho |  .94179599   (fraction of variance due to u_i)
------------------------------------------------------------------------------
F test that all u_i=0:   F(97, 161) =    46.51    Prob > F = 0.0000


The lack of any statistically significant beneficial effect of any of
the interventions on BMI does not surprise me, given the generally
intractable nature of obesity.  The general futility of simple
office-based exhortations to lose weight is part of my point.

But what is the power of this study?  I don't know how to calculate
that.  Am I failing to see statistically significant beneficial effects
on BMI because of inadequate power?

Hence my question about resources or references about power and sample
size calculations in panel analysis.

Thanks for any advice or pointers to references.

-- 
Christopher W. Ryan, MD
SUNY Upstate Medical University Clinical Campus at Binghamton
and Wilson Family Practice Residency, Johnson City, NY
cryanatbinghamtondotedu
GnuPG and PGP public keys available at http://pgp.mit.edu

"If you want to build a ship, don't drum up the men to gather wood,
divide the work and give orders. Instead, teach them to yearn for the
vast and endless sea."  [Antoine de St. Exupery]


Rodrigo A. Alfaro wrote:
> I am assuming that power is computed for fixed T. 
> Rodrigo.
> 
> ----- Original Message ----- 
> From: "Christopher W. Ryan" <[email protected]>
> To: "Statalist" <[email protected]>
> Sent: Monday, October 16, 2006 10:08 PM
> Subject: st: power analysis for panel data
> 
> 
> Can anyone direct me to any resources that would help me understand
> power analysis in the context of panel data?  Thanks.
> 
> --Chris
*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index