Model Regresi Logit Multinomial
Nofrida Elly
6/10/2021
Menggunakan data GSS 2006 untuk mengetahui ketertarikan responden akan politik (polint) menurut jenis kelamin (sex) dan ras/warna kulit (race) responden serta interaksinya.
Tujuan : 1. Mengetahui ketertarikan responden akan politik (polint) menurut jenis kelamin (sex) dan ras/warna kulit (race) 2. Menghitung peluang laki-laki ras kulit hitam sangat tertarik pada politik
Keterangan: Reference level polint = Not at All Interseted Reference level sex = female Reference level race =others
setwd("D:/PASCASARJANA/Semester 3/1. ADK/02. Praktikum/13")Input data
sex <- factor(rep(c("Male", "Female"), 15))
polint <- factor(rep(c("VI", "FI", "SI", "NVI", "NAI"), rep(6,5)))
race <- rep(factor(rep(c("White", "Black", "Other"), rep(2,3))),5)
counts <- c(134,86,22,15,12,9,
139,138,18,32,17,17,
123,182,21,46,26,29,
64,97,10,21,25,17,
64,78,14,18,18,24)
#### Referensi
polint <- relevel(polint, ref = "NAI")
sex <- relevel(sex, ref = "Female")
race <- relevel(race, ref = "Other")
#### Membuat data frame
dataku <- data.frame(race, sex, polint, counts)
datapolint <- dataku[rep(row.names(dataku), counts), 1:3]
#### Cek Struktur dan Ukuran data
head(datapolint)## race sex polint
## 1 White Male VI
## 1.1 White Male VI
## 1.2 White Male VI
## 1.3 White Male VI
## 1.4 White Male VI
## 1.5 White Male VIdim(datapolint)## [1] 1516 3#install.packages("foreign")
#install.packages("nnet")
library("foreign")
library("nnet")## Warning: package 'nnet' was built under R version 3.6.3model1 <- multinom (polint ~ race + sex, data=datapolint)## # weights: 25 (16 variable)
## initial value 2439.907875
## iter 10 value 2361.773697
## iter 20 value 2359.061081
## final value 2359.059069
## convergedsummary(model1)## Call:
## multinom(formula = polint ~ race + sex, data = datapolint)
##
## Coefficients:
## (Intercept) raceBlack raceWhite sexMale
## FI -0.29944790 0.6794937 0.8883678 0.16985529
## NVI 0.04339509 -0.0431156 0.1209414 -0.08916677
## SI 0.35266179 0.4472863 0.4858761 -0.17452473
## VI -1.08119659 0.9245297 1.1649180 0.66953332
##
## Std. Errors:
## (Intercept) raceBlack raceWhite sexMale
## FI 0.2476602 0.3241398 0.2530080 0.1739927
## NVI 0.2371914 0.3342628 0.2469414 0.1909547
## SI 0.2206358 0.2977949 0.2289730 0.1696374
## VI 0.2894266 0.3625649 0.2895054 0.1853683
##
## Residual Deviance: 4718.118
## AIC: 4750.118- Diatas adalah hasil output program R pemodelan regresi logistik multinomial pada data di atas dengan peubah responnya adalah ketertarikan responden akan politik dan peubah bebasnya adalah Jenis kelamin (Sex) dan warna kulit (Race)
#### Interpretasi odds
exp(summary(model1)$coefficients)## (Intercept) raceBlack raceWhite sexMale
## FI 0.7412273 1.9728785 2.431158 1.1851333
## NVI 1.0443504 0.9578007 1.128559 0.9146930
## SI 1.4228498 1.5640620 1.625599 0.8398561
## VI 0.3391894 2.5206824 3.205660 1.95332551. Persamaan prediksi odds untuk responden FI (Fairly Interested) dibanding dengan responden yang NAI (Not At All Interested) terhadap politik adalah :
Interpretasi pada tiap nilai dugaan parameter model
expa(-0.2994) = 0.7412 Untuk responden berjenis kelamin Female (perempuan) dengan other (ras lainnya), odds responden Fairly Interested (FI)) politik dibanding yang Not at All Interested (NAI) politik adalah 0.7412
expa(0.1698) = 1.185 Untuk responden berjenis kelamin Male (laki-laki) dengan other (ras lainnya), odds responden Fairly Interested (FI) politik dibanding yang Not at All Interested (NAI) politik adalah 1.1851 kali odds responden berjenis kelamin perempuan (Female)
expa(0.6795) = 1.9729 Untuk responden dengan ras black, odds responden Fairly Interested (FI) politik dibanding yang Not at All Interested (NAI) politik adalah 1.9729 kali odds ras lainnya
expa(0.8884) = 2.4311 Untuk responden dengan ras putih odds responden Fairly Interested (FI) politik dibanding yang Not at All Interested (NAI) politik adalah 2.4311 kali odds ras lainnya.
2. Persamaan prediksi odds untuk responden NVI (Not Very Interested) dibanding dengan responden yang NAI (Not At All Interested) terhadap politik adalah :
Interpretasi pada tiap nilai dugaan parameter model
expa(0.0434) = 1.0434 Untuk responden berjenis kelamin Female (perempuan) dengan Other (ras lainnya), odds responden Not Very Interested (NVI) politik dibanding yang Not at All Interested (NAI) terhadap politik adalah 1.0443
expa(-0.0892) = 0.9174 Untuk responden berjenis kelamin Male (laki-laki) dengan Other (ras lainnya), odds responden Not Very Interested (NVI) politik dibanding yang Not at All Interested politik adalah 0.9147 kali odds responden berjenis kelamin Female (perempuan)
expa(-0.0431) = 0.9587 Untuk responden dengan ras Black, odds responden Not Very Interested (NVI) politik dibanding yang Not at All Interested (NAI) politik adalah 0.9578 kali odds ras lainnya (Other)
expa(0.1209) = 1.1268 Untuk responden dengan ras putih odds responden Not Very Interested (NVI) politik dibanding yang Not at All Interested politik adalah 1.1286 kali odds ras lainnya (Other)
3. Persamaan prediksi odds untuk responden SI (Somewhat Interested) dibanding dengan responden yang NAI (Not At All Interested) terhadap Politik adalah :
Interpretasi pada tiap nilai dugaan parameter model
expa(0.3527) = 1.4223 Untuk responden berjenis kelamin Female (perempuan) dengan ras lainnya (other), odds responden Somewhat Interested politik (SAI) dibanding yang Not at All Interested terhadap politik adalah 1.4223
expa(-0.1745) = 0.839 Untuk responden berjenis kelamin male (laki-laki) dengan ras lainnya (other), odds responden Somewhat Interested (SI) politik dibanding yang Not at All Interested (NAI) politik adalah 0.8398 kali odds responden berjenis kelamin perempuan (female)
expa(0.4473) = 1.56411 Untuk responden dengan ras kulit hitam, odds responden Somewhat Interested (SI) politik dibanding yang Not at All Interested politik adalah 1.5641 kali odds ras lainnya (Other)
expa(0.4859) = 1.6256 Untuk responden dengan ras kulit putih (Race White) odds responden Somewhat Interested (SI) politik dibanding yang Not at All Interested politik adalah 1.6256 kali odds ras lainnya (Other)
4. Persamaan prediksi odds untuk responden VI (Very Interested) dibanding dengan responden yang NAI (Not At All Interested) pada Politik adalah :
Interpretasi pada tiap nilai dugaan parameter model
expa(-1.0812) = 0.33922 Untuk responden berjenis kelamin Female (perempuan) dengan ras lainnya (other), odds responden Very Interested (VI) politik dibanding yang Not at All Interested politik adalah 0.3392
expa(0.6695) = 1.9533 Untuk responden berjenis kelamin male (laki-laki) dengan ras lainnya (other), odds responden Very Interested (VI) politik dibanding yang Not at All Interested politik adalah 1.9533 kali odds responden berjenis kelamin Female (perempuan)
expa(0.9245) = 2.520 Untuk responden dengan ras Hitam (black race), odds responden Very Interested (VI) politik dibanding yang Not at All Interested politik adalah 1.2.5207 kali odds ras lainnya (other)
expa(1.1649) = 3.20577 Untuk responden dengan ras putih (white race) odds responden Very Interested (VI) politik dibanding yang Not at All Interested (NAI) politik adalah 3.2057 kali odds ras lainnya (Other)
Jika diketahui seseorang berjenis kelamin laki-laki dan ras kulit hitam. Berdasarkan model pada poin a maka ketertarikan tersebut dalam politik adalah :
#predict reglog
dataku2=data.frame(sex=as.factor("Male"),race=as.factor("Black"))
predict(model1,newdata=dataku2,"probs")## NAI FI NVI SI VI
## 0.1391374 0.2411365 0.1273036 0.2600530 0.2323694Dari hasil output diatas dapat terlihat bahwa berdasarkan model pada poin a, jika diketahui seseorang berjenis kelamin laki-laki (Male) dengan ras hitam (Black Race) maka ketertarikan seseorang terhadap politik untuk Not at All Interested adalah 0.14, untuk Fairly Interseted (FI) sebesar 0.24, untuk Not Very Interested (NVI) sebesar 0.13, untuk Somewhat Interested sebesar 0.26 dan untuk Very Interested sebesar 0.23