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

Meaning of PosteriorPs() #75

Open
malisas opened this issue Dec 1, 2020 · 3 comments
Open

Meaning of PosteriorPs() #75

malisas opened this issue Dec 1, 2020 · 3 comments

Comments

@malisas
Copy link
Contributor

malisas commented Dec 1, 2020

Hi, I am wondering whether PosteriorPs() returns
(1) the posterior probabilities of response in the stim condition or
(2) the posterior mean of proportions for the stim condition.

The description of PosteriorPs() says it returns
The posterior probability that the samples subjected to the 'treatment', or 'stimulated', condition responded.

However, in the plot.COMPASSResult() description, the following example is given, which makes it seem like PosteriorPs() returns proportions:

## visualize the proportion of cells belonging to a category
plot(CR, measure=PosteriorPs(CR))

Which description is correct (or are they both correct and I am simply interpreting it incorrectly)?


Also, Supplementary Figure 11 of the COMPASS paper plots a “Heatmap of background corrected posterior mean of proportions (equation 11)”, which is defined as:

Background corrected posterior mean of proportions (Pi):
max(0, ¯psi − ¯pui)

I am curious as to which which function call would be used to obtain the quantity plotted in that figure. My guess is that it is a transformed call to PosteriorDiff().

Thank you for your help!

@gfinak
Copy link
Member

gfinak commented Dec 1, 2020

Thanks, Malisa
Indeed, the description is wrong, it returns the posterior proportions of cells in the stimulated samples. The posterior probability of stimulation is the mean_gamma matrix.

@gfinak gfinak closed this as completed in 1de78e4 Dec 1, 2020
@gfinak gfinak reopened this Dec 1, 2020
@gfinak
Copy link
Member

gfinak commented Dec 1, 2020

Still need to regenerate the doc.

@malisas
Copy link
Contributor Author

malisas commented Dec 1, 2020

Thanks for clarifying, Greg!

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