Robert A. LaBudde
2007-Jul-23 02:33 UTC
[R] Conditional logistic regression on n:m matched "cohort" data
I am designing an interlaboratory validation study for a presence/absence alternative method test kit vs. a presence/absence reference method test kit. There will be 10 laboratories conducting tests using both methods. In each laboratory, there will be 5 specimens tested, each of the 5 specimens twice by both methods (alternative, standard). The total number of data are 10 x 5 x 4 = 200. The general structure is as follows: id: sequential result label, 1 to 200. lab: a factor which has 10 levels, one for each participating lab. specimen: 1 to 5 in each lab, no connection between labs. method: A or R. result: 0 or 1 (absence or presence). This experiment appears to be equivalent to a matched cohort study. The matching is done on specimens, and there are 2 alternative and 2 reference results for each specimen. The sketchy description for clogit() of package "survival" presumes a matched case-control study is being analyzed. I am looking for the correct syntax to use in clogit() for this experiment. Is library('survival') clogit(result ~ lab + strata(specimen), method="exact", data=mydata) all that is required? Thanks. ===============================================================Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"
Seemingly Similar Threads
- Conditional logistic regression on n:m matched "cohort" data [Corrected]
- Confidence interval for binomial variance
- S-Plus "resample" package and associated functions
- R function which finds confidence interval for binomial variance
- Singularity in simple ANCOVA problem