data(andes)
andes <- andes$ages
# remove missing data
andes <- andes[complete.cases(andes)]
# remove outliers
andes <- sort(andes[which(andes < 10)])
stratCI(andes, method="Strauss-Sadler89",
nparams="one.par", C=0.95, endpoint="first")
stratCI(andes, method="Strauss-Sadler89",
nparams="one.par", C=0.95, endpoint="last")
stratCI(andes, method="Strauss-Sadler89",
nparams="two.par", C=0.95)
stratCI(andes, method="Marshall94", confidence = 0.95,
quantile = 0.8)
stratCI(andes, method="Marshall94", confidence = 0.95,
quantile = 0.95)
Run the code above in your browser using DataLab