#' Aggregates data and performs poisson regression using GLM.
#' @title Poisson regression by aggregated data using glm
#' @param formula The formula as specified in glm.
#' @param data dataset.
#' @param aggvars vector of variables in data to be aggregated (summarized).
#' @param byvars vector of variables in data to aggregate by.
#' @export
poissonregression <- function(formula,data,aggvars,byvars){
.SD=NULL
if (class(data)[1]!="data.table") data <- data.table(data)
aggdata <- data[,lapply(.SD,sum),by=byvars,.SDcols=aggvars]
fit <- stats::glm(formula,data=aggdata,family=stats::poisson())
fit
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.