Fit a constant only linear regression model using 'fitlm'

7 vues (au cours des 30 derniers jours)
M
M le 1 Mai 2015
Commenté : M le 1 Mai 2015
I have the following simple regression model
y(t) = B + u(t)
where B=1 and u(t) are random drawings from the standard normal distribution. Also n = 100. I would like to fit a constant only linear regression model but am unsure how to do so. I imagine I have to use 'fitlm', but for some reason cannot specify that there are no predictor variables . My code so far is simply:
b = 1
u = randn(100,1)
y = b + u

Réponse acceptée

Brendan Hamm
Brendan Hamm le 1 Mai 2015
If you wanted to fit this using fitlm you could do the following:
fit1 = fitlm(ones(size(y)),y,'y~1');

Plus de réponses (1)

John D'Errico
John D'Errico le 1 Mai 2015
Well, you CAN use a tool like fitlm to do this. But that would be like using a Mack truck to take a pea to Boston.
The linear regression estimator of a constant model is just the mean. So...
B = mean(u);
Of course it fails to give you statistics on the model. But it is the mean for god sakes! It won't fit very well, unless your data is constant. So plot your data, and look at the residuals. If there is a pattern in the data, then a constant model might be inadequate for your purposes.
  2 commentaires
M
M le 1 Mai 2015
Thanks for the reply. I understand what you're saying, however I was required to show that calculating the OLS estimate of beta would yield the same result as mean(u) so this is why I asked for help with fitlm. Cheers
M
M le 1 Mai 2015
mean(y), rather!

Connectez-vous pour commenter.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by