Skip to content

Conversation

@chrisvittal
Copy link
Collaborator

@chrisvittal chrisvittal commented Mar 2, 2022

the actual expression for x/y ploidy is:

    per_sample = coverage.transmute_cols(autosomal_mean_dp=auto_dp,
                                         x_mean_dp=x_dp,
                                         x_ploidy=2 * x_dp / auto_dp,
                                         y_mean_dp=y_dp,
                                         y_ploidy=2 * y_dp / auto_dp)

Comment on lines 374 to 376
- ``x_ploidy`` (*float64*): Estimated ploidy on X chromosome. Equal to ``x_mean_dp / autosomal_mean_dp * 2``.
- ``y_mean_dp`` (*float64*): Mean depth on calling intervals on chromosome.
- ``y_ploidy`` (*float64*): Estimated ploidy on Y chromosome. Equal to ``2 * autosomal_mean_dp / y_mean_dp``.
- ``y_ploidy`` (*float64*): Estimated ploidy on Y chromosome. Equal to ``y_mean_db / autosomal_mean_dp * 2``.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This associates as (x_mean_dp / autosomal_mean_dp) * 2, right? Could you write it 2 * x_mean_dp / autosomal_mean_dp so there's no ambiguity?

the actual expression for x/y ploidy is:

```python3
    per_sample = coverage.transmute_cols(autosomal_mean_dp=auto_dp,
                                         x_mean_dp=x_dp,
                                         x_ploidy=2 * x_dp / auto_dp,
                                         y_mean_dp=y_dp,
                                         y_ploidy=2 * y_dp / auto_dp)
```
@chrisvittal chrisvittal force-pushed the vds/fix-impute-sex-docs branch from eabd145 to 1fa9f01 Compare March 3, 2022 17:55
@danking danking merged commit 7f96d4e into hail-is:main Mar 4, 2022
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

Successfully merging this pull request may close these issues.

3 participants