$ smcsmc -Np 10
scrm 2 1 -t 8000.000000 -r 800.000000 20000000
Indirectly called
EM step 0
Particle filtering step 100% completed.
Segment data is beyond loci length
Inference step completed.
MODEL IS RESET at base 2e+07
set recombination rate 1.00507e-09(782.4 / 7.78454e+11)
popsize is equal to 9899.92 ( 1.55112e+07 / 783.4/2)
check_model_updated_Ne
Updated Ne at time 0 : | 9899.92 ( 783.4)
Forest state was created 10 times
Forest state destructor was called 0 times
Actual recombination 7824
Actual recomb op is 7.78454e+12
End of EM_step 0
###########################
# smcsmc log #
###########################
Program was compiled on: Fri-Mar-27-16:00:31-UTC-2015
smcsmc version: b98e8df8e4c84f6615c17bbc48048611906c8635
scrm version: 47202c9266494341766205e69be1381d272ced26
smcsmc parameters:
HIST saved in file: pfARGHIST
Ne saved in file: pfARGNe
Segment Data file: empty
EM steps = 0
N = 10
ESS = 0.5 (by default)
scrm model parameters:
Extract window = 0
Random seed =18446744073709551615
Sample size = 2
Seq length = 20000000
mutation rate = 1e-08
recomb rate = 1e-09
inferred recomb rate = 1.00507e-09
Pop size (at Generation):
( 0 ) | 9899.92
Actual recombination 7824
and
$ smcsmc -v
Program was compiled on: Fri-Mar-27-16:00:31-UTC-2015
smcsmc version: b98e8df8e4c84f6615c17bbc48048611906c8635
scrm version: 47202c9266494341766205e69be1381d272ced26
The output is now, something look like this:
and