Questions tagged [sandwich]
Sandwich, or sandwich variance estimation, refers to a method of estimating standard errors from estimating equations that is robust to many model based assumptions. The preferred tag is "robust-standard-error"
45 questions
5
votes
1
answer
241
views
Does it mean that we don't need a normal assumption for using sandwich estimator in normal linear regression?
According to this post,
the blogger uses the theory of estimating equations to construct the robust sandwich variance estimator.
In this post, it said that:
Now we ...
- 344
4
votes
2
answers
493
views
Sandwich variance estimator or bootstrap-based variance for stabilized inverse probability weighting (IPW)
Multiple published papers describe IPW as akin to having population with multiply copies of the same individuals. Hence, the correlation should be accounted and corrected using sandwich variance ...
- 905
0
votes
0
answers
15
views
QMLE and Confidence Intervals
In the setting of Huber's paper, the true density is $y\sim f_0$, and the density under consideration is $y\sim f(\cdot|\theta)$. If we choose to maximize $y\sim f(\cdot|\theta)$, we will find that $\...
- 387
0
votes
0
answers
17
views
Diagnostic checks before and after adjusting for standard errors in R
Currently, I'm fitting a regression model for my dataset. As there are clusters, I refitted the model using the coeftest() function from the ...
- 227
3
votes
1
answer
131
views
Can a subject appear multiple times in a risk set in recurrent event analyses?
To fit a Cox model to recurrent event data (Andersen & Gill), using the R survival package, requires the user to cast the data into counting process format (see [1]).
For recurrent event analyses, ...
- 143
3
votes
1
answer
30
views
How can I tell if a clutser-randomised crossover trial has made a unit of analysis error?
I am studying the following paper:
https://jamanetwork.com/journals/jama/fullarticle/2698491
This is a cluster randomised control trial with crossover. I want to ensure they have not made a unit of ...
- 503
1
vote
0
answers
71
views
Sandwich Package is not Working Properly [closed]
The sandwich package in R is not working properly. All of the following methods are giving the same results:
Regular standard error
Robust (HC1)Standard
One way ...
- 11
2
votes
0
answers
113
views
Does the sandwich-estimate eliminate heteroskedasticity in a model residuals versus fitted plot, or simply make the estimation robust to heteroskedas? [closed]
Does the sandwich-estimator/Huber-White/GEE eliminate heteroskedasticity in a residuals versus fitted plot, or simply make the model robust to it?
user318514
1
vote
0
answers
132
views
If the normality assumption in the for the GLS estimation fails, would you switch to GEE?
I want a marginal model, ideally fit via GLS. But the normality of residuals doesn't hold. It isn't much skewed, I don't want any transformations. It's just non-normal in shape. Yet still reporting ...
- 41
1
vote
0
answers
93
views
M-estimator: There is no "of something" in the definition
I see that when talking about estimator, we have "of something", where "something" refers to a fixed parameter. For example, we say that the sample mean is an estimator of the ...
- 862
1
vote
1
answer
1k
views
what does generalized estimating equations with robust standard error mean
I am currently working on logistic regression and came across some articles stating "generalized estimating equations with robust standard error" or "with robust sandwich estimators for ...
- 361
0
votes
1
answer
761
views
How to compute sandwich estimator with QMLE and Poisson regression in R (using glmrob package)
I am trying to applying a sandwich estimator to a Poisson regression with QMLE in R, where I used the glmrob function. The code I am using looks like:
...
- 1
3
votes
2
answers
860
views
How to compute the sandwich variance ML estimator in R
I'm currently estimating a DCC-type model by maximum likelihood. Im using the command solnp and it return an object where I can compute the Hessian H evaluated at ...
1
vote
1
answer
242
views
Decision between vcovPC and vcovPL (sandwich)
I want to do a linear probability model with clustered errors. The data also has a panel structure.
In the R package “sandwich,” there are two functions: vcovPC() ...
- 125
6
votes
0
answers
823
views
r quantreg - quantile regression with clustered standard errors
I fit a quantile regression using quantreg:::rq on clustered data. I use the Huber sandwich estimator to obtain cluster-corrected standard errors, which is ...
- 225
2
votes
1
answer
597
views
Robust SE clustered GLM Gamma Log Link to match GEE Robust SE
How do I get the robust standard errors/sandwich variance estimators for GLM using a Gamma family with a log-link to match the robust standard errors from the GEE output?
...
1
vote
1
answer
51
views
How to determine if correlated errors are a result of using an incorrect functional form?
Assume you have identified that the errors of a regression model are correlated. How should you determine if this is a sign of using incorrect an functional form? That is, a situation that should be ...
1
vote
0
answers
145
views
Why cannot the sandwich SE be used when the Kenward-Roger denominator df are in use?
This question is about mixed models, and the use of both the Kenward-Roger correction (K-R) for small sample sizes and robust (sandwich) estimators for standard errors in the same model. Is it ...
- 205
2
votes
1
answer
124
views
Variance Estimation for Least Squares with Probability Weights
I'm running a simulation study and finding that the nominal SEs of the estimated coefficients when using weights in lm in R are an underestimate of the simulation SE. I have confirmed that $\hat{\beta}...
- 497
2
votes
0
answers
176
views
Calculating sandwich estimator
Considering design matrix $X \in \mathbb{R}^{n\times p}$ $(n>p)$ and response $y\in \mathbb{R}^{n}$. The sandwich estimator can be calculated directly using
$$(X^TX)^{-1}X^T diag(r^2) X (X^TX)^{-...
- 75
1
vote
0
answers
464
views
Robust regression with Sandwich estimator
I understand that rlm (robust regression) addresses issues of outliers and influential observations, but does not address heteroskedasticity. I have come to learn ...
- 715
3
votes
2
answers
1k
views
How to estimate robust sandwich standard errors when estimating parameters using optim() in R?
Currently I am using numerical optimization in R via the optim() function to estimate some parameters in a complicated ...
- 369
3
votes
1
answer
2k
views
How to calculate sandwich standard errors for generalized least squares models?
Dependent data can be modeled using covariance structures like compound symmetry, spherical, AR-1, and other. Using generalized least squares, inference can be made on the regression coefficients ...
- 64.8k
1
vote
1
answer
2k
views
Wondering what type of covariance correction for standard errors is better: Hansen-Hodrick or Newey-West?
I am wondering what type of covariance correction for standard errors is better: Hansen-Hodrick or Newey-West?
Also, does someone know if StatsModels package that ...
- 361
1
vote
0
answers
992
views
Different optimal bandwidths of Newey West (1994) in R and STATA
R and STATA gave very different optimal bandwidths for the same data set. It will be greatly appreciated if someone can give me any hint why this happens. Here are two sample codes from R and STATA ...
- 121
6
votes
1
answer
9k
views
Standard Errors with Weighted Least Squares Regression
For OLS, $\hat{\beta} = (X'X)^{-1}X'y$, and $\text{var}(\hat{\beta}) = (X'X)^{-1} X' \sigma^2 I X (X'X)^{-1}$. I can reproduce these "by hand".
For WLS, with heteroskedastic errors and weights in ...
- 1,197
1
vote
1
answer
2k
views
What does setting vcov=sandwich result in? I was expecting that it will result in HC consistent standard errors [closed]
However, when I checked the results using, vcovHC=sandwich, I found different results. Can anyone help me to figure out the difference?
...
- 11
2
votes
0
answers
851
views
Multinomial tests of heterogeneity using robust standard errors?
I am looking at differences in frequencies of categorical variables collected at different sites, adjusting for stratification factors of age and sex. I would like to use robust standard error ...
- 64.8k
5
votes
1
answer
1k
views
Difference between Quasi-Poisson and Sandwich Covariance
I understand that both methods can be utilized to obtain correct inference in overdispersed Poisson data. What I don't understand is the difference between them: why the analyst would choose one over ...
4
votes
1
answer
347
views
How to implement the sandwich estimator in a semi-parametric situation?
I am trying to implement a sandwich estimator described in Zhang et al. (2012, p. 1012) in very brief terms. The information they give is not enough for me to understand what has been actually done, ...
- 6,704
0
votes
1
answer
122
views
Reference for incremental sandwich covariance from biglm?
I am working on some similar methods to Lumley's biglm wrapper around Miller's AS274 algorithm, and I can't seem to find a reference for his incremental Huber/White ...
- 155
1
vote
0
answers
79
views
I can't correct the OLS model with heteroskedasticity by the lmtest means
i'm using a selvaggio model to explain the behavior of deposits in a bank's data, and i need to use the estimated parameters,
the problem is the heteroskedasticity that i detectect with breusch-pagan ...
4
votes
1
answer
2k
views
Clustered (grouped) standard errors MLE in R
I'm doing the following maximum likelihood estimation using mle2 function from bbmle package:
...
- 305
1
vote
0
answers
366
views
Application of Huber-White Variance Estimates in GLMER
I'm currently working on an analysis in R using GLMER mixed-effects model with a logistic regression framework under the lme4 package. I would like to include empirical (Huber–White sandwich) variance ...
- 369
2
votes
1
answer
213
views
effect estimate without constant variance (r sandwich vcovHC)
For the sandwich package in R, when using vcovHC(model) where model <- lm(Y~X), what gives the effect estimate regarding the ...
4
votes
0
answers
2k
views
Sandwich Estimator in Maximum Likelihood Estimation of Logit
I am estimating a discrete choice model using mixed logit using Halton Draws. So everything is effectively done with MCMC. The code is written in MATLAB. I am using MATLAB's ...
- 69
3
votes
0
answers
804
views
How to calculate robust standard error with offset?
I think the default vcovHC in R's sandwich package does not handle offsets in Poisson models. We see this because robust (...
- 64.8k
2
votes
0
answers
237
views
Why sandwich estimators aren't always used in OLS regression?
I asked before what is the intuition behind sandwich estimators. I must still missing something because I don't understand why sandwich estimators are not always applied to OLS residuals.
Can you ...
- 4,588
3
votes
1
answer
3k
views
Robust OLS versus ML with sandwich estimator
If you compare the standard errors of the OLS coefficients with the White correction, versus the ML estimates with the variance estimated with the sandwich estimator, which standard errors do you ...
- 3,489
2
votes
0
answers
629
views
Do robust standard errors protect you from proportional odds assumptions?
Cox Proportional Hazards models are traditionally taught alongside proportional hazards assumptions. There is a corresponding test of proportionality. However, if standard errors are calculated from ...
- 64.8k
1
vote
2
answers
1k
views
Clustered data WITHOUT multilevel / GEE model?
I have a data-set with around 700 observations from 12 centres. Although the clustering effect as tested in a random intercept model didn't seem significant, it seems more appropriate to use a ...
- 125
0
votes
1
answer
525
views
Huber sandwich estimator in quantile regression
I need the description of Huber sandwich estimate method for quantile regression.
I found this "a Huber sandwich estimate using a local estimate of the sparsity function". Sparsity function looks ...
- 9
6
votes
1
answer
3k
views
Pairwise comparisons for a regression with sandwich estimates (in R)
The question in short
I run a regression in R and made a boxplot of the response variable with grouping by one of the predictor variables. On this boxplot I'd like to add some information about the ...
- 5,182
11
votes
2
answers
15k
views
Robust regression inference and Sandwich estimators
Can you give me an example of the use of sandwich estimators in order to perform robust regression inference?
I can see the example in ?sandwich, but I don't quite ...
- 5,182
5
votes
1
answer
1k
views
Robust Residual standard error (in R)
I have a question regarding to the concept of robust standard errors. What I found about that topic is, that one can estimate the robust standard error for regression coefficients to eliminate ...
- 193