-
Notifications
You must be signed in to change notification settings - Fork 13
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
6 changed files
with
178 additions
and
8 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,8 +1,8 @@ | ||
| Type: Package | ||
| Package: mrvalidation | ||
| Title: madrat data preparation for validation purposes | ||
| Version: 2.54.2 | ||
| Date: 2024-01-10 | ||
| Version: 2.55.0 | ||
| Date: 2024-01-25 | ||
| Authors@R: c( | ||
| person("Benjamin Leon", "Bodirsky", , "[email protected]", role = c("aut", "cre")), | ||
| person("Stephen", "Wirth", role = "aut"), | ||
|
|
||
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,166 @@ | ||
| #' @title calcValidFeedConversion | ||
| #' @description calculates various feed indicators | ||
| #' | ||
| #' @param livestockSystem if TRUE, ruminant meat and milk are aggregated, and poultry meat and egg are aggregated | ||
| #' | ||
| #' @return List of magpie objects with results on country level, weight on country level, unit and description. | ||
| #' @author Benjamin Leon Bodirsky, github Copilot | ||
| #' @seealso | ||
| #' \code{\link{calcFoodSupplyPast}}, | ||
| #' \code{\link{calcValidLivestockShare}} | ||
| #' @examples | ||
| #' \dontrun{ | ||
| #' calcOutput("ValidFeed") | ||
| #' } | ||
| #' | ||
| #' @importFrom magpiesets reporthelper | ||
| #' @importFrom magclass dimOrder | ||
| calcValidFeedConversion <- function(livestockSystem = TRUE) { | ||
| ### calculate product specific feed conversion efficiency as quotient between | ||
| ### feed and animal products | ||
|
|
||
| mb <- calcOutput("FAOmassbalance", aggregate = FALSE) | ||
|
|
||
| mb <- mb[, , c("ge", "nr")] | ||
| x <- NULL | ||
| weight <- NULL | ||
|
|
||
| ap <- mb[, , "production"][, , findset("kli")] | ||
| quotientProductspecfic <- collapseNames(ap) | ||
| # extract feed items from mb | ||
| feed <- mb[, , c("feed", "feed_livst_chick", "feed_livst_egg", | ||
| "feed_livst_milk", "feed_livst_pig", | ||
| "feed_livst_rum")] | ||
| feedGrouped <- mbind( | ||
| add_dimension( | ||
| dimSums(feed[, , c("tece", "trce", "rice_pro", "maiz", "brans")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Cereal Intensity"), | ||
| add_dimension( | ||
| dimSums(feed[, , c("soybean", "groundnut", "puls_pro", "rapeseed", "oils", "oilcakes", | ||
| "cassav_sp", "oilpalm", "sunflower")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Oilcrop intensity"), | ||
| add_dimension( | ||
| dimSums(feed[, , c("potato", "cassav_sp", "others", | ||
| "sugr_beet", "sugr_cane", "molasses")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Roots and other crops"), | ||
| add_dimension( | ||
| dimSums(feed[, , c("pasture")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Pasture intensity"), | ||
| add_dimension( | ||
| dimSums(feed[, , c("foddr", "res_cereals", "res_fibrous", "res_nonfibrous")], | ||
| dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Other roughage intensity"), | ||
| add_dimension( | ||
| dimSums(feed[, , c("livst_chick", "livst_egg", "livst_milk", "livst_pig", | ||
| "livst_rum", "fish")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Animal product intensity"), | ||
| add_dimension( | ||
| dimSums(feed[, , c("distillers_grain", "alcohol", "scp")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Other processed items intensity") | ||
| ) | ||
|
|
||
|
|
||
| feedProductspecific <- collapseNames(feedGrouped[, , c("feed_livst_chick", "feed_livst_egg", | ||
| "feed_livst_milk", "feed_livst_pig", | ||
| "feed_livst_rum")]) | ||
| getNames(feedProductspecific, dim = 2) <- substring(getNames(feedProductspecific, dim = 2), 6) | ||
|
|
||
| if (livestockSystem == TRUE) { | ||
| feedProductspecific <- mbind( | ||
| add_dimension( | ||
| dimSums(feedProductspecific[, , c("livst_milk", "livst_rum")], dim = c("ElementShort")), | ||
| dim = 3.2, add = "ElementShort", nm = "Ruminant meat and dairy"), | ||
| add_dimension( | ||
| dimSums(feedProductspecific[, , c("livst_chick", "livst_egg")], dim = c("ElementShort")), | ||
| dim = 3.2, add = "ElementShort", nm = "Poultry meat and eggs"), | ||
| add_dimension( | ||
| dimSums(feedProductspecific[, , c("livst_pig")], dim = c("ElementShort")), | ||
| dim = 3.2, add = "ElementShort", nm = "Monogastric meat") | ||
| ) | ||
| quotientProductspecfic <- mbind( | ||
| add_dimension( | ||
| dimSums(quotientProductspecfic[, , c("livst_milk", "livst_rum")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Ruminant meat and dairy"), | ||
| add_dimension( | ||
| dimSums(quotientProductspecfic[, , c("livst_chick", "livst_egg")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Poultry meat and eggs"), | ||
| add_dimension( | ||
| dimSums(quotientProductspecfic[, , c("livst_pig")], dim = c("ItemCodeItem")), | ||
| dim = 3.1, add = "ItemCodeItem", nm = "Monogastric meat") | ||
| ) | ||
| } else { | ||
| getNames(feedProductspecific, dim = 2) <- reportingnames(getNames(feedProductspecific, dim = 2)) | ||
| getNames(quotientProductspecfic, dim = 1) <- reportingnames(getNames(quotientProductspecfic, dim = 1)) | ||
| } | ||
|
|
||
| # Calculate feed conversion efficiency total | ||
| quotientTmp <- collapseNames(dimSums(quotientProductspecfic, dim = c("ItemCodeItem"))[, , "ge"]) | ||
| indicatorTmp <- collapseNames(dimSums(feedProductspecific, | ||
| dim = c("ItemCodeItem", "ElementShort"))[, , "ge"]) / quotientTmp | ||
| nameIndicator <- "Productivity|Feed conversion efficiency" | ||
| indicatorTmp <- setNames(collapseNames(indicatorTmp), paste0(nameIndicator, " (", "GE per GE", ")")) | ||
| quotientTmp <- setNames(collapseNames(quotientTmp), paste0(nameIndicator, " (", "GE per GE", ")")) | ||
| x <- mbind(x, indicatorTmp) | ||
| weight <- mbind(weight, quotientTmp) | ||
|
|
||
| # calculate feed conversion efficiency Livestock specific | ||
| quotientTmp <- quotientProductspecfic[, , c("ge", "nr")] | ||
| indicatorTmp <- dimSums(feedProductspecific, dim = c("ItemCodeItem"))[, , c("ge", "nr")] / quotientTmp | ||
| prefix <- "Productivity|Feed conversion efficiency|" | ||
| nameIndicator <- paste0(prefix, getNames(indicatorTmp, dim = 1), " (", "GE per GE", ")") | ||
| x <- mbind(x, setNames(collapseNames(indicatorTmp[, , "ge"]), nameIndicator)) | ||
| weight <- mbind(weight, setNames(collapseNames(quotientTmp[, , "ge"]), nameIndicator)) | ||
| nameIndicator <- paste0(prefix, getNames(indicatorTmp, dim = 1), " (", "Nr per Nr", ")") | ||
| x <- mbind(x, setNames(collapseNames(indicatorTmp[, , "nr"]), nameIndicator)) | ||
| weight <- mbind(weight, setNames(collapseNames(quotientTmp[, , "nr"]), nameIndicator)) | ||
|
|
||
| # calculate feed conversion efficiency Livestock and product specific | ||
| quotientTmp <- collapseNames(quotientProductspecfic[, , "ge"]) | ||
| indicatorTmp <- collapseNames(feedProductspecific[, , "ge"]) / quotientTmp | ||
| prefix <- "Productivity|Feed conversion efficiency|" | ||
| for (item in getNames(feedProductspecific, dim = 2)) { | ||
| nameIndicator <- paste0(prefix, item, "|+|", | ||
| getNames(indicatorTmp, dim = 1), " (", "GE per GE", ")") | ||
| x <- mbind(x, setNames(collapseNames(indicatorTmp)[, , item], nameIndicator)) | ||
| weightTmp <- collapseNames(indicatorTmp[, , item]) * NA | ||
| weightTmp[, , ] <- collapseNames(quotientTmp[, , item]) | ||
| weight <- mbind(weight, setNames(weightTmp, nameIndicator)) | ||
| } | ||
|
|
||
| # calculate feed conversion efficiency Product specific | ||
| quotientTmp <- collapseNames(dimSums(quotientProductspecfic, dim = c("ItemCodeItem"))[, , "ge"]) | ||
| indicatorTmp <- collapseNames(dimSums(feedProductspecific, dim = c("ElementShort"))[, , "ge"]) / quotientTmp | ||
| prefix <- "Productivity|Feed conversion efficiency|+|" | ||
| nameIndicator <- paste0(prefix, getNames(indicatorTmp, dim = 1), " (", "GE per GE", ")") | ||
| x <- mbind(x, setNames(collapseNames(indicatorTmp), nameIndicator)) | ||
| weightTmp <- collapseNames(indicatorTmp) * NA | ||
| weightTmp[, , ] <- collapseNames(quotientTmp) | ||
| weight <- mbind(weight, setNames(weightTmp, nameIndicator)) | ||
|
|
||
| # roughage share for ruminants as quotient of roughage and total feed | ||
| quotientTmp <- dimSums(feedProductspecific, dim = "ItemCodeItem") | ||
| indicatorTmp <- dimSums( | ||
| feedProductspecific[, , c("Other roughage intensity", "Pasture intensity")], | ||
| dim = "ItemCodeItem") / quotientTmp | ||
| if (livestockSystem == TRUE) { | ||
| quotientTmp <- quotientTmp[, , "Ruminant meat and dairy"] | ||
| indicatorTmp <- indicatorTmp[, , "Ruminant meat and dairy"] | ||
| } | ||
| prefix <- "Productivity|Roughage share|" | ||
| nameIndicator <- paste0(prefix, getNames(indicatorTmp, dim = 1), " (", "GE per GE", ")") | ||
| x <- mbind(x, setNames(collapseNames(indicatorTmp[, , "ge"]), nameIndicator)) | ||
| weight <- mbind(weight, setNames(quotientTmp[, , "ge"], nameIndicator)) | ||
|
|
||
| x <- add_dimension(x, dim = 3.1, add = "scenario", nm = "historical") | ||
| x <- add_dimension(x, dim = 3.2, add = "model", nm = "Weindl_et_al2017") | ||
|
|
||
| getNames(x) <- sub("\\|$", "", getNames(x)) | ||
| getNames(weight) <- sub("\\|$", "", getNames(weight)) | ||
|
|
||
|
|
||
| return(list(x = x, | ||
| weight = weight, | ||
| unit = "GE feed per GE product", | ||
| description = "Agricultural Demand") | ||
| ) | ||
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters