Displaying 1 result from an estimated 1 matches for "varlm".
Did you mean:
carlm
2012 Sep 18
1
Contradictory results between different heteroskedasticity tests
Hi all,
I'm getting contradictory results from bptest and ncvTest on a model
calculated by GLS as:
olslm = lm(log(rr)~log(aloi)*reg*inv, data)
varlm = lm(I(residuals(olslm)^2)~log(aloi)*reg*inv, data)
glslm = lm(log(rr)~log(aloi)*reg*inv, data, weights=1/fitted(varlm))
Testing both olslm and glslm with both ncvTest and bptest gives:
> ncvTest(olslm)
Non-constant Variance Score Test
Variance formula: ~ fitted.values
Chisquare = 46.88206...