data(AL)
group <- AL$group
table(group)
##take the negated kfront marker measurements
kfront <- -AL$kfront
x <- kfront[group=="D-"]
y <- kfront[group=="D0"]
z <- kfront[group=="D+"]
vus.cut <- VUS.CutPoint(x=x,y=y,z=z,cut.seq=NULL)
Run the code above in your browser using DataLab