Problems running non-censored test

69 views
Skip to first unread message

Alejandro Gonzalez

unread,
Jan 11, 2011, 7:50:34 AM1/11/11
to browni...@googlegroups.com
Hello,

I am trying to run a non-censored analysis in Brownie to test whether there is support for different rate of evolution of a continuous character depending on the state of a discrete trait (presence - absence). I used SIMMAP to create about 50 Stochastic character mappings of my discrete trait on the phylogeny and created a Nexus file with the continuous trait in a Data block and the SIMMAP mutation mappings in the Trees block. It appears that at least the configuration of the file is correct as I can load it into Brownie and run a One rate Brownian motion model. However when I try to run the BM model with different rate parameters for each state on the tree I get the following Error:

|gsl: lu.c:262: ERROR: matrix is singular
Default GSL error handler invoked.
Abort trap

Does this have to do with the SIMMAP mutation mappings or is this a Brownie problem?
Is there any other means of obtaining the mutation mappings other than SIMMAP? I would like to obtain a number of mappings, if possible over a subsample of trees from the posterior distribution of a BEAST analysis.

Thank you for any assistance!

Cheers

Alejandro

__________________________________

Alejandro Gonzalez Voyer

Post-doc

NEW ADDRESS

Estación Biológica de Doñana
Consejo Superior de Investigaciones Científicas (CSIC)
Av Américo Vespucio s/n
41092 Sevilla
Spain




Alejandro Gonzalez

unread,
Jan 11, 2011, 9:24:56 AM1/11/11
to browni...@googlegroups.com
Hi,

Sorry I forgot to include the file. Here it is as attachment.

Cheers

Alejandro

TestBrownie.nex
Reply all
Reply to author
Forward
0 new messages