Re: [R-sig-phylo] Multivariate ASR with Discrete Characters

2020-03-03 Thread Donald Miles
Hello William,

You can use the threshBayes function in the package phytools to implement
Felsensteins threshold model. One can combine continuous and discrete
traits in the same analysis.

Cheers,

Don Miles

On Tue, Mar 3, 2020 at 11:50 AM Julien Clavel 
wrote:

> Hi William,
>
> If you want to perform a joint reconstruction of the ancestral states for
> both continuous and discrete traits (e.g. to look at the evolutionary
> correlations) I think that you should have a look at the threshold model of
> Felsenstein (a bivariate version is available in phytools; I also remember
> an Evolution meeting where E. Goolsby was working on an efficient
> implementation of the model). The model assumes that discrete traits are
> governed by a latent process which is continuous which makes possible to
> fit a multivariate model with other continuous traits (although maybe with
> a reduced power).
>
> All the best,
>
> Julien
>
> ___
> R-sig-phylo mailing list - R-sig-phylo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
> Searchable archive at
> http://www.mail-archive.com/r-sig-phylo@r-project.org/
>

[[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://www.mail-archive.com/r-sig-phylo@r-project.org/


Re: [R-sig-phylo] to log or not to log

2017-01-26 Thread Donald Miles
I agree with Ted. The comparison should be how the transformation affects
the distribution of the response variable. As Ted mentioned, you should
examine the behavior of the residuals to determine whether they conform to
the assumptions of the statistical analysis after transformation.

Donald Miles

On Thu, Jan 26, 2017 at 7:49 PM, Theodore Garland 
wrote:

> I don't think you can compare models like this for different transforms of
> the dependent variable.  The likelihood, etc., values are not comparable,
> as your results suggest.  But I am sure someone will correct me if I am
> wrong!
>
> Assuming this is some sort of regression model (i.e., you have one or more
> independent variables), then you CAN look at residuals to see if they are
> better/worse behaved (e.g., approximately normal, no nasty outliers,
> homoscedasticity).
>
> Cheers,
> Ted
>
>
> On Thu, Jan 26, 2017 at 2:25 PM, Jacob Berv  com>
> wrote:
>
> > Dear R-sig-phylo,
> >
> > When analyzing comparative data in a PCM framework, is it be appropriate
> > to use an AIC score to advocate for log transforming input data? ie
> > model(data) vs model(log(data))
> >
> > I’m running a few OUwie models and it’s not entirely clear from the
> > biology whether or not the ‘standard’ log transformation makes sense. In
> my
> > case, the model(log(data)) results have much lower AIC values relative to
> > the model(data) results.
> >
> > Cheers,
> > Jake Berv
> > ___
> > R-sig-phylo mailing list - R-sig-phylo@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
> > Searchable archive at http://www.mail-archive.com/r-
> > sig-ph...@r-project.org/
>
> [[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://www.mail-archive.com/r-
> sig-ph...@r-project.org/

[[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://www.mail-archive.com/r-sig-phylo@r-project.org/

Re: [R-sig-phylo] Package phytools load failed

2016-01-09 Thread Donald Miles
Dear All,

I receive the same error as Gabriela when I try to load phytools in R or R
Studio.

Here is the text:

library("phytools",
lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck
= vI[[j]]) :
  there is no package called ‘Biostrings’
Error: package or namespace load failed for ‘phytools’

I am running R v 3.2.3 and phytools v 0.5-10.

If I try to load Biostrings I receive an error message that Biostrings is
not compatible with R 3.2.3.

Cheers,

Don Miles
Ohio University



On Sat, Jan 9, 2016 at 12:50 PM, Vojtěch Zeisek  wrote:

> Hello
>
> Dne So 9. ledna 2016 13:11:45, Gabriela Wofkova napsal(a):
> > ​Hello,
> >
> > I have an problem with loading package "phytools". I have the newest
> > version of R studio​ and it cannot load package "Biostrings" which is
> > reqired to package "phytools". Also it not works in R program, too. Can
> you
> > have any idea, how to evade this?
> >
> > I need to use "read.newick" and load phylogenetic tree.
>
> What exactly are You doing? phytools does not require Biostrings and
> Biostrings does not require phytools. :-)
> Which operating system are You running? Which error did You get? Errors are
> keys to solutions... Do You have newest versions of the packages (I'd try
> to
> reinstall them first to see if they weren't corrupted)?
> Both load fine for me (on Linux):
>
> > library(phytools)
> Loading required package: ape
> Loading required package: maps
>
>  # ATTENTION: maps v3.0 has an updated 'world' map.#
>  # Many country borders and names have changed since 1990. #
>  # Type '?world' or 'news(package="maps")'. See README_v3. #
>
> > library(Biostrings)
> Loading required package: BiocGenerics
> Loading required package: parallel
>
> Attaching package: ‘BiocGenerics’
>
> The following objects are masked from ‘package:parallel’:
>
> clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport,
> clusterMap, parApply, parCapply, parLapply, parLapplyLB,
> parRapply, parSapply, parSapplyLB
>
> The following objects are masked from ‘package:stats’:
>
> IQR, mad, xtabs
>
> The following objects are masked from ‘package:base’:
>
> anyDuplicated, append, as.data.frame, as.vector, cbind, colnames,
> do.call,
> duplicated, eval, evalq, Filter, Find, get, grep, grepl,
> intersect, is.unsorted, lapply, lengths, Map, mapply, match, mget,
> order,
> paste, pmax, pmax.int, pmin, pmin.int, Position, rank,
> rbind, Reduce, rownames, sapply, setdiff, sort, table, tapply, union,
> unique, unlist, unsplit
>
> Loading required package: S4Vectors
> Loading required package: stats4
> Loading required package: IRanges
> Loading required package: XVector
>
> You do not need Biostrings to use read.newick.
> Note there is also read.tree in ape package reading trees in Newick format.
>
> > ​Thanks a lot for any useful answer!​
> >
> > Gabriela, Charles University in Prague​
>
> Sincerely,
> Vojtěch
>
> --
> Vojtěch Zeisek
> http://trapa.cz/en/
>
> Department of Botany, Faculty of Science
> Charles University in Prague
> Benátská 2, Prague, 12801, CZ
> http://botany.natur.cuni.cz/en/
>
> Institute of Botany, Academy of Science
> Zámek 1, Průhonice, 25243, CZ
> http://www.ibot.cas.cz/en/
>
> Czech Republic
>
> ___
> R-sig-phylo mailing list - R-sig-phylo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
> Searchable archive at
> http://www.mail-archive.com/r-sig-phylo@r-project.org/
>

[[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://www.mail-archive.com/r-sig-phylo@r-project.org/