# NOT RUN {
# Do not run
# Let's say high emitter is 5 times the normal ef.
co_efhe <- ef_cetesb(p = "COd", "PC_G") * 5
# Let's say that the perfil of high emitters increases linearly
# till 30 years and after that percentage is constant
perc <- c(seq(0.01, 0.3, 0.01), rep(0.3, 20))
# Now, lets use our ef with normal deterioration
co_ef_normal <- ef_cetesb(p = "COd", "PC_G")
efd <- ef_whe(efhe = co_efhe, phe = perc, ef = co_ef_normal)
# now, we can plot the three ef
plot(co_efhe)
lines(co_ef_normal, pch = 16, col = "red" )
lines(efd, pch = 16, col = "blue")
# }
Run the code above in your browser using DataLab