pbLRT: Parametric bootstrap likelihood ratio test
In corncob: Count Regression for Correlated Observations with the Beta-Binomial

pbLRT

R Documentation

Parametric bootstrap likelihood ratio test

Description

Parametric bootstrap likelihood ratio test

Usage

pbLRT(mod, mod_null, B = 1000)

Arguments

`mod`	an object of class `bbdml`
`mod_null`	an object of class `bbdml`, should be nested within `mod`
`B`	Integer. Defaults to `1000`. Number of bootstrap iterations.

Value

P-value from parametric bootstrap likelihood ratio test.

Examples

data(soil_phylum_small_otu1)
mod1 <- bbdml(formula = cbind(W, M - W) ~ DayAmdmt,
phi.formula = ~ DayAmdmt,
data = soil_phylum_small_otu1)

mod2 <- bbdml(formula = cbind(W, M - W) ~ 1,
phi.formula = ~ 1,
data = soil_phylum_small_otu1)
pbLRT(mod1, mod2, B = 50)

corncob documentation built on May 29, 2024, 7:15 a.m.