Storey and Tibshirani 2003 method for multitest.multipletests

[Travis Wrightsman]

Dear statsmodels devs,

I wanted to correct some p-values into q-values according to the algorithm described in Remark B of Storey and Tibshirani 2003 in Python. I didn't find any implementations so I wrote one based off of WGCNA's implementation in R and the paper remark. Is this something that you all would be interested in adding as a new method to multitest.multipletests? I don't see it in the current list of methods but maybe it is very similar to an existing one.

[josef...@gmail.com]

On Mon, Jun 28, 2021 at 12:16 PM Travis Wrightsman <tw...@cornell.edu> wrote:

Dear statsmodels devs,

I wanted to correct some p-values into q-values according to the algorithm described in Remark B of Storey and Tibshirani 2003 in Python. I didn't find any implementations so I wrote one based off of WGCNA's implementation in R and the paper remark. Is this something that you all would be interested in adding as a new method to multitest.multipletests? I don't see it in the current list of methods but maybe it is very similar to an existing one.

We cannot use code that has been based on or is a translation of R code.

The R package is GPL licensed which is not compatible with our BSD license.

I cannot say anything about the method yet. I will have to read up or it. 

Maybe Kerby is familiar with it, he implemented related functions in multipletests.

Josef

 

--
You received this message because you are subscribed to the Google Groups "pystatsmodels" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pystatsmodel...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/pystatsmodels/602bf939-0e7c-4fea-9e16-15e3077d6950n%40googlegroups.com.

[Travis Wrightsman]

I see. It's not a direct translation, but even so it seems like my code
wouldn't be compatible with your license.

Jun 28, 2021 1:22:52 PM josef...@gmail.com:

>
>
> On Mon, Jun 28, 2021 at 12:16 PM Travis Wrightsman <tw...@cornell.edu>
> wrote:
>> Dear statsmodels devs,
>>
>> I wanted to correct some p-values into q-values according to the
>> algorithm described in Remark B of Storey and Tibshirani

>> 2003[https://doi.org/10.1073/pnas.1530509100] in Python. I didn't find

>> any implementations so I wrote one based off of WGCNA's

>> implementation[https://github.com/cran/WGCNA/blob/master/R/qvalue.R]

>> in R and the paper remark. Is this something that you all would be
>> interested in adding as a new method to multitest.multipletests? I
>> don't see it in the current list of methods but maybe it is very
>> similar to an existing one.
>
> We cannot use code that has been based on or is a translation of R
> code.
> The R package is GPL licensed which is not compatible with our BSD
> license.
>
> I cannot say anything about the method yet. I will have to read up or
> it. 
> Maybe Kerby is familiar with it, he implemented related functions in
> multipletests.
>
> Josef
>  
>> --
>> You received this message because you are subscribed to the Google
>> Groups "pystatsmodels" group.
>> To unsubscribe from this group and stop receiving emails from it, send
>> an email to pystatsmodel...@googlegroups.com.
>> To view this discussion on the web visit

>> https://groups.google.com/d/msgid/pystatsmodels/602bf939-0e7c-4fea-9e16-15e3077d6950n%40googlegroups.com[https://groups.google.com/d/msgid/pystatsmodels/602bf939-0e7c-4fea-9e16-15e3077d6950n%40googlegroups.com?utm_medium=email&utm_source=footer].
> --
> You received this message because you are subscribed to a topic in the
> Google Groups "pystatsmodels" group.
> To unsubscribe from this topic, visit
> https://groups.google.com/d/topic/pystatsmodels/zOiJ_8VK3NE/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to

> pystatsmodel...@googlegroups.com.
> To view this discussion on the web visit

> https://groups.google.com/d/msgid/pystatsmodels/CAMMTP%2BC-r3U1%3Di7G83jsnBgRihvNDg5Q78T0nNKLWz5VoAHQsQ%40mail.gmail.com[https://groups.google.com/d/msgid/pystatsmodels/CAMMTP%2BC-r3U1%3Di7G83jsnBgRihvNDg5Q78T0nNKLWz5VoAHQsQ%40mail.gmail.com?utm_medium=email&utm_source=footer].