Nothing
#' Summarize the Multivariate Regression Discontinuity
#'
#' \code{summary.mrd} is a \code{summary} method for class \code{"mrd"}.
#' It is based on \code{summary.RD} function in the "rdd" package.
#'
#' @method summary mrd
#'
#' @param object An object of class \code{"mrd"}, usually a result of a call to \code{\link{mrd_est}}.
#' @param level A numeric value between 0 and 1 specifying the confidence level for confidence intervals. The default is 0.95.
#' @param digits A non-negative integer specifying the number of digits to display.
#' The default is \code{max(3, getOption("digits") - 3)}.
#' @param ... Additional arguments passed to \code{summary}.
#'
#' @return \code{summary.mrd} returns a list which has the following components depending on methods
#' implemented in the \code{"mrd"} object:
#' \item{center_coefficients}{A matrix containing bandwidths, number of observations, estimates,
#' SEs, confidence intervals, z-values and p-values for each estimated bandwidth and/or parametric model.}
#' \item{univR_coefficients}{A matrix containing bandwidths, number of observations, estimates,
#' SEs, confidence intervals, z-values and p-values for each estimated bandwidth and/or parametric model.}
#' \item{univM_coefficients}{A matrix containing bandwidths, number of observations, estimates,
#' SEs, confidence intervals, z-values and p-values for each estimated bandwidth and/or parametric model.}
#' \item{front_coefficients}{A matrix containing estimates and confidence intervals (if applicable)
#' for the complete model.}
#' \item{front_ht_coefficients}{A matrix containing estimates and confidence intervals (if applicable)
#' for the heterogeneous treatment model.}
#' \item{front_t_coefficients}{A matrix containing estimates and confidence intervals (if applicable)
#' for the treatment only model.}
#'
#' @references Drew Dimmery (2016). rdd: Regression Discontinuity Estimation. R package
#' version 0.57. https://CRAN.R-project.org/package=rdd
#'
#' @importFrom stats residuals
#'
#' @include mrd_est.R
#'
#' @export
summary.mrd <- function(object, level = 0.95, digits = max(3, getOption("digits") - 3), ...) {
call.copy <- object$call
if ("data" %in% names(call.copy) && length(call.copy$data) > 1) {
call.copy$data <- "(.)"
}
if ("subset" %in% names(call.copy) && length(call.copy$subset) > 1) {
call.copy$subset <- "(.)"
}
cat("\n")
cat("Call:\n")
# print(object$call)
print(call.copy)
cat("\n")
out = list()
if ('center' %in% names(object)){
cat("Centering Method:\n")
center_obj = object$center$tau_MRD
center_summary = summary(center_obj, level, digits)
center_out = center_summary$coefficients
out[['center_coefficients']] = center_out
cat("\n")
}
if ('univ' %in% names(object)){
cat("Univariate Method on Frontier R:\n")
univR_obj = object$univ$tau_R
univR_summary = summary(univR_obj, level, digits)
univR_out = univR_summary$coefficients
out[['univR_coefficients']] = univR_out
cat("Univariate Method on Frontier M:\n")
univM_obj = object$univ$tau_M
univM_summary = summary(univM_obj, level, digits)
univM_out = univM_summary$coefficients
out[['univM_coefficients']] = univM_out
cat("\n")
}
if ('front' %in% names(object)){
cat("Frontier Method:\n")
front_obj = object$front$tau_MRD
front_summary = summary(front_obj, level, digits)
front_out = front_summary$coefficients
front_ht_out = front_summary$ht_coefficients
front_t_out = front_summary$t_coefficients
out[['front_coefficients']] = front_out
out[['front_ht_coefficients']] = front_ht_out
out[['front_t_coefficients']] = front_t_out
cat("\n")
}
return(invisible(out))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.