st: R: Wilcoxon signed-rank test for clustered data - flag: Stata 9/2 SE

["Carlo Lazzaro" <carlo.lazzaro@tin.it>]

Hasha may want to consider a bootstrap procedure (that allows for
clusters)like in the following examples:
--------- code begins-----------------------------------
drop _all
set obs 10
input id time
            id       time
  1 0
  1 1
  2 0
  2 1
  3 0
  3 1
  4 0
  4 1
  5 0
  5 1
g acu_sight=2*uniform() if  time==0
replace acu_sight=5*uniform() if time==1
qui sum  acu_sight, d
bootstrap r(p50), reps(1000) cluster(id) : qui sum  acu_sight, d
--------- code ends-------------------------------------

In your query, the issue with clustering seems to be double, since you have
two operated eyes for each principal sample unit (ie: patient) and two
different points in time for measurement (before and after ocular surgery).
Please note that the previous example takes into account only one source of
clustering, that is visual acuity before and after the surgery. 
Last but absolutely not least, I would strongly recommend you to take a look
at Roger Newson's webpage (http://www.imperial.ac.uk/nhli/r.newson/). 
Roger has addressed a lot of  his research activity in creating Stata
user-written programmes, some of those focusing on non parametric median
statistics.

HTH and Kind Regards,
Carlo

-----Messaggio originale-----
Da: owner-statalist@hsphsun2.harvard.edu
[mailto:owner-statalist@hsphsun2.harvard.edu] Per conto di Harsha
Inviato: martedì 5 ottobre 2010 6.46
A: statalist@hsphsun2.harvard.edu
Oggetto: st: Wilcoxon signed-rank test for clustered data

Dear All,

I have a data comparing a few parameters in a group of patients before and
after a particular surgery. These surgeries are done on both eyes of a few
patients resulting in a correlated data. A few variables which I am trying
to compare are non-parametric. I am trying to apply Wilcoxon signed rank
test but am faced with the problem of the observations being not independent
between the 2 eyes of the patients. There is a program on SAS by Rosner's
group to use Wilcoxon for clustered data. Is there something similar in
stata? Is there any other ways to compare the medians between groups in
correlated data?

Kindly suggest.

Harsha

*
*   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/