Extracting the Bonferroni limits from the Studentized Residual Plot

[#1]

Hi,

I would like to get the numerical value of the Bonferroni limit from a existing Studentized Residuals plot and maybe see it change when changing the number of parameters/observations. I already tried to ask the Assisitant to write a JSL code, but that did not work. I also tried to implement the script from this post to my own data and that also did not work (I am a newbie at scripting) Solved: Automate Outlier identification with Bonferroni Limits - User Community

Thanks in advance!!

Best regards, Robin

[#2 · Accepted solution]

Hi @ConvolutionOx19 

I'll take the JSL from the post you've put and annotate it just to make it clearer what you need to do (essentially just copy paste your fit model details):

Names Default To Here( 1 );

// example
dt = Current Data Table();


///////////////////////***********CHANGE THIS BIT **************////////////////

// fit least squares with studentized residuals
obj = dt << 

//////COPY PASTE YOUR FIT MODEL TEXT HERE///////////////
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 )}
	)
);

//////////////////END OF COPY PASTE//////////////////////

// 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 );

To paste the right details, fit your model, go to the red triangle options on the top of the model page and get 'save script > to clipboard'.

Thanks,
Ben

[#3]

Hi @Ben_BarrIngh,

This worked. Thank you very much!

Regards, Robin