enviro <- environmental
## make an ordered factor (so it will not be reordered)
enviro$smell <- cut(enviro$ozone, breaks=c(0,30,50,Inf),
labels=c("ok","hmmm","yuck"), ordered=TRUE)
marginals(enviro)
Run the code above in your browser using DataLab