regression line in quadratic form
14 ビュー (過去 30 日間)
古いコメントを表示
Hi, all
During my exercise, I'd like to draw a regression line fitted with the following quadratic equation.
The followings are the code I wrote down.
y = log(wage(index));
x = [ones(length(y), 1) experience(index) experience(index).^2];
b = (x'*x)^(-1)*x'*y;
yhat = x*b;
figure;
plot(x(:, 2), y, '.', x(:, 2), yhat);
title('exercise');
xlabel('experience');
ylabel('log(wage)');
print('-dpdf', 'excercise');
close;
And this is the result.
I firstly guess the single curve will be shown on the figure, but... why is there a bunch of lines, and how can I fix this with minimization of fixing?
Thank you in advance, guys.
0 件のコメント
回答 (1 件)
J. Alex Lee
2020 年 9 月 9 日
編集済み: J. Alex Lee
2020 年 9 月 12 日
[edited] Alternatively, just create a new variable to hold an already sorted "model" experience
N = 200; % choose an arbitrary number of points to define model
ExperienceMdl = linspace(0,70,N)';
xMdl = [ones(N,1),ExperienceMdl,ExperienceMdl.^2]
yhat = ExperienceMdl * b;
And also, your compuation b looks inefficient, you can simply do
b = x\y
9 件のコメント
J. Alex Lee
2020 年 9 月 13 日
You saw multiple regression lines because you plotted the results of your model based on the data, which is not ordered in "experience" variable.
One of the previous answeres suggested you to sort your data in order of your "experience" variable. But you can realize that there is no need to plot your model based on the data points alone.
linspace(a,b,N) will generate a vector of N-values evenly spaced (and sorted) from a to b. I chose N=200 points because that should be more than enough to resolve the shape of your model over the domain 0~70, which is roughly the extent of your data.
参考
カテゴリ
Help Center および File Exchange で Linear Regression についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!