gammacap_ols_hc_qcap_generic: Leverage Adjustment - Generic

In jeksterslab/gammaMatrix: Asymptotic Covariance Matrix of Covariances

Description

Leverage Adjustment - Generic

Usage

gammacap_ols_hc_qcap_generic(
  h,
  k,
  type = "hc3",
  g1 = 1,
  g2 = 1.5,
  constant = 0.7
)

Arguments

h	Numeric vector. Leverage values.
k	Positive integer. p number of regressors plus 1.
type	Character string. Correction type. Possible values are "hc0", "hc1", "hc2", "hc3", "hc4", "hc4m", and "hc5".
g1	Numeric. g1 value for type = "hc4m" or type = "hc5".
g2	Numeric. g2 value for type = "hc4m".
constant	Numeric. Constant for type = "hc5"

Value

A vector.

Author(s)

Ivan Jacob Agaloos Pesigan

References

Dudgeon, P. (2017). Some improvements in confidence intervals for standardized regression coefficients. Psychometrika, 82, 928-951. doi:10.1007/s11336-017-9563-z.

See Also

Other Gamma Matrix Functions: gammacap_adfnb(), gammacap_adf(), gammacap_gen(), gammacap_mvnadj1(), gammacap_mvnadj2(), gammacap_mvn(), gammacap_nb(), gammacap_ols_generic(), gammacap_ols_hc_generic(), gammacap_ols_hc_qcap(), gammacap_ols_hc(), gammacap_ols(), gammacapnames(), gammacap()

Examples

set.seed(42)
n <- 1000
k <- 2
z <- matrix(
  data = rnorm(n = n * k), nrow = n, ncol = k
)
q <- chol(
  matrix(
    data = c(1.0, 0.5, 0.5, 1.0),
    nrow = k, ncol = k
  )
)
x <- as.data.frame(z %*% q)
colnames(x) <- c("y", "x")
obj <- lm(y ~ x, data = x)
h <- hatvalues(obj)
head(1 / ((1 - h)^2))

head(gammacap_ols_hc_qcap_generic(h = h, k = 2))

jeksterslab/gammaMatrix documentation built on Dec. 20, 2021, 10:10 p.m.