Subscribe Bookmark RSS Feed

Confidence and prediction intervals for orthogonal regression?

JesperJohansen

Community Trekker

Joined:

Apr 15, 2013

I can tell that this question was posted back in 2012, but never answered, so here it goes again:

Is there a way to plot confidence and prediction intervals for orthogonal regression in the Bivariate platform. Also I would like to save the interval equations. Both options exist for "Fit line", i.e. OLS linear regression.

If the option does not exist, I would like to request it be implementet. 

BR
Jesper
4 REPLIES
JesperJohansen

Community Trekker

Joined:

Apr 15, 2013

Noone from JMP reading this?

BR
Jesper
ian_jmp

Staff

Joined:

Jun 23, 2011

With JMP Pro, you could try bootstrapping:

NamesDefaultToHere(1);

// Sample data
dt = NewTable("Orthogonal Regression",
		NewColumn("x1", Formula(RandomUniform(-1, 1))),
		NewColumn("x2", Formula(RandomUniform(-1, 1))),
		AddRows(50)
);

// Fit orthogonal
fo = dt << Bivariate(Y( :x2 ), X( :x1 ), Fit Orthogonal( Equal Variances ));

// Bootstrap fitted parameters
dt2 = Report(fo)[TableBox(2)] << Bootstrap(1000);

// Look at the result
dt2 << Distribution(
				Continuous Distribution( Column( :Intercept ) ),
				Continuous Distribution( Column( :Slope ) ),
				SendToReport(
					Dispatch( {"Intercept"}, "Quantiles", OutlineBox, {Close( 1 )} ),
					Dispatch( {"Intercept"}, "Summary Statistics", OutlineBox, {Close( 1 )} ),
					Dispatch( {"Slope"}, "Quantiles", OutlineBox, {Close( 1 )} ),
					Dispatch( {"Slope"}, "Summary Statistics", OutlineBox, {Close( 1 )} )
				)
			);
JesperJohansen

Community Trekker

Joined:

Apr 15, 2013

Alas, I don't have Pro :'-(
BR
Jesper
akrantasiwakoti

New Contributor

Joined:

Aug 10, 2017

Hi Ian,

 

Thank you for this JSL script. However, the script only boostraps estimates of slope and intercept. I was curious how would I then go about estimating prediction interval in deming regression using this output? I would appreciate any assistance.

 

Thank you