smshuai / U1-snRNA

GNU General Public License v3.0
1 stars 2 forks source link

Typo in beta_binomial_model.py? #1

Open stephan-hutter opened 4 years ago

stephan-hutter commented 4 years ago

Line 191 looks like this: Ns_n = normal_panel.loc[(row.chrom, row.pos, row.alt), 'dpP'].values

I assume you are looking at total depth on the negative strand. Shouldn't the 'dpP' then be 'dpN'?

smshuai commented 4 years ago

You are right about this. This shall affect both normal and tumour scores in the negative strand the same way. I need to check what exactly is affected. Thanks again for pointing this out!