Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula pada Data Keuangan
Alfi Khairiati 1 , Retno Budiarti 2* , I Gusti Putu Purnaba 3
1,2,3 Departemen Matematika, Fakultas MIPA, Institut Pertanian Bogor, Jalan Meranti, Kampus IPB Dramaga 16680, Jawa Barat, Indonesia
* Corresponding author. Email: retno.budiarti@gmail.com

## JAMBURA JOURNAL OF MATHEMATICS
Jambura J. Math. Vol. 4, No. 2, pp. 209–219, July 2022 Journal Homepage: http://ejurnal.ung.ac.id/index.php/jjom DOI: https://doi.org/10.34312/jjom.v4i2.13829

## ABSTRAK
Analisis regresi adalah analisis statistika untuk memprediksi atau menjelaskan hubungan antara peubah respon dengan satu atau beberapa peubah penjelas. Analisis regresi yang paling sederhana dan umum digunakan adalah regresi linear. Regresi copula sering dijadikan sebagai metode alternatif untuk mengatasi masalah dilanggarnya asumsi-asumsi dari model regresi linear. Keuntungan penggunaan regresi copula adalah peubah respon tidak harus mengikuti sebaran tertentu dan regresi copula dapat menjelaskan hubungan tak linear. Pada penelitian ini, copula yang digunakan adalah copula Gaussian dan copula student’s t . Tujuan utama dari penelitian ini adalah membandingkan hasil analisis regresi linear dengan regresi copula pada data keuangan. Pada metode regresi linear, tujuannya adalah menentukan nilai dugaan peubah respon dan menganalisis pengaruh faktor makroekonomi terhadap harga saham BMRI. Sementara itu, regresi copula digunakan untuk menduga parameter copula dengan metode Maximum Likelihood Estimation dan menentukan nilai dugaan peubah respon dengan pendekatan metode Monte Carlo. Ukuran keakuratan model yang digunakan adalah MAPE ( Mean Absolute Persentage Error ). Pada penelitian ini menggunakan data keuangan yang terdiri dari data harga saham BMRI sebagai peubah respon, serta IHSG dan nilai tukar rupiah sebagai peubah penjelas. Hasil penelitian menunjukkan bahwa nilai MAPE regresi linear dan regresi copula kecil dan tidak berbeda secara signifikan, artinya kedua regresi cukup baik dalam memodelkan data keuangan.

## Kata Kunci:
Regresi Linear; Copula Gaussian; Copula Student’s t; Regresi Copula

## ABSTRACT
Regression analysis is a statistical analysis to predict or explain the relationship between the response variable and one or more explanatory variables. The simplest and most commonly used regression analysis is linear regression. Copula regression is often used as an alternative method to overcome the problem of violating the assumptions of the linear regression model. The advantage of using copula regression is that the response variable does not have to follow a certain distribution and copula regression can explain nonlinear relationships. In this study, the copula used is the Gaussian copula and the student’s t copula. The main objective of this study is to compare the results of linear regression analysis with copula regression on financial data. In the linear regression method, the objective is to determine the estimated value of the response variable and to analyze the effect of macroeconomic factors on BMRI’s stock price. Meanwhile, copula regression was used to estimate the copula parameters using the Maximum Likelihood Estimation method and to determine the estimated value of the response variable using the Monte Carlo method. The measure of the accuracy of the model used is
e-ISSN: 2656-1344 c 2022 A. Khairiati, R. Budiarti, I.G.P. Purnaba | Under the license CC BY-NC 4.0 Received: 17 March 2022 | Accepted: 20 April 2022 | Online: 1 June 2022
JJoM | Jambura J. Math. 209 Volume 4 | Issue 2 | July 2022
MAPE (Mean Absolute Percentage Error). This study uses financial data consisting of BMRI stock price data as response variables, as well as IHSG and the rupiah exchange rate as explanatory variables. The results showed that the MAPE values for linear regression and copula regression were small and not significantly different, meaning that both regressions were quite good in modeling financial data.

## Keywords:
Linear Regression; Gaussian Copula; Student’s t Copula; Copula Regression

## Format Sitasi:
A. Khairiati, R. Budiarti, and I.G.P. Purnaba, ”Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula pada Data Keuangan”, Jambura J. Math. , vol. 4, No. 2, pp. 209–219, 2022, doi: https://doi.org/10.34312/jjom.v4i2.13829

