Hi all, I was wondering if anyone knew if it is possible to extract the lambda value from a compar.gee model?
Here is the pet model used in the compar.gee help file require(ape) tr <- "((((Homo:0.21,Pongo:0.21):0.28,Macaca:0.49):0.13,Ateles:0.62):0.38,Galago:1.00);" tree.primates <- read.tree(text = tr) X <- c(4.09434, 3.61092, 2.37024, 2.02815, -1.46968) Y <- c(4.74493, 3.33220, 3.36730, 2.89037, 2.30259) M<-compar.gee(X ~ Y, corStruct=corPagel(value=0.5, phy=tree.primates)) If I am correct, corPagel uses the ‘value’ argument as the initial value of lambda, but then estimates a different value during model fitting. I would like to report that estimated value. I have tried searching for lambda in str(M) and it doesn’t appear to be there. I am not sure where else to look. Many thanks for any help provided, _______________________________ Cody Dey PhD. Candidate McMaster University Department of Biology [email protected] +1 905 525 9140 ext. 24032 _______________________________________________ R-sig-phylo mailing list - [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/[email protected]/
