Description Usage Arguments Value Note Examples
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 user-friendly 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 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.