data(HELPrct)
# calculate the observed difference
mean(age ~ sex, data=HELPrct)
obs <- compareMean(age ~ sex, data=HELPrct); obs
# calculate the permutation distribution
nulldist <- do(100) * compareMean(age ~ shuffle(sex),
data=HELPrct)
xhistogram(~ result, groups=(result >= obs), nulldist,
xlab="difference in means")
Run the code above in your browser using DataLab