Re: [R] implementing user defined covariance matirx

From: Thilo Kellermann <tkellermann_at_ukaachen.de>
Date: Thu 27 Jul 2006 - 01:58:53 EST

maybe the help of one of the corrStruct classes is of interest to you:

?corCompSymm
?corSymm
?corAR1
?corCar1
?corARMA
?corExp
?corGaus
?corLin
?corRatio
?corSphere

Good luck,
Thilo

On Wednesday 26 July 2006 17:34, Jonathan Smith wrote:
> 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

-- 
________________________
Thilo Kellermann
Department of Psychiatry und Psychotherapy
RWTH Aachen University
Pauwelstr. 30
52074 Aachen
Tel.: +49 (0)241 / 8089977
Fax.: +49 (0)241 / 8082401
E-Mail: tkellermann@ukaachen.de



______________________________________________ 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 03:21:39 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 - 04:19:06 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.