pcosnew$Q1_Y n missing unique Mean 1016 5 6 2009 1940 1994 1995 2004 2008 2009 Frequency 1 1 1 1 5 1007 % 0 0 0 0 0 99
1 = Watchful waiting 2 = Rad prost + XRT + hormone 3 = Rad prost + XRT 4 = XRT + prost 5 = Rad prost + XRT - unable to sequence 6 = Rad prost + horm 7 = Hormone + rad prost 8 = Rad prost + hormone - unable to sequence 9 = XRT + hormone 10 = hormone + XRT 11 = XRT + hormone - unable to sequence 12 = Rad prost only 13 = XRT only 14 = Hormone only
> with(trts, table(trtment, primtrts, exclude = NULL)) primtrts trtment 1 2 3 4 5 6 7 8 9 10 99 <NA> 1 561 0 0 0 0 0 0 0 0 7 110 0 2 1 14 2 2 0 12 0 13 0 0 0 0 3 0 2 74 0 0 0 0 19 1 0 0 0 6 0 0 0 21 0 0 0 5 0 0 0 0 7 0 0 0 0 0 130 0 14 0 5 1 0 8 2 0 0 0 0 1 0 4 0 0 0 0 9 0 0 0 0 23 0 0 0 5 1 0 0 10 3 0 0 0 1 0 112 0 10 15 4 0 11 3 0 0 0 0 0 0 0 2 5 6 0 12 33 0 0 0 0 0 0 1392 0 0 4 0 13 35 0 0 0 2 0 0 0 607 1 23 0 14 30 0 0 0 0 0 0 0 0 485 25 0 <NA> 0 0 0 0 0 0 0 0 0 0 0 1593 * Double check inclusion criteria "clinically localized." Frank used metastasis (E6) and clinstg from abstract: allg$Stage <- ifelse(allg$metastasis==1,"Metastatic", ifelse(allg$clinstg%in%CL , "Clinically Localized", ifelse(allg$clinstg%in%RA, "Regionally Advanced", "NA" ))) _Answer: high psa is generally greater than 20. First throw out people with psa over 50. Then use clinstg to only include 11-15. Can ignore e6, which is metastasis._correct!