# NOT RUN {
x <- climate_data
phiDelta <- phiDelta.stats(x[,-1],x[,1])
phiDelta.plot(phiDelta$phi, phiDelta$delta)
phiDelta.plot(phiDelta$phi, phiDelta$delta,
ratio = phiDelta$ratio,
border = "green",
iso_neg_predictive_value = TRUE,
crossing = FALSE)
# }
Run the code above in your browser using DataLab