Julia的Bootstraping加权函数

Question

我试图使用Bootstrap.jl函数来获得加权函数(例如加权中值)的标准误差(SE)。

请参阅下面的Bootstrap.bootstrap代码，以获得未加权中值的SE。

using StatsBase, DataFrames, Bootstrap

v = collect(1:1:20)

bootstrap(median, v, BasicSampling(100))

现在，我需要向上面的median传递第二个参数，以获得加权median的SE。在bootstrap函数之外，如下所示：

w = collect(0.1:0.1:2)
median(v, Weights(w))          

如何将第二个参数传递给median函数在bootstrap中包含权重？请注意，引导带重采样应该应用于这两个向量，为它们绘制相同的索引。

Answer 1

您可以将包含两个向量的DataFrame传递给bootstrap的第二个参数。然后编写一个匿名函数来使用median中的每个列。例如。

df = DataFrame(v = collect(1:1:20),
          w = collect(0.1:0.1:2))

bootstrap(d -> median(d[!,:v], Weights(d[!,:w])), df, BasicSampling(100))