R/dexr_output_statistics.R

In UniK-INES/dexR: Analysing and Reporting DEX Experiments

#' Calculate aggregate figures of traded energy.
#' @param dexpa 
#' @param cinfos 
#' @param requestdata
#' @return markdown table
#' 
#' @author Sascha Holzhauer
#' @export
output_statistics_comp_energy <- function(dexpa, cinfos, requestdata) {

	requestdata <- dplyr::filter(requestdata, status %in% c(1,2,3))
	requestTable <- plyr::ddply(requestdata, c("id"), function(data) {
				# data <- requestdata[requestdata$id == "enavi_08-01",]
				result <- data.frame(
				id = unique(data$id),
				energy_accepted_supply = sum(data[data$energy_accepted < 0, "energy_accepted"], na.rm=T),
				energy_accepted_demand = sum(data[data$energy_accepted > 0, "energy_accepted"], na.rm=T),
				energy_requested_supply = sum(data[data$energy_requested < 0, "energy_requested"], na.rm=T),
				energy_requested_demand = sum(data[data$energy_requested > 0, "energy_requested"], na.rm=T))
				result
			})
	
	cinfosTable <- plyr::ddply(cinfos, c("id"), function(data) {
				result <- data.frame(
						id = unique(data$id),
						energy_cleared = sum(data$energy_cleared, na.rm=T))
				result
			})
	
	if (nrow(requestTable) == 0) {
		futile.logger::flog.warn("Request data has no rows!", 
				"dexr.output.stats.comp.energy")
	} else if (nrow(cinfosTable) == 0) {
		futile.logger::flog.warn("Clearing information has no rows!", 
				"dexr.output.stats.comp.energy")
	} else {
		aggTable <- merge(requestTable, cinfosTable, by="id", all=T)
		
		knitr::kable(aggTable, format="markdown", caption="Aggregated energy information",
				digits = 3,
				col.names = c(	"ID",
						"Energy accepted (asks/supply)",
						"Energy accepted (bids/demand)",
						"Requested energy (asks/supply)",
						"Requested energy (bids/demand)",
						"Energy cleared"))
	}
}
#' Calculate aggregate figures of costs.
#' @param dexpa 
#' @param cinfos 
#' @return markdown table
#' 
#' @author Sascha Holzhauer
#' @export
output_statistics_comp_costs <- function(dexpa, cinfos, requestdata) {
	
	requestdata <- dplyr::filter(requestdata, status %in% c(1,2,3))
	requestTable <- plyr::ddply(requestdata, c("id"), function(data) {
				# data <- requestdata[requestdata$id == "enavi_08-01",]
				result <- data.frame(
						id = unique(data$id),
						price_requested_supply = sum(data[data$energy_requested < 0, "price_requested"], na.rm=T),
						price_requested_demand = sum(data[data$energy_requested > 0, "price_requested"], na.rm=T),
						price_accepted_supply = sum(data[data$energy_requested < 0, "price_accepted"], na.rm=T),
						price_accepted_demand = sum(data[data$energy_requested > 0, "price_accepted"], na.rm=T))
				result
			})
	
	cinfosTable <- plyr::ddply(cinfos, c("id"), function(data) {
				result <- data.frame(
						id = unique(data$id),
						energy_cleared = sum(data$energy_cleared, na.rm=T))
				result
			})
	
	if (nrow(requestTable) == 0) {
		futile.logger::flog.warn("Request data has no rows!", 
				"dexr.output.stats.comp.energy")
	} else if (nrow(cinfosTable) == 0) {
		futile.logger::flog.warn("Clearing information has no rows!", 
				"dexr.output.stats.comp.energy")
	} else {
		knitr::kable(requestTable, format="markdown", caption="Aggregated energy information",
				digits = 3,
				col.names = c(	"ID",
						"Requested price (asks/supply)",
						"Requested price (bids/demand)",
						"Accepted price (asks/supply)",
						"Accepted price (bids/demand)"))
	}
}
#' Calculate aggregate figures of number of clients.
#' @param dexpa 
#' @param requestdata 
#' @return markdown table
#' 
#' @author Sascha Holzhauer
#' @export
output_statistics_comp_numclients <- function(dexpa, requestdata) {
	
	requestdata <- dplyr::filter(requestdata, status %in% c(1,2,3))
	requestTable <- plyr::ddply(requestdata, c("id"), function(data) {
				# data <- requestdata[requestdata$id == "enavi_08-01",]
				result <- data.frame(
						id = unique(data$id),
						numclients = length(unique(data$username)))
				result
			})
	
	if (nrow(requestTable) == 0) {
		futile.logger::flog.warn("Request data has no rows!", 
				"dexr.output.stats.comp.clients")
	} else {
		knitr::kable(requestTable, format="markdown", caption="Aggregated client information",
				digits = 3,
				col.names = c(	"ID",
						"Number of clients"))
	}
}
#' Show aggregate gini figures in table per ID.
#' @param dexpa 
#' @param data data.frame with column "id' and individual data point in column 'gini' 
#' @return markdown table
#' 
#' @author Sascha Holzhauer
#' @export
output_statistics_comp_gini <- function(dexpas, data) {
	tableData <- plyr::ddply(data, c("Kind", "id"), function(df) {
				# data <- requestdata[requestdata$id == "enavi_08-01",]
				result <- data.frame(
						Mean = mean(df$values, na.rm = T),
						Sum = sum(df$values, na.rm = T))
				result
			})
	
	if (nrow(tableData) == 0) {
		futile.logger::flog.warn("Gini data has no rows!", 
				"dexr.output.stats.comp.clients")
	} else {
		knitr::kable(tableData, format="markdown", caption="Aggregated client information",
				digits = 3,
				col.names = c(	"Type",
						"ID", "Gini mean", "Gini sum"))
	}
}