Exporting summary plm results to latex

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

David Winsemius, MD
Alameda, CA, USA