[R] nlme: creating pdMats

Scott Rifkin scott.rifkin at yale.edu
Sun Apr 11 05:18:33 CEST 2004

I have an experiment with an interaction (VT) between Varieties and Time.  
There are 12 Varieties and 2 Times (Early and Late), and Time is a fixed 

I have been estimating a single variance from this interaction using 

I would like to test 2 more complicated models.

1) Estimate 2 variances, one for Variety*Time(Early) and one for 
Variety*Time(Late), again with a diagonal matrix.

2) The above with a covariance term between the Times across Varieties.  
To do this, I think I would want a variance covariance matrix like (1) but
with one more parameter wherever a row and a column share a Variety but 
differ in Times.  The rest would be zero.

I would appreciate if someone would explain how I can do this using the 
standard pdMat classes, or, more likely, how to create a new one to do 

Thanks much,
Scott Rifkin

More information about the R-help mailing list