# NOT RUN {
# }
# NOT RUN {
data(saliva)
data(throat)
data(tonsils)
### Combine the data sets into a single list
group.data <- list(saliva, throat, tonsils)
### Get PI using MLE with CI
mle <- Est.PI(group.data)$MLE
### Plot with Error Bars
Plot.PI(mle)
### Plot without Error Bars
Plot.PI(mle, FALSE)
### Plot with Error Bars and scaling
Plot.PI(mle, TRUE, TRUE)
# }
Run the code above in your browser using DataLab