WOLFRAM

Wolfram 语言与系统参考资料中心

将非参数和参数分布拟合于加权数据

    
随机生成位于到的实数，并利用 WeightedData 根据 NormalDistribution 的 PDF 对其进行加权：
Wolfram Language code: data = RandomReal[{-5, 5}, 10 ^ 4]; wts = PDF[NormalDistribution[], data];
创建一个 WeightedData 对象并计算其均值和方差：
Wolfram Language code: 𝒲 = WeightedData[data, wts]; {Mean[𝒲], Variance[𝒲]}
用最大似然估计法将数据拟合为 NormalDistribution：
Wolfram Language code: 𝒟 = EstimatedDistribution[𝒲, NormalDistribution[μ, σ]]; {Mean[𝒟], Variance[𝒟]}
用 HistogramDistribution 生成一个加权直方图：
Wolfram Language code: 𝒟NP = HistogramDistribution[𝒲]; {Mean[𝒟], Variance[𝒟]}
Wolfram Language code: Labeled[Plot[Evaluate@{PDF[𝒟, t], PDF[𝒟NP, t]}, {t, -4, 4}, Exclusions -> None, Filling -> Axis, Axes -> {True, False}, ImageSize -> 450, PlotStyle -> Thick, PlotLegends -> Placed[Style[#, Bold, FontFamily -> "Helvetica"]& /@ {"NormalDistribution", "HistogramDistribution"}, Bottom], Frame -> True, GridLines -> Automatic, GridLinesStyle -> Directive[Gray, Dotted]], Column[{Style["Parametric and Nonparametric Estimates", Bold, FontFamily -> "Helvetica", FontSize -> 18], Style["Weighted Data", Bold, FontFamily -> "Helvetica", FontColor -> Gray, FontSize -> 14]}], {{Top, Left}}]