Filter nan before qvalue calculations about fsharp.stats

[bvenn]

quickfix for reference:

let getQvalueFromNanCorrupted (inputPvals:float[]) = 

    let pvalValid =
        inputPvals 
        |> Array.indexed   
        |> Array.filter (fun (i,p) -> not (nan.Equals p))

    let indicesOfNan = 
        inputPvals 
        |> Array.indexed   
        |> Array.filter (fun (i,p) -> (nan.Equals p))
	// if all pvalues are none, just return the nans
    if pvalValid.Length = 0 then 
        nan,indicesOfNan |> Array.map snd

    else 

        let pi0 = Testing.MultipleTesting.Qvalues.pi0Bootstrap (pvalValid |> Array.map snd)
        let qv = Testing.MultipleTesting.Qvalues.ofPValuesBy pi0 snd pvalValid

        let aggChunks =
            Array.map2 (fun (i,pval) qval -> i,qval) pvalValid qv
            |> Array.append indicesOfNan
            |> Array.sortBy fst
            |> Array.map snd
        pi0,aggChunks 

[bvenn]

closed by 5c95836