get_lr_tbl: Get LR test results

In MiguelRodo/gamlssutils:

Description

Get the GAMLSS LR test results into a neatly formatted table.

Usage

get_lr_tbl(mod, mod_null, mod_null_name, silent = FALSE)

Arguments

mod, mod_null	objects of class 'gamlss'. The model and a corresponding null model, respectively.
silent	logical. If FALSE, then warning messages are suppressed. Default is FALSE.
mod_null_null	character. Name by which to refer to null model.

Value

A table with columns Comp. model, LR DF, LR Chi^2 and LR P-value.

Examples

set.seed(1)
data_mod <- data.frame(x = rnorm(20, 5) + rgamma(20, 0.01, 0.0001))
data_mod[['y']] <- 3 * data_mod[['x']] + rnorm(20, sd = 5)
mod_mu <- gamlss::gamlss(formula = y ~ -1 +  x,
                         family = "NO",
                         data = data_mod,
                         control = gamlss::gamlss.control(trace = FALSE))
mod_sigma <- gamlss::gamlss(formula = y ~ -1 +  x,
                            sigma.formula =  ~ x,
                            family = "NO",
                            data = data_mod,
                            control = gamlss::gamlss.control(trace = FALSE))
lr_tbl <- get_lr_tbl(mod = mod_sigma, mod_null = mod_mu, mod_null_name = "sigma")

MiguelRodo/gamlssutils documentation built on July 9, 2020, 12:48 a.m.