-
Hello, big fan of sceptre so far. What is practically different from using the gRNA integration strategy "union" vs "bonferroni". I couldn't find the answer in the e-book, I apologize if it is in there and I missed it. I have this question because when I run the same sceptre analysis with all three gRNA integration strategies (data frame with 1699 discovery pairs), I get very different numbers of negative control pairs that turn up as significant. union: 45/978 negative control pairs as significant I was hoping to learn what is happening in the bonferroni integration strategy. Additionally, is bonferroni clearly the gRNA integration strategy that I should use because it produces the least amount of false positives? |
Beta Was this translation helpful? Give feedback.
Replies: 3 comments 2 replies
-
For more information on gRNA integration strategies, please read Section 2.4 of the e-book. In short, the union strategy treats cells with any gRNA targeting a given element as all the same and carries out one test comparing cells with any gRNA targeting a given element against a set of control cells, whereas the Bonferroni strategy tests each gRNA targeting a given element separately, and then combines those p-values via the Bonferroni global test. We expect the union strategy to have better power if gRNAs targeting the same element have roughly similar effects, whereas we expect the Bonferroni strategy to have better power if gRNAs have very different effects (e.g. due to some of your gRNAs not working). In our experience, we have found the union strategy to have better power more often than the Bonferroni strategy, which is why it is our default. Having said that, it appears you have some non-negligible miscalibration present in your union analysis. If you want to get a reliable answer as quickly as possible, you could choose the Bonferroni strategy since it appears to have much better calibration. Otherwise, I would recommend investing some effort in improving the calibration of the union strategy by following the suggestions in Section 5.4 of the e-book. If you were successful, you could then compare the power of the union and Bonferroni strategies if you have any positive control pairs in your data, choosing the one with the better power. (Note that comparing power in the presence of miscalibration can be misleading, because one method may have more true positives simply by virtue of having more false positives.) Otherwise, if you have no positive control pairs, I would recommend using the union strategy if you can improve its calibration. Finally, if you are unable to improve the calibration of the union strategy, then you can go ahead and use the Bonferroni strategy. |
Beta Was this translation helpful? Give feedback.
-
Thanks for the kind words. I just wanted to add that section 10.4 of the book also contains information about gRNA integration strategies that may be relevant. |
Beta Was this translation helpful? Give feedback.
-
Thank you both for the responses. I will look into fixing some of the mis-calibration issues and re-run with the union parameter |
Beta Was this translation helpful? Give feedback.
For more information on gRNA integration strategies, please read Section 2.4 of the e-book. In short, the union strategy treats cells with any gRNA targeting a given element as all the same and carries out one test comparing cells with any gRNA targeting a given element against a set of control cells, whereas the Bonferroni strategy tests each gRNA targeting a given element separately, and then combines those p-values via the Bonferroni global test. We expect the union strategy to have better power if gRNAs targeting the same element have roughly similar effects, whereas we expect the Bonferroni strategy to have better power if gRNAs have very different effects (e.g. due to some of your g…