Linear Regression in Matlab Programming

x = [1 2 3 4 5 6];
y = [2.2 2.8 3.6 4.5 5.1 5.9];

% Fit 1st-degree (linear) polynomial coeff = polyfit(x, y, 1);
slope = coeff(1);
intercept = coeff(2);
fprintf('Slope (m): %.4f, Intercept (c): %.4f ', slope, intercept);

% Construct design matrix with intercept column X = [x(:) ones(length(x), 1)];
Y = y(:);
% Solve normal equations: theta = (X^T * X)^(-1) * X^T * Y theta = (X' * X)  (X' * Y);
slope = theta(1);
intercept = theta(2);
fprintf('Slope (m): %.4f, Intercept (c): %.4f ', slope, intercept);

% Fit LinearModel object mdl = fitlm(x, y);
disp(mdl);

y_fit = slope * x + intercept;
figure;
scatter(x, y, 'filled', 'MarkerFaceColor', 'b');
hold on;
plot(x, y_fit, 'r-', 'LineWidth', 2);
title('Linear Regression Model Fit');
xlabel('Independent Variable (x)');
ylabel('Dependent Variable (y)');
legend('Observed Data', 'Fitted Line', 'Location', 'northwest');
grid on;
hold off;

SS_res = sum((y - y_fit).^2);
SS_tot = sum((y - mean(y)).^2);
R_squared = 1 - (SS_res / SS_tot);
fprintf('R-squared (R^2): %.4f ', R_squared);