2$sigma$Pinv) # partial correlations
Best wishes,
Julien
De : R-sig-phylo de la part de
Oliver Betz
Envoyé : dimanche 13 juin 2021 14:45
À : r-sig-phylo@r-project.org
Objet : [R-sig-phylo] phylogenetic correlation analysis
Dear all:
I would like to perform a phylogenetic corr
Rphylopars is an option too (https://github.com/ericgoolsby/Rphylopars)
[[alternative HTML version deleted]]
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http:
Whether one gets them from PGLS directly or from contrasts, one can
get correlations by just inferring the covariance matrix and then
calculating
r(x,y) = Cov(x,y) / (Var(x) Var(y))^(1/2)
where of course Var(x) is also Cov(x,x), and so on.
You would not need a separate run to get correlations.
J
nificant”
> partial correlations from the model fit:
>
> fit_bm2 <- mvgls(data~1, tree=tree, model="BM", method="PL",
> penalty="LASSO")
> cov2cor(fit_bm2$sigma$Pinv) # marginal correlations
> cor2pcor(fit_bm2$sigma$Pinv) # partial correlations
>
Envoyé : dimanche 13 juin 2021 14:45
À : r-sig-phylo@r-project.org
Objet : [R-sig-phylo] phylogenetic correlation analysis
Dear all:
I would like to perform a phylogenetic correlation analysis (simlar to
PGLS, but correlation instead of regression), so that I get a
correlation matrix, whe
Dear all:
I would like to perform a phylogenetic correlation analysis (simlar to
PGLS, but correlation instead of regression), so that I get a
correlation matrix, where all the Pearson or Spearman correlation
coefficients between all of my variables are listed. One solution
might be to ca