[R] implementing user defined covariance matirx

From: Jonathan Smith <jon3smith_at_hotmail.com>
Date: Thu 27 Jul 2006 - 01:34:06 EST


I am trying to implement my own covariance matrix into R. Then be able to use it in gls or lme or nlme to analyze some data. I simply want to use corr=mymodel(form ~tim/peep) same as can use corr= corAR1 and many others. Having a terrible time trying to figure out how to do this. I have found documentation saying that you can do this but cant find out how. Any suggestings would be greatly appreciated.

My covariance matrix is:

| 1 s1^g s2^g s3^g |
| s1^g 1 s1^g s2^g |
| s2^g s1^g 1 s1^g ... |
| s3^g s2^g s1^g 1 |
| : : : : : :: |
| : : : : :::: |

Thankyou kindly for your time
Jonathan Smith



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 and provide commented, minimal, self-contained, reproducible code. Received on Thu Jul 27 05:04:41 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 27 Jul 2006 - 06:16:17 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.