## Not run:
# bystats(sex==2, county, city)
# bystats(death, race)
# bystats(death, cut2(age,g=5), race)
# bystats(cholesterol, cut2(age,g=4), sex, fun=median)
# bystats(cholesterol, sex, fun=quantile)
# bystats(cholesterol, sex, fun=function(x)c(Mean=mean(x),Median=median(x)))
# latex(bystats(death,race,nmiss=FALSE,subset=sex=="female"), digits=2)
# f <- function(y) c(Hazard=sum(y[,2])/sum(y[,1]))
# # f() gets the hazard estimate for right-censored data from exponential dist.
# bystats(cbind(d.time, death), race, sex, fun=f)
# bystats(cbind(pressure, cholesterol), age.decile,
# fun=function(y) c(Median.pressure =median(y[,1]),
# Median.cholesterol=median(y[,2])))
# y <- cbind(pressure, cholesterol)
# bystats(y, age.decile,
# fun=function(y) apply(y, 2, median)) # same result as last one
# bystats(y, age.decile, fun=function(y) apply(y, 2, quantile, c(.25,.75)))
# # The last one computes separately the 0.25 and 0.75 quantiles of 2 vars.
# latex(bystats2(death, race, sex, fun=table))
# ## End(Not run)
Run the code above in your browser using DataCamp Workspace