# NOT RUN {
if (require(mosaicData)) {
data(HELPrct)
# calculate the observed difference
mean(age ~ sex, data=HELPrct)
obs <- diffmean(age ~ sex, data=HELPrct); obs
# calculate the permutation distribution
nulldist <- do(100) * diffmean(age ~ shuffle(sex),
data=HELPrct)
histogram(~ diffmean, groups=(diffmean >= obs), nulldist,
xlab="difference in means")
}
# }
Run the code above in your browser using DataLab