>Ah, I think I see what you want. Try this on each pair of exclusive sets:
>
n_total<-dim(mydataset)[1]
under65<-mydataset$age <= 65
n_under65<-sum(under65)
under65row<-c(sum(mydataset$dose[under65] == "FD"),
sum(mydataset$dose[under65] == "RD"),
sum(mydataset$vitalstatus[under65] == "dead" &
mydataset$dose[under65] == "FD"),
sum(mydataset$vitalstatus[under65] == "dead" &
mydataset$dose[under65] == "RD"))
over65row<-c(sum(mydataset$dose[!under65] == "FD"),
sum(mydataset$dose[!under65] == "RD"),
sum(mydataset$vitalstatus[!under65] == "dead" &
mydataset$dose[!under65] == "FD"),
sum(mydataset$vitalstatus[!under65] == "dead" &
mydataset$dose[!under65] == "RD"))
>
>Then under65row and over65row should be the first two rows of your result.
>Can't test this at the moment, but I don't think it's too far wrong.
Thanks Jim.
Yes it looks like that code should do the job. I was really hoping for a code
like "SummariseForSubsetAnalysis(mydataset, by=mydatatset$dose, subsets=c(age,
renal, sex, toxicity), event=survival )" which would magically do it for me ;-)
I guess if this is something I start having to do lots I might have to write
one.
Surprised one doesn't seem to exist - perhaps the number of variations in what
people want would be too complex.
Calum
********************************************************************************************************************
This message may contain confidential information. If yo...{{dropped:21}}
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
