similar to: Results of vcovCL (sandwich) and of cluster() in Stata

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

Dear useRs, here is the announcement for the next "new" package: sandwich 0.1-3. sandwich provides heteroskedasticity (and autocorrelation) consistent covariance matrix estimators (also called HC and HAC estimators). The former are implemented in the function vcovHC() (which was available in strucchange before - and independently in hccm() in John Fox's car package). And the

Dear useRs, here is the announcement for the next "new" package: sandwich 0.1-3. sandwich provides heteroskedasticity (and autocorrelation) consistent covariance matrix estimators (also called HC and HAC estimators). The former are implemented in the function vcovHC() (which was available in strucchange before - and independently in hccm() in John Fox's car package). And the

Is it possible to get the empirical (sandwich) S.E. estimates for the fixed effects in lme () (thus allowing possibly correlated errors within the group)? In SAS you can get it by the 'empirical' option to PROC MIXED. Cheers, Michael -- Na (Michael) Li, Ph.D. Division of Biostatistics A443 Mayo Building, MMC 303 School of Public Health 420 Delaware

Dear useRs, a new version 2.0-0 of the sandwich package for estimating sandwich covariance matrices is available from the CRAN mirrors. The tools for computing heteroskedasticity (and autocorrelation) consistent covariance matrix estimators (also called HC and HAC estimators, including the Eicker-Huber-White estimator) have been generalized over the last releases from linear regression to

Dear R users: I'm happy to announce the first CRAN release of the clubSandwich package: https://cran.r-project.org/web/packages/clubSandwich clubSandwich provides several variants of the cluster-robust variance estimator for ordinary and weighted least squares linear regression models, including the bias-reduced linearization estimator of Bell and McCaffrey (2002). The package includes

Dear R users: I'm happy to announce the first CRAN release of the clubSandwich package: https://cran.r-project.org/web/packages/clubSandwich clubSandwich provides several variants of the cluster-robust variance estimator for ordinary and weighted least squares linear regression models, including the bias-reduced linearization estimator of Bell and McCaffrey (2002). The package includes

Hi All, I have a dataset, longitudinal in nature, each row is a 'visit' to a clinic, which has numerous data fields and a count variable for the number of 'events' that occurred since the previous visit. ~50k rows, ~2k unique subjects so ~25 rows/visits per subject, some have 50 some have 3 or 4. In STATA there is an adjustment for the fact that you have multiple rows per

Hi, I would like to estimate coefficients using poisson regression and then get standard errors that are adjusted for heteroskedasticity, using a complex sample survey data. Then I will calculate prevalence ratio and confidence intervals. Can sandwich estimator of variance be used when observations aren?t independent? In my case, observations are independent across groups (clusters), but

Dear R-help, I have a four-part question about regression, matrices, and sandwich package. 1) In the sandwich package, I would like to better understand the meat() function. >From the bread() documentation, for a simple OLS regression, bread() returns (1/n * X'X)^(-1) That is, for a simple regression (per the documentation on bread()): MyLM <- lm(y ~ x) bread(MyLM)

Dear List, I'm just teaching myself semi-parametric techniques. Apologies in advance for the long post. I've got observational data and a longitudinal, semi-parametric model that I want to fit in GAM (or potentially something equivalent), and I'm not sure how to do it. I'm posting this to ask whether it is possible to do what I want to do using "canned" commands

Dear R-Users: I am trying to find the robust (or sandwich) estimates of the standard error of fixed effects parameter estimates using the package "lmer2". In model-1, I used "robust=TRUE" on the other, in model-2, I used "robust=FALSE". Both models giving me the same estimates. So my question is, does the robust option works in lmer2 to get the robust estimates of

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

Hi all, I installed SuSe 9.1 with Samba 3.0.4-suse. I can not get the smbd to start from swat, but nmbd is running, When I start smdb and nmbd with the SWAT restat then ps -A shows that the processes are there and it seems to run, but not according to the 'status' of swat. Copying files from samba and printing with the shared printers work though. The printers show "access denied, no