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

drmfit 
a fitted DRM object (an output from the

Summaries of the fitted DRM object, including the basic information about the fitted DRM, maximum empirical likelihood estimates (MELE) 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.
For more detailed but less userfriendly formatted
information about a fitted DRM object, one should directly
check the output from the drmdel
function.
See drmdel
for details.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21  # 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 builtin 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)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.