2014-11-24 46 views

回答

4

什麼使它「未縮放」是它不會被sigma^2縮放,即:solve(t(X) %*% X)與係數的(縮放)方差相反:solve(t(X) %*% X)*sigma^2

如果您想要縮放的方差,即:solve(t(X) %*% X)*sigma^2,那麼您可以簡單地縮放它或僅使用vcov()來代替。一個小例子如下:

x <- 1:100 
y <- x + rnorm(100, 2) 
fit <- lm(y ~ x) 
summary(fit)$cov*summary(fit)$sigma^2 # One way 
vcov(fit) # Another way 
相關問題