How do I get the constant term when performing multiple linear regression using STEPWISEFIT?

I'm using STEPWISEFIT to perform multiple linear regression of the response values in the n-by-1 vector y on the predictive terms in the n-by-p matrix X as follows:

[b,SE,PVAL,INMODEL,STATS,NEXTSTEP,HISTORY] = stepwisefit(X,y);

However, the size of b contains only p-1 coefficients i.e., the constant term is not included.

Answer 1

[MathWorks Support Team]

The constant term of the model corresponds to the "intercept" field of the output structure STATS and can be obtained as follows:

b0 = STATS.intercept