## 1. Pendahuluan
Analisis regresi adalah analisis statistika untuk memprediksi atau menjelaskan hubungan antara peubah respon dengan satu atau beberapa peubah penjelas. Analisis regresi yang paling sederhana dan umum digunakan adalah regresi linear. Menurut Pindyck dan Rubinfeld [1] dan dan Nurrizqi [2] regresi linear tersebut memiliki beberapa asumsi, yaitu: (1) Peubah Y dan peubah X memiliki hubungan linear, (2) Variabel X 1 , X 2 , . . . , X n adalah variabel nonstokastik dan saling bebas, (3) Nilai harapan galat adalah nol, (4) Varian galat konstan untuk semua pengamatan, (5) Galat bersifat saling bebas satu dengan lainnya, dan (6) Galat berdistribusi normal. Namun begitu, asumsi-asumsi tersebut dalam penerapannya sering tidak terpenuhi dan oleh karena itu menurut Hakim [3] ini diperlukan model yang lebih fleksibel dalam memodelkan peubah respon yang dikenal dengan Generalized Linear Models (GLM). GLM dapat dikatakan memperumum regresi linear karena melonggarkan beberapa asumsi pada regresi linear. Namun ternyata ada yang lebih fleksibel dibandingkan dengan GLM, yaitu regresi copula. Ada beberapa keunggulan copula, antara lain kelonggaran asumsi terhadap sebaran khususnya sebaran normal, dapat menjelaskan hubungan yang tidak linear, mudah membangun sebaran bersamanya karena sebaran marginal dari peubah acak bisa berbeda atau bahkan sebaran marginalnya tidak diketahui [4, 5]. Oleh karena itu, regresi copula mempunyai cakupan analisis antar peubah menjadi lebih luas seperti, dimungkinkannya peubah penjelas merupakan peubah diskret maupun kontinu.
Penggunaan copula dalam bidang keuangan menawarkan metode yang lebih sederhana dan fleksibel untuk memodelkan struktur kebergantungan peubah banyak. Kelebihan copula dalam data keuangan adalah imbal hasil aset keuangan dibatasi oleh asumsi kenormalan. Menurut Darwis [4] ada beberapa penelitian menggunakan pendekatan copula diantaranya untuk memodelkan struktur hubungan pada klimatologi dan meteorologi. Hasil penelitian menunjukkan metode copula memiliki kinerja yang lebih baik dalam kondisi asumsi-asumsi kenormalan terlanggar. Dalam manajemen risiko, regresi copula telah ditunjukkan oleh Parsa and Klugman [6] untuk menyediakan pemodelan yang lebih baik daripada OLS dan GLM. Menurut Hikmah [7] penggunaan analisis copula pada kasus multivariat sangat sulit dilakukan untuk mendapatkan fungsi copula sehingga dikembangkan pendekatan lain dengan menggunakan vine copula. Sedangkan menurut Bramantya [8], pada ilmu studi konvensional aset-aset keuangan yang ada diasumsikan saling bebas padahal dalam kenyataannya di zaman pasar global seperti ini aset-aset keuangan tersebut tidak sepenuhnya bebas karena setiap bursa saham yang ada hampir dipastikan saling berinteraksi, baik secara langsung maupun tidak langsung. Pada studi kontemporer hubungan antar satu aset
JJoM | Jambura J. Math. 209 Volume 4 | Issue 2 | July 2022
Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula. . .
dengan yang lainnya dapat dimodelkan dengan sebuah fungsi yang bernama copula. Pada umumnya, kasus-kasus dibidang keuangan seringkali dijumpai data yang sebarannya tidak normal atau bahkan dipaksakan dengan asumsi sebaran normal. Salah satu kondisi dibidang keuangan yang tidak pasti yaitu saham. Saham adalah salah satu instrumen investasi yang nilai jualnya berdasarkan kinerja perusahaan penerbit saham. Pergerakan pasar yang tidak pasti memunculkan resiko yang berasal dari internal perusahaan penerbit saham seperti laporan keuangan atau kondisi keuangan perusahaan [4].
Data yang digunakan pada penelitian ini adalah data keuangan. Analisis yang dilakukan yaitu analisis hubungan antara harga saham BMRI (Bank Mandiri) sebagai peubah respon dengan faktor-faktor makroekonomi sebagai peubah penjelas menggunakan regresi linear dan regresi copula. Kemudian keakuratan kedua metode diukur menggunakan Mean Absolute Percentage Error (MAPE). Semakin kecil nilai MAPE, maka semakin akurat metode tersebut. Tujuan utama dari penelitian ini adalah membandingkan analisis regresi linear dengan regresi copula pada data keuangan. Kebaruan yang digunakan penulis pada karya ilmiah ini adalah mengkaji model regresi copula untuk kasus satu peubah respon dan dua peubah penjelas yang diterapkan pada bidang keuangan. Penelitian ini menggunakan bantuan software open-source R Studio.

## 2. Metode
Pada penelitian ini data keuangan yang digunakan adalah data harga saham BMRI sebagai peubah respon, IHSG (Imbal Harga Saham Gabungan) dan nilai tukar rupiah terhadap USD sebagai peubah penjelas. Periode data yang digunakan adalah periode bulanan dari Januari 2009 sampai Desember 2015 untuk setiap data.
Langkah-langkah analisis data yang dilakukan pada penelitian ini adalah :
1. Melakukan analisis regresi linear
(a) Plot data harga penutupan saham terhadap IHSG dan nilai tukar.
(b) Memodelkan hubungan antar variabel dengan model regresi linear.
(c) Estimasi parameter model regresi linear.
(d) Melakukan uji asumsi model regresi linear.
(e) Menentukan nilai dugaan perubah respon dan menganalisis pengaruh faktor makroekonomi terhadap harga saham BMRI.
(f) Menghitung Mean Absolute Percentage Error (MAPE) sebagai ukuran keakuratan model.
2. Melakukan analisis hubungan antar peubah dengan regresi copula (a) Mengidentifikasi sebaran untuk masing-masing peubah respon dan peubah penjelas.
(b) Transformasi peluang.
(c) Menggunakan copula Gaussian dan copula Student’s t untuk analisis regresi copula.
(d) Menduga parameter copula dengan metode Maximum Likelihood Estimation (MLE).
(e) Menentukan nilai dugaan peubah respon dengan pendekatan metode Monte Carlo [9].
(f) Menghitung Mean Absolute Percentage Error (MAPE) sebagai ukuran keakuratan model.
JJoM | Jambura J. Math. 210 Volume 4 | Issue 2 | July 2022

## 2.1. Analisis Regresi
Analisis regresi merupakan suatu analisis statistika yang digunakan untuk memprediksi atau menjelaskan hubungan suatu variabel dependen Y dengan menggunakan satu atau lebih variabel independen X 1 , . . . , X n . Jika n > 1, regresi yang terbentuk disebut regresi berganda dengan persamaan sebegai berikut :
Y t = δ 0 + δ 1 X 1 t + δ 2 X 2 t + · · · + δ n X nt + ε t ; t = 1, 2, . . . , N
dengan ε t merupakan galat pada waktu t . Model ditulis dalam bentuk matrik :
Y = X δ + ε .
Asumsi yang mendasari model regresi menurut Pindyck dan Rubinfeld [1] adalah:
1. Peubah Y dan peubah X memiliki hubungan linear
2. Variabel X 1 , . . . , X n adalah variabel nonstokastik dan saling bebas. 3. Nilai harapan galat adalah nol : E ( ε t ) = 0. 4. Varian galat konstan untuk semua pengamatan : E ε 2 t = σ 2 .
5. Galat bersifat saling bebas satu dengan lainnya.
6. Galat berdistribusi normal.

## 2.2. Ordinary Least Square (OLS)
Metode Ordinary Least Square (OLS) bertujuan untuk menemukan penduga parameter regresi dengan meminimumkan jumlah kuadrat galat. Akan dicari vektor parameter b δ yang meminimumkan fungsi berikut :
Q = N ∑ t = 1 ε 2 t = ε 0 ε
dengan,
ε = Y − b Y
dan
b Y = X b δ ε merupakan vektor galat berukutan N × 1 dengan b Y merupakan vektor penduga dari Y berukuran N × 1 [1].

## 2.3. Copula
Menurut Nelsen [10], fungsi C : [ 0, 1 ] 2 → [ 0, 1 ] adalah copula berdimensi dua, jika memenuhi sifat-sifat berikut :
1. C ( u 1 , 1 ) = u 1 dan C ( 1, u 2 ) = u 2 , 2. C grounded artinya C ( u 1 , 0 ) = C ( 0, u 2 ) = 0, 3. C two-increasing artinya ∀ u 1 , u 2 , v 1 , v 2 ∈ [ 0, 1 ] , u 1 < u 2 , v 1 < v 2 , maka berlaku C ( u 2 , v 2 ) − C ( u 1 , v 2 ) − C ( u 2 , v 1 ) + C ( u 1 , v 1 ) ≥ 0.
Berikut ini akan diperlihatkan suatu teorema yang sangat penting yaitu Teorema Sklar, yang merupakan teori tentang copula dan merupakan dasar dari banyak aplikasi teori statistika. Teorema Sklar ini menjelaskan tentang peranan copula terhadap fungsi sebaran multivariat dengan sebaran-sebaran marjinal univariat [10].
JJoM | Jambura J. Math. 211 Volume 4 | Issue 2 | July 2022
Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula. . .

## 2.3.1. Teorema Sklar
Misalkan F dan G adalah sebaran marjinal berturut-turut bagi peubah acak X dan Y maka fungsi H adalah fungsi sebaran bersama bagi ( X , Y ) . Ada copula C sehingga untuk setiap x , y ∈ R berlaku,
H ( x , y ) = C ( F ( x ) , G ( y )) = C ( u , v ) .

## 2.3.2. Copula Elliptical
Ada dua keluarga copula yang sering digunakan yaitu copula Elliptical dan copula Archimedean. Copula yang paling penting pada keluarga copula Elliptical adalah copula Gaussian dan copula Student’s t.

## 2.3.3. Copula Gaussian
Copula Gaussian 3 dimensi merupakan copula yang memenuhi persamaan berikut:
C Gauss R ( u , v , w ) = Φ R ( Φ − 1 ( u ) , Φ − 1 ( v ) , Φ − 1 ( w ))
dengan Φ adalah fungsi sebaran normal kumulatif dan Φ R adalah fungsi sebaran bersama ( joint distribution function ) sebaran normal dengan mean 0 dan matriks korelasi R [11].

## 2.3.4. Copula Student’s t
Copula Student’s t 3 dimensi merupakan copula yang memenuhi persamaan berikut :
C t V ( u , v , w ) = t V , R ( t − 1 V ( u ) , t − 1 V ( v ) , t − 1 V ( w ))
dimana, t V , R adalah sebaran t multivariat, t V adalah sebaran t dengan derajat kebebasan V dan R adalah matriks korelasi [11].

## 2.4. Regresi Copula
Pada penelitian ini dibatasi satu peubah respon dan dua peubah penjelas. Misalkan tiga peubah acak X 1 , X 2 , dan Y memiliki fungsi sebaran berturut-turut F 1 , F 2 , dan G serta H fungsi sebaran bersama antara X 1 , X 2 , dan Y . Misalkan peubah acak baru U = F 1 ( X 1 ) , V = F 2 ( X 2 ) dan W = G ( Y ) . Berdasarkan Teorema Sklar, fungsi sebarang tiga peubah X 1 , X 2 , dan Y didefinisikan dengan,
H ( x 1 , x 2 , y ) = P ( X 1 ≤ x 1 , X 2 ≤ x 2 , Y ≤ y ) = P ( F 1 ( X 1 ) ≤ F 1 ( x 1 ) , F 2 ( X 2 ) ≤ F 2 ( x 2 ) , G ( Y ) ≤ G ( y )) = P ( U ≤ u , V ≤ v , W ≤ w ) = C ( u , v , w )
Misalkan X 1 , X 2 , dan Y peubah acak kontinu dengan fungsi sebaran berturut-turut F 1 , F 2 , dan G , fungsi kepekatan marginal f 1 , f 2 , dan g , fungsi sebaran bersama H dan fungsi kepekatan peluang bersama h . Dari transformasi U = F 1 ( X 1 ) , V = F 2 ( X 2 ) dan W = G ( Y ) , didapat X 1 = F − 1 1 ( U ) , X 2 = F − 1 2 ( V ) dan Y = G − 1 ( W ) , maka dapat dinyatakan
JJoM | Jambura J. Math. 212 Volume 4 | Issue 2 | July 2022
c ( u , v , w ) = h F − 1 1 ( u ) , F − 1 2 ( v ) , G − 1 ( w ) | J | (1)
dengan
J = ∂ x 1 ∂ u ∂ x 1 ∂ v ∂ x 1 ∂ w ∂ x 2 ∂ u ∂ x 2 ∂ v ∂ x 2 ∂ w ∂ y ∂ u ∂ y ∂ v ∂ y ∂ w = 1 f 1 ( x 1 ) 0 0 0 1 f 2 ( x 2 ) 0 0 0 1 g ( y )
dan
c ( u , v , w ) = ∂ 3 C ( u , v , w ) ∂ u ∂ v ∂ w .
Sehingga persamaan (1) menjadi
c ( u , v , w ) = h F − 1 1 ( u ) , F − 1 2 ( v ) , G − 1 ( w ) 1 f 1 ( x 1 ) f 2 ( x 2 ) g ( y )
atau, persamaan fungsi h dapat ditulis
( x 1 , x 2 , y ) = C ( F 1 ( x 1 ) , F 2 ( x 2 ) , G ( y )) f 1 ( x 1 ) f 2 ( x 2 ) g ( y ) . (2)
Berikut ini, diperlihatkan kronologi penentuan regresi copula dengan memanfaatkan persamaan (2). Fungsi kepekatan peluang dari peubah acak Y bersyarat X 1 = x 1 , X 2 = x 2 adalah
m ( y | x 1 , x 2 ) = h ( x 1 , x 2, y ) f 1,2 ( x 1 , x 2 )
dengan asumsi X 1 dan X 2 peubah acak yang saling bebas, artinya f 1,2 ( x 1 , x 2 ) = f 1 ( x 1 ) f 2 ( x 2 ) sehingga didapatkan
m ( y | x 1 , x 2 ) = c ( F 1 ( x 1 ) , F 2 ( x 2 ) , G ( y )) f 1 ( x 1 ) f 2 ( x 2 ) g ( y ) f 1 ( x 1 ) f 2 ( x 2 ) (3) = c ( F 1 ( x 1 ) , F 2 ( x 2 ) , G ( y )) g ( y )
Berdasarkan persamaan (3), didapat nilai harapan peubah acak Y bersyarat X 1 = x 1 , X 2 = x 2 sebagai berikut:
E ( Y | X 1 = x 1 , X 2 = x 2 ) = Z ∞ − ∞ y m ( y | x 1 , x 2 ) dy (4) = Z ∞ − ∞ yc ( F 1 ( x 1 ) , F 2 ( x 2 ) , G ( y )) g ( y ) dy
Karena w = G ( y ) maka y = G − 1 ( w ) dan dw = g ( y ) dy , Persamaan (4) dapat ditulis
E ( Y | X 1 = x 1 , X 2 = x 2 ) = Z 1 0 G − 1 ( w ) c ( u , v , w ) dw (5)
JJoM | Jambura J. Math. 213 Volume 4 | Issue 2 | July 2022
Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula. . .
Persamaan (5) dikenal sebagai regresi copula, seperti berikut :
ˆ y = E ( Y | X 1 = x 1 , X 2 = x 2 ) = Z ∞ − ∞ yc ( F 1 ( x 1 ) , F 2 ( x 2 ) , G ( y )) g ( y ) dy

## 2.5. Maximum Likelihood Estimation
Maximum Likelihood Estimation (MLE) adalah metode untuk mengestimasi parameter model. Metode maximum likelihood merupakan metode pendugaan parameter suatu distribusi peluang dengan cara memaksimalkan fungsi likelihood . Misalkan Z merupakan peubah acak dengan fungsi kepekatan peluang f ( z | θ ) dengan θ merupakan parameter yang belum diketahui. Jika z 1 , . . . , z N adalah nilai sampel acak, maka fungsi likelihood dapat didefinisikan sebagai:
L ( θ ) = f ( z 1 , . . . , z N | θ ) · · · = N ∏ t = 1 f ( z t | θ )
MLE dari θ adalah nilai θ yang memaksimumkan L ( θ ) . MLE dari θ adalah nilai penyelesaian dari persamaan,
∂ ∂θ L ( θ ) = 0 [ 12 ] .

## 2.6. Akaike Information Criterian (AIC)
Metode Akaike Information Criterian (AIC) yang dikenalkan oleh Hirotugu Akaike (1974) merupakan salah satu ukuran mengevaluasi model. Pemilihan model terbaik dengan nilai AIC terkecil, dengan persamaan :
AIC = − 2log L + 2 k
dengan,
L adalah logaritma fungsi kemungkinan; k adalah banyaknya parameter [13].

## 2.7. Bayesian Information Criterian (BIC)
Metode Bayesian Infromation Criterian (BIC) yang dikenalkan oleh Gideon E.Schwarz merupakan salah satu ukuran untuk mengevaluasi model. Pemilihan model terbaik dengan nilai BIC terkecil, dengan persamaan :
BIC = − 2log L + k log n
dengan,
L adalah logaritma fungsi kemungkinan;
k adalah banyaknya parameter;
n adalah banyaknya data [13].

