# NOT RUN {
data(smokerlead)
attach(smokerlead)
# The following example reproduces parts of the first column of Table 3 in Rosenbaum (2012).
w<-rank(abs(lead))
wd<-rank(abs(lead))*dose
wl<-rank(abs(llead))
H<-cbind(w,wd,wl)
tt(lead,H,gamma=2.8)
tt(lead,H,gamma=3.2)
rm(w,wd,wl,H)
detach(smokerlead)
# }
Run the code above in your browser using DataLab