pres <- c(1, 12.5, 25, 50, 100, 200, 400, 800, 1600)
VR <- c(0.99, 0.97, 0.95, 0.89, 0.78, 0.61, 0.43, 0.23, 0.01)
plot(VR ~ log10(pres), type = "b") # find the 'mcp'
sigmaP(VR, pres, method = "casagrande", mcp = 1.6)
sigmaP(VR, pres, method = "pacheco")
sigmaP(VR, pres, method = "reg3")
# End (not run)
Run the code above in your browser using DataLab