Menu

Strange error when conducting species delimitation for two populations

BP&P Q&A
Anonymous
2017-02-21
2017-02-21
  • Anonymous

    Anonymous - 2017-02-21

    Hi everyone,

    I'm attempting to conduct species delimitation for two species (labeled as 'east' and west'), but I continue to get the following error:
    """
    [Note: Ancestral nodes in order: 3 eastwest]

    MCMC settings: 10000 burnin, sampling every 10, 10000 samples
    Approximating posterior, using sequence data
    (Settings: cleandata=1 print=1 saveconP=1 moveinnode=1)

    Starting rjMCMC...
    PrSplit = 0.500000
    rj algorithm 0: new theta from sliding window with c = 2.00

    Starting species-delimitation model: 1

    Error: should not be here?
    """"

    I've attempted with changing the priors/etc, but I can't seem to figure this one out. Below I give my Bpp control file:

    CONTROLFILE:
    """"
    seed = -1

       seqfile = sampleX_SeqBpp.txt
      Imapfile = sampleX.Imap.txt
       outfile = sampleX_SeqBpp.txt.out
      mcmcfile = sampleX_SeqBpp.txt.mcmc
    

    speciesdelimitation = 1 0 2
    speciesmodelprior = 1
    speciestree = 0

    species&tree = 2 east west
    48 54
    (east,west);

       usedata = 1
         nloci = 2600
     cleandata = 0
    
    thetaprior =  1 100 
      tauprior = 1 10000.0
    
       finetune = 1: .01 .01 .01 .01 .01 .01 .01 .01
    
         print = 1 0 0 0
        burnin = 10000
      sampfreq = 10
       nsample = 10000
    

    """
    Thank you for the help!
    -HA

     
  • Ziheng Yang

    Ziheng Yang - 2017-04-07

    This is caused by a bug. you can check whether this is fixed in ver 3.3 or 3.3a.
    ziheng

     

Anonymous
Anonymous

Add attachments
Cancel





Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.