[R] lme/gls heteroscedastic correlation structure

From: Steve Roberts <steve.roberts_at_manchester.ac.uk>
Date: Thu 16 Jun 2005 - 00:55:11 EST


We want to fit lme/gls models with a compound symmetry variance structure, but where the ICC takes different values across different groups of clusters. (Specifically a two arm clustrered trial with different ICC in each arm) Equivalent models can be fitted with heteroscedastic variance functions, but these have positivity constraints which we want to avoid and we would prefer to use the correlation structure directly. Is there a way to do this? Writing a new variance structure seems rather non-trivial - but maybe I haven't tried hard enough..

Any ideas?


  Dr Steve Roberts

Senior Lecturer in Medical Statistics,
CMMCH NHS Trust and University of Manchester Biostatistics Group, 0161 275 5192/5764 / 0161 276 5785

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Jun 16 00:58:14 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:42 EST