load_la

close all;
clear; % Limpar Variaveis
clc; % Limpar console
BASE = load("housing.data"); % Carregar base de dados

% Separar valores da base de dados
X = BASE(:,1:end-1);
y = BASE(:, end);

plot(X,y,"xb"); hold on; % Plot X*y

X = [ones(size(X),1) X];

W = pinv(X'*X)*X'*y; % Definir valores da regressão linear

X_novos = [3:0.01:9]';
X_novos = [ones(size(X_novos,1),1) X_novos];

y_h = X_novos*W;

plot(X_novos(:,2),y_h,"om"); % Plot reta de regressão

y_erro = X*W;
Erro = (y_erro - y).^2;
SME = (1/size(X,1)*sum(Erro));