567691_zh-TW

liliconcom

使用 Bonferroni 限值自動識別離群值

各位，

 

我對 Bonferroni 限制有疑問。 我的目標是開發一個數據評估流程，其中可以使用 Bonferroni 極限檢測異常值。 我的數據集由發布測試組成，定期向其中添加額外的數據。 當然，我也可以使用螺柱。 擬合模型平台的殘差圖，但我想通過 Bonferroni Limits 自動執行異常值檢測，而不必使用該圖手動執行。 有沒有辦法提取螺柱的 Bonferroni 極限值。 殘差圖？ 一個小腳本在這裡也對我有用。 我已經和他討論過這個問題@喬納斯_林內他建議我在社區中解決這個問題。

 

謝謝它提前！

 

最佳約翰

 

 

 

基礎數據分析和建模數據混合和清理回复：使用 Bonferroni 限值自動識別離群值

非常感謝您的解決方案，對我來說，該腳本完美運行！

 

最佳約翰

回复：使用 Bonferroni 限值自動識別離群值

請使用此示例來說明如何訪問用於繪製 Bonferroni 極限線的值。

 

Names Default To Here( 1 );

// example
dt = Open( "$SAMPLE_DATA/Tiretread.jmp" );

// fit least squares with studentized residuals
obj = dt << Fit Model(
 Y( :ABRASION ),
 Effects( :SILICA, :SILANE, :SULFUR ),
 Personality( "Standard Least Squares" ),
 Emphasis( "Minimal Report" ),
 Run(
  :ABRASION << {Summary of Fit( 1 ), Analysis of Variance( 1 ),
  Parameter Estimates( 1 ), Scaled Estimates( 0 ),
  Plot Actual by Predicted( 0 ), Plot Regression( 0 ),
  Plot Residual by Predicted( 0 ), Plot Studentized Residuals( 1 ),
  Plot Effect Leverage( 1 ), Plot Residual by Normal Quantiles( 0 ),
  Box Cox Y Transformation( 0 )}
 )
);

// access report layer
rpt = obj << Report;

// access plot frame
frame = rpt["Studentized Residuals"][FrameBox(1)];

// access line segments in plot frame
line = frame << Find Segs;

// ask for value of line representing Bonferroni limits (only positive necessary)
bonferroni = (line[2] << Get Y Values)[1];

// display limit in Log
Show( bonferroni );