unconditionalStdErr: Get the unconditional standard error for each variable from...

Description Usage Arguments Details Value See Also

View source: R/oadaAICtable.R

Description

Calculates the unconditional standard error for each variable fitted in a set of models fitted using oadaAICtable or codetadaAICtable.

Usage

1
2
3
unconditionalStdErr(nbdaAICtable, typeFilter = NULL, netFilter = NULL,
  baselineFilter = NULL, includeAsocial = TRUE, includeNoILVs = TRUE,
  nanReplace = FALSE)

Arguments

nbdaAICtable

an object of class oadaAICtable or tadaAICtable.

typeFilter

an optional string allowing the user to get the support for variables within a specified type of models. e.g. typeFilter="additive" gets the variable support in the subset of additive models. See typeSupport for an explanation of the different model types.

netFilter

an optional string allowing the user to get the support for variables from models with a specific combination of network effects, e.g. netFilter = "0:1:0" gets the variable support in the subset of models containing only network 2. See networksSupport for an explanation of network combinations.

baselineFilter

an optional string allowing the user to get the support for variables within the subset of models with a specific baseline function, e.g. typeFilter="gamma" gets the variable support in the subset of models with a gamma baseline rate function. e.g. typeFilter="additive" gets the variable support in the subset of additive models.

includeAsocial

logical indicating whether asocial models should be included.

includeNoILVs

logical indicating whether models with no ILVs should be included.

nanReplace

logical indicating whether standard errors recorded as NaNs should be replaced with the model averaged mean SE across all models for which the SE for that variable could be calculated.

Details

A standard error (SE) from an individual model is conditional on that model. The unconditional standard error (USE) is not really "unconditional" but is rather conditional on the whole set of models fitted rather than an individual model. It takes into account the SEs within each model and the variation in maximum likelihood estimates among models. SEs cannot be always be calulated in an NBDA, meaning we are left with a situation in which USEs cannot be calculated if 1 or more models have SE=NaN- even if there are only a few models of low weight for which this is the case. Here we offer a pragmatic solution to such cases- by replacing the SEs in such models with a model-averaged SE across all models in which the SE can be calculated, we are able to estimate a USE that is likely to be approximately correct. However, we recommend that this solution only be used if SEs are absent for only a few models of low weight. In other cases we recommend that the user use confidence intervals from profLikCI conditional on the best model as a measure of uncertainty in parameter estimates. For s parameters multiModelLowerLimits and multiModelPropST then provide a means to assess the robustness of the findings to model selection uncertainty.

Value

numeric vector giving the model averaged estimate for each variable.

See Also

oadaAICtable, tadaAICtable, typeSupport, typeByNetworksSupport, networksSupport, modelAverageEstimates, unconditionalStdErr.


whoppitt/NBDA documentation built on April 25, 2021, 7:55 a.m.