Nothing
phylogeneticMean <- function(traits, phy, lambda = 1){
if (!is.null(names(traits)))
traits <- traits[phy$tip.label];
if (inherits(phy, 'phylo')) {
vmat <- vcv.phylo(phy);
}else{
vmat <- phy;
}
dd <- diag(vmat);
vmat <- vmat * lambda;
diag(vmat) <- dd;
onev <- matrix(rep(1, length(traits)), nrow=length(traits), ncol = 1);
anc <- as.vector(solve(t(onev)%*% solve(vmat) %*% onev) %*% (t(onev)%*% solve(vmat) %*% traits));
beta <- as.vector((t(traits-anc) %*% solve(vmat) %*% (traits-anc))/length(traits));
return(list(anc = as.vector(anc), beta = as.vector(beta)));
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.