spaMM (version 3.2.0)

confint.HLfit: Confidence intervals for fixed-effect parameters

Description

This computes confidence intervals for a given fixed effect parameter, based on the p_v-based approximation of the profile likelihood ratio for this parameter. The profiling is other all other fitted parameters: other fixed effects, as well as variances of random effects and spatial correlations if these were fitted.

Usage

# S3 method for HLfit
confint(object, parm, level=0.95, verbose=TRUE,...)

Arguments

object

An object of class HLfit, as returned by the fitting functions in spaMM.

parm

The name of a parameter to be fitted, or its position in the the object's $fixef vector. Valid names are those of the object's $fixef;

level

The coverage of the interval;

verbose

whether to print the interval or not. As the function returns its more extensive results invisibly, this printing is the only visible output;

…

Additional arguments (maybe not used, but conforming to the generic definition of confint).

Value

A list including the confidence interval for the target parameter, and the fits lowerfit and upperfit giving the profile fits at the confidence bounds. This is returned invisibly.

Examples

# NOT RUN {
<!-- % checked in test-confint.R -->
# }
# NOT RUN {
 data("wafers")
 wfit <- HLfit(y ~X1+(1|batch),family=Gamma(log),data=wafers,HLmethod="ML")
 confint(wfit,"X1")
# }