labels.eoar: labels.eoar - Extract parameter labels from 'eoar' objects.
In tmcd82070/EoAR: Evidence of Absence Regression (EoAR)

Description Usage Arguments Details Author(s) See Also Examples

labels.eoar - Extract parameter labels from eoar objects.

1 2	## S3 method for class 'eoar' labels(obj, type = "coef")

`obj`	An `eoar` model object. See `eoar`.
`type`	The type of parameter label requred or a regular expression. Parameter type possibilities are "coef" or "derived". Regular expressions are used to match parameter labels using `grep`. If `type` is not "coef" and not "derived" and the regular expression fails to match anything, all parameter labels.

Coefficient labels are for variables in the log-linear model for lambda. Derived parameter labels are non-coefficient parameters link M, Mtot, and lambda.

Trent McDonald

coef, eoar

# A 3 year study of 7 sites. 21 "cells". lambda change = 20/year
set.seed(9430834) # fixes Y and g of this example, but not the RNG's used in chains
ns <- 3

ny <- 7
g <- data.frame(
 alpha = rnorm(ns*ny,70,2),
 beta = rnorm(ns*ny,700,25)
)
Y <- rbinom(ns*ny, c(rep(20,ny), rep(40,ny), rep(60,ny)), g$alpha/(g$alpha+g$beta))

df <- data.frame(year=factor(c(rep("2015",ny),rep("2016",ny),rep("2017",ny))),
   Year=c(rep(1,ny),rep(2,ny),rep(3,ny)))

# Uninformed eoar (use low number of iterations because it's and example)
eoar.1 <- eoar(Y~year, g, df, nburn = 1000, niters= 50*10, nthins = 10 )

labels(eoar.1)
labels(eoar.1,"derived")
labels(eoar.1,"^M")  # all M parameters
labels(eoar.1,"\\[3\\]$") # M[3] and lambda[3]
labels(eoar.1,".") # all parameter labels

plot(eoar.1$out[,labels(eoa.1)])  # trace plot of coefficients.