[R-sig-ME] FW: Estimating an "unbalanced" SUR with gls

Lehner, Stephan stephan.lehner at wu.ac.at
Wed Mar 28 17:09:18 CEST 2018

Dear Sir or Madam!

I was asked to forward my message to this mailing list (see below). For replication purposes I can share the dataset, if needed.

All the best,
Stephan Lehner

From: Douglas Bates <bates at stat.wisc.edu>
Sent: Dienstag, 27. März 2018 17:36
To: Lehner, Stephan <stephan.lehner at wu.ac.at>
Cc: Peer, Stefanie <stefanie.peer at wu.ac.at>
Subject: Re: Estimating an "unbalanced" SUR with gls

It is best to send questions like this to the R-SIG-Mixed-Models at R-project.org<mailto:R-SIG-Mixed-Models at R-project.org> mailing list. See https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models for details.

On Tue, Mar 27, 2018 at 6:44 AM Lehner, Stephan <stephan.lehner at wu.ac.at<mailto:stephan.lehner at wu.ac.at>> wrote:
Dear Prof. Bates,

I am writing you concerning some issues, I have with the correlation option in gls. For a meta-analysis about parking prices, I am trying to estimate an “unbalanced” SUR model using the gls function in R.  I estimated the SUR-model with Stata using the following approach: http://ageconsearch.umn.edu/bitstream/116272/2/sjart_st0079.pdf.

I tried the same with R using the following code lines. The data is the same. This means that the data is already “prepaired” i. e.  (just the fifth step “Fitting the SUR model” is different) .
id … the time covariate
id2 … the group covariate
Code in R:
cs1CompSymm <- corSymm(form = ~ id | id2)
cs1CompSymm <- Initialize(cs1CompSymm, data = MySUR)
gls_meta.fit.mC1 <- gls(elasticity ~ 0 + consO + consD + consQ + … , data = MySUR, method="REML", correlation=cs1CompSymm)
Code in Stata:
tsset  id2 id
xtgee  elasticity conso consd consq …, family(gaussian) link(identity) corr(unstructured) noconstant

Problem 1:
The estimated correlation matrix within-id2 in R is completely different to the correlation matrix in Stata. Hence, I am asking myself if the two methods are comparable or if I am missing something.

Estimated within-id2 correlation matrix R:

         c1       c2       c3
r1   1.0000
r2   0.4941   1.0000
r3  -0.0289  -0.1557   1.0000

Correlation Structure: General
Formula: ~id | id2
 Parameter estimate(s):
                1              2
2             0.647
3             0.655     0.066

Problem 2:
The results in R change, whether I am using id2 (an integer variable) or id22 (a String variable). Id22 is the integer variable plus a prefix “ID:”.
cs1CompSymm <- corSymm(form = ~ id | id2)
cs1CompSymm <- corSymm(form = ~ id | id22)

My data is highly unbalanced (just 8 out of 199 observations have a complete group). If is of any help, I am happy to share my dataset with you. I can also send you more information If needed.


Stephan Lehner, MSc
Research Assistant | PhD candidate
Department for Multi-Level Governance and Development
Vienna University of Economics and Business
Welthandelsplatz 1, D4, Room 2.216
1020 Vienna

+43-1-313 36-5938<tel:+43%201%20313365938>

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list