# NOT RUN {
mod1 <- cnfa(x = climdat.hist, s.dat = ABPR, field = "CODE")
dep <- departure(x = climdat.hist, y = climdat.fut, s.dat = ABPR)
vuln <- vulnerability(cnfa = mod1, dep = dep)
vuln.map <- vulnerability_map(vuln)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab