INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS. VOL. NO. AUGUST 2024 Analyzing Factors Contributing to Gender Inequality in Indonesia using the Spatial Geographically Weighted Logistic Ordinal Regression Model Hani Khaulasari and Yuniar Farida AbstractAiGender inequality is a condition of discrimination caused by social systems and structures. The main objective of this research is to identify factors that influence gender inequality in each province in Indonesia and obtain classification accuracy values using Geographically Weighted Ordinal Logistic Regression (GWOLR). The dataset used in this research consists of a response variable, namely the gender inequality index where the index value is divided into ordinal categories . ow, medium, and hig. and four predictor variables from the dimensions of health, education, human empowerment, social-culture, and work. The results of this study show that the classification accuracy of the GWOLR model is 85%. The mapping of provinces in Indonesia based on influential variables forms three groups. The first group . is influenced by the percentage of women who give birth with the assistance of health workers (X1 ) and the female Human Development Index (HDI) (X3 ). The second group . is influenced by the ratio of womenAos Pure Participation Rate (APM) (X2 ) and the percentage of rape crimes against women (X4 ). The third group . is influenced by the percentage of women who give birth with the assistance of health workers (X1 ), the ratio of womenAos Pure Participation Rate (APM) (X2 ), the percentage of womenAos Human Development Index (HDI) ratio (X3 ), and the percentage of womenAos rape crimes (X4 ). INTRODUCTION ENDER includes inherent characteristics, behavior, differentiation of roles, positions, responsibilities, and division of labor between men and women as a result of the formation of culture or social environment in which humans grow and are raised. Gender inequality is an unfair condition caused by social systems and structures so that women and men become victims of the system . Gender justice will be realized if conditions are created where the portions and social cycles of women and men are equal, harmonious, balanced, and harmonious . Gender inequality is a condition of discrimination between men and women due to social systems and structures or conditions where there is a gap between men and women in family, community, national, and state life . Gender justice and gender equality can be achieved by providing equal opportunities to contribute to development. Forms of gender inequality caused by gender differences are First, marginalization of access to resources, for example, information and technology, education, and employment, which results in poverty . and can impact both men and women. Second, subordination is an attitude of lowering social position or status. Third, stereotypes are labels Khaulasari and Y. Farida are with Mathematics. Sunan Ampel State Islamic University Surabaya. Indonesia email: hani. khaulasari@uinsby. Manuscript received July 18, 2023. accepted September 14, 2023 or markings of a certain group that are often detrimental and cause injustice. Fourth, violence against a personAos physical and mental psychological integrity. Fifth, women have a longer and larger workload . , . Gender inequality is detrimental to all areas of development . ocial, economic, education, defense, and securit. Gender inequality has a strong relationship with poverty and unequal access to education, health services, and access to finance . The average condition of IndonesiaAos Gender Inequality Index (IKG) in 2022 has decreased com-pared to previous years. However, the condition of gender inequality in each province is still above average and many provinces still experience gender inequality towards women. The issue of gender inequality became the National Medium Term Development Plan (RPJMN) for 2020-2024 and became the Sustainable Development Goals (SDG. for 2030 so the Indonesian Government issued INPRES Number 9 of 2000 . , . According to . Indicators of gender inequality include risk factors in the dimensions of health . aternal mortality and adolescent fertilit. , empowerment . ducation and parliamen. , and the labor market. Based on BPS research in 2020, there is a negative relationship between the Gender Empowerment Index (IDG). Human Development Index (IPM), and Gender Development Index (IPG) on gender gaps in Indonesia. IndonesiaAos gender inequality index varies, and some provinces still have prominent levels of gender inequality. The Research conducted . focus-es on analyzing gender gaps using spatial The findings of this research reveal that analysis using spatial methods was not applied in this research. The spatial data used is only a description of the location of gender inequality without dis-cussing spatial statistical analysis. The results of the OLS global regression model show that important factors that influence gender inequality are the percentage of individuals aged ten years and over who have never received a formal education, the number of women aged 20-24 years who married or lived together before the age of 18, and the percentage of women aged 20-24 years. 24-year-olds who married or lived together before the age of 18. births that occur outside health facilities. HDI of the female population, and IDG influence gender inequality in Indonesia. This research only produced a goodness of fit model of 74. this is because there are still alternative statistical methods that can provide more accurate results, namely by using spatial statistical analysis methods. Each province in Indonesia has unique characteristics and conditions that can influence gender inequality. The Central Statistics Agency (Badan Pusat INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS. VOL. NO. AUGUST 2024 Statistica. categorizes the gender inequality index into four categories ranging from low to high . The appropriate method for analyzing spatial data with an ordinal categorical response variable is Geographically Weighted Ordinal Logistic Regression (GWOLR). , . The GWOLR model in the case of human index development using exponential kernel function weighting shows superior classification accuracy compared to ordinal logistic regression, as measured by the Apparent Error Rate (APER) which reaches 86. 84% . Other research using the GWOLR method can map very well, 7% . GWOLR method is better than ordinal logistic regression [?]. The main objective of this research is to identify factors that influence gender inequality in each region in Indonesia and obtain classification accuracy values in mapping gender gaps in each province in Indonesia using Geographically Weighted Logistic Ordinal Regression. II. METHODOLOGY Logistic Ordinal Regression The ordinal logistic regression method is a statistical technique used to ana-lyze response variables with an ordinal scale, comprising three or more categories. In this model, predictor variables can include both categorical and continuous data, involving two or more variables . The logistic regression model can be referred to as the logit model. If the response variable employs an ordinal scale with A LA cate-gories and A X A refers to the vector of independent variables in AonAo observations, the OLR (Ordinal Logistic Regressio. model can be formulated as follows equation . P(Y O l | xi ) = 0l xi T k logit[[P(Y O l | xi )]] = ln 1 Oe P(Y O l | xi ) . With l = 1, 2, . L Oe 1, i = 1, 2, . , n, k = 1, 2, . , p. P (Y O l | xi ) describes the cumulative probability that is less than or equal to the l-category on xk . Geographically Weighted Logistic Ordinal Regression The GWOLR model is a fusion of the Geographically Weighted Regression (GWR) model and the ordinal logistic regression model . The purpose of the GWOLR model is to depict the relationship between an ordinal-scaled response variable and predictor variables, where each regression coefficient depends on the location of the observed data . If the response variable consists of L categories, then the GWOLR model will be used for site I. This can be expressed as Equation . The GWOLR model utilizes the Weighted Maximum Likelihood Estimation (MLE) method. Subsequently, a likelihood function in the form of ln is constructed by transforming the original likelihood function. The weighting factors in the GWOLR model are associated with geographic locations, with each location having distinct values that reflect local characteristics within the GWOLR model. Consequently, a weighting scheme is applied to the likelihood of the GWOLR Parameter estimation is performed by calculating the first partial derivative of the weighted ln-likelihood equation with respect to the parameter to be estimated, followed by equating the derivative to zero. Due to the nonlinearity of this initial partial derivative, the solution employs the NewtonRaphson iteration method Hypothesis testing on the GWOLR model includes testing the suitability between the GWOLR model and the ordinal logistic regression model, testing parameters simultaneously, and testing parameters individually. The goodness of fit test of the GWOLR model and Ordinal Logistic Regression model. This test aims to test the significance of geographical factors. The hypothesis used is . H0 : k . i , vi ) = k . here was no difference between the GWOLR model and the ordinal logistic regression H1 : At least one k . i , vi ) = k [There are differences between the GWOLR model and the ordinal logistic regression mode. Test statistics as equation . D . C) /d f 1 D . C O ) /d f 2 Reject H0 if Fount > F(,d f 1,d f . In the study by . , weighting functions generate different parameter estimations for each According to . the Kernel Exponential function dan Gaussian function can be a weighting function. Fixed Gaussian kernel weighting function can be written as follows equation. ! 1 di j 2 w j . i , vi ) = exp Oe Fixed exponential kernel weighting function can be written as follows equation . di j 2 w j . i , vi ) = exp Oe P(Y O l | xi ) logit[[P(Y O l | xi )]] = ln = 0l . i , vi ) xi T k . i , vi ) With di j = . i Oe u j ) . i Oe v j ) , which represents the dis1 Oe P(Y O l | xi ) . tance between locations . i , vi ) and . j , v j ) And b as a known The cumulative probability of the g-the response category can non-negative parameter commonly called the smoothing pabe expressed as equation . Selecting an optimal bandwidth is crucial in ensuring the modelAos accuracy towards the data. Bandwidth exp. i , vi ) xi k . i , vi )) P(Y O l | xi ) = , l = 1, 2, . L Oe 1 can be achieved by utilizing Cross Validation (CV) as follows 1 Oe exp. i , vi ) xi k . i , vi )) . The probability of the response variable at the location-I belonging to category l can be expressed as equation . CV . = Oc yi Oe yC=i . i ) . i , vi ) xi T k . i , vi )) Al . i ) = , l = 1, 2, . L Oe 1 A simultaneous test is utilized to assess the collective signifi1 Oe exp. i , vi ) xi T k . i , vi )) . cance of the variable parameters in the GWOLR model, which INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS. VOL. NO. AUGUST 2024 does not appear to be modifying the subject. The hypothesis is as follows . H0 : 1 . 1 , v1 ) = 2 . 2 , v2 ) = . = k . i , vi ) = 0 (No variables that affect the mode. H1 : at least one of k . i , vi ) = 0 . t least there is one predictor variable that affects the mode. Test Statistics as equation . d Oe ln L(O) G2 = Oe2 ln L(E) . Source Data and Variables This research uses a combination of secondary data sourced from the Central Statistics Agency (BPS) and the Ministry of WomenAos Empowerment and Child Pro-tection (KemenpA). Data was obtained through the website w. id, the gender gap publication book by BPS and the gender publication book by Kemenpp-pa. The experimental units are locations in 34 provinces in Indonesia. as depicted in Fig. d maximized value of the likelihood function Where L(E) below H0 d maximized value of the likelihood function under L(O) the population Reject H0 if G2 > N(,d f ) where d f = trace(S ) where S is a matrix with the i-th row and j-th column elements. Z j . i , v j ) Si j = Ri j Z j . i , v j ) where R is a matrix with the i-th row element Ri j = xi X T W . i , v j )dA. i , v j )X Oe1X T W . i , v j )dA. i , v j ) Afterward, the significance of the parameters is tested individually or partially with the following hypotheses . H0 : k . i , vi ) = 0, i = 1, 2, . , n. k = 1, 2, . , p (No variables that affect the respons. H1 : k . i , vi ) = 0 (Variable predictor affects the response variabl. Test Statistics as equation . Zcount = Ck . i , vi ) SE Ck . i , vi ) . H0 rejected if |Zcount | > ZOy/2 According to . the best regression model is determined by comparing the ordinal logistic regression model and GWOLR. One of the criteria used to select the best model is the model that achieves the highest classification accuracy. In this study, the classification accuracy is measured using APER (Apparent Error Rat. The corresponding error rate can be obtained by calculating the APER value. Hence, to determine the classification accuracy, one can use 1-APER. Fig. 1: Research data collection locations in 34 provinces in Indonesia The dataset in this study consisted of response variables, namely the gender inequality index . , where the index values were divided into three categories with ordinal scale, . low (IKG < 0. medium (IKG 0. 449 and . high (IKG > 0. while the predictor variable consisted of The percentage of women aged 15-49 giving birth with assistance from healthcare personnel (X1 ), the Ratio of net enrollment rate (APM) of women in Higher Education (X2 ), the Percentage of Human Development Index (IPM) Women (X3 ). Percentage of rape crimes against women (X4 ). per-centage of womenAos involvement in parliament (X5 ) and location factors point of Latitude and Point of Longitude. The analysis steps in this study, as illustrated in the flowchart Fig. TABLE I: Confusion matrix table of the classification Observation Prediction Class =1 Class=2 Class=3 Class=1 F11 F21 F31 Actual Class Class=2 Class=3 F12 F13 F22 F23 F32 F33 Referring to Table 1I, mistakes made in categorizing objects can be computed through APPER . , which is described as Equation . APER = F12 F13 F21 F23 F31 F33 F11 F12 F13 F21 F22 F23 F31 F32 F33 Meanwhile, to assess the accuracy of the classification, the following Equation . can be used: Total Accuracy Rate = 1 Oe APER Fig. 2: Flow chart of research analysis steps i. R ESULT The gender inequality index is categorized into three, namely . low (IKG < 0. medium (IKG 0. INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS. VOL. NO. AUGUST 2024 high (IKG > 0. , as presented in Fig. Based on Fig 3. The gender gap situation in every province Fig. 3: Description of the gender inequality index in Indonesia in the year 2022. in Indonesia has shown improvement from year to year. 2022, the highest gender inequality index was 0. 492, which marked a decrease compared to the 2021 figure of 0. This highest index was observed in Southeast Sulawesi. On the other hand. Bali Province consist-ently had the lowest gender inequality index from 2021 to 2022. Several provinces in Indonesia exhibit low gender inequality, including Aceh. North Sumatra. Riau. West Sumatra. Lampung. Banten. Java Island. Bali. East Nusa Tenggara (NTT). West Nusa Tenggara (NTB). Gorontalo. Yogyakarta Special Region (DIY), and Jakarta. However, there are provinces that require immediate attention from the government due to their high gender inequality These provinces are primarily located in the Eastern Indonesia (WIT) region, comprising Maluku. North Maluku. West Pa-pua, and Papua, as well as the Central Indonesia region (WITA), including Jambi. Bengkulu. West Kalimantan. Central Kalimantan. Southeast Sulawesi, and West Sulawesi. Additionally, provinces with moderate gender gaps are closely monitored to prevent them from falling into the high inequality These provinces in-clude South Sumatra. West Java. Banten. West Nusa Tenggara. South Kalimantan. North Kalimantan, and Southeast Sulawesi Based on the examination of multicollinearity, it was found that the values of the Variance Inflation Factor (VIF) for each predictor variable are below There-fore, it can be concluded that there is no significant correlation between the predic-tor variables. Ordinal Logistic Regression (OLR) model on the Gender Inequality Index in Indonesia as in equation . logit [P (Y O . ] logit [P (Y O . ] Oe66, 551 Oe 0, 401 Oe 0, 016X2 Oe 0, 366X3 Oe 0, 620X4 Oe 0, 046X5 Oe66, 551 Oe 0, 401 Oe 0, 016X2 Oe 0, 366X3 Oe 0, 620X4 Oe 0, 046X5 In the OLR model, the factors that influence global gender inequality in Indonesia are obtained The percentage of women aged 15-49 giving birth with assistance from healthcare personnel (X1 ), the Percentage of Human Development Index (IPM) Women (X3 ), and Percentage of rape crimes against women (X4 ). Next. The initial step in modeling the gender inequality index in each province of Indonesia is to determine the geographical locations based on Latitude and longitude in each region. Based on the results of the Breusch-Pagan test, it can be concluded that spatial heterogeneity exists in each province in Indonesia. Then, the Euclidean distance is calculated from locations i and j. Subsequently, the optimal bandwidth is determined using the Cross Validation (CV) Using an optimal bandwidth of 1. 091, the best weight matrix is obtained by incorporating the Euclidean distance and the optimal bandwidth into the exponential kernel weighting The weights obtained for each research location are then used to estimate the GWOLR parameters in each province of Indonesia using Newton Raphson iteration. After getting the GWOLR model estimation, tests are conducted to determine the similarity between the GWOLR model and the Ordinal Logistic Regression model regarding overall and individual parameter levels involved. The following are the hypotheses used to test the equivalence between the GWOLR model and the Ordinal Logistic Regression model can be presented Table II. TABLE II: test of equality between the GWOLR model and the Ordinal Logistic Regression model. Model OLR GWOLR Devians Deviant/Df The computed F test value for testing the equivalence between the GWOLR and Ordinal Logistic Regression models The added F test value is greater than the critical F value of F . , which is 1. Therefore, it can be concluded that there is a significant difference between the GWOLR and Ordinal Logistic Regression (OLR) models. The modeling of the Gender Inequality Index in the provinces of Indonesia in 2022 yields significantly different results between the GWOLR and Ordinal Logistic Regression models. The following is the simultaneous testing of the parameters of the GWOLR model in this study is obtained. G2 statistic value is 49. 661, compared to the N 2 . value of 11,889. Therefore, the G2 statistic value is greater than the N 2 . Hence, it can be concluded that at least one predictor variable significantly influences the Gender Inequality Index in the provinces of Indonesia. The partial testing of the parameters of the GWOLR model is conducted using the following hypotheses: H0 : k . i , vi ) = 0, i = 1, 2, . , 34. k = 1, 2, . , 5 H1 : k . i , vi ) = 0 Each province in Indonesia has a different model, resulting in other significant variables. For example, testing will be conducted on the parameter k that influences variable H1 for each i, k with i = 1, 2, . , 34 and k = 1, 2, . , 5. Assuming that k at the eleventh location . 11 , v11 ) representing East Java province has a significant effect, the Z-test values can be observed in Table i. GWOLR model on the Gender Inequality Index in East Java Province . th locatio. as in Equation . logit [P (Y O . ] logit [P (Y O . ] 85, 14 Oe 0, 6783X1 Oe 0, 1783X2 Oe 0, 2007X3 0, 0494X4 Oe 0, 006X5 80, 32 Oe 0, 6783X1 Oe 0, 1783X2 Oe 0, 2007X3 0, 0494X4 Oe 0, 006X5 Based on the model and Table 3, the variables that influence the GWOLR model in East Java Province are The predictor variables that influence the Gender Inequality Index in East Java Province are the percentage of women aged 15-49 years who gave birth with the assistance of health workers (X1 ) and the percentage of women on the Human Development Index INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS. VOL. NO. AUGUST 2024 TABLE i: Partial Parameter Significance Test in Location of East Java Province Parameter Estimate 01 . 11 , v11 ) 85,14 02 . 11 , v11 ) 80,32 1 . 11 , v11 ) -0,6783 2 . 11 , v11 ) -0,1783 3 . 11 , v11 ) -0,2007 4 . 11 , v11 ) 0,0494 5 . 11 , v11 ) -0,006 *Significant of = 5% 24,456 23,45 0,2345 0,1732 0,1003 0,0426 0,1153 Zccount 0,905 1,208 -2,893 -1,029 -2,001 1,160 -0,052 P-Value 0,817 0,886 0,001* 0,152 0,002* 0,877 0,479 Odd ratio 1,971 1,195 1,222 1,051 1,006 (HDI) (X3 ). A one percent decrease in the percentage of evermarried women aged 15-49 years whose childbirth is assisted by health workers (X1 ) is associated with East Java being 971 times more likely to become a high gender inequality Similarly, a one percent decrease in the percentage of women in the Human Development Index (HDI) (X3 ) is associated with the East Java region being 1. 222 times more likely to become a high gender inequality area. Mapping of provinces based on significant variables in the GWOLR model can be seen in Fig 4. The comparison of model logistic regression Fig. 4: Mapping of provinces based on significant variables in the GWOLR model using the Expo-nential Kernel Function and GWOLR is conducted by comparing the accuracy rate and the error value (APER) in classification as present-ed by Table IV. TABLE IV: shows that in modeling the Gender Inequality Index, the use of the GWOLR model yields a higher classification accuracy than the ordinal logistic regression model. The GWOLR model is more suitable for modeling the Gender Inequality Index in Indonesian provinces than the ordinal logistic regression model. Model Ordinal Logistic Regression GWOLR APER Accuracy IV. C ONCLUSION The Factors that influence the Gender Inequality Index according to the GWOLR model can be categorized into several dimensions, such as Health. Educa-tion. Community Empowerment, and Socio-Culture. The factors that influence the Gender Inequality Index in each province are different. The mapping of provinces in Indonesia based on influential variables forms three groups. The first group . is influenced by the percentage of women aged 15-49 years who gave birth with the assistance of health workers (X. and the share of womenAos Human Development Index (HDI) (X3 ). Provinces included in group one are Aceh. Bali. DI Yogyakarta. DKI Jakarta. Gorontalo. Central Java. East Java. East Kalimantan. Bangka Belitung Islands. Riau Islands. Lampung. East Nusa Tenggara. Riau. South Sulawesi. North Sulawesi. West Sumatra. South Sumatra, and North Sumatra. The second group . is influenced by the ratio of pure enrollment rates (APM) of women in higher education (X2 ) and the percentage of rape crimes against women (X4 ). Provinces in the second group are West Java. South Kalimantan. North Kalimantan. West Nusa Tenggara, and Central Sulawesi. The third group . is influenced by the percentage of women aged 15-49 years who gave birth with the assistance of health workers (X1 ), the ratio of the Pure Participation Rate (APM) of women in higher education (X2 ), the percentage of the Human Development Index ratio (HDI) of women (X3 ), and the percentage of rape crimes against women (X4 ). Provinces in the third group are Banten. Bengkulu. Jambi. West Kalimantan. Central Kalimantan. Maluku. North Maluku. Papua. West Pa-pua. West Sulawesi, and Southeast Sulawesi. Modeling the factors influencing the 2022 Gender Inequality Index using the GWOLR method with exponential kernel weighting has proven to be superior to using ordinal logistic regression based on classification accuracy values. These findings are consistent with previous research . , . The results obtained show that the GWOLR model is superior to the global ordinal logistic regression model. How-ever, the level of accuracy obtained in this study does not exceed previous The research produced a classification accuracy of R EFERENCES . Judiasih. AuImplementasi Kesetaraan Gender Dalam Beberapa Aspek Kehidupan Bermasyarakat Di Indonesia,Ay ACTA DIURNAL Jurnal Ilmu Hukum Kenotariatan, vol. 5, no. Art. Jun. 2022, doi: 23920/acta. Rahim. AuGender Dalam Perspektif Islam,Ay SOSIORELIGIUS: JURNAL ILMIAH SOSIOLOGI AGAMA, vol. 3, no. Art. Oct. 2018, doi: 10. 24252/sosioreligius. Ponthieux and D. Meurs. AuChapter 12 - Gender Inequality,Ay in Handbook of Income Distribution. Atkinson and F. Bourguignon. Eds. , vol. Elsevier, 2015, pp. 981Ae1146, doi: 10. 1016/B978-0-44459428-0. Valsecchi. Iacoviello. Jacques. Borinca, and J. Pichastor. AuMenAos Gender Norms and Gender-Hierarchy-Legitimizing Ideologies: The Effect of Priming Traditional Masculinity Versus a Feminization of MenAos Norms,Ay Gender Issues, vol. 40, no. 1, pp. 465Ae480, 2023, doi: 1007/s12147-022-09308-8. Fadila. Riyanto. Avisya. Irianti, and D. Radianto. AuKesetaraan Gender,Ay Humantech: Jurnal Ilmiah Multidisiplin Indonesia, vol. 2, no. Art. Jun. Aprilia and M. Triani. AuAnalisis Pengaruh Ketimpangan Gender. Rasio Ketergantungan Dan Kesehatan Terhadap Kemiskinan Di Indonesia,Ay Jurnal Kajian Ekonomi dan Pembangunan, vol. 4, no. Art. Sep. 2022, doi: 10. 24036/jkep. BPS RI. AuIndeks Ketimpangan Gender (IKG) 2022,Ay w. INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS. VOL. NO. AUGUST 2024 Razi. Rencana Strategis Kementerian Agama Tahun 2020-2024. Jakarta: Kementerian Agama Republik Indonesia, 2020. Tim Kemen pA. Pembangunan Manusia Berbasis Gender. Jakarta: Kementerian Pemberdayaan Perempuan dan Perlindungan Anak Republik Indonesia, 2022. [Onlin. Available: https://w. BPS. Kajian Penghitungan Indeks Ketimpangan Gender 2021. Jakarta: BPS Republik Indonesia, 2021. Marsono. AuDeteksi Spasial Pada Model Indeks Ketimpangan Gender Indonesia,Ay Buana Gender, vol. 6, no. 1, pp. 49Ae66, 2021. TIM BPS RI. Kajian Penghitungan Indeks Ketimpangan Gender Jakarta: BPS-Statistics Indonesia, 2022. [Onlin. Available: https://w. id/publication/ . Zuhdi. Saputro, and P. Widyaningsih. AuParameters Estimation of Geographically Weighted Ordinal Logistic Regression (GWOLR) Model,Ay J. Phys. : Conf. Ser. , vol. 855, p. Jun. 2017, doi: 1088/1742-6596/855/1/012064. Amalah. Jaya, and N. Sirajang. AuPemodelan Geographically Weighted Logistic Regression dengan Metode Ridge,Ay ESTIMASI: Journal of Statistics and Its Application, pp. 130Ae143. Aug. 2023, doi: 20956/ejsa. Pradita. Yasin, and D. Safitri. AuPemodelan Faktor-Faktor Yang Mempengaruhi Indeks Pembangunan Manusia Kabupaten/Kota Di Jawa Timur Menggunakan Geographically Weighted Ordinal Logistic Regression,Ay Jurnal Gaussian, vol. 4, no. Art. Jul. 2015, doi: 14710/j. Rochmah and V. Ratnasari. AuPemodelan Ketahanan Pangan di Jawa Timur Menggunakan Metode Geographically Weighted Ordinal Logistic Regression (GWOLR),Ay JSSITS, vol. 8, no. 2, pp. D397AeD404. Feb. 2020, doi: 10. 12962/j23373520. Fathurahman. AuHypothesis testing of Geographically weighted bivariate logistic regression,Ay Journal of Physics: Conference Series, vol. 1417, pp. 1Ae8, 2019, doi: 10. 1088/1742-6596/1417/1/012008. Wardani. Gusriani, and D. Kusuma. AuPemetaan Zonasi Resiko Covid-19 Di Provinsi Jawa Barat Menggunakan Model Geographically Weighted Ordinal Logistic Regression (GWOLR),Ay Teorema: Teori dan Riset Matematika, vol. 7, no. 1, pp. 193-204, 2022. Sari and P. Purhadi. AuPemodelan Indeks Pembangunan Manusia Provinsi Jawa Barat. Jawa Timur dan Jawa Tengah Tahun 2019 dengan Menggunakan Metode Regresi Logistik Ordinal,Ay Jurnal Gaussian, vol. 10, no. Art. Feb. 2021, doi: 10. 14710/j. Mishra. Kumar. Prasad, and M. Punia. AuGeographically Weighted Method Integrated with Logistic Regression for Analyzing Spatially Varying Accuracy Measures of Remote Sensing Image Classification,Ay Journal of the Indian Society of Remote Sensing, vol. 49, no. 5, pp. 1189Ae1199, 2021, doi: 10. 1007/s12524-020-01286-2. Dong. Nakaya, and C. Brunsdon. AuGeographically weighted regression models for ordinal categorical response variables: An application to geo-referenced life satisfaction data,Ay Computers. Environment and Urban Systems, vol. 70, pp. 35Ae42. Jul. 2018, doi: 1016/j. Fotheringham and T. Oshan. AuGeographically weighted regression and multicollinearity: Dispelling the myth,Ay Journal of Geographical Systems, vol. 18, pp. 303Ae329, 2016. Syinfi. Pemodelan Rata-rata Umur Kawin Pertama (UKP) Wanita Di Provinsi Jawa Timur Tahun 2012 Dengan Pendekatan Geographically Weighted Ordinal Logistic Regression (GWOLR). Thesis Program Statistika. Surabaya: ITS, 2015. Fadliana. Pramoedyo, and R. Fitriani. AuImplementation Of Locally Compensated Ridge-Geographically Weighte Regression Model In Spatial Data With Multicollinearity Problems(Case Study: Stunting among Children Aged under Five Years in East Nusa Tenggara Provinc. ,Ay MEDIA STATISTIKA, vol. 13, no. 2, pp. 125Ae135. Dec. Hasanah. AuPerbandingan Metode Klasifikasi Artificial Neural Network Backprogation dan Regresi Logistik (Studi Kasus: Bank Internasional Indonesi. ,Ay Jurnal Statistika dan Matematika, vol. 1, no. Feb. 2019, doi: 10. 32493/sm. Anselin. Spatial Econometrics: Methods and Models. Netherlands: Kluwer Academic Publishers, 2013. Aldila. Cahya, and K. Suryadi. AuSpatio-Temporal Analysis for Detecting COVID-19 Spreading Pattern in Indonesia,Ay Journal of Mathematical and Fundamental Sciences, vol. 52, no. 3, pp. 250Ae268. Aug. 2020, doi: 10. 5614/j. Long. Regression Models for Categorical and Limited Dependent Variables. California: SAGE Publications. Inc. , 1997.