R/VegfruitShare.R

In pik-piam/magpie4: MAgPIE outputs R package for MAgPIE version 4.x

#' @title VegfruitShare
#' @description Calculates the share of fruits, vegetables and nuts in total food supply from the food demand model
#' 
#' @export
#' 
#' @param gdx GDX file
#' @param file a file name the output should be written to using write.magpie
#' @param level Level of regional aggregation; "iso" ISO country codes, "reg" (regional), "glo" (global)
#' @param after_shock FALSE is using the exogenous real income and the prices before a shock, TRUE is using the endogeenous real income that takes into account food price change on real income
#' @param calibrated if FALSE, the true regression outputs are used, if TRUE the values calibrated to the start years are used
#' @param attributes unit: kilocalories per day ("kcal"), g protein per day ("protein"). Mt reactive nitrogen ("nr").
#' 
#' @return magpie object with the livestock share in a region or country. Unit is dimensionsless, but value depends on chosen attribute
#' @author Benjamin Bodirsky
#' @importFrom magclass dimSums
#' @examples
#' 
#'   \dontrun{
#'     x <- VegfruitShare(gdx)
#'   }
#' 


VegfruitShare<- function(gdx, 
                          file=NULL, 
                          level="reg", 
                          after_shock=TRUE, 
                          calibrated=TRUE,
                          attributes="kcal"){
  
  kcal<-Kcal(gdx, level=level, products="kfo", 
                   product_aggr=FALSE, 
                   after_shock=after_shock, 
                   calibrated=calibrated,
                   attributes=attributes,
                   per_capita=TRUE)
  out=dimSums(kcal[,,"others"],dim=3.1)/dimSums(kcal,dim=3.1)
  
  out(out,file)
}