similar to: Robust standard errors in logistic regression

Displaying 20 results from an estimated 2000 matches similar to: "Robust standard errors in logistic regression"

2009 Apr 13
3
Clustered data with Design package--bootcov() vs. robcov()
Hi, I am trying to figure out exactly what the bootcov() function in the Design package is doing within the context of clustered data. From reading the documentation/source code it appears that using bootcov() with the cluster argument constructs standard errors by resampling whole clusters of observations with replacement rather than resampling individual observations. Is that right, and is
2004 Sep 06
4
Cox regression for prevalence estimates
Hello, I'm an MD working in an eye clinic. I'm learning by myself to use R for use in my research works and for implementation in a software project. There are some authors who recomends the use of Cox regression as a substitute for Logistic regression (<a href="http://www.biomedcentral.com/1471-2288/3/21.pdf"> Barros AJD, Hirakata VN. BMCMedical Research Methodology, 2003;
2005 Jan 17
2
Omitting constant in ols() from Design
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
2006 Jul 05
2
p-values
Dear All, When I run rlm to obtain robust standard errors, my output does not include p-values. Is there any reason p-values should not be used in this case? Is there an argument I could use in rlm so that the output does include p-values? Thanks in advance, Celso [[alternative HTML version deleted]]
2006 Jul 07
2
Diverging results with SPSS
Dear List, I apologize in advance if this is silly. I tried to replicate an analysis I did previously in SPSS using R, and was surprised to find different results. So my question is: shouldn't the following SPSS syntax REGRESSION DEPENDENT INC89 /METHOD=ENTER hiedyrs experien SE93rec. Yeld the same results of the following R command modelB<-lm(INC89~HIEDYRS+EXPERIEN+SE93REC) I
2003 Nov 29
1
problems with xlab
Dear R-list members, I'm using the 'effects' package to plot a graph with fixed values. For some reason, I cannot change the label for the 'x' axis using 'xlab=""'. The usual commands work for the y axis and for the main title. My last attempt used the following syntax: model.effects<-(all.effects(model.1)) plot (model.effects, main="Figure
2010 May 10
2
Robust SE & Heteroskedasticity-consistent estimation
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
2009 Dec 02
1
Incorporating the results of White's HCCM into a linear regression:
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
2007 Nov 09
1
White's test again
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
2010 Dec 06
1
waldtest and nested models - poolability (parameter stability)
Dear All, I'm trying to use waldtest to test poolability (parameter stability) between two logistic regressions. Because I need to use robust standard errors (using sandwich), I cannot use anova. anova has no problems running the test, but waldtest does, indipendently of specifying vcov or not. waldtest does not appear to see that my models are nested. H0 in my case is the the vector of
2011 Jul 11
1
Robust vce for heckman estimators
When using function heckit() from package ‘sampleSelection’, is there anyway to make t-tests for the coefficients using robust covariance matrix estimator? By “robust” I mean something like if a had an object ‘lm’ called “reg” and then used: > coeftest(reg, vcov = vcovHC(reg)). I’m asking this because in Stata we could use function heckman and then use vce option “robust”. We could do the
2008 Sep 04
2
Correct for heteroscedasticity using car package
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
2008 Sep 16
3
How to do Clustered Standard Errors for Regression in R?
I can't seem to find the right set of commands to enable me to do perform a regression with cluster-adjusted standard-errors. There seems to be nothing in the archives about this -- so this thread could help generate some useful content. I've searched everywhere. Can anyone point me to the right set of commands? An example would be most helpful as well. Bo [[alternative HTML version
2008 May 08
2
poisson regression with robust error variance ('eyestudy
Ted Harding said: > I can get the estimated RRs from > RRs <- exp(summary(GLM)$coef[,1]) > but do not see how to implement confidence intervals based > on "robust error variances" using the output in GLM. Thanks for the link to the data. Here's my best guess. If you use the following approach, with the HC0 type of robust standard errors in the
2011 Jul 25
1
predict() and heteroskedasticity-robust standard errors
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
2011 Jan 01
2
robust standard error of an estimator
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]]
2006 Apr 04
1
F test for clustered data regression ?
I am using the Design library and robcov to compute variance-covariance matrices for clustered data regression. Is there an easy way to compute the F-test (i.e. linear hypothesis) for clustered data regression ? Thanks in advance! Benn
2006 Oct 06
1
Goodness of fit with robust regression
Dear list members, I have been doing robust regressions in R, using the MASS package for rlm and robustbase for logistic regressions. I must be doing something wrong, because my output does not include r-squares (or adjusted r-squares), or, in the case of glmrob, -2log likelihoods. Does anyone know how to get an output that includes these? Thanks so much for the help
2010 Jun 08
2
how to ignore rows missing arguments of a function when creating a function?
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
2013 Apr 05
1
white heteroskedasticity standard errors NLS
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