similar to: Latex outputs of multilevel models

Hello, I've got a question concerning the usage of robust standard errors in regression using lm() and exporting the summaries to LaTeX using the memisc-packages function mtable(): Is there any possibility to use robust errors which are obtained by vcovHC() when generating the LateX-output by mtable()? I tried to manipulate the lm-object by appending the "new" covariance

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 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

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)

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

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

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 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]]

Hola. Estoy utilizando lyx con el módulo de knitr y tengo un problemilla con la función toLatex del paquete memisc. Pongo un ejemplo mejor. En Rstudio lo puedo hacer como sigue en un fichero Rnw. y la tabla en el pdf aparece alineada en el pdf. \documentclass{article} \usepackage{booktabs} \usepackage{dcolumn} \begin{document} <<>>= library(memisc) X1 <- rnorm(1000) X2 <-

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 <-

It is not a pure R question,but I hope some one can give me advices. I want to use analysis my data with the multilevel model.The data has 2 levels---- the second level has 52 units and each second level unit has 19-23 units.I think the sample size is quite small,but just now I can't make the sample size much bigger.So I want to ask if I use the multilevel model to analysis the data set,will

Dear all, I am struggling with a problem which I have been reading on the forums about and it did not seem to me that there is a precise answer to my question. However, I still hope there is one. I am working with http://timss.bc.edu/ PIRLS data and trying to conduct multilevel analysis. There are different weights for each level of analysis in the PIRLS dataset (e.g. there is a school

Dear All, What is the preferred way to get Bayesian analysis results (such as those from MCMCpacki, MCMCglmm, and DPpackage) into LaTeX table automatically? I have been using the "apsrtable" package and similar functions in "memisc" package, but neither seems to handle MCMC output directly. Many thanks. Shige

Dear ListeRs, I am trying to fit a heteroscedastic multilevel model using lmer{lme4- package). Take, for instance, the (fictive) model below. lmer(test.result ~ homework + Sex -1 + (1 | School)) Suppose that I suspect the error terms in the predicted values to differ between men and women (so, on the first level). In order to model this, I want the 'Sex'-variable to be random on

Goodmorning everybody, i'm an italian statistician and i'm using R for research. Could someone tell me some indices to see the goodness of fit in multilevel modelling? I'm using the lmer function, and I want to know if my model fit well my data. I actually want to justify the use of multilevel model instead the classical one. Hope someone can help me. Thank you. Greetings

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) ~

hello everybody i was seeking a converter beetween R regression output (eg with summary) and the conventional way to present regression output in paper: every model as a vertical vector with \beta, t beetween parenthesis below the first, and other statistics (R^2 etc) . I've seen hmisc and xtable, and if I didn't miss something, they don't include something like that. Thank you

Hi, I'm trying to reshape my data set from wide to tall format for multilevel modeling. Unfortunately, the function I typically use (make.univ from the multilevel package) does not appear to work with unbalanced data frames, which is what I'm dealing with. Below is an example of the columns of a data frame similar to what I'm working with: ID a1 a2 a4 b2 b3 b4 b5 b6 Below

I am using a multilevel modeling approach to model change in a person's symptom score over time (i.e., longitudinal individual growth models). I have been using the lme function in the multilevel package for the analyses, but my problem is that my outcome (symptoms) and one of my predictors (events) are count data, and are non-normal. Do you have any suggestions for how to deal with them?

Hi, does anyone out there have a recommendation for multilevel / random effects and longitudinal analysis? My dream book would be something that's both accessible to a non-statistician but rigorous (because I seem to be slowly turning into a statistician) and ideally would use R. Peter