similar to: memisc-Tables with robost standard errors

Hi Folks, I've got a question regarding the 'quantreg' package maintained by Roger Koenker: I tried to produce LaTeX tables using the longtable and dcolumn options as decribed in the manual, but the function latex() doesn't seem to react on _any_ other options than 'digits' and 'transpose'. To reproduce these results the following minimal example may be used:

Hello Everyone (Quick) question: Does anyone know a package/command or simply a way of creating a table of different panel data estimations (estimated using /*plm()*/ ) just as *mtable()* does for models estimated with /*lm()*/? It seems *mtable* (and *apsrtable* equally) only support /*lm*/ and some other classes but unfortunately not /*plm*/. I am pretty sure others must have encountered this

Hi Folks, I've got a question concerning the calculation of the Schwarz-Criterion (BIC) done by summary.regsubsets() of the leaps-package: Using regsubsets() to perform subset-selection I receive an regsubsets object that can be summarized by summary.regsubsets(). After this operation the resulting summary contains a vector of BIC-values representing models of size i=1,...,K. My problem

Hi there, I have a dataset from a longitudinal study with a lot of drop-out. I want to implement the inverse probability weighting method by Robins 1995 JASA paper "Analysis of semiparametric regression models for repeated outcomes in the presence of missing data". Does anyone know if there is a package to do it in R (or other software)? Thanks a lot! Lei

Hello, first of all: I'm new to R and have only used SPSS befor this (which can't do this at all...). I'm trying to output some regression results to latex. The regressions are normal OLS and I'm trying to output the results with standardized regression coefficients and t-statistics based on "Huber-White sandwich estimator for variance". The final result should be

Hi, I have a number of multilevel models and I would like some Latex outputs of them. I usually use the "apsrtable" package, but it does not accept "lme" outputs. Neither does the "mtable" function in the "memisc" package. Is there any good alternative that I am missing? Thanks, Jonas

Hi - I have a series of 9 GAM regressions with about 5 parametric effects and three non-parametric effects in each. What is a good library or command for turning GAM outputs into publication-ready regression tables? I tried apsrtable and the mtable command in memisc but neither seemed to work with the gam output. I'd be okay with two separate tables - one for the parametric effects and

Hello, I'm using memisc to output regression results to tables and latex. My problem is that the output that Latex needs must be in between $ $ so that it is read as formula but memisc does not output the result between $ $. For example, latex needs: $0.05^{***}$ and memisc outputs 0.05^{***} in an entry. I am new to Latex and I imagine it could also be a latex 'problem' and not

Dear R-Users, I want to use mtable from package "memisc" to produce Latex-style estimation output. However, mtable() only gives me a p-value and not the corresponding test-statistic. Does anyone know how to extract it, either from a glm/anova object or mtable? Here is a short example: # Run this #################### install.packages("memisc") library(memisc) set.seed(1)

Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG <- glm (low ~ X, family = "binomial", data = DATOS) Ejemplo: library(MASS) data(birthwt, package="MASS") birthwt$low <- factor(birthwt$low) birthwt$race <- factor(birthwt$smoke) REG_LOG <- glm (low ~ smoke, family =

Greetings, Is there a way to control the number of digits after the decimal from linear regression output using the memisc package? I have tried the following code, but it does not work: fm <- lm(y ~ X) mtable(fm, digits=9) The default seems to be 3 digits after the decimal. Also, is there a way in general to control the number of digits after the decimal? Thank you! Geoff

Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Rejected)~Gender,data=berkeley,family="binomial") berk2 <-

Hi, I want to export the result of glht in R into a LaTeX table, such as that result: Linear Hypotheses: Estimate Std. Error z value Pr(>|z|) Group1 - Group2 == 0 -0.14007 0.01589 -8.813 <0.001 "***" Group1 - Group3 == 0 -0.09396 0.01575 -5.965 <0.001 *** --- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05

Hi all, Does anyone know a way to change the significant stars in mtable (package memisc)? The default is Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1, however I need it to be Signif. codes: 0 '***' 0.01 '**' 0.05 '*' 0.1 ' ' 1 Kind regards, Erich [[alternative HTML version deleted]]

Hi all, I have used the apsrtable package to generate tables (using LaTeX and Sweave) for quite some. However, suddenly the notes option in the package appears to have stopped working (or I am doing something wrong). So, when I try to run the following commands (in the R console): > attach(NatExp.df) > model1.lm <- lm((1-schoolspend) ~ jadif) > model2.lm <- lm((1-nohousdisc) ~

Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Rejected)~Gender,data=berkeley,family="binomial") berk2 <-

Dear All, I am trying to recode a variable using the functions provided by "memisc" package. Actually I am following the examples on page 9-10 of the vignette: -------------------------------------------------------------------------- d.fig <- within(d.fig,{ sev <- recode(sev, 1 -> 0.9, 2 -> 1.0, 3 -> 1.1, 4 ->

Dear all, how can I use vcovHC() to get robust/corrected standard errors from an arima() output? I ran an arima model with AR(1) and got the estimate, se, zvalue and p-value using coeftest(arima.output). However, I cannot use vcovHC(arima.output) to get corrected standard errors. It seems vcovHC works only with lm and plm objects? Is there another way I can get robust/corrected

Dear useRs, I am pleased to announce the availability of version 0.95 of package 'memisc' on CRAN (http://cran.r-project.org/web/packages/memisc/). 'memisc' does not implement any new estimators, but focuses on providing an infrastructure for data analysis especially of survey data. It may be of interest to all those who have had to use a commercial package like SPSS or Stata

Dear useRs, I am pleased to announce the availability of version 0.95 of package 'memisc' on CRAN (http://cran.r-project.org/web/packages/memisc/). 'memisc' does not implement any new estimators, but focuses on providing an infrastructure for data analysis especially of survey data. It may be of interest to all those who have had to use a commercial package like SPSS or Stata