将多个正态混合拟合到数据

[#1]

你好，

我想获得与我的数据相似的多个正态混合拟合，类似于此处显示的内容

arXiv:1306.5856

想过如何在 JMP 或 JSL 中完成吗？ 我可以为每个峰值位置设置下限/上限以换取可能更差的拟合吗？

[#2]

如果您使用 Shadowgram 选项，则分发平台中有类似的功能，然后您可以在 Shadowgram 上放置 2 和 3 混合分布

 

[#3]

我认为这可以在分发平台中完成 -> 适合普通 3 混合，请参见屏幕截图。

可能您必须将目标列转换为整数，以便分发平台中的频率角色起作用，请参见下面的示例脚本。 该脚本将生成一个表，您可以执行图形构建器和分发平台。

如果您需要拟合更复杂的曲线，您可以查看可靠性平台或自定义拟合非线性平台，后者非常强大，但很难拟合。

还有其他帖子处理这个问题：

已解决：如何在直方图上绘制单个峰拟合曲线（拟合 2 个正态混合物） - JMP 用户社区

BR

New Table( "20210331_NormalMixture",
 Add Rows( 100 ),
 New Script(
  "density vs. x",
  Graph Builder(
   Variables( X( :x ), Y( :density ) ),
   Elements( Points( X, Y, Legend( 6 ) ), Smoother( X, Y, Legend( 7 ) ) )
  )
 ),
 New Script(
  "Distribution of x",
  Distribution(
   Freq( :density ),
   Continuous Distribution( Column( :x ), Fit Normal 3 Mixture ),
   SendToReport(
    Dispatch( {"x"}, "Compare Distributions", OutlineBox, {Close( 1 )} ),
    Dispatch( {"x"}, "Quantiles", OutlineBox, {Close( 1 )} ),
    Dispatch( {"x"}, "Summary Statistics", OutlineBox, {Close( 1 )} )
   )
  )
 ),
 New Column( "x",
  Numeric,
  "Continuous",
  Format( "Best", 12 ),
  Formula( Row() ),
  Set Display Width( 53 )
 ),
 New Column( "density",
  Numeric,
  "Continuous",
  Format( "Best", 12 ),
  Formula(
   Round(
    1000000000 * Normal Mixture Density(
     :x,
     [20, 40, 70],
     [10, 10, 10],
     [0.2, 0.6, 0.2]
    )
   )
  ),
  Set Selected,
  Set Display Width( 91 )
 )
);