## 2.8. Mean Absolute Persentage Error
Menurut Rinaldi [14], menentukan kesalahan prediksi dapat diukur menggunakan berbagai ukuran error , dimana diketahui bahwa ˆ y t adalah data prediksi dari
JJoM | Jambura J. Math. 214 Volume 4 | Issue 2 | July 2022
perhitungan model pada waktu t . Dalam penelitian ini ukuran error yang digunakan dalam menghitung nilai ketepatan prediksi yaitu MAPE ( Mean Absolute Persentage Error ). MAPE merupakan ukuran ketepatan relatif yang digunakan untuk mengetahui persentase penyimpangan hasil peramalan [15], dengan persamaan sebagai berikut:
MAPE = 1 n n ∑ i = 1 | PE i |
Galat persentase ( percentage error ) :
PE i = y t − ˆ y t y t × 100%.

## 3. Hasil dan Pembahasan

## 3.1. Analisis Regresi Linear
Model regresi linear yang akan digunakan untuk data keuangan harga saham BMRI sebagai berikut:
y t = δ 0 + δ 1 I HSG t + δ 2 NT t + ε t
dengan:
y t : Harga saham BMRI pada waktu t δ 0 , δ 1 , δ 2 : Koefisien model I HSG t : Harga saham IHSG pada waktu t NT t : Nilai tukar pada waktu t ε t : galat pada waktu t .
Langkah berikutnya, menentukan pendugaan parameter model menggunakan metode Ordinary Least Square . Hasil pendugaan disajikan pada Tabel 1.
Tabel 1. Hasil metode Ordinary Least Square
Variabel Intercept IHSG NT Parameter -1,126 1,101 0,053 p-value 0,000 0,001
Berdasarkan hasil estimasi model pada Tabel 1, persamaan regresi linear dinyatakan dengan,
ˆ y t = − 1, 126 + 1, 101 I HSG + 0, 053 NT .
Nilai p-value dari perubahan IHSG dan nilai tukar < 5% yang berarti perubahan IHSG dan nilai tukar dapat dianggap berpengaruh terhadap harga saham BMRI. Dari persamaan di atas menunjukkan bahwa koefisien IHSG sebesar 1,101. Hal ini berarti bahwa setiap kenaikan pada IHSG sebesar satu-satuan, maka akan menyebabkan harga saham BMRI naik sebesar 1,101 satuan dengan asumsi variabel lain bernilai tetap. Selain itu, koefisien nilai tukar adalah sebesar 0,053. Hal ini berarti bahwa setiap kenaikan nilai tukar sebesar satu satuan, maka akan menyebabkan harga penutupan saham BMRI naik sebesar 0,053 satuan dengan asumsi variabel lain bernilai tetap. Selanjutnya asumsi model regresi linear diuji dan hasilnya yang ditunjukkan pada Tabel 2.
Tabel 2 menunjukkan bahwa semua asumsi model regresi linear terpenuhi.
JJoM | Jambura J. Math. 215 Volume 4 | Issue 2 | July 2022
Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula. . .
Tabel 2. Hasil uji asumsi regresi linear
Uji Metode Hasil Kesimpulan Kenormalan Jarque Bera Nilai p = 0,835 Normal Homoskedastisitas Breush Pagan Nilai p = 0,189 Homoskedastisitas Autokorelasi Durbin Watson Nilai p = 0,8513 Tidak ada Autokorelasi Multikoliniearitas VIF VIF X 1 =1,198404 Tidak terjadi VIF X 2 =1,198404 Multikolinieritas

