avar | R Documentation |
avar
calculated the asymptotic variance of the regression estimation
avar(X, e, group = NULL, J = NULL)
X |
A matrix or vector of independent variable(s) |
e |
The residuals from the regression |
group |
(optional) The cluster identifier (hhid from |
J |
(optional) This is assumed to be X'X, and can be input if pre-calculated |
A matrix of the covariates
hhid <- c("a","b","c","a","b","c" ,"a","b","c" ,"a","b","c" ,"a","b","c") tid <- c("1","1" ,"1" ,"2","2" ,"3","3","3" ,"4","4","5" ,"5","6","6" ,"6") w <- rep(1, 15) x1 <- rnorm(15, mean=50, sd=10) x2 <- rnorm(15, mean=50, sd=10) y <- x1 + rnorm(15, mean=50, sd=10) struc <- projdummies(hhid, tid, w) projvar_list <- projvar(x1, struc) x1p <- projvar_list$var projvar_list <- projvar(x2, struc) x2p <- projvar_list$var projvar_list <- projvar(y, struc) yp <- projvar_list$var reg <- regress1(yp, data.frame(x1p,x2p)) matCov <- avar(data.frame(x1p, x2p), reg$res, struc$hhid, reg$XX)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.