search

gchen98 / macs

Automatically exported from code.google.com/p/macs
16 stars 6 forks source link

Question about -ej option in macs #33

Open cdquinto opened 9 years ago

cdquinto commented 9 years ago

Hello again,

When the option -ej is used, what happens to the size of population j? Is it the same as specified in the command? In particular, I want to know if these 2 commands will have the same behavior. Command 1: macs 10 100000 -t 0.004 -r 0.004 -h 1e5 -I 2 5 5 -n 1 2 -n 2 2 -ej 0.0003 2 1 -en 0.000300001 1 2 -en 0.00005 1 1 Command 2: macs 10 100000 -t 0.004 -r 0.004 -h 1e5 -I 2 5 5 -n 1 2 -n 2 2 -ej 0.0003 2 1 -en 0.00005 1 1

Thank you

gchen98 commented 9 years ago

For population join events all chromosomes from the population specified in the first argument are migrated to the population specfied in the second argument.

On 07/08/2015 06:08 PM, cdquinto wrote:

Hello again,

When the option -ej is used, what happens to the size of population j? Is it the same as specified in the command? In particular, I want to know if these 2 commands will have the same behavior. Command 1: macs 10 100000 -t 0.004 -r 0.004 -h 1e5 -I 2 5 5 -n 1 2 -n 2 2 -ej 0.0003 2 1 -en 0.000301 1 2 -en 0.00005 1 1 Command 2: macs 10 100000 -t 0.004 -r 0.004 -h 1e5 -I 2 5 5 -n 1 2 -n 2 2 -ej 0.0003 2 1 -en 0.00005 1 1

Thank you

— Reply to this email directly or view it on GitHub https://github.com/gchen98/macs/issues/33.

cdquinto commented 9 years ago

In my example,

macs 10 100000 -t 0.004 -r 0.004 -h 1e5 -I 2 5 5 -n 1 2 -n 2 2 -ej 0.0003 2 1

Population 1 has size 2 before receiving the chromosomes from population 2.

After these populations are joined, does population 1 still has size 2?

gchen98 commented 9 years ago

Yes. It receives migrants from 2, but the population size is an independent parameter.

On 07/08/2015 10:03 PM, cdquinto wrote:

In my example,

macs 10 100000 -t 0.004 -r 0.004 -h 1e5 -I 2 5 5 -n 1 2 -n 2 2 -ej 0.0003 2 1

Population 1 has size 2 before receiving the chromosomes from population 2.

After these populations are joined, does population 1 still has size 2?

— Reply to this email directly or view it on GitHub https://github.com/gchen98/macs/issues/33#issuecomment-119819826.