[R] Exporting summary plm results to latex

Sebastian Barfort sb3730 at nyu.edu
Thu Oct 11 04:18:53 CEST 2012 

Hi David, 
this looks very promising, but I am afraid I can't see clearly how this would work. 
If possible, would you mind a short explanation? perhaps using the attached exampled?

Thank you for your advice,
Sebastian



On Oct 10, 2012, at 9:16 PM, David Winsemius <dwinsemius at comcast.net> wrote:

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