## 3.2. Analisis Regresi Copula
Regresi copula yang digunakan dalam penelitian ini adalah persamaan (5). Model regresi copula yang terbentuk bergantung pada fungsi sebaran copula. Nilai y dugaan didapatkan dengan mencari nilai harapan dari fungsi peluang bersyarat seperti persamaan berikut,
ˆ y = E ( Y | X 1 , X 2 ) + ( y − E ( Y | X 1 , X 2 )) = E ( Y | X 1 , X 2 ) . (6)
Dengan memanfaatkan persamaan (6) perhitungan ˆ y =
E ( Y | X 1 , X 2 ) dengan pendekatan metode Monte Carlo, yaitu:
ˆ y = Z 1 0 G − 1 ( w ) c ( u , v , w ) dw ≈ 1 n n ∑ t = 1 G − 1 ( w t ) c ( u t , v t , w t ) dw .
Langkah pertama dalam analisis regresi copula adalah melakukan identifikasi sebaran marjinal dengan mencocokkan grafik sebaran empirik dengan beberapa sebaran khususnya [16], sebagaimana ditampilkan pada Gambar 1.
Gambar 1. Grafik sebaran empiris dengan sebaran marginal peubah respon
Pada Gambar 1 sebaran empiris peubah respon diduga memiliki sebaran Weibull, normal, dan lognormal. Oleh karena itu, untuk memutuskan sebaran peubah respon yang paling sesuai ditentukan menggunakan nilai AIC dan BIC. Hal yang sama dilakukan untuk peubah penjelas dan hasilnya dirangkum pada Tabel 3.
Pada Tabel 3 menunjukkan hasil AIC dan BIC terkecil dari peubah harga saham BMRI
JJoM | Jambura J. Math. 216 Volume 4 | Issue 2 | July 2022
Tabel 3. Sebaran yang sesuai untuk peubah respon dan peubah penjelas
Peubah Sebaran yang dipilih Nilai AIC Nilai BIC BMRI Weibull* 271.834 276.696 Normal 272.678 277.539 Lognormal 300.623 305.485 IHSG Weibull* 1403.396 1408.258 Normal 1409.036 1413.897 Lognormal 1436.249 1441.110 Nilai Tukar Weibull 1494.534 1499.396 Normal 1483.685 1488.546 Lognormal* 1475.321 1480.183
Keterangan: *) menunjukkan nilai AIC dan BIC terkecil
dan IHSG yaitu sebaran Weibull, sedangkan peubah nilai tukar yaitu sebaran lognormal. Hasil penduga parameter sebaran yang terpilih pada Tabel 3 ditunjukkan pada Tabel 4 dengan memilih hasil nilai AIC dan BIC terkecil dari tiga peubah yang mendekati sebaran empirisnya.
Tabel 4. Penduga parameter sebaran yang terpilih
Peubah Sebaran Parameter BMRI Weibull b µ = 3,570478 b σ = 4,146079 IHSG Weibull b µ = 4,666251 b σ = 4290,993 Nilai Tukar Lognormal
b µ = 9,246011
b σ = 0,148512

## 3.3. Pendugaan Parameter Copula
Parameter copula Gaussian dan Student’s t berturut-turut adalah matriks korelasi R 1 dan R 2 yang diestimasi menggunakan metode Maximum Likelihood Estimation [17, 18]. Adapun nilai estimasi dari matriks R 1 dan R 2 adalah sebagai berikut:
R 1 =   1 0.9888 0.4755 0.9888 1 0.4828 0.4755 0.4828 1   R 2 =   1 0.9890 0.4752 0.9890 1 0.4829 0.4752 0.4829 1 
 .

## 3.4. Evaluasi Model Linear dan Copula
Setelah dilakukan pendugaan peubah respon kedua model regresi linear dan regresi copula, selanjutnya dilakukan pengukuran keakuratan model menggunakapan MAPE. Evaluasi model tersebut diukur ketepatan nilai dugaan menggunakan Mean Absolute Percentage Error (MAPE), dengan nilai MAPE disajikan pada Tabel 5.

## Tabel 5. Nilai MAPE model linear dan copula
Kode Saham MAPE model linear MAPE model copula BMRI 0,046 0,077
Untuk memperjelaskan hasil pengukuran keakuratan kedua model, dilakukan plot antara harga saham BMRI aktual dengan harga saham BMRI dugaan berdasarkan regresi linear dan regresi copula pada Gambar 2.
Gambar 2 mendukung nilai MAPE yang kecil pada kedua model yaitu 4,6% untuk model
JJoM | Jambura J. Math. 217 Volume 4 | Issue 2 | July 2022
Perbandingan Analisis Regresi Linear dengan Analisis Regresi Copula. . .
Gambar 2. Plot Y aktual dengan ˆ Y regresi linear dan ˆ Y regresi copula
linear dan 7,7% untuk model copula, yaitu harga saham BMRI hasil dugaan kedua model terlihat mendektai harga saham BMRI aktual.

## 4. Kesimpulan
Berdasarkan hasil analisis regresi linear dan analisis regresi copula, penelitian ini telah mengaplikasikan harga saham BMRI dengan faktor makroekonominya pada data yang digunakan adalah periode bulanan dari Januari 2009 – Desember 2015, dapat disimpulkan bahwa: (1) Pada analisis regresi linear dengan membandingkan Y aktual dan Y dugaan diperoleh nilai MAPE sebesar 4,6%; (2) Pada analisis regresi copula dengan membandingkan Y aktual dan Y dugaan diperoleh nilai MAPE sebesar 7,7%; (3) Dengan membandingkan nilai MAPE kedua model, dikatakan bahwa kedua regresi cukup baik dalam memodelkan data keuangan ini. Hal tersebut ditunjukkan dengan nilai MAPE kedua model kurang dari 10%. Pada penelitian ini, hanya dibatasi penggunaan peubah penjelasnya saling bebas, sehingga disarankan untuk penelitian selanjutnya mengkaji dan mengaplikasikan pemodelan dengan menggunakan peubah penjelas yang tidak saling bebas. Kemudian disarankan dalam penggunaan copula selanjutnya menggunakan keluarga copula yang lebih kompleks.

