langkah 1 menginstal packages

library(plm)
library(lmtest)

langkah 2 membaca data

library(readxl)
Data_Panel_uppi <- read_excel("C:/Users/ASUS/OneDrive/Documents/Data Panel uppi.xlsx")
View(Data_Panel_uppi)

langkah 3 membuat deskripsi data yang di gunakan

str(Data_Panel_uppi)
summary(Data_Panel_uppi)
names(Data_Panel_uppi)

langkah 4 membuat estimasi modal regresi data paneldaei keiga dari ketiga modal yang ada

  1. common efect modal
common=plm(Miskin~TPT+IPM, data = Data_Panel_Miskin, model = "pooling")
summary(common)
  1. model fixed effect
FE=plm(Miskin~TPT+IPM,data=Data_Panel_uppi,model="within","effect"="time",index=c("Kabupaten","Tahun"))
summary(FE)
  1. model random effect
RE=plm(Miskin~TPT+IPM,data=Data_Panel_uppi,model="random",index=c("Kabupaten","Tahun"))
summary(RE)

langkah 5 menentukan model terbaik yang akan di gunakan

phtest(FE,RE)

Langkah 6: Melakukan Pengujian Breusch-Pangan

m=plm(Miskin~TPT+IPM, data = Data_Panel_uppi, model = "within", index = c("Kabupaten","Tahun"))
plmtest(m, effect = "twoways", type = "bp")#uji efek individu maupun waktu
plmtest(m,effect = "individual") # uji efek individu
plmtest(m,effect = "time",type = "bp")#uji efek waktu

Langkah 6: Pembuatan Model

model1=plm(Miskin~TPT+IPM, data = Data_Panel_uppi, model = "within", effect = "time", index = c("Kabupaten","Tahun"))
summary(model1)
fixef(model1,type="level")

Langkah 7: Pengujian Diagnostik

  1. Uji Korelasi Serial (Uji Multikolineritas)
pbgtest(FE)
  1. Uji Heteroskedasitas
bptest(FE)