Julia中的等效np.nanquantile？答案

【问题标题】：Equivalent np.nanquantile in Julia?Julia中的等效np.nanquantile？
【发布时间】：2021-07-07 07:26:21
【问题描述】：

在 Julia 中是否有 np.nanquantile 等价物？我有一个二维数组并沿一个轴计算分位数，但该数组包含NaN-values。我当前的代码块：

quantiles = Array{Float32}(undef, size(array, 1), 2)
p=0.1
quantiles[:, 1] = mapslices(x -> quantile(x, p), array, dims = 2)
quantiles[:, 2] = mapslices(x -> quantile(x, 1 - p), array, dims = 2)

【问题讨论】：

请确认您的数据中是否有NaN 或missing 值（因为这会影响答案）。您是如何在 array 中获得 NaN 值的？
@BogumiłKamiński 我有 NaN 值，它们来自 Python 计算。

标签： arrays numpy statistics julia

【解决方案1】：

最简单的做法是使用以下内容：

x -> quantile(filter(!isnan, x), p)

例如

julia> array = [1   NaN 3 4
                NaN 2   3 4]
2×4 Matrix{Float64}:
   1.0  NaN    3.0  4.0
 NaN      2.0  3.0  4.0

julia> mapslices(x -> quantile(filter(!isnan, x), 0.5), array, dims = 2)
2×1 Matrix{Float64}:
 3.0
 3.0

【讨论】：