## Referensi
[1] R. S. Pindyck and D. L. Rubinfeld, Econometrics Models and Economic Forecast . Singapore: MC Graw Hill, 1998.
[2] S. T. Nurrizqi, Analisis Pengaruh Faktor Makroekonomi terhadap Imbal Hasil Saham Sektor Keuangan . Bogor: IPB University, 2020.
[3] L. Hakim, Studi Simulasi Pembandingan Regresi Generalized Linear Model dan Regresi Berbasis Copula . Bogor: IPB University, 2021.
[4] D. Darwis, Analisis Hubungan dan Prediksi Indeks Harga Saham Gabungan dengan Faktor Makroekonomi melalui Pendekatan Copula . Bogor: IPB University, 2016.
[5] G. Masarotto and C. Varin, “Gaussian copula marginal regression,” Electronic Journal of Statistics , vol. 6, no. none, jan 2012, doi: http://dx.doi.org/10.1214/12-EJS721.
[6] A. R. Parsa and S. A. Klugman, “Copula Regression,” Casualty Actuarial Society , vol. 5, no. 1, pp. 46–50, 2011.
[7] I. R. Hikmah, Identifikasi Struktur Dependensi dan Prediksi Indeks Harga Saham Gabungan dengan Pendekatan C-V Vine Copula . Bogor: IPB University, 2017.
[8] I. S. Bramantya, R. Budiarti, and I. G. P. Purnaba, “Pemodelan Indeks Harga Saham Gabungan dan Penentuan Rank Correlation dengan Menggunakan Copula,” in Prosiding
JJoM | Jambura J. Math. 218 Volume 4 | Issue 2 | July 2022
Seminar Nasional Matematika, Statistika, Pendidikan Matematika, dan Komputasi , 2014, pp. 502– 512.
[9] N. Metropolis and S. Ulam, “The Monte Carlo Method,” Journal of the American Statistical Association , vol. 44, no. 247, pp. 335–341, sep 1949, doi: http://dx.doi.org/10.1080/01621459. 1949.10483310.
[10] R. B. Nelsen, An Introduction to Copulas , ser. Springer Series in Statistics. New York, NY: Springer New York, 2006, doi: http://dx.doi.org/10.1007/0-387-28678-0.
[11] J. Rank, “Copulas: From theory to application in Finance,” 2006.
[12] H. Joe, Multivariate Models and Multivariate Dependence Concepts . London: Crc Press, 1997.
[13] K. P. Burnham and D. R. Anderson, “Multimodel Inference,” Sociological Methods & Research , vol. 33, no. 2, pp. 261–304, nov 2004, doi: http://dx.doi.org/10.1177/0049124104268644.
[14] G. A. Rinadi, L. R. Sasongko, and B. Susanto, “Regresi Median Pada Copula Bivariat,” JTAM — Jurnal Teori dan Aplikasi Matematika , vol. 3, no. 1, pp. 07–14, apr 2019, doi: http://dx.doi.org/10.31764/jtam.v3i1.728.
[15] W. Wei, Time Series Analysis: Univariate and Multivariate Methods . New York: Addison- Wesley Publishing Company, Inc., 1990.
[16] S. A. Klugman, H. H. Panjer, and G. E. Willmot, Loss Models: Further Topics . Canada: John Willey & Sons, Inc, 2013.
[17] L. Zhang and V. P. Singh, “Non-Archimedean Copulas,” in Copulas and their Applications in Water Resources Engineering . Cambridge University Press, jan 2019, pp. 261–303, doi: http://dx.doi.org/10.1017/9781108565103.008.
[18] Q. Zhang, V. P. Singh, J. Li, F. Jiang, and Y. Bai, “Spatio-temporal variations of precipitation extremes in Xinjiang, China,” Journal of Hydrology , vol. 434-435, no. 7, pp. 7–18, apr 2012, doi: http://dx.doi.org/10.1016/j.jhydrol.2012.02.038.
This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution-NonCommercial 4.0 International License . Editorial of JJoM: Department of Mathematics, Universitas Negeri Gorontalo, Jln. Prof. Dr. Ing. B.J. Habibie, Moutong, Tilongkabila, Kabupaten Bone Bolango, Provinsi Gorontalo 96119, Indonesia.
219 Volume 4 | Issue 2 | July 2022