Displaying 2 results from an estimated 2 matches for "corphyl".
2013 May 17
2
peering inside functions in a package?
...quot;phylo"))
stop("object \"phy\" is not of class \"phylo\"")
attr(value, "formula") <- form
attr(value, "fixed") <- TRUE
attr(value, "tree") <- phy
class(value) <- c("corBrownian", "corPhyl", "corStruct")
value
}
<environment: namespace:ape>
[[alternative HTML version deleted]]
2006 May 10
1
ape comparative analysis query
..., :
NA/NaN/Inf in foreign function call (arg 14)
(NB using the 'data=' argument seems to be necessary in the gls case- if
I don't, there is a further problem: 'Row names in dataframe do not
match tree tip names. data taken to be in the same order as in tree. in:
Initialize.corPhyl(X[[1]], ...)')
This seems to go away if I remove the root taxon (A1):
tree2<-drop.tip(tree, "A1")
x2<-c(B1=47,D1=43,C1=45,B2=50,B3=47,B4=48)
y2<-c(B1=2.9,D1=5.4,C1=2.8,B2=3.5,B3=3.2,B4=3.5)
dat2<-as.data.frame(cbind(x2,y2))
model2<-gls(y2~x2,data=dat2, correlation=corB...