## How to calculate slope and intercept

Community Trekker

Joined:

May 20, 2016

Hello,

I am facing 2 issues when trying tocalculate Slope and Intercept values:

1. I did not find an easy way to get the Slope and intercept except the use of "Bivariate + fit line"
2. Using the Bivariatein a loop seems to result in a slow execution of the code: I think that creating a table and the associated fit lile  several times is very slow

My question: is there a way to get the slope and intercept without using the bivariate to get faster ?

Here is the code I use in a loop to extract the slopes and intercepts:

``````biv1 = Bivariate( Y( MyTab:X ), X( MyTab:Y ), Fit Line( {Line Color( "Red" )} ) );
colBox = Report( biv1 )[Number Col Box( 7 )];
beta = colBox << GetAsMatrix;
My_Slope = beta[2];   // slope
My_Intercept = beta[1];   // intercept
Close( MyTab, nosave );
);``````

Thanks.

2 ACCEPTED SOLUTIONS

Accepted Solutions

Super User

Joined:

Jun 23, 2011

Solution

for starters, have a look at this option:

http://www.pega-analytics.co.uk/blog/linear-regression-matrix-form/

it is a solution offered by David Burnham in his blog. using matrix calculations he achieves high speeds for calculating the regression coefficients.

best,

ron

Staff

Joined:

Jun 23, 2011

Solution

Rather than having a loop, you could also consider making a big table with a 'By' variable, sticking with the Bivariate platform, then using 'Make Combined Data Table'. I expect the matrix approach will be faster, though.

5 REPLIES

Super User

Joined:

Jun 23, 2011

Solution

for starters, have a look at this option:

http://www.pega-analytics.co.uk/blog/linear-regression-matrix-form/

it is a solution offered by David Burnham in his blog. using matrix calculations he achieves high speeds for calculating the regression coefficients.

best,

ron

Community Trekker

Joined:

May 20, 2016

Thanks Ron,

That's what I was looking for: I will implement and test the spead :)

Staff

Joined:

Jun 23, 2011

Solution

Rather than having a loop, you could also consider making a big table with a 'By' variable, sticking with the Bivariate platform, then using 'Make Combined Data Table'. I expect the matrix approach will be faster, though.

Community Trekker

Joined:

Nov 16, 2016

Hi Ian,

I have followed your approach using the 'By' variable and sticking with the bivariate platform. I am not sure of the correct syntax for 'Make Combined Data Table' to work for multiple Fit Lines.  Can you help?

``````kXPlot = dt << Bivariate( Y( :LMS_dX ), X( :X_Nominal ), By(Column(::bcol)), Fit Line());
kXPlotDT = kXPlot < ???``````

Thank you,

Ry

Community Trekker

Joined:

Nov 16, 2016

I figured it out.

``kXPlotDT = Report( kXPlot )["Parameter Estimates"][Table Box(1)] << Make Combined Data Table;``

Regards,

Ry