search for: corbrownian

Let's say I would like to look inside the function corBrownian in library (ape). When I type in the function name I get the following, which is not nearly the detail that goes into this function. I am wondering how to begin cracking this function open (and others) so I can learn more about it and perhaps code my own corClass one day. Thanks. > corBrowni...

...1=48) y<-c(B1=2.9,D1=5.4,C1=2.8,B2=3.5,B3=3.2,B4=3.5,A1=1.8) picx<-pic(x,phy=tree) picy<-pic(y,phy=tree) Which seems to make sense. However, if I try to use anything more sophisticated I run into problems: >dat<-as.data.frame(cbind(x,y)) >model1<-gls(y~x,data=dat, correlation=corBrownian(phy=tree)) Error in corFactor.corStruct(object) : NA/NaN/Inf in foreign function call (arg 1) or equally: > model2<-compar.gee(y~x,data=dat, phy=tree) [1] "Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27" [1] "running glm to get initial regression estimate" [1...

I have spatial data on a sphere (the Earth) for which I would like to run an gls model assuming that the errors are autcorrelated, i.e. including a corSpatial correlation in the model specification. In this case the distance metric should be calculated on the sphere, therefore metric = "euclidean" in (for example) corSpher would be incorrect. I would be grateful for help on how to