# NOT RUN {
data(HCD)
HCDcounts<-HCD[,-1]
HCDf<-HCD[,1]
SimpsonS<-Simpsonci(X=HCDcounts, f=HCDf, type = "Sequen",
alternative = "greater", conf.level = 0.95, dist = "MVN")
summary(SimpsonS)
# }
Run the code above in your browser using DataCamp Workspace