similar to: After heteroskedasticity correction, how can I get new confidential interval?

Hello there, I have a linear regression model for which I estimated heteroskedasticity-robust (Huber-White) standard errors using the coeftest function in the lmtest-package. Now I would like to inspect the predicted values of the dependent variable for particular groups and include a confidence interval for this prediction. My question: is it possible to estimate confidence intervals for the

Hi all, Does anyone know how to report heteroskedasticity and autocorrelation-consistent standard errors when using the "spm" command in SemiPar package? Suppose the original command is sp1<-spm(y~x1+x2+f(x3), random=~1,group=id) Any suggestion would be greatly appreciated. Thanks, Susan [[alternative HTML version deleted]]

Hi, I'm dealing with time series. I usually use stl() to estimate trend, stagionality and residuals. I test for normality of residuals using shapiro.test(), but I can't test for autocorrelation and heteroskedasticity. Is there a way to perform Durbin-Watson test and Breusch-Pagan test (or other simalar tests) for time series? I find dwtest() and bptest() in the package lmtest, but it

Hi, I have a simple linear regression for which I want to obtain HAC corrected standard errors, since I have significant serial/auto correlation in my residuals, and also potential heteroskedasticity. Would anyone be able to direct me to the function that implements this in R? It's a basic question and I'm sure I'm missing something obvious here. I looked up this post:

Hi all, It seems that I can get White's (HC3) test using MASS. The syntax I used for the particular problem is anova(scireg3, white.adjust="hc3") where scireg3 is an object from the lm function. But, the anova summary table is all I get. I don't get the new estimates or standard errors correcting for heteroskedasticity. Is there a way to get that information? Thanks

Hello Is there any function to calculate White's standard errors in R in an NLS regression. The sandwich and car package do it but they need an lm object to calculate the error's. Does anyone have idea how to do it for an NLS object ? Regards The woods are lovely, dark and deep But I have promises to keep And miles before I go to sleep And miles before I go to sleep ----- [[alternative

Hi all, I just discovered R a couple of days ago and I must say it rocks. I've been looking for heteroskedasticity tests and couldn't find any, however. Particularly, I've been told in one of my courses on econometrics of White's method (>< white.test()). The test's statistic is beta / sqrt(W), where W is Var(beta) "? la White", that is the beta(i) matrix is

Dear all, Trying to test for heteroskedasticity I tried several test from the car package respectively lmtest. Now that they produce rather different results i am somewhat clueless how to deal with it. Here is what I did: 1. I plotted fitted.values vs residuals and somewhat intuitively believe, it isn't really increasing... 2. further I ran the following tests bptest (studentized

Hey I am trying to run /coeftest()/ using a heteroskedasticity-consistent estimation of the covariance matrix and i get this error: # packages >library(lmtest) >library(sandwich) #test > coeftest(*GSm_inc.pool*, vcov = vcovHC(*GSm_inc.pool*, method="arellano", > type="HC3")) /Fehler in 1 - diaghat : nicht-numerisches Argument f?r bin?ren Operator/ something like:

Hi, R core team I am using the function {aidsEst} in package [micEcon] to do an AIDS model now. So far, everything is good. But I want to test the auto correlation and heteroskedasticity of the individual equation from AIDS demand system. How can I return the individual equation? PS: serial correlation test is {bgtest} in package [lmtest] and heteroskedasticity is {bptest} in package

Hi, I have ove the robust standard error of an estimator but I don't know how to do this. The code for my regression is the following: reg<-lm(fsn~lctot) But then what do I need to do? -- Charlène Lisa Cosandier [[alternative HTML version deleted]]

thank you, achim. I will try chol2inv. sandwich is a very nice package, but let me make some short suggestions. I am not a good econometrician, so I do not know what prewhitening is, and the vignette did not explain it. "?coeftest" did not work after I loaded the library. automatic bandwidth selection can be a good thing, but is not always. as to my own little function, I like the

Dear R users, I have a question here library(AER) library(plm) library(sandwich) ## take the following data data("Gasoline", package="plm") Gasoline$f.year=as.factor(Gasoline$year) Now I run the following regression rhs <- "-1 + f.year + lincomep+lrpmg+lcarpcap" m1<- lm(as.formula(paste("lgaspcar ~", rhs)), data=Gasoline) ###Now I want to find the

Dear all, Sorry if this is too obvious. I am trying to fit my multiple regression model using lm() Before starting model simplification using step() I checked whether the model presented heteroscedasticity with ncv.test() from the CAR package. It presents it. I want to correct for it, I used hccm() from the CAR package as well and got the Heteroscedasticity-Corrected Covariance Matrix. I am not

Hi all, Within a foreach loop with doSNOW, we cant call functions which come from the non-default package. We need to load(require/library) the package once more within the foreach loop. Anyone knows why would happen like this? Is it caused by the snow package and something happened when "snow" parallelize the job? Other than load the package once more with in the foreach loop, is

Hello everyone, I'm working on a current linear model Y = a0 + a1* X1 + ... + a7*X7 + residuals. And I know that this model presents both heteroskedasticity (tried Breusch-Pagan test and White test) and residuals autocorrelation (using Durbin Watson test). Ultimately, this model being meant to be used for predictions, I would like to be able to remove this heteroskedasticity and residuals

I am trying to compute the white heteroskedasticity-robust standard errors (also called the Huber standard errors) in a linear model, but I can't seem to find a function to do it. I know that the design library in S+ has something like this (robcov?), but I have not yet seen this library ported to R. Anyone know if there is already a function built into R to do this relatively simple job?

Hi, I am relatively new to R; when creating functions, I run into problems with missing values. I would like my functions to ignore rows with missing values for arguments of my function) in the analysis (as for example is the case in STATA). Note that I don't want my function to drop rows if there are missing arguments elsewhere in a row, ie for variables that are not arguments of my

Can anyone please tell me how to get the r-squared output from a piecewise (segmented) regression using the strucchange package? Here is the R code I have tried thus far. library(lmtest) library(strucchange) data <- ts(c(rnorm(30), runif(30)), frequency = 12, start = c(2005, 01)) bpts <- breakpoints(data ~ 1) print(bpts) summary(bpts) coeftest(bpts) [[alternative HTML version

Hello, I've had no luck finding an R package that has the ability to estimate a Tobit model allowing for heteroskedasticity (multiplicative, for example). Am I missing something in survReg? Is there another package that I'm unaware of? Is there an add-on package that will test for heteroskedasticity? Thanks for your help. Cheers, Alan Spearot -- Alan Spearot Department of Economics