520 likes | 1.05k Views
Pertemuan 5 Metnum 2011 Bilqis. PENCOCOKAN KURVA: ANALISA REGRESI. Pertemuan VI. Materi Minggu Ini. Pencocokan Kurva Regresi Kuadrat Terkecil (RKT) - RKT untuk Kurva Linier - RKT untuk Kurva Non-Linier Regresi Polynomial Tugas V. Tujuan. Pencocokan Kurva (1).
E N D
PENCOCOKAN KURVA:ANALISA REGRESI Pertemuan VI KomNum
Materi Minggu Ini • Pencocokan Kurva • Regresi Kuadrat Terkecil (RKT) - RKT untuk Kurva Linier - RKT untuk Kurva Non-Linier • Regresi Polynomial • Tugas V KomNum
Pencocokan Kurva(1) Seringkali data tersajikan dalam bentuk rangkaian nilai diskrit (deretan angka2 dalam urutan yang kontinu), tanpa disertai bentuk fungsi yang menghasilkan data tsb. Dalam kasus di atas, kita dapat men-”generate” fungsi sederhana untuk mengaproksimasi bentuk fungsi sebenarnya dengan memanfaatkan rangkaian data yang ada. Mengapa bentuk fungsi begitu penting bagi kita?... Selain untuk memenuhi kebutuhan proses numeris, seperti integrasi atau mendapatkan solusi pendekatan dari persamaan differensial, seringkali kita harus menganalisa tren atau melakukan pengujian hipotesa terhadap nilai2 diskrit yang dihasilkan oleh fungsi tsb. KomNum
Pencocokan Kurva(2) Terhadap keberadaan rangkaian data diskrit yang tidak diketahui fungsi asalnya, terdapat 2 hal yang dapat dilakukan : Pertama, berusaha mencari bentuk kurva (fungsi) yang mewakili rangkaian data diskrit tersebut. Kedua, berusaha meng-estimasi/memperkirakan nilai data pada titik2 yang belum diketahui. Keduanya dikenal sebagai teknik Pencocokan Kurva (curve fitting). KomNum
Pencocokan Kurva(3) Pendekatan2 yang lazim digunakan untuk melakukan Pencocokan Kurva antara lain adalah : Regresi Kuadrat Terkecil (least-square regresion) Metode ini digunakan apabila data yang tersaji memiliki tingkat kesalahan berarti (akurasi rendah). Anda hanya perlu membuat sebuah garis lurus yang merepresentasikan tren umum dari data2 tsb. KomNum
Pencocokan Kurva(3) Interpolasi Jika tingkat ketelitian data yang disajikan lebih baik, maka metode Interpolasi dapat dipakai. Anda dapat menggunakan segmen2 garis lurus untuk menghubungkan titik2 data (interpolasi linier). Atau, dengan menggunakan kurva (interpolasi polynomial) untuk menghubungkan titik2 data. Analisis regresi menggunakan sedikit notasi dan perhitungan statistik. Ini artinya, ada sedikit yang perlu anda ingat kembali… KomNum
Sedikit Notasi Statistik Rerata/rata-rata data ( y ) adalah jumlah nilai data (y)dibagi jumlah data (n) : yi y = n Deviasi Standar ( ) atau penyebaran nilai data adalah : ( yi – y )2 = n – 1 Atau dapat pula dinyatakan dalam bentuk Varians ( 2) : ( yi – y )2 2 = n – 1 √ KomNum
Regresi Kuadrat Terkecil(1) Jika data yang tersaji memiliki tingkat kesalahan yang cukup siginifikan, maka penggunaan interpolasi bukanlah pilihan yang bijaksana. Karena (kemungkinan besar) hasil pendekatannya akan kurang memuaskan. Ada cara yang lebih mudah, yaitu dengan membuat kurva yang dapat mewakili titik2 data tersebut. Katakan kurva ini adalah kurva dari fungsi g(x) yang ‘mirip’ dengan fungsi sebenarnya. Tetapi jika penyebaran titik datanya sangat besar? … Rasanya koq sulit cara di atas bisa berhasil dengan baik Kecuali jika, … kita dapat membuat kurva g(x) yang mampu meminimalkan perbedaan (selisih) antara titik2 data dengan kurva g(x). Teknik ini yang disebut dengan metode Regresi Kuadrat Terkecil. KomNum
Regresi Kuadrat Terkecil(2) Secara umum prosedur untuk mengaplikasikan metode Regresi Kuadrat Terkecil (RKT) ini adalah sbb : • Titik2 data diplotkedalamkoordinatcartesian. Dari poladatanyabisadilihat, apakahkurvapendekatan yang akankitabuatberbentukgarislurusataugarislengkung; • Tentukanlahsebuahfungsi g(x) yang dptmewakilifungsi titik2 data f(x) : • g(x) = a0 + a1x + a2x2 + … + arxr • g(x) fungsi buatan kita, fungsi aproximasi, bukan fungsi sebenarnya • Jika a0, a1, …, aradalah parameter fungsi g(x), makatentukannilai parameter2tsbsdmkhingga g(x) dptmendekati titik2 data; • Jikakoordinat titik2 data adalah M(xi, yi), makaselisihdenganfungsi g(x) adalah : • Ei = MiGi • = yi – g(xi; a0, a1, …, ar) • = yi – (a0 + a1xi + a2xi2 + a3xi3 + … + arxir) • Yi data sebenarnya KomNum
Regresi Kuadrat Terkecil(3) • Tingkat kesalahanfungsi g(x) diukurmelaluirumusberikut : • n n • D2 = Ei2 = { yi – g(xi)}2 • i=1 i=1 • Carinilai parameter a0, a1, …, arsdmkhingga D2dapatseminimummungkin. D2bernilai minimum jikaturunanpertamanyaterhadap a0, a1, …, ar = 0. • ∂D2 ∂D2 ∂D2 • = 0 = 0 … = 0 • ∂a0 ∂a1 ∂ar • 7. Persamaan no 6 diatasakanmemberikannilai parameter a0, a1, …, ar. Sehinggapersamaankurvaterbaik yang mewakili titik2 data dapatdiperoleh. Sebenarnya aproximasi KomNum
RKT untuk Kurva Linier(1) Permasalahan kita akan menjadi sederhana apabila kurva (yang diwakili oleh titik2 data) berbentuk garis lurus. Sepertiandaketahui, bahwapersamaanumumsebuahgarislurusadalah : g(x) = a + bx Sekarang kita mencari a dan b Dan jumlah-kuadrat-kesalahandapatdihitungmelaluipersamaan : nn D2 = ∑ Ei2 = ∑ { yi – a – bxi }2 i=1i=1 Agar nilai D2 dapat seminimal mungkin, maka persamaan jumlah- kuadrat-kesalahan harus diturunkan terhadap parameter a dan b (untuk kemudian disamadengankan 0). KomNum
RKT untuk Kurva Linier(2) Turunan pertama thd parameter b : ∂D2 = 0 ∂b ∂ n ( ∑ yi – a – b xi )2 = 0 ∂b i=1 n -2 ∑ [( yi – a – b xi ) xi ] = 0 i=1 ∑ yi xi - ∑ a xi - ∑ b xi2 = 0 … (2) Turunan pertama thd parameter a : ∂D2 = 0 ∂a ∂ n ( ∑ yi – a – b xi )2 = 0 ∂a i=1 n -2 ∑ ( yi – a – b xi ) = 0 i=1 ∑ yi - ∑ a - ∑ b xi = 0 … (1) KomNum
RKT untuk Kurva Linier(3) Jika ∑ a dapatdiasumsikansenilaidengan n a (sebagaiakibatpenjumlahanakumulatif suku-1 sampaisuku-n), makapersamaan(1)dapatditulissebagai : ∑ yi - ∑ a - ∑ b xi = 0 … (1) n a + ∑ b xi = ∑ yi n a = ∑ yi - ∑ b xi a = 1/n (∑ yi - ∑ b xi) … (3) a = y – b x… (4) KomNum
RKT untuk Kurva Linier(4) Sementara, persamaan(2)dapatditulis : ∑ yi xi - ∑ a xi - ∑ b xi2 = 0 … (2) ∑ a xi + ∑ b xi2 = ∑ yi xi… (5) Interpolasipersamaan(3)kepersamaan(5)akanmenghasilkan : ∑ xi 1/n ( ∑ yi - ∑ b xi )+ ∑ b xi2 = ∑ yi xi ∑ xi ∑ yi – ( ∑ xi )2 b + n ∑ b xi2 = n ∑ yi xi b [ n ∑ xi2 – ( ∑ xi )2 ] = n ∑ yi xi - ∑ yi ∑ xi atau, n ∑ yi xi - ∑ yi ∑ xi b = … (6) n ∑ xi2 – (∑ xi)2 KomNum
RKT untuk Kurva Linier(5) Persamaan (4) dan (6) dapat dimanfaatkan untuk mendapatkan koefisien a dan b, sehingga fungsi g(x) dapat diperoleh. Sedangkan untuk mengetahui derajat kesesuaian dari persamaan yang dicari, dapat dihitung melalui koefisien korelasi yang berbentuk : Dt2 – D2 r2 = … (7) Dt2 dengan, r = koefisien korelasi nn Dt2 = ∑ ( yi – y )2 D2 = ∑ ( yi – a - bxi )2 i=1i=1 Jika r = 1, maka fungsi g(x) memiliki tingkat kesesuaian yang tinggi dengan persamaan aslinya (f(x)). Atau r = 0 jika sebaliknya. KomNum
RKT untuk Kurva Linier(6) contoh : tentukan persamaan garis yang mewakili data berikut : x 1 2 3 4 5 6 7 8 9 10 y 4 6 8 10 14 16 20 22 24 28 KomNum
RKT untuk Kurva Linier(7) Nilaireratauntuk x dan y adalah : x = ∑ x / n = 55/10 = 5,5 y = ∑ y / n = 152/10 = 15,2 Jikapersamaanumumgarisdinyatakansebagai : y = a + bx, dan, n ∑ xiyi - ∑ xi ∑ yi 10 . 1058 – 55 . 152 2220 b = = = = 2,690909 n ∑ xi2 – (∑ xi)2 10 . 385 – (55)2 825 a = y – bx = 15,2 - 2,690909 . 5,5 = 0,4 Jadipersamaangaris yang mendekatirangkaian data tersebutadalah : y = 0,4 + 2,69.x KomNum
RKT untuk Kurva Non-Linier(1) Di dalam praktek akan sering kita jumpai kasus dimana plotting titik2 data memiliki tren berupa kurva lengkung. Sehingga persamaan yang sudah dikenalkan sebelumnya (RKT utk Kurva Linier) tidak dapat langsung digunakan. Dan lagi, kurva lengkung yang didekati dengan sebarang garis lurus akan menimbulkan kesalahan yang berarti. KomNum
RKT untuk Kurva Non-Linier(2) Kecuali untuk beberapa bentuk fungsi yang memang dapat didekati dengan metode Linierisasi Kurva Non-Linier. Fungsi2 tersebut antara lain : • Fungsi Eksponensial • y = a ebx dengan a1 dan b1 adalah konstanta • persamaan di atas dapat dilinierkan dengan logaritma-natural spt berikut : • ln y = ln a + b x ln e • jika ln e = 1, maka ln y = ln a + b x • persamaan di atas berbentuk garis lurus dengan kemiringan b, dan memotong sumbu ln y di ln a. KomNum
RKT untuk Kurva Non-Linier(3) 2. Fungsi Berpangkat fungsi berpangkat adalah contoh lain fungsi dengan kurvanya yang non-linier. y = a xb dengan a dan b adalah konstanta me-linier-kan fungsi di atas juga dapat dilakukan menggunakan persamaan logaritmik spt berikut : log y = b log x + log a persamaan di atas berbentuk garis lurus dengan kemiringan b dan memotong sumbu log y di log a. KomNum
RKT untuk Kurva Non-Linier(4) Contoh : Tentukan persamaan kurva lengkung yang diwakili serangkaian data berikut : x 1 2 3 4 5 y 0,5 1,7 3,4 5,7 8,4 Penyelesaian masalah di atas dilakukan melalui 2 fashion, transformasi log dan ln. KomNum
RKT untuk Kurva Non-Linier(5) Transformasi log (fungsi asli diamsusikan sebagai fungsi berpangkat) Misal persamaan yang dicari adalah : y = axb Persamaan tersebut dapat dilinierisasi melalui fungsi logaritmik sbb : log y = log axb atau dapat pula dinyatakan sebagai : log y = b log x + log a Atau jika p = log y; A = log a; B = b; q = log x; Maka persamaan log di atas dapat ditulis menjadi lebih sederhana : p = Bq + A KomNum
RKT untuk Kurva Non-Linier(6) Dari tabel tersebut dapat diperoleh beberapa parameter penting, seperti : y = ∑ yi / n = 19,7 / 5 = 3,94 q = ∑ log xi / n = 2,0791 / 5 = 0,4158 p = ∑ log yi / n = 2,1411 / 5 = 0,42822 Sedangkan koefisien A dan B dihitung melalui persamaan (4) dan (6) : n ∑ qi pi - ∑ qi ∑ pi 5 (1,4240) – 2,0791 (2,1411) 2,6684 B = = = = 1,7572 n ∑ qi2 – (∑ qi)2 5 (1,1692) – 2,0791 (2,0791) 1,5233 A = p – B q = 0,42822 – 1,7572 . 0,4158 = - 0,3024 Karena A = log a maka a = 0,4984 Karena B = b maka b = 1,7572 Dengan demikian fungsi yang dicari adalah : y = 0,4984 x1,7572 KomNum
RKT untuk Kurva Non-Linier(7) Transformasi ln (fungsi asli diasumsikan sebagai fungsi eksponensial) Misal persamaan yang dicari adalah : y = a ebx Persamaan tersebut dapat dilinierisasi melalui fungsi logaritmik sbb : ln y = ln axbx atau dapat pula dinyatakan sebagai : ln y = ln a + ln ebxatau ln y = ln a + bx Atau jika p = ln y; A = ln a; B = b; q = x; Maka persamaan log di atas dapat ditulis menjadi lebih sederhana : p = A + Bq KomNum
RKT untuk Kurva Non-Linier(8) Dari tabeltersebutdapatdiperolehbeberapa parameter penting, seperti : y = ∑ yi / n = 19,7 / 5 = 3,94 q = ∑ qi / n = 15 / 5 = 3 p = ∑ pi / n = 4,93/ 5 = 0,99 Sedangkankoefisien A dan B dihitungmelaluipersamaan (4) dan (6) : n ∑ qi pi - ∑ qi ∑ pi 5 (21,6425) – 15 (4,93) 34,2625 B = = = = 0,6852 n ∑ qi2 – (∑ qi)2 5 (55) – (15)2 50 A = p – B q = 0,99 – 0,68525 . 3,0 = - 1,06575 Karena A = ln a maka a = 0,34447 Karena B = b maka b = 0,68525 Dengandemikianfungsi yang dicariadalah : y = 0,34447 e0,68525x KomNum
RKT untuk Kurva Non-Linier(9) Sekarang waktunya memilih. Mana di antara 2 pendekatan yang memberikan akurasi lebih bagus. Caranya adalah dengan menghitung koefisien korelasi (7) : Dt2 – D2 r2 = Dt2 dengan n n n Dt2 = ∑ (yi – y)2 D2 = ∑ (yi – axb)2 D2 = ∑ (yi – aebx)2 i=1 i=1 i=1 KomNum
RKT untuk Kurva Non-Linier(9) Dari tabel tersebut dapat dicari nilai r untuk transformasi log : Dt2 – D2 ½ 40,132 – 0,00238 ½ r = = = 0,99997 Dt2 40,132 Sedangkan untuk transformasi ln : Dt2 – D2 ½ 40,132 – 5,60746 ½ r = = = 0,92751 Dt2 40,132 Dapat dilihat bahwa koefisien korelasi r untuktransformasi log lebih mendekati nilai 1 dibanding transformasi ln. Sehingga bisa disimpulkan bahwa persamaan yang didapat dari transformasi log memberikan pendekatan yang lebih baik. KomNum
Regresi Polynomial(1) Persamaan garis lurus dapat didekati dg Metode Kuadrat Terkecil. Sementara untuk kurva lengkung, pendekatan yang cukup logis adalah melalui transformasi data aslinya ke bentuk lain yang sesuai. Atau, untuk kurva lengkung dapat pula didekati dengan Regresi Polynomial. Persamaan polynomial orde-r mempunyai bentuk : y = a0 + a1x + a2x2 + … + arxr Jumlah kuadrat dari kesalahan adalah : n D2 = ∑ (yi – a0 – a1xi – a2xi2 - … - arxir)2 … (9) i=1 KomNum
Regresi Polynomial(2) Jika didiferensialkan terhadap setiap koefisiennya, maka persamaan (9) akan menjadi : ∂D2 n = - 2 ∑ (a0 - a1x - a2x2 - … - arxr) ∂a0 i=1 ∂D2 n = - 2 ∑ xi (a0 - a1x - a2x2 - … - arxr) ∂a1 i=1 ∂D2 n = - 2 ∑ xi2 (a0 - a1x - a2x2 - … - arxr) ∂a2 i=1 . . . . ∂D2 n = - 2 ∑ xir (a0 - a1x - a2x2 - … - arxr) ∂ar i=1 Bisa anda lihat, bentuk persamaan2 di samping menyerupai sistem persamaan KomNum
Regresi Polynomial(3) n ∑xi ∑xi2 … ∑xir a0 ∑yi ∑xi ∑xi2 ∑xi3 … ∑xir+1 a1 ∑xiyi ∑xi2 ∑xi3 ∑xi4 … ∑xir+2 a2 = ∑xi2yi . . . . . . . . . . . . . . . . . . ∑xir ∑xir+1 ∑xir+2 … ∑xir+r ar ∑xiryi Berarti kita bisa menyatakan himpunan persamaan turunan tersebut menjadi persamaan matriks AX = B. Dan selanjutnya kita dapat menggunakan metode eliminasi Gauss, Gauss-Jordan, dll untuk mencari nilai a0, a1, …, ar. KomNum
Regresi Polynomial(4) contoh : carilah persamaan kurva polynomial orde-2 yang mewakili data berikut : x 0 1 2 3 4 5 y 2,1 7,7 13,6 27,2 40,9 61,1 Persamaan polynomial orde-2 memiliki bentuk : g(x) = a0 + a1x + a2x2 Ei = yi – g(x) kesalahan untuk setiap nilai g(xi) thd yi Ei2 = ∑ (yi – a0 – a1x – a2x2)2 jumlah-kuadrat-kesalahan dr seluruh rangkaian nilai D2 = Ei2 Diferensialisasi D2 terhadap setiap koefisiennya akan menghasilkan bentuk : n ∑xi ∑xi2 a0 ∑yi ∑xi∑xi2 ∑xi3 a1 = ∑xiyi ∑xi2∑xi3 ∑xi4 a2 ∑xi2yi KomNum
Regresi Polynomial(5) Denganmenggunakan data tabel, makadiperoleh : 6a0 + 15a1 + 55a2 = 152,6 a2 = 1,86071 15a0 + 55a2 + 225a3= 585,6 a1 = 2,35929 55a0 + 225a2+ 979a3 = 2488,8 a0 = 2,47857 Jadipersamaankurva yang dicariadalah : y = 2,47857+2,35929 x + 1,86071 x2 KomNum
Latihan (1) • Tentukan: (a) rerata; (b) deviasistandar; dan (c) varian; dari data-data berikut • 0,95 1,42 1,54 1,55 1,63 • 1,32 1,15 1,47 1,95 1,25 • 1,46 1,47 1,92 1,35 1,05 • 1,85 1,74 1,65 1,78 1,71 • 2,39 1,82 2,06 2,14 2,27 • Gunakanregresikuadratterkeciluntukmenaksirfungsigarislurusdari data berikut : • x 1 3 5 7 10 12 13 16 18 20 • y 3 2 6 5 8 7 10 9 12 10 • 3. Gunakanregresikuadratterkeciluntukmenaksirfungsigarislurusdari data berikut : • x 4 6 8 10 14 16 20 22 24 28 28 34 36 38 • y 30 18 22 28 14 22 16 8 20 8 14 14 0 8 • 4. GunakanregresikuadratterkeciluntukmenaksirfungsikurvaLog daridata berikut : • x 1 2 2,5 4 6 8 8,5 • y 0,4 0,7 0,8 1,0 1,2 1,3 1,4 KomNum