drmdel (version 1.0)

summaryDRM: Print simple summaries of a DRM fit object

Description

Print simple summaries of a density ratio model (DRM) fit object, i.e. an output from drmdel function.

Usage

summaryDRM(drmfit)

Arguments

drmfit

a fitted DRM object (an output from the drmdel function). See drmdel for details.

Value

• Summaries of the fitted DRM object, including the basic information about the fitted DRM, maximum dual empirical likelihood estimates (MDELE) of the DRM parameters, results of the dual-empirical-likelihood ratio (DELR) test for hypothesis about the DRM parameter $\beta$ and some summary statistics of the estimated distribution functions $F_k$, $k = 0, 1, ..., m$.

Examples

# Data generation
set.seed(25)
n_samples <- c(100, 200, 180, 150, 175)  # sample sizes
x0 <- rgamma(n_samples[1], shape=5, rate=1.8)
x1 <- rgamma(n_samples[2], shape=12, rate=1.2)
x2 <- rgamma(n_samples[3], shape=12, rate=1.2)
x3 <- rgamma(n_samples[4], shape=18, rate=5)
x4 <- rgamma(n_samples[5], shape=25, rate=2.6)
x <- c(x0, x1, x2, x3, x4)

# Fit a DRM with the basis function q(x) = (x, log(abs(x))), which
# is the basis function for gamma family. This basis function is
# the built-in basis function 6.
drmfit <- drmdel(x=x, n_samples=n_samples, basis_func=6)

# A brief summary of the DRM fit
summaryDRM(drmfit)

# List the details of the DRM fit.  See help(drmdel) for details.
names(drmfit)