Dummy Tradeoff Data

a <- sample(1000:8000, size = 100, replace = TRUE) / 1000  # yield dummy
b <- sample(9000:40000, size = 100, replace = TRUE) / 1000  # carbon dummy
cl <- b / a   # carbon loss:yield
astd <- (a - min(a)) / diff(range(a))   # standardized yield
bstd <- 1 - (b - min(b)) / diff(range(b))  # standardized carbon

Distribution of data

par(mfrow = c(3, 3))
hist(a, main = "Dummy yield")
hist(b, main = "Dummy carbon")
hist(astd, main = "Standardized yield")
hist(bstd, main = "Standardized carbon")
hist(cl, main = "C/yld")
hist(bstd * astd, main = "C std * Y std")
hist(bstd + astd, main = "C std + Y std")

Comparing algorithms to one another

par(mfrow = c(2, 2))
plot(cl, astd + bstd, main = "C/yld vs C std + Y std")
plot(cl, astd * bstd, main = "C/yld vs C std * Y std")
plot(astd * bstd, astd + bstd, main = "C std * Y std vs C std + Y std")


PrincetonUniversity/agroEcoTradeoff documentation built on May 8, 2019, 3:12 p.m.