[R] Exporting summary plm results to latex
David Winsemius
dwinsemius at comcast.net
Thu Oct 11 03:16:59 CEST 2012
On Oct 10, 2012, at 4:45 PM, Sebastian Barfort wrote:
> I am also interested in the standard errors, but beneath not next to the point estimates which is standard in the xtable package.
Last year Mark Difford offered code to do that.
From: Mark Difford <mark_difford at yahoo.co.uk>
Subject: Re: [R] Latex Table Help on R
Date: July 21, 2011 12:55:02 AM PDT
To: r-help at R-project.org help <r-help at r-project.org>
#------------------
> However, I would like the standard deviations under the means in brackets.
> Can anyone check this code to see how this can be adjusted?
Jim,
You need to use "underset," a LaTeX command. The bare-bones call is
$\underset{}{}$, where the underset value goes in the first curly and your
main value goes in the second curly (i.e. is typeset above the underset).
I don't use xtable but rather use [Ron:sic] Harrell's functions in Hmisc package,
then pass it through his latex() function, so can't take you further.
##
paste('$\\underset','{',data$SDs,'}','{',data$means,'}$', sep="")
#------------------
I sent Mark a correction at the time because it is Frank Harrell's package ... not Ron.
Frank's web page on literate program can be very useful.
http://biostat.mc.vanderbilt.edu/wiki/Main/StatReport
> If you by any chance remember the name of the package or how to do it that would be much appreciated!
>
> Cheers,
> Sebastian
>
>
> On Oct 10, 2012, at 7:10 PM, Duncan Mackay <mackay at northnet.com.au> wrote:
>
>> Hi
>>
>> If you just want the coefficients.
>>
>> xtable(summary(fe)$coef)
>> % latex table generated in R 2.15.1 by xtable 1.7-0 package
>> % Thu Oct 11 09:04:59 2012
>> \begin{table}[ht]
>> \begin{center}
>> \begin{tabular}{rrrrr}
>> \hline
>> & Estimate & Std. Error & t-value & Pr($>$$|$t$|$) \\
>> \hline
>> x & 0.12 & 0.07 & 1.78 & 0.08 \\
>> \hline
>> \end{tabular}
>> \end{center}
>> \end{table}
>>
>> There is another package whose name eludes me which may help for tables which have different outputs to the output of lm etc
>>
>> HTH
>>
>> Duncan
>>
>> Duncan Mackay
>> Department of Agronomy and Soil Science
>> University of New England
>> Armidale NSW 2351
>> Email: home: mackay at northnet.com.au
>>
>>
>>
>> At 05:09 11/10/2012, you wrote:
>>> HI,
>>>
>>> May be you can use library(texreg):
>>>
>>> library(plm)
>>>
>>> #generating some data
>>> x <- rnorm(270)
>>> y <- rnorm(270)
>>> t <- rep(1:3,30)
>>> i <- rep(1:90, each=3)
>>>
>>> data <- data.frame(i,t,x,y)
>>>
>>> fe <- plm(y~x,data=data,model="within")
>>> summary(fe)
>>> library(texreg)
>>> fe1<-extract.plm(fe) #extract the plm object
>>>
>>> library(xtable)
>>>
>>> xtable(do.call(rbind,lapply(fe1,function(x) data.frame(x))))
>>> % latex table generated in R 2.15.0 by xtable 1.7-0 package
>>> % Wed Oct 10 14:59:10 2012
>>> \begin{table}[ht]
>>> \begin{center}
>>> \begin{tabular}{rr}
>>> \hline
>>> & x \\
>>> \hline
>>> Estimate & -0.03 \\
>>> Std. Error & 0.08 \\
>>> Pr($>$$|$t$|$) & 0.68 \\
>>> R\$\verb|^|2\$ & 0.00 \\
>>> Adj. R\$\verb|^|2\$ & 0.00 \\
>>> Num. obs. & 270.00 \\
>>> \hline
>>> \end{tabular}
>>> \end{center}
>>> \end{table}
>>> #Another example. In this case, you can create two tables from the zz1 list
>>> data("Produc", package = "plm")
>>> zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, data = Produc, index = c("state","year"))
>>> zz1<-extract.plm(zz)
>>>
>>>
>>> lapply(lapply(zz1,function(x) data.frame(x)),xtable)
>>> [[1]]
>>> % latex table generated in R 2.15.0 by xtable 1.7-0 package
>>> % Wed Oct 10 15:08:02 2012
>>> \begin{table}[ht]
>>> \begin{center}
>>> \begin{tabular}{rrrr}
>>> \hline
>>> & Estimate & Std..Error & Pr...t.. \\
>>> \hline
>>> log(pcap) & -0.03 & 0.03 & 0.37 \\
>>> log(pc) & 0.29 & 0.03 & 0.00 \\
>>> log(emp) & 0.77 & 0.03 & 0.00 \\
>>> unemp & -0.01 & 0.00 & 0.00 \\
>>> \hline
>>> \end{tabular}
>>> \end{center}
>>> \end{table}
>>>
>>> [[2]]
>>> % latex table generated in R 2.15.0 by xtable 1.7-0 package
>>> % Wed Oct 10 15:08:02 2012
>>> \begin{table}[ht]
>>> \begin{center}
>>> \begin{tabular}{rr}
>>> \hline
>>> & x \\
>>> \hline
>>> R\$\verb|^|2\$ & 0.94 \\
>>> Adj. R\$\verb|^|2\$ & 0.88 \\
>>> Num. obs. & 816.00 \\
>>> \hline
>>> \end{tabular}
>>> \end{center}
>>> \end{table}
>>>
>>>
>>> Hope it helps.
>>>
>>> A.K.
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> ----- Original Message -----
>>> From: Sebastian Barfort <sb3730 at nyu.edu>
>>> To: r-help at r-project.org
>>> Cc:
>>> Sent: Wednesday, October 10, 2012 1:07 PM
>>> Subject: [R] Exporting summary plm results to latex
>>>
>>> Dear all,
>>>
>>> I am trying to export my fixed effect results to Latex. I am using the plm package with the summary function. However, it does not look like apsrtable, stargazer, or any other package can accompany using the plm package.
>>>
>>> I am interested in a classic table with the coefficient in one row followed by the standard error in paranthesis in the next row and stars by the coefficient to show relevant coefficient level.
>>>
>>> coefficient 1 xxx**
>>> (xxx)
>>>
>>> Here is a reproducible example:
>>>
>>> library(plm)
>>>
>>> #generating some data
>>> x <- rnorm(270)
>>> y <- rnorm(270)
>>> t <- rep(1:3,30)
>>> i <- rep(1:90, each=3)
>>>
>>> data <- data.frame(i,t,x,y)
>>>
>>> fe <- plm(y~x,data=data,model="within")
>>> summary(fe)
>>>
>>> If there is an alternative to using the plm package that works with any of the export to latex packages, I would be very interested to know. Otherwise, any ideas of how to solve this problem are very welcome. I almost exclusively use fixed effect panel models, and the problem of exporting results to Latex is one of the things preventing me from switching entirely from Stata to R.
>>>
>>>
>>> Kind regards,
>>> Sebastian
>>>
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
David Winsemius, MD
Alameda, CA, USA
More information about the R-help
mailing list