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