similar to: Robust standard errors

I am trying to use gee() to calculate the robust standard errors for a logit model. My dataset (zol) has 195019 observations; winner, racebl, raceas, racehi are all binary variables. ID is saved as a vector of length 195019 with alternating 0's and 1's. I get the following error message. I also tried the same command with corstr set to "independence" and got the same

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

This is a question about how to calculate similarities/distances among items that are classified by hierarchical attributes for the purpose of visualizing the relations among items by means of clustering, MDS, self-organizing maps, and so forth. I have a set of ~260 items that have been classified using two sets of hierarchically-organized codes on the basis of form and content. The data looks

I would like to know if there is a function that will run a linear regression and report the White (heteroscedasticity consistent) std. errors. I've found the hccm() function in the car library, but that just gives me the White covariance matrix. I'd like to be able to see the White std. errors without having to do much more work, if possible. Thanks, Brian

Using hccm() I got a heteroscedasticity correction factor on the diagonal of the return matrix, but I don't know how to incorporate this into my linear model: METHOD 1: > OLS1 <- lm(formula=uer92~uer+low2+mlo+spec+degree+hit) Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -0.0623377 0.0323461 -1.927 0.057217 . uer 0.2274742 0.0758720

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'm using maxlik with functions specified (L, his gradient & hessian). Now I would like determine some robust standard errors of my estimators. So I 'm try to use vcovHC, or hccm or robcov for example but in use one of them with my result of maxlik, I've a the following error message : Erreur dans terms.default(object) : no terms component Is there some attributes

dear R experts: I wonder whether it is possible to extend the summary method for the lm function, so that it uses an option "hccm" (well, model "hc0"). In my line of work, it is pretty much required in reporting of almost all linear regressions these days, which means that it would be very nice not to have to manually library car, then sqrt the diagonal, and recompute

Does anybody have any suggestions regarding applying standard regression packages lm(), hccm(), and others within a parallel environment? Most of the packages I've found only deal with iterative processes (bootstrap) or simple linear algebra. While the latter might help, I'd rather not program the estimation code. I'm currently using a IA-64 Teragrid system through UC San Diego.

Hi! I need to run ols regressions with Huber-White sandwich estimators and the correponding standard errors, without an intercept. What I'm trying to do is create an ols object and then use the robcov() function, on the order of: f <- ols(depvar ~ ind1 + ind2, x=TRUE) robcov(f) However, when I go f <- ols(depvar ~ ind1 + ind2 -1, x=TRUE) I get the following error: Error in

Dear R wizards: I am trying (finally) to build a function that might be useful to others. In particular, I want to create a summary.lme (extended lm) method that [a] adds normalized coefficients and [b] white heteroskedasticity adjusted se's and T's. I believe I already know how to do the programming to do these two, at least in simple unweighted cases. Now my challenges are just [1]

Hi everyone, I have tried to solve a simple problem for days but I can't figure out how to run it properly. If someone could give me a hint, this would be really great. Basically, I want to run a standard economist's fixed, and random effects regression (corresponds to xtreg in STATA) but with _variable_ weights (they correspond to changing industry shares in the market). Here is

Readers of this list might be interested in the following commenta about R. In a recent report, by Michael N. Mitchell http://www.ats.ucla.edu/stat/technicalreports/ says about R: "Perhaps the most notable exception to this discussion is R, a language for statistical computing and graphics. R is free to download under the terms of the GNU General Public License (see http://www.r-project.

luxury watches, + cheap and delivery. - &#104;ttp&#58;/&#47;&#103;oo.&#103;l&#47;7&#99;5A&#83;a &#104;ttp&#58;&#47;/go&#111;.gl&#47;&#68;8zJ&#103;4 h&#116;tp:&#47;/&#103;&#111;o&#46;g&#108;/Fm&#119;&#122;&#111;2 ozzko ky jpmuf uievj cx c edhs vdcb o itj ky jxin nisbd awtbi nid nuyzv el a beuj