Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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


From   "Harsha" <harshabl.76@gmail.com>
To   "'Carlo Lazzaro'" <carlo.lazzaro@tin.it>, <statalist@hsphsun2.harvard.edu>
Subject   st: RE: Wilcoxon signed-rank test for clustered data - flag: Stata 9/2 SE
Date   Tue, 5 Oct 2010 14:30:39 +0530

Thanks Carlo. I will look into Roger Newson's webpage. I will also run the
bootstrap and get back if I have some queries.
Harsha


-----Original Message-----
From: Carlo Lazzaro [mailto:carlo.lazzaro@tin.it] 
Sent: 05 October 2010 13:02
To: statalist@hsphsun2.harvard.edu
Cc: 'Harsha'
Subject: R: Wilcoxon signed-rank test for clustered data - flag: Stata 9/2
SE

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/


© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index