Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Is logR and logOR corrected for purity? #73

Open
dariober opened this issue Mar 26, 2018 · 1 comment
Open

Is logR and logOR corrected for purity? #73

dariober opened this issue Mar 26, 2018 · 1 comment

Comments

@dariober
Copy link

Hi- As per title, I'm wondering whether the logR and logOR estimates are corrected by the purity estimate. I.e. as if the reported values are those one would see if the purity were 100%.

Let's make an example... You have two tumor samples (A and B) from the same tumor tissue. Assume A and B are the same in terms of genomic profile, including copy number variations, and there are no technical artifacts in library prep & sequencing.

HOWEVER, sample A has purity 40% and sample B has purity 80%. Assuming Facets correctly estimates these purities, what would be the logR (and logOR) profiles of these two samples? Would they be the same in terms of profile and absolute values? Or sample A would have half the log values of sample B since its purity is half?

Thank you!

@veseshan
Copy link
Collaborator

Hi Dario,

logR and logOR are locus level estimates from the data. They are not corrected for purity. The only correction applied is GC-correction for logR. So the sample with higher purity will show more pronounced differences in logR levels between segments and higher separation in logR. For instance, a single copy gain will show a difference in logR of log2(1.2) and log2(1.4) for 40% and 80% purity and logOR levels of log(1.4) and log(1.8) respectively. The noise within segments can be similar since that is a function of coverage. Only the copy number states (tcn and lcn) are corrected for purity.

Venkat

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants