[Octave] Cara Menyelesaikan Multi Regresi Linear


Saya mengambil data dari ilmu statistik yang saya simpan di x.dat
1 30
1 35
2 42
2 47
3 50
3 62
4 64
5 71
6 79
7 85

Sedangkan data y, disimpan di y.dat
40
44
49
53
60
65
69
78
85
92

Saya membuat script kecil di octave untuk menyelesaikan kasus diatas. Script ini menggunakan persamaan matriks transpose biasa
\Theta = (X^T X)^{-1}X^T y
dimana theta adalah koeficien yang kita cari.

Berikut script yang saya buat

x = load('x.dat');
y = load('y.dat');

figure;
plot(x,y,'o');
xlabel('x1 and x2');
ylabel('y');
legend('x1','x2');
m = length(y);  
x = [ones(m, 1), x] ;
xt = transpose(x);
disp(xt);

thetabaru = pinv(xt*x)*xt*y;
disp(thetabaru);
hold on;
plot(x(:, 2), x*thetabaru, '-',x(:,3),x*thetabaru,'-');

Maka akan menghasilkan theta : 20.63777 3.74223 0.53343
dan gambar grafik :

  1. assalamua’laikum…
    maaf ukhti mau nanya itu slide number yang di bawah untuk penunjuk ke tulisan terdahulu bikin manual atau gimana..???
    apa sudah setingan themanya..???
    terima kasih…

    • yogi
    • November 29th, 2013

    nice post , please visit back😀 http://yosmantri.student.ipb.ac.id/ thanks😀

  1. No trackbacks yet.

Tinggalkan Balasan

Isikan data di bawah atau klik salah satu ikon untuk log in:

Logo WordPress.com

You are commenting using your WordPress.com account. Logout / Ubah )

Gambar Twitter

You are commenting using your Twitter account. Logout / Ubah )

Foto Facebook

You are commenting using your Facebook account. Logout / Ubah )

Foto Google+

You are commenting using your Google+ account. Logout / Ubah )

Connecting to %s