Wages Index Final - only employed
Amra Fetahovic
2023-11-24
Ismir mail 03.01.
Ovaj wage model sadrzi samo one osobe koje su zaposlene, a dodacu i novu dummy varijablu full time dummy ili part time dummy.
To znaci da cu da gledam samo Active 1 i 2 tj. full-time i part time, ali pod vrstom prihoda ovih uzimam sve one koje sam koristila i u prethodnom modelu.
library(haven)
library(dplyr)##
## Attaching package: 'dplyr'## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionsetwd("C:/Users/Amra/Documents/HBS database inspecting")
HBS_2015_HH_short <- read_sav("~/HBS database inspecting/HBS_2015_HH - short.sav")
hbs_short <- HBS_2015_HH_short01 OZNACAVA DA SE RADI O PRVOM CLANU DOMACINSTVA KOJI JE PRIMAO ZARADE YN Oznacava DA-NE tj. odnosi se na da li su ili nisu primali naknade IN Mjesecni iznos odredjene naknade NM Broj mjeseci za koji je primljen taj iznos
Sabracu za svakog clana domacinstva koji su rekli da su primili neku od nakanda i vidjeti koliko mi to pokriva opstina
hbs_short$amr.NEWINCOM <- rowSums(hbs_short[,c("Q02_YN_01_S12","Q02_YN_02_S12","Q02_YN_03_S12","Q02_YN_04_S12","Q02_YN_05_S12","Q02_YN_06_S12", "Q02_YN_07_S12", "Q02_YN_08_S12", "Q02_YN_09_S12","Q02_YN_10_S12","Q02_YN_11_S12","Q02_YN_12_S12", "Q02_YN_13_S12")],na.rm = T)
table(hbs_short$amr.NEWINCOM)##
## 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
## 44 592 1583 1297 1598 900 697 500 231 134 63 39 11 7 1 3
## 18
## 2Prema ovoma samo u 44 domacinstava ni jedan clan ne dobiva nsita sto je navedeno u piatanjama QA01 - QA23.
Plan za dalje:
- Za svaku varijablu koja ima neki inkom izracunati mjesecne izdatke i za svakog clana. jer imaju i iznos a ima i koliko mjesecno je primao, tako da to treba mnoziti sa kolonom broj mjeseci pa dijeli sa brojem 12.
- Nakon toga cu da izracunam mjesecna primanja za svakog clana domacinstva
- Necu racunati allowance jer to nije zarada niti cu stavljati iznamljivanje zemlje, garaze i slicen rente jer i to mi nije zarada.
- Onda cu to sve sabrati i imacu mjesecnu zaradu po domacinstvu.
- Zadrzavam samo one koje su Active 1 i Active 2
- Formiram full time dummy
ZA PLATE UZIMAM SVAKOG CLANA PORODICE POSEBNO
- NECU BRISATI _S2, _S6, _S7, _S10, COICOP, jer cu mozda samo za ove likove racunati transport.
- Samo prvih osam clanova razmatram, jer npr za 9-og clana QA01_YN_09_S12 samo je jedan primio remitencices, 10 nije nisa,
- _S11 brisati to su penzije, QBizbaciti to su penzije, a ostaviti QA, to je income, izbaciti QC, QD,
library(dplyr)
hbs_wage <- hbs_short %>% select(contains(c( "SifraDom","Entity","Canton","naziv_opcine","RurUrb","ncomp","finweight","Relat","Sex","Age", "Educ","Activ", "HH_Monthly", "QA", "_S2", "_S6", "_S7", "_S10", "COICOP")))#remove additional unnecessary variables
hbs_wage <- hbs_wage %>% select(-c(M_Entity, M_Canton_Region, Relat_9, Relat_10, Relat_11, Relat_12, Relat_13, Sex_9, Sex_10, Sex_11, Sex_12, Sex_13, Age_9, Age_10, Age_11, Age_12, Age_13, Educ_9, Educ_10, Educ_11, Educ_12,Educ_13,Activ_9, Activ_10, Activ_11, Activ_12, Activ_13))#izbacumjem nepotrebne varijable iz Educ (Class, now, Level)
hbs_wage <- hbs_wage %>% select(!contains(c("Class","now","Level","QA16","QA17","QA18", "QA19", "QA20", "QA22", "QA21","_09_S12", "_10_S12")))FAZA 2
SADA PRVO FORMIRAM WAGE ZA SVAKOG CLANA DOMACINSTVA ZA INCOME
Mjesecno all wage for 1st household
library(haven)
hbs_wage$QA01_01_mj <- (hbs_wage$QA01_IN_01_S12*hbs_wage$QA01_NM_01_S12)/12
hbs_wage$QA02_01_mj <- (hbs_wage$QA02_IN_01_S12*hbs_wage$QA02_NM_01_S12)/12
hbs_wage$QA03_01_mj <- (hbs_wage$QA03_IN_01_S12*hbs_wage$QA03_NM_01_S12)/12
hbs_wage$QA04_01_mj <- (hbs_wage$QA04_IN_01_S12*hbs_wage$QA04_NM_01_S12)/12
hbs_wage$QA05_01_mj <- (hbs_wage$QA05_IN_01_S12*hbs_wage$QA05_NM_01_S12)/12
hbs_wage$QA06_01_mj <- (hbs_wage$QA06_IN_01_S12*hbs_wage$QA06_NM_01_S12)/12
hbs_wage$QA07_01_mj <- (hbs_wage$QA07_IN_01_S12*hbs_wage$QA07_NM_01_S12)/12
hbs_wage$QA08_01_mj <- (hbs_wage$QA08_IN_01_S12*hbs_wage$QA08_NM_01_S12)/12
hbs_wage$QA09_01_mj <- (hbs_wage$QA09_IN_01_S12*hbs_wage$QA09_NM_01_S12)/12
hbs_wage$QA10_01_mj <- (hbs_wage$QA10_IN_01_S12*hbs_wage$QA10_NM_01_S12)/12
hbs_wage$QA11_01_mj <- (hbs_wage$QA11_IN_01_S12*hbs_wage$QA11_NM_01_S12)/12
hbs_wage$QA12_01_mj <- (hbs_wage$QA12_IN_01_S12*hbs_wage$QA12_NM_01_S12)/12
hbs_wage$QA13_01_mj <- (hbs_wage$QA13_IN_01_S12*hbs_wage$QA13_NM_01_S12)/12
hbs_wage$QA14_01_mj <- (hbs_wage$QA14_IN_01_S12*hbs_wage$QA14_NM_01_S12)/12
hbs_wage$QA15_01_mj <- (hbs_wage$QA15_IN_01_S12*hbs_wage$QA15_NM_01_S12)/12
hbs_wage$QA23_01_mj <- (hbs_wage$QA23_IN_01_S12*hbs_wage$QA23_NM_01_S12)/12Mjesecno sva primanja za prvog householda
hbs_wage$QA_01_all_mj <- rowSums(hbs_wage[,c("QA01_01_mj","QA02_01_mj","QA03_01_mj","QA04_01_mj","QA05_01_mj","QA06_01_mj", "QA07_01_mj", "QA08_01_mj", "QA09_01_mj", "QA10_01_mj", "QA11_01_mj", "QA12_01_mj","QA13_01_mj", "QA14_01_mj", "QA15_01_mj", "QA23_01_mj")],na.rm = T)
summary(hbs_wage$QA_01_all_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 0.0 0.0 290.0 466.7 18350.0plot(density(hbs_wage$QA_01_all_mj))
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za prvog clana
hbs_wage$QAnotall_01_mj <- rowSums(hbs_wage[,c("QA01_01_mj","QA02_01_mj","QA07_01_mj", "QA08_01_mj", "QA09_01_mj", "QA14_01_mj", "QA15_01_mj")],na.rm = T)
summary(hbs_wage$QAnotall_01_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 0.0 0.0 244.1 369.2 18350.0plot(density(hbs_wage$QAnotall_01_mj))
mjesecno all wage for 2nd household
Mjesecno sva primanja za drugog householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za drugog clana
hbs_wage$QAnotall_02_mj <- rowSums(hbs_wage[,c("QA01_02_mj","QA02_02_mj","QA07_02_mj", "QA08_02_mj", "QA09_02_mj", "QA14_02_mj", "QA15_02_mj")],na.rm = T)
summary(hbs_wage$QAnotall_02_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 0.0 0.0 127.5 0.0 12833.3plot(density(hbs_wage$QAnotall_02_mj))
mjesecno all wage for 3nd household
Mjesecno sva primanja za treci householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za treceg clana
hbs_wage$QAnotall_03_mj <- rowSums(hbs_wage[,c("QA01_03_mj","QA02_03_mj","QA07_03_mj", "QA08_03_mj", "QA09_03_mj", "QA14_03_mj", "QA15_03_mj")],na.rm = T)
summary(hbs_wage$QAnotall_03_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 67.53 0.00 12000.00plot(density(hbs_wage$QAnotall_03_mj))
mjesecno all wage for 4th household
Mjesecno sva primanja za cetvrta householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za cetvrtog clana
hbs_wage$QAnotall_04_mj <- rowSums(hbs_wage[,c("QA01_04_mj","QA02_04_mj","QA07_04_mj", "QA08_04_mj", "QA09_04_mj", "QA14_04_mj", "QA15_04_mj")],na.rm = T)
hbs_wage$QAnotall_04_mj <- as.numeric(hbs_wage$QAnotall_04_mj)
summary(hbs_wage$QAnotall_04_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 18.41 0.00 2750.00str(hbs_wage$QAnotall_04_mj)## num [1:7702] 0 0 0 0 0 0 0 0 0 0 ...plot(density(hbs_wage$QAnotall_04_mj))
mjesecno all wage for 5th household
Mjesecno sva primanja za petog householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za petog clana
hbs_wage$QAnotall_05_mj <- rowSums(hbs_wage[,c("QA01_05_mj","QA02_05_mj","QA07_05_mj", "QA08_05_mj", "QA09_05_mj", "QA14_05_mj", "QA15_05_mj")],na.rm = T)
hbs_wage$QAnotall_05_mj <- as.numeric(hbs_wage$QAnotall_05_mj)
summary(hbs_wage$QAnotall_05_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 4.682 0.000 1500.000str(hbs_wage$QAnotall_05_mj)## num [1:7702] 0 0 0 0 0 0 0 0 0 0 ...plot(density(hbs_wage$QAnotall_05_mj))
mjesecno all wage for 6th household
Mjesecno sva primanja za sesti householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za sestog clana
hbs_wage$QAnotall_06_mj <- rowSums(hbs_wage[,c("QA01_06_mj","QA02_06_mj","QA07_06_mj", "QA08_06_mj", "QA09_06_mj", "QA14_06_mj", "QA15_06_mj")],na.rm = T)
hbs_wage$QAnotall_06_mj <- as.numeric(hbs_wage$QAnotall_06_mj)
summary(hbs_wage$QAnotall_06_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 1.657 0.000 1200.000str(hbs_wage$QAnotall_06_mj)## num [1:7702] 0 0 0 0 0 0 0 0 0 0 ...plot(density(hbs_wage$QAnotall_06_mj))
mjesecno all wage for 7th household
Mjesecno sva primanja za sedmog householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za petog clana
hbs_wage$QAnotall_07_mj <- rowSums(hbs_wage[,c("QA01_07_mj","QA02_07_mj","QA07_07_mj", "QA08_07_mj", "QA09_07_mj", "QA14_07_mj", "QA15_07_mj")],na.rm = T)
hbs_wage$QAnotall_07_mj <- as.numeric(hbs_wage$QAnotall_07_mj)
summary(hbs_wage$QAnotall_07_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0000 0.0000 0.0000 0.4068 0.0000 1200.0000str(hbs_wage$QAnotall_07_mj)## num [1:7702] 0 0 0 0 0 0 0 0 0 0 ...plot(density(hbs_wage$QAnotall_07_mj))
mjesecno all wage for 8th household
Mjesecno sva primanja za osmog householda
Mjesecno samo QA01_, QA02, QA_07, QA_08, QA_09, QA14, QA15 za osmog clana
hbs_wage$QAnotall_08_mj <- rowSums(hbs_wage[,c("QA01_08_mj","QA02_08_mj","QA07_08_mj", "QA08_08_mj", "QA09_08_mj", "QA14_08_mj", "QA15_08_mj")],na.rm = T)
hbs_wage$QAnotall_08_mj <- as.numeric(hbs_wage$QAnotall_08_mj)
summary(hbs_wage$QAnotall_08_mj)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00000 0.00000 0.00000 0.02597 0.00000 200.00000str(hbs_wage$QAnotall_08_mj)## num [1:7702] 0 0 0 0 0 0 0 0 0 0 ...plot(density(hbs_wage$QAnotall_08_mj))
Samo topli obrok i prevoz za sve
Ovo racunam da vidim da li i u kojoj mjeri razlikuju prevoz i topli obrok
hbs_wage$QA07_mj_all <- rowSums(hbs_wage[,c("QA07_01_mj","QA07_02_mj","QA07_03_mj","QA07_04_mj","QA07_05_mj","QA07_06_mj", "QA07_07_mj", "QA07_08_mj")],na.rm = T) #topli obrok
hbs_wage$QA08_mj_all <- rowSums(hbs_wage[,c("QA08_01_mj","QA08_02_mj","QA08_03_mj","QA08_04_mj","QA08_05_mj","QA08_06_mj", "QA08_07_mj", "QA08_08_mj")],na.rm = T) #prevoztapply (hbs_wage$QA07_mj_all, hbs_wage$Entity, summary) #topl ob## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 44.89 0.00 1000.00
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 8.143 0.000 800.000
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0000 0.0000 0.0000 0.7301 0.0000 100.0000tapply (hbs_wage$QA08_mj_all, hbs_wage$Entity, summary) # prev## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 6.261 0.000 500.000
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 2.239 0.000 400.000
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0000 0.0000 0.0000 0.2212 0.0000 100.0000hbs_wage %>% group_by(Entity) %>% summarise(mean07= mean(QA07_mj_all, na.rm=T), mean08 = mean(QA08_mj_all,na.rm=T) , n())## # A tibble: 3 × 4
## Entity mean07 mean08 `n()`
## <dbl> <dbl> <dbl> <int>
## 1 1 44.9 6.26 4643
## 2 2 8.14 2.24 2607
## 3 3 0.730 0.221 452NOVI FRAME SAMO SA _mj VARIJABLAMA I SOCIO-EK
hbs_wage.NEW <- hbs_wage %>% select(contains(c("SifraDom", "Entity", "naziv_Opcine", "Relat", "Age", "Activ", "Educ","Rur","Sex", "all_mj","notall_","QA07_0","QA08_0")))#all_mj dodajem radi ukupnih prihoda po clanu a Q07_i Q08_ radi toplog obr i prevoza po svakom clanu; notall_mj su mjesecni incomi od rada po osnovu zaposlenosti, rada u bih ali strani poslodavac, usluge(npr.pravnici),poljoprivreda, topli obrok prevoz i smjestaj
names(hbs_wage.NEW)## [1] "SifraDom" "Entity" "naziv_Opcine" "Relat_1"
## [5] "Relat_2" "Relat_3" "Relat_4" "Relat_5"
## [9] "Relat_6" "Relat_7" "Relat_8" "Age_1"
## [13] "Age_2" "Age_3" "Age_4" "Age_5"
## [17] "Age_6" "Age_7" "Age_8" "Activ_1"
## [21] "Activ_2" "Activ_3" "Activ_4" "Activ_5"
## [25] "Activ_6" "Activ_7" "Activ_8" "Educ_1"
## [29] "Educ_2" "Educ_3" "Educ_4" "Educ_5"
## [33] "Educ_6" "Educ_7" "Educ_8" "RurUrb"
## [37] "Sex_1" "Sex_2" "Sex_3" "Sex_4"
## [41] "Sex_5" "Sex_6" "Sex_7" "Sex_8"
## [45] "QA_01_all_mj" "QA_02_all_mj" "QA_03_all_mj" "QA_04_all_mj"
## [49] "QA_05_all_mj" "QA_06_all_mj" "QA_07_all_mj" "QA_08_all_mj"
## [53] "QAnotall_01_mj" "QAnotall_02_mj" "QAnotall_03_mj" "QAnotall_04_mj"
## [57] "QAnotall_05_mj" "QAnotall_06_mj" "QAnotall_07_mj" "QAnotall_08_mj"
## [61] "QA07_01_mj" "QA07_02_mj" "QA07_03_mj" "QA07_04_mj"
## [65] "QA07_05_mj" "QA07_06_mj" "QA07_07_mj" "QA07_08_mj"
## [69] "QA08_01_mj" "QA08_02_mj" "QA08_03_mj" "QA08_04_mj"
## [73] "QA08_05_mj" "QA08_06_mj" "QA08_07_mj" "QA08_08_mj"FAZA 3- FROM WIDE TO LONG
to summarize specific variables - code
library(tidyr)
hbs_wage.NEW.long <- hbs_wage.NEW %>% pivot_longer(cols = c(starts_with("Age"), starts_with("Relat"), starts_with("Activ"),starts_with("Educ"),starts_with ("Sex") ,contains("_mj")),names_to = c(".value", "Relat"), names_pattern = "(\\w+).*?(\\d)")dim(hbs_wage.NEW.long)## [1] 61616 14names(hbs_wage.NEW.long) #QA07_0 je topli obrok za svakog clana QA_0 je ukupni wage za svakog clana i QA08_0 je prevoz za svakog clana QAnotall wage ali redukovani## [1] "SifraDom" "Entity" "naziv_Opcine" "RurUrb" "Relat"
## [6] "Age_" "Relat_" "Activ_" "Educ_" "Sex_"
## [11] "QA_0" "QAnotall_0" "QA07_0" "QA08_0"prvo izbaci godine NA i sve ispod 15 godina i preko 80
hbs_wage.NEW.long.1 <- hbs_wage.NEW.long %>% dplyr::filter(14 < Age_ & Age_<80)
dim(hbs_wage.NEW.long.1)## [1] 18548 14names(hbs_wage.NEW.long.1)## [1] "SifraDom" "Entity" "naziv_Opcine" "RurUrb" "Relat"
## [6] "Age_" "Relat_" "Activ_" "Educ_" "Sex_"
## [11] "QA_0" "QAnotall_0" "QA07_0" "QA08_0"#18548 obseryvacijahead(hbs_wage.NEW.long.1)## # A tibble: 6 × 14
## SifraDom Entity naziv_Opcine RurUrb Relat Age_ Relat_ Activ_ Educ_ Sex_
## <dbl> <dbl> <chr> <dbl+lbl> <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 2258 1 B.Grahovo 1 [Urban] 1 75 1 8 5 1
## 2 2258 1 B.Grahovo 1 [Urban] 2 79 2 8 5 2
## 3 2258 1 B.Grahovo 1 [Urban] 3 47 3 4 5 1
## 4 2273 1 B.Grahovo 1 [Urban] 2 54 3 4 5 1
## 5 2299 1 B.Grahovo 1 [Urban] 1 71 1 8 5 1
## 6 2299 1 B.Grahovo 1 [Urban] 2 64 2 5 4 2
## # ℹ 4 more variables: QA_0 <dbl>, QAnotall_0 <dbl>, QA07_0 <dbl>, QA08_0 <dbl>#QA_O je wage tj.sva posmatrana primanja
table (hbs_wage.NEW.long.1$Activ_) #koji su employed 6429##
## 1 2 3 4 5 6 7 8 9
## 5676 753 1189 1600 3729 1918 433 3066 184tapply(hbs_wage.NEW.long.1$QA_0, hbs_wage.NEW.long.1$Activ_, summary)## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 250.0 600.0 656.7 920.0 18350.0
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 50.0 175.0 250.9 333.3 3000.0
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 69.44 25.00 2193.33
##
## $`4`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 0.0 0.0 10.3 0.0 1000.0
##
## $`5`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 11.33 0.00 3000.00
##
## $`6`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 2.678 0.000 1200.000
##
## $`7`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 1.447 0.000 266.667
##
## $`8`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 13.36 0.00 3333.33
##
## $`9`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 59.46 0.00 1500.00tapply(hbs_wage.NEW.long.1$QAnotall_0, hbs_wage.NEW.long.1$Activ_, summary)## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 50.0 520.0 580.6 850.0 18350.0
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0 0.0 100.0 192.8 285.0 3000.0
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 41.38 0.00 2193.33
##
## $`4`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 5.628 0.000 1000.000
##
## $`5`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 8.912 0.000 1800.000
##
## $`6`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 1.521 0.000 1200.000
##
## $`7`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 1.186 0.000 266.667
##
## $`8`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 10.35 0.00 3333.33
##
## $`9`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 28.38 0.00 1500.00# Zadrzi samo Activ_==1 i/ili ==2 tj. samo one koji su zaposleni
hbs_wage.NEW.long.1 <- hbs_wage.NEW.long.1 %>% dplyr::filter(Activ_==1 | Activ_==2)
table (hbs_wage.NEW.long.1$Activ_)##
## 1 2
## 5676 753dim(hbs_wage.NEW.long.1) #dakle 6429 je ili bilo puno zaposleno ili pola radnog vremena## [1] 6429 14table (hbs_wage.NEW.long.1$QA_0>0)##
## FALSE TRUE
## 1185 5244#od svih koji su radili 6034 ih je koji su imali prihode vece od 0
table (hbs_wage.NEW.long.1$QAnotall_0 >0)##
## FALSE TRUE
## 1650 4779#od svih koji su radili 5280 ih je koji su imali prihode vece od 0#FAZA 4 SREDI QA_0 ### Remove zeros
hbs_wage.NEW.long.2 <- hbs_wage.NEW.long.1 %>% dplyr::filter(QA_0>0)
summary(hbs_wage.NEW.long.2$QA_0)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.083 400.000 630.000 746.782 977.000 18350.000#FAZA 4a SREDI QAnotall_0 ### Remove zeros
hbs_wage.NEW.long.2a <- hbs_wage.NEW.long.1 %>% dplyr::filter(QAnotall_0>0)
summary(hbs_wage.NEW.long.2a$QAnotall_0)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 4.0 400.0 616.7 719.9 950.0 18350.0BRISATI PERCENTILE
dim(hbs_wage.NEW.long.2)#imam samo 5163 varijabli## [1] 5244 14hbs_wage.NEW.long.2 %>%
select(naziv_Opcine) %>%
group_by(naziv_Opcine) %>%
summarise (
n = n()
)%>%
arrange (desc(n))%>%
print(n=137) #izlistalo 28 opstina koje imaju manje od 10 obs## # A tibble: 130 × 2
## naziv_Opcine n
## <chr> <int>
## 1 Brcko 312
## 2 Bijeljina 245
## 3 Banja Luka 236
## 4 N.G.Sarajevo 217
## 5 Grad Mostar 190
## 6 Zenica 174
## 7 Travnik 136
## 8 Tuzla 104
## 9 Bihac 97
## 10 Bugojno 93
## 11 Sanski Most 90
## 12 Samac 85
## 13 Prijedor 84
## 14 Vitez 81
## 15 N.Travnik 80
## 16 Doboj 78
## 17 Tesanj 74
## 18 Cazin 72
## 19 Zavidovici 71
## 20 Zivinice 70
## 21 Konjic 69
## 22 C.Sarajevo 68
## 23 N.Sarajevo 67
## 24 Zvornik 67
## 25 V.Kladusa 62
## 26 Ilidza 61
## 27 Novi Grad 61
## 28 Modrica 60
## 29 Visoko 60
## 30 S.Brijeg 59
## 31 Prnjavor 55
## 32 D.Vakuf 52
## 33 Gradacac 51
## 34 S.G.Sarajevo 51
## 35 Gracanica 50
## 36 Lopare 48
## 37 Teslic 47
## 38 Laktasi 45
## 39 Maglaj 45
## 40 Derventa 44
## 41 G.Vakuf 44
## 42 Kakanj 42
## 43 Kalesija 42
## 44 Ljubuski 42
## 45 Trebinje 41
## 46 Gorazde 40
## 47 Lukavac 39
## 48 Brod 38
## 49 Kiseljak 38
## 50 B.Krupa 36
## 51 Livno 36
## 52 Busovaca 33
## 53 Orasje 33
## 54 Grude 32
## 55 Foca - RS 31
## 56 Citluk 30
## 57 Capljina 28
## 58 Gacko 28
## 59 Banovici 26
## 60 Vogosca 26
## 61 Ugljevik 25
## 62 Kresevo 24
## 63 Pale 24
## 64 Gradiska 22
## 65 Rogatica 22
## 66 Rudo 22
## 67 Jajce 21
## 68 Posusje 21
## 69 Tomislavgrad 21
## 70 Drvar 20
## 71 Ilijas 19
## 72 Ribnik 19
## 73 Srbac 19
## 74 I.Ilidza 18
## 75 Kljuc 18
## 76 Visegrad 18
## 77 Vlasenica 18
## 78 Breza 17
## 79 Prozor 17
## 80 Stolac 17
## 81 Buzim 16
## 82 Cajnice 16
## 83 Kupres 16
## 84 N.Gorazde 16
## 85 K.Dubica 15
## 86 Kladanj 15
## 87 Fojnica 14
## 88 Zepce 14
## 89 Bileca 13
## 90 Bratunac 13
## 91 Glamoc 13
## 92 Kotor Varos 13
## 93 B.Petrovac 12
## 94 Knezevo 12
## 95 Nevesinje 12
## 96 Sokolac 12
## 97 Celinac 11
## 98 Doboj-Jug 11
## 99 Mrkonjic 11
## 100 Srebrenica 11
## 101 Srebrenik 11
## 102 Doboj-Istok 9
## 103 Hadzici 9
## 104 Milici 8
## 105 Odzak 8
## 106 Pelagicevo 8
## 107 Teocak 8
## 108 B.Grahovo 7
## 109 Han Pijesak 7
## 110 Kostajnica 7
## 111 Celic 6
## 112 D-Samac 6
## 113 Olovo 6
## 114 Sekovici 6
## 115 Vares 6
## 116 Vukosavlje 6
## 117 Ostra Luka 5
## 118 Usora 5
## 119 Jablanica 4
## 120 Petrovo 4
## 121 Trnovo 4
## 122 Berkovici 3
## 123 Donji Zabar 3
## 124 Kalinovik 3
## 125 Sapna 3
## 126 Sipovo 3
## 127 Osmaci 2
## 128 Krupa na Uni 1
## 129 Neum 1
## 130 Ravno 1dim(hbs_wage.NEW.long.2a)#imam samo 5163 varijabli## [1] 4654 14hbs_wage.NEW.long.2a %>%
select(naziv_Opcine) %>%
group_by(naziv_Opcine) %>%
summarise (
n = n()
)%>%
arrange (desc(n))%>%
print(n=137) #izlistalo 28 opstina koje imaju manje od 10 obs## # A tibble: 130 × 2
## naziv_Opcine n
## <chr> <int>
## 1 Brcko 267
## 2 Banja Luka 218
## 3 Bijeljina 213
## 4 N.G.Sarajevo 201
## 5 Grad Mostar 176
## 6 Zenica 163
## 7 Travnik 130
## 8 Bihac 87
## 9 Tuzla 82
## 10 Prijedor 81
## 11 Samac 80
## 12 N.Travnik 78
## 13 Sanski Most 76
## 14 Vitez 76
## 15 Doboj 75
## 16 Bugojno 73
## 17 Cazin 68
## 18 Konjic 64
## 19 Zvornik 63
## 20 Visoko 58
## 21 C.Sarajevo 57
## 22 N.Sarajevo 56
## 23 S.Brijeg 56
## 24 Tesanj 56
## 25 Modrica 55
## 26 Novi Grad 55
## 27 Zavidovici 54
## 28 Zivinice 51
## 29 Ilidza 50
## 30 S.G.Sarajevo 50
## 31 Gradacac 45
## 32 Laktasi 45
## 33 Prnjavor 45
## 34 V.Kladusa 44
## 35 D.Vakuf 43
## 36 Lopare 43
## 37 Gracanica 42
## 38 Trebinje 40
## 39 Kakanj 39
## 40 Gorazde 38
## 41 Teslic 38
## 42 Derventa 37
## 43 Kiseljak 35
## 44 G.Vakuf 34
## 45 B.Krupa 33
## 46 Livno 33
## 47 Ljubuski 32
## 48 Lukavac 32
## 49 Busovaca 31
## 50 Foca - RS 31
## 51 Grude 31
## 52 Maglaj 29
## 53 Orasje 29
## 54 Citluk 28
## 55 Brod 27
## 56 Capljina 27
## 57 Gacko 27
## 58 Banovici 25
## 59 Kalesija 25
## 60 Ugljevik 25
## 61 Vogosca 25
## 62 Pale 24
## 63 Kresevo 23
## 64 Gradiska 21
## 65 Drvar 20
## 66 Rogatica 20
## 67 Jajce 19
## 68 Posusje 19
## 69 Srbac 18
## 70 Vlasenica 18
## 71 Breza 17
## 72 I.Ilidza 17
## 73 Ribnik 17
## 74 Rudo 17
## 75 Ilijas 16
## 76 Kljuc 16
## 77 Kupres 16
## 78 Stolac 16
## 79 Cajnice 15
## 80 Visegrad 15
## 81 Fojnica 14
## 82 K.Dubica 14
## 83 Bileca 13
## 84 Bratunac 13
## 85 Buzim 13
## 86 Glamoc 13
## 87 Kladanj 13
## 88 Kotor Varos 13
## 89 Prozor 13
## 90 Knezevo 12
## 91 Nevesinje 12
## 92 Sokolac 12
## 93 Tomislavgrad 12
## 94 B.Petrovac 11
## 95 Celinac 11
## 96 Srebrenik 11
## 97 Doboj-Jug 10
## 98 Mrkonjic 10
## 99 N.Gorazde 10
## 100 Srebrenica 10
## 101 Doboj-Istok 9
## 102 Hadzici 9
## 103 B.Grahovo 7
## 104 Han Pijesak 7
## 105 Kostajnica 7
## 106 Milici 7
## 107 Odzak 7
## 108 Teocak 7
## 109 Zepce 7
## 110 Celic 6
## 111 D-Samac 6
## 112 Pelagicevo 6
## 113 Vares 6
## 114 Vukosavlje 6
## 115 Olovo 5
## 116 Ostra Luka 5
## 117 Sekovici 4
## 118 Trnovo 4
## 119 Usora 4
## 120 Berkovici 3
## 121 Donji Zabar 3
## 122 Jablanica 3
## 123 Kalinovik 3
## 124 Petrovo 3
## 125 Sipovo 3
## 126 Osmaci 2
## 127 Krupa na Uni 1
## 128 Neum 1
## 129 Ravno 1
## 130 Sapna 1unique(hbs_wage$naziv_Opcine) #137 opcina## [1] "B.Grahovo" "B.Krupa" "B.Petrovac" "Banovici" "Banja Luka"
## [6] "Berkovici" "Bihac" "Bijeljina" "Bileca" "Bratunac"
## [11] "Brcko" "Breza" "Brod" "Bugojno" "Busovaca"
## [16] "Buzim" "C.Sarajevo" "Cajnice" "Capljina" "Cazin"
## [21] "Celic" "Celinac" "Citluk" "D-Samac" "D.Vakuf"
## [26] "Derventa" "Doboj" "Doboj-Istok" "Doboj-Jug" "Dobretici"
## [31] "Donji Zabar" "Drvar" "Foca" "Foca - RS" "Fojnica"
## [36] "G.Vakuf" "Gacko" "Glamoc" "Gorazde" "Gracanica"
## [41] "Grad Mostar" "Gradacac" "Gradiska" "Grude" "Hadzici"
## [46] "Han Pijesak" "I.Ilidza" "I.N.Sarajevo" "I.Stari Grad" "Ilidza"
## [51] "Ilijas" "Jablanica" "Jajce" "Jezero" "K.Dubica"
## [56] "Kakanj" "Kalesija" "Kalinovik" "Kiseljak" "Kladanj"
## [61] "Kljuc" "Knezevo" "Konjic" "Kostajnica" "Kotor Varos"
## [66] "Kresevo" "Krupa na Uni" "Kupres" "Laktasi" "Livno"
## [71] "Lopare" "Lukavac" "Ljubinje" "Ljubuski" "Maglaj"
## [76] "Milici" "Modrica" "Mrkonjic" "N.G.Sarajevo" "N.Gorazde"
## [81] "N.Sarajevo" "N.Travnik" "Neum" "Nevesinje" "Novi Grad"
## [86] "Odzak" "Olovo" "Orasje" "Osmaci" "Ostra Luka"
## [91] "Pale" "Pale-FBiH" "Pelagicevo" "Petrovo" "Posusje"
## [96] "Prijedor" "Prnjavor" "Prozor" "Ravno" "Ribnik"
## [101] "Rogatica" "Rudo" "S.Brijeg" "S.G.Sarajevo" "Samac"
## [106] "Sanski Most" "Sapna" "Sekovici" "Sipovo" "Sokolac"
## [111] "Srbac" "Srebrenica" "Srebrenik" "Stolac" "Teocak"
## [116] "Tesanj" "Teslic" "Tomislavgrad" "Travnik" "Trebinje"
## [121] "Trnovo" "Tuzla" "Ugljevik" "Usora" "V.Kladusa"
## [126] "Vares" "Visegrad" "Visoko" "Vitez" "Vlasenica"
## [131] "Vogosca" "Vukosavlje" "Zavidovici" "Zenica" "Zepce"
## [136] "Zivinice" "Zvornik"unique(hbs_wage.NEW.long.2$naziv_Opcine) #130 ## [1] "B.Grahovo" "B.Krupa" "B.Petrovac" "Banovici" "Banja Luka"
## [6] "Berkovici" "Bihac" "Bijeljina" "Bileca" "Bratunac"
## [11] "Brcko" "Breza" "Brod" "Bugojno" "Busovaca"
## [16] "Buzim" "C.Sarajevo" "Cajnice" "Capljina" "Cazin"
## [21] "Celic" "Celinac" "Citluk" "D-Samac" "D.Vakuf"
## [26] "Derventa" "Doboj" "Doboj-Istok" "Doboj-Jug" "Donji Zabar"
## [31] "Drvar" "Foca - RS" "Fojnica" "G.Vakuf" "Gacko"
## [36] "Glamoc" "Gorazde" "Gracanica" "Grad Mostar" "Gradacac"
## [41] "Gradiska" "Grude" "Hadzici" "Han Pijesak" "I.Ilidza"
## [46] "Ilidza" "Ilijas" "Jablanica" "Jajce" "K.Dubica"
## [51] "Kakanj" "Kalesija" "Kalinovik" "Kiseljak" "Kladanj"
## [56] "Kljuc" "Knezevo" "Konjic" "Kostajnica" "Kotor Varos"
## [61] "Kresevo" "Krupa na Uni" "Kupres" "Laktasi" "Livno"
## [66] "Lopare" "Lukavac" "Ljubuski" "Maglaj" "Milici"
## [71] "Modrica" "Mrkonjic" "N.G.Sarajevo" "N.Gorazde" "N.Sarajevo"
## [76] "N.Travnik" "Neum" "Nevesinje" "Novi Grad" "Odzak"
## [81] "Olovo" "Orasje" "Osmaci" "Ostra Luka" "Pale"
## [86] "Pelagicevo" "Petrovo" "Posusje" "Prijedor" "Prnjavor"
## [91] "Prozor" "Ravno" "Ribnik" "Rogatica" "Rudo"
## [96] "S.Brijeg" "S.G.Sarajevo" "Samac" "Sanski Most" "Sapna"
## [101] "Sekovici" "Sipovo" "Sokolac" "Srbac" "Srebrenica"
## [106] "Srebrenik" "Stolac" "Teocak" "Tesanj" "Teslic"
## [111] "Tomislavgrad" "Travnik" "Trebinje" "Trnovo" "Tuzla"
## [116] "Ugljevik" "Usora" "V.Kladusa" "Vares" "Visegrad"
## [121] "Visoko" "Vitez" "Vlasenica" "Vogosca" "Vukosavlje"
## [126] "Zavidovici" "Zenica" "Zepce" "Zivinice" "Zvornik"unique(hbs_wage.NEW.long.2a$naziv_Opcine) #131## [1] "B.Grahovo" "B.Krupa" "B.Petrovac" "Banovici" "Banja Luka"
## [6] "Berkovici" "Bihac" "Bijeljina" "Bileca" "Bratunac"
## [11] "Brcko" "Breza" "Brod" "Bugojno" "Busovaca"
## [16] "Buzim" "C.Sarajevo" "Cajnice" "Capljina" "Cazin"
## [21] "Celic" "Celinac" "Citluk" "D-Samac" "D.Vakuf"
## [26] "Derventa" "Doboj" "Doboj-Istok" "Doboj-Jug" "Donji Zabar"
## [31] "Drvar" "Foca - RS" "Fojnica" "G.Vakuf" "Gacko"
## [36] "Glamoc" "Gorazde" "Gracanica" "Grad Mostar" "Gradacac"
## [41] "Gradiska" "Grude" "Hadzici" "Han Pijesak" "I.Ilidza"
## [46] "Ilidza" "Ilijas" "Jablanica" "Jajce" "K.Dubica"
## [51] "Kakanj" "Kalesija" "Kalinovik" "Kiseljak" "Kladanj"
## [56] "Kljuc" "Knezevo" "Konjic" "Kostajnica" "Kotor Varos"
## [61] "Kresevo" "Krupa na Uni" "Kupres" "Laktasi" "Livno"
## [66] "Lopare" "Lukavac" "Ljubuski" "Maglaj" "Milici"
## [71] "Modrica" "Mrkonjic" "N.G.Sarajevo" "N.Gorazde" "N.Sarajevo"
## [76] "N.Travnik" "Neum" "Nevesinje" "Novi Grad" "Odzak"
## [81] "Olovo" "Orasje" "Osmaci" "Ostra Luka" "Pale"
## [86] "Pelagicevo" "Petrovo" "Posusje" "Prijedor" "Prnjavor"
## [91] "Prozor" "Ravno" "Ribnik" "Rogatica" "Rudo"
## [96] "S.Brijeg" "S.G.Sarajevo" "Samac" "Sanski Most" "Sapna"
## [101] "Sekovici" "Sipovo" "Sokolac" "Srbac" "Srebrenica"
## [106] "Srebrenik" "Stolac" "Teocak" "Tesanj" "Teslic"
## [111] "Tomislavgrad" "Travnik" "Trebinje" "Trnovo" "Tuzla"
## [116] "Ugljevik" "Usora" "V.Kladusa" "Vares" "Visegrad"
## [121] "Visoko" "Vitez" "Vlasenica" "Vogosca" "Vukosavlje"
## [126] "Zavidovici" "Zenica" "Zepce" "Zivinice" "Zvornik"#7 se ne nalaye u hbs_wage.NEW.long.2 jer vec su izbacene dakle 29+7 = 36 opstine manje###share of topli obrok i prevoz u FBiH
summary(hbs_wage.NEW.long.2$QA07_0, na.rm=T)#mean TO = 175## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 8.333 137.500 166.667 175.588 200.000 800.000 3954tapply(hbs_wage.NEW.long.2$QA07_0,hbs_wage.NEW.long.2$Entity, summary)## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 8.333 146.667 170.000 179.769 200.000 800.000 2070
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 30.0 110.0 129.2 144.5 180.0 800.0 1575
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 100 100 100 100 100 100 309summary(hbs_wage.NEW.long.2$QA08_0) #prevoz mean je 76.54## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 4.167 51.000 53.000 76.877 91.667 400.000 4799tapply(hbs_wage.NEW.long.2$QA08_0,hbs_wage.NEW.long.2$Entity, summary)## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 4.417 50.000 53.000 73.638 82.500 300.000 2825
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 4.167 60.000 85.000 98.032 117.500 400.000 1663
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 100 100 100 100 100 100 311#ovdje ima dosta NAova a 0 nema tako da to bas nije u potpunosti tacno. Ima ljudi i u FBiH i u RS koji ne primaju topli ob i prevoz a rade###prebacicu NA u 0 i opet izracunati to i prevoz
hbs_wage.NEW.long.2$QA07_0 <- ifelse (is.na(hbs_wage.NEW.long.2$QA07_0==T),0,paste(hbs_wage.NEW.long.2$QA07_0))
table (is.na (hbs_wage.NEW.long.2$QA07_0))##
## FALSE
## 5244hbs_wage.NEW.long.2$QA07_0 <- as.numeric(hbs_wage.NEW.long.2$QA07_0)
summary(hbs_wage.NEW.long.2$QA07_0)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 43.19 0.00 800.00tapply(hbs_wage.NEW.long.2$QA07_0, hbs_wage.NEW.long.2$Entity, summary)## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 63.88 150.00 800.00
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.00 0.00 0.00 12.26 0.00 800.00
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0000 0.0000 0.0000 0.9615 0.0000 100.0000hbs_wage.NEW.long.2$QA08_0 <- ifelse (is.na(hbs_wage.NEW.long.2$QA08_0==T),0,paste(hbs_wage.NEW.long.2$QA08_0))
table (is.na (hbs_wage.NEW.long.2$QA08_0))##
## FALSE
## 5244hbs_wage.NEW.long.2$QA08_0 <- as.numeric(hbs_wage.NEW.long.2$QA08_0)
summary(hbs_wage.NEW.long.2$QA08_0)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 6.524 0.000 400.000tapply(hbs_wage.NEW.long.2$QA08_0, hbs_wage.NEW.long.2$Entity, summary)## $`1`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 8.852 0.000 300.000
##
## $`2`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.000 0.000 0.000 3.304 0.000 400.000
##
## $`3`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.0000 0.0000 0.0000 0.3205 0.0000 100.0000#FAZA 5 BRISATI OPSTINE SA MANJE OD 10 OBZERVACIJA=hbs_wage.NEW.long.2
Delete low frequency obs (<10) ne mozemo ih stavljati u model
Provjeri da li su izbrisane sve ispod 10
hbs_wage.NEW.long.3 %>%
select(naziv_Opcine) %>%
group_by(naziv_Opcine) %>%
summarise (
n = n()
)%>%
arrange ((n))%>%
print(n=10) #sve su izbrisane## # A tibble: 101 × 2
## naziv_Opcine n
## <chr> <int>
## 1 Celinac 11
## 2 Doboj-Jug 11
## 3 Mrkonjic 11
## 4 Srebrenica 11
## 5 Srebrenik 11
## 6 B.Petrovac 12
## 7 Knezevo 12
## 8 Nevesinje 12
## 9 Sokolac 12
## 10 Bileca 13
## # ℹ 91 more rowsunique(hbs_wage.NEW.long.3$naziv_Opcine)## [1] "B.Krupa" "B.Petrovac" "Banovici" "Banja Luka" "Bihac"
## [6] "Bijeljina" "Bileca" "Bratunac" "Brcko" "Breza"
## [11] "Brod" "Bugojno" "Busovaca" "Buzim" "C.Sarajevo"
## [16] "Cajnice" "Capljina" "Cazin" "Celinac" "Citluk"
## [21] "D.Vakuf" "Derventa" "Doboj" "Doboj-Jug" "Drvar"
## [26] "Foca - RS" "Fojnica" "G.Vakuf" "Gacko" "Glamoc"
## [31] "Gorazde" "Gracanica" "Grad Mostar" "Gradacac" "Gradiska"
## [36] "Grude" "I.Ilidza" "Ilidza" "Ilijas" "Jajce"
## [41] "K.Dubica" "Kakanj" "Kalesija" "Kiseljak" "Kladanj"
## [46] "Kljuc" "Knezevo" "Konjic" "Kotor Varos" "Kresevo"
## [51] "Kupres" "Laktasi" "Livno" "Lopare" "Lukavac"
## [56] "Ljubuski" "Maglaj" "Modrica" "Mrkonjic" "N.G.Sarajevo"
## [61] "N.Gorazde" "N.Sarajevo" "N.Travnik" "Nevesinje" "Novi Grad"
## [66] "Orasje" "Pale" "Posusje" "Prijedor" "Prnjavor"
## [71] "Prozor" "Ribnik" "Rogatica" "Rudo" "S.Brijeg"
## [76] "S.G.Sarajevo" "Samac" "Sanski Most" "Sokolac" "Srbac"
## [81] "Srebrenica" "Srebrenik" "Stolac" "Tesanj" "Teslic"
## [86] "Tomislavgrad" "Travnik" "Trebinje" "Tuzla" "Ugljevik"
## [91] "V.Kladusa" "Visegrad" "Visoko" "Vitez" "Vlasenica"
## [96] "Vogosca" "Zavidovici" "Zenica" "Zepce" "Zivinice"
## [101] "Zvornik"#101 opcina od ukupno 137 koje su u HBS, dakle u wage- imam manje 36 opstina#FAZA 5 BRISATI OPSTINE SA MANJE OD 10 OBZERVACIJA=hbs_wage.NEW.long.2a
Delete low frequency obs (<10) ne mozemo ih stavljati u model
Provjeri da li su izbrisane sve ispod 10
hbs_wage.NEW.long.3a %>%
select(naziv_Opcine) %>%
group_by(naziv_Opcine) %>%
summarise (
n = n()
)%>%
arrange ((n))%>%
print(n=10) #sve su izbrisane## # A tibble: 100 × 2
## naziv_Opcine n
## <chr> <int>
## 1 Doboj-Jug 10
## 2 Mrkonjic 10
## 3 N.Gorazde 10
## 4 Srebrenica 10
## 5 B.Petrovac 11
## 6 Celinac 11
## 7 Srebrenik 11
## 8 Knezevo 12
## 9 Nevesinje 12
## 10 Sokolac 12
## # ℹ 90 more rowsunique(hbs_wage.NEW.long.3a$naziv_Opcine)## [1] "B.Krupa" "B.Petrovac" "Banovici" "Banja Luka" "Bihac"
## [6] "Bijeljina" "Bileca" "Bratunac" "Brcko" "Breza"
## [11] "Brod" "Bugojno" "Busovaca" "Buzim" "C.Sarajevo"
## [16] "Cajnice" "Capljina" "Cazin" "Celinac" "Citluk"
## [21] "D.Vakuf" "Derventa" "Doboj" "Doboj-Jug" "Drvar"
## [26] "Foca - RS" "Fojnica" "G.Vakuf" "Gacko" "Glamoc"
## [31] "Gorazde" "Gracanica" "Grad Mostar" "Gradacac" "Gradiska"
## [36] "Grude" "I.Ilidza" "Ilidza" "Ilijas" "Jajce"
## [41] "K.Dubica" "Kakanj" "Kalesija" "Kiseljak" "Kladanj"
## [46] "Kljuc" "Knezevo" "Konjic" "Kotor Varos" "Kresevo"
## [51] "Kupres" "Laktasi" "Livno" "Lopare" "Lukavac"
## [56] "Ljubuski" "Maglaj" "Modrica" "Mrkonjic" "N.G.Sarajevo"
## [61] "N.Gorazde" "N.Sarajevo" "N.Travnik" "Nevesinje" "Novi Grad"
## [66] "Orasje" "Pale" "Posusje" "Prijedor" "Prnjavor"
## [71] "Prozor" "Ribnik" "Rogatica" "Rudo" "S.Brijeg"
## [76] "S.G.Sarajevo" "Samac" "Sanski Most" "Sokolac" "Srbac"
## [81] "Srebrenica" "Srebrenik" "Stolac" "Tesanj" "Teslic"
## [86] "Tomislavgrad" "Travnik" "Trebinje" "Tuzla" "Ugljevik"
## [91] "V.Kladusa" "Visegrad" "Visoko" "Vitez" "Vlasenica"
## [96] "Vogosca" "Zavidovici" "Zenica" "Zivinice" "Zvornik"#103 opcina od ukupno 137 koje su u HBS, dakle u wage- imam manje 34 opstinaFAZA 6 ISPITAJ VARIJABLE
Rename varijables
library(dplyr)
WI.df <- hbs_wage.NEW.long.3
names(WI.df)## [1] "SifraDom" "Entity" "naziv_Opcine" "RurUrb" "Relat"
## [6] "Age_" "Relat_" "Activ_" "Educ_" "Sex_"
## [11] "QA_0" "QAnotall_0" "QA07_0" "QA08_0"WI.df <- rename(WI.df, c(wages = QA_0, Age = Age_, Active = Activ_, Educ = Educ_, Munic = naziv_Opcine))library(dplyr)
WI.dfa <- hbs_wage.NEW.long.3a
names(WI.dfa)## [1] "SifraDom" "Entity" "naziv_Opcine" "RurUrb" "Relat"
## [6] "Age_" "Relat_" "Activ_" "Educ_" "Sex_"
## [11] "QA_0" "QAnotall_0" "QA07_0" "QA08_0"WI.dfa <- rename(WI.dfa, c(wages = QAnotall_0, Age = Age_, Active = Activ_, Educ = Educ_, Munic = naziv_Opcine))CHECK NA
table(is.na(WI.df$wages), WI.df$Munic)##
## B.Krupa B.Petrovac Banovici Banja Luka Bihac Bijeljina Bileca Bratunac
## FALSE 36 12 26 236 97 245 13 13
##
## Brcko Breza Brod Bugojno Busovaca Buzim C.Sarajevo Cajnice Capljina
## FALSE 312 17 38 93 33 16 68 16 28
##
## Cazin Celinac Citluk D.Vakuf Derventa Doboj Doboj-Jug Drvar Foca - RS
## FALSE 72 11 30 52 44 78 11 20 31
##
## Fojnica G.Vakuf Gacko Glamoc Gorazde Gracanica Grad Mostar Gradacac
## FALSE 14 44 28 13 40 50 190 51
##
## Gradiska Grude I.Ilidza Ilidza Ilijas Jajce K.Dubica Kakanj Kalesija
## FALSE 22 32 18 61 19 21 15 42 42
##
## Kiseljak Kladanj Kljuc Knezevo Konjic Kotor Varos Kresevo Kupres
## FALSE 38 15 18 12 69 13 24 16
##
## Laktasi Livno Lopare Lukavac Ljubuski Maglaj Modrica Mrkonjic
## FALSE 45 36 48 39 42 45 60 11
##
## N.G.Sarajevo N.Gorazde N.Sarajevo N.Travnik Nevesinje Novi Grad Orasje
## FALSE 217 16 67 80 12 61 33
##
## Pale Posusje Prijedor Prnjavor Prozor Ribnik Rogatica Rudo S.Brijeg
## FALSE 24 21 84 55 17 19 22 22 59
##
## S.G.Sarajevo Samac Sanski Most Sokolac Srbac Srebrenica Srebrenik
## FALSE 51 85 90 12 19 11 11
##
## Stolac Tesanj Teslic Tomislavgrad Travnik Trebinje Tuzla Ugljevik
## FALSE 17 74 47 21 136 41 104 25
##
## V.Kladusa Visegrad Visoko Vitez Vlasenica Vogosca Zavidovici Zenica
## FALSE 62 18 60 81 18 26 71 174
##
## Zepce Zivinice Zvornik
## FALSE 14 70 67table(is.na(WI.dfa$wages), WI.dfa$Munic)##
## B.Krupa B.Petrovac Banovici Banja Luka Bihac Bijeljina Bileca Bratunac
## FALSE 33 11 25 218 87 213 13 13
##
## Brcko Breza Brod Bugojno Busovaca Buzim C.Sarajevo Cajnice Capljina
## FALSE 267 17 27 73 31 13 57 15 27
##
## Cazin Celinac Citluk D.Vakuf Derventa Doboj Doboj-Jug Drvar Foca - RS
## FALSE 68 11 28 43 37 75 10 20 31
##
## Fojnica G.Vakuf Gacko Glamoc Gorazde Gracanica Grad Mostar Gradacac
## FALSE 14 34 27 13 38 42 176 45
##
## Gradiska Grude I.Ilidza Ilidza Ilijas Jajce K.Dubica Kakanj Kalesija
## FALSE 21 31 17 50 16 19 14 39 25
##
## Kiseljak Kladanj Kljuc Knezevo Konjic Kotor Varos Kresevo Kupres
## FALSE 35 13 16 12 64 13 23 16
##
## Laktasi Livno Lopare Lukavac Ljubuski Maglaj Modrica Mrkonjic
## FALSE 45 33 43 32 32 29 55 10
##
## N.G.Sarajevo N.Gorazde N.Sarajevo N.Travnik Nevesinje Novi Grad Orasje
## FALSE 201 10 56 78 12 55 29
##
## Pale Posusje Prijedor Prnjavor Prozor Ribnik Rogatica Rudo S.Brijeg
## FALSE 24 19 81 45 13 17 20 17 56
##
## S.G.Sarajevo Samac Sanski Most Sokolac Srbac Srebrenica Srebrenik
## FALSE 50 80 76 12 18 10 11
##
## Stolac Tesanj Teslic Tomislavgrad Travnik Trebinje Tuzla Ugljevik
## FALSE 16 56 38 12 130 40 82 25
##
## V.Kladusa Visegrad Visoko Vitez Vlasenica Vogosca Zavidovici Zenica
## FALSE 44 15 58 76 18 25 54 163
##
## Zivinice Zvornik
## FALSE 51 63wages
plot(density(WI.df$wages))
summary(WI.df$wages)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.083 400.000 630.000 748.772 981.667 18350.000sd (WI.df$wages)## [1] 640.5685Age
plot(density(WI.df$Age))
summary(WI.df$Age)## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 17.00 33.00 43.00 42.56 52.00 79.00sd (WI.df$Age)## [1] 11.6293Active
Status tekuce zaposlenosti: 1. puno radno vrijeme 2. pola radnog vremena Od 3-9 je izbaceno 3. nezaposlen 4 - trazi prvo zaposlenje 5 - domacica 6 - student 7 - nesposoban za rad 8 - penzioner 9 - ostalo
table(WI.df$Active)##
## 1 2
## 4515 580Prema instrukciji Ismira ovu cu samo preimenovati u Full time i Part time i stavicu Full. time da mi je referentna
Educ
Steceno obrazovanje: 1 - bez skole; 2- nepotpuna osnovna 8.g.skola 3 - nepotpuna osnovna 9.g. skola; 4-osnovna skola, 5-srednja skola, 6-specijalizacia poslije s.s; 7 - Visa skola; 8 - Fakultet; 9 - Master; 10 - Doktorat
table(WI.df$Educ)##
## 1 2 3 4 5 6 7 8 9 10
## 39 60 4 658 3370 30 184 666 76 8Od ove varijable napravicu bez osnovne skole (1,2,3), osnovna skola (4), srednja skola (5,6), fakultet (7,8,9,10)
WI.df <- WI.df %>% mutate(Full.time = case_when(Active == 1 ~"full time", Active == 2 ~"part time"),
Educ.1 = case_when(Educ %in% c(1,2,3) ~ "no elementary", Educ == 4 ~ "elementary", Educ %in% c(5,6) ~ "high schl", Educ %in% c(7,8,9,10) ~ "faculty")) WI.dfa <- WI.dfa %>% mutate(Full.time = case_when(Active == 1 ~"full time", Active == 2 ~"part time"),
Educ.1 = case_when(Educ %in% c(1,2,3) ~ "no elementary", Educ == 4 ~ "elementary", Educ %in% c(5,6) ~ "high schl", Educ %in% c(7,8,9,10) ~ "faculty")) Check if it is ok
Active
table (WI.df$Full.time)##
## full time part time
## 4515 580table (WI.dfa$Full.time)##
## full time part time
## 4078 433round (prop.table(table(WI.df$Activ))*100,1)## Warning: Unknown or uninitialised column: `Activ`.## numeric(0)round (prop.table(table(WI.dfa$Active))*100,1)##
## 1 2
## 90.4 9.6##Educ
table (WI.df$Educ)##
## 1 2 3 4 5 6 7 8 9 10
## 39 60 4 658 3370 30 184 666 76 8table (WI.df$Educ.1)##
## elementary faculty high schl no elementary
## 658 934 3400 103round (prop.table(table(WI.df$Educ.1))*100,1)##
## elementary faculty high schl no elementary
## 12.9 18.3 66.7 2.0table (WI.dfa$Educ)##
## 1 2 3 4 5 6 7 8 9 10
## 32 47 4 562 2989 24 169 608 69 7table (WI.dfa$Educ.1)##
## elementary faculty high schl no elementary
## 562 853 3013 83round (prop.table(table(WI.dfa$Educ.1))*100,1)##
## elementary faculty high schl no elementary
## 12.5 18.9 66.8 1.8Relevel
Full.time
Sada treba relevel varijablu da full time bude
referentna vrijednost
WI.dfa$Full.time <- fct_relevel (WI.dfa$Full.time, "full time", "part time")
table(WI.dfa$Full.time)##
## full time part time
## 4078 433str(WI.dfa$Full.time)## Factor w/ 2 levels "full time","part time": 1 1 1 1 1 1 1 1 1 1 ...Educ
Sada treba relevel varijablu da no elementary bude
referentna vrijednost
Sex
round (prop.table(table(WI.df$Sex_))*100,1)##
## 1 2
## 65.8 34.2table(WI.df$Munic)##
## B.Krupa B.Petrovac Banovici Banja Luka Bihac Bijeljina
## 36 12 26 236 97 245
## Bileca Bratunac Brcko Breza Brod Bugojno
## 13 13 312 17 38 93
## Busovaca Buzim C.Sarajevo Cajnice Capljina Cazin
## 33 16 68 16 28 72
## Celinac Citluk D.Vakuf Derventa Doboj Doboj-Jug
## 11 30 52 44 78 11
## Drvar Foca - RS Fojnica G.Vakuf Gacko Glamoc
## 20 31 14 44 28 13
## Gorazde Gracanica Grad Mostar Gradacac Gradiska Grude
## 40 50 190 51 22 32
## I.Ilidza Ilidza Ilijas Jajce K.Dubica Kakanj
## 18 61 19 21 15 42
## Kalesija Kiseljak Kladanj Kljuc Knezevo Konjic
## 42 38 15 18 12 69
## Kotor Varos Kresevo Kupres Laktasi Livno Lopare
## 13 24 16 45 36 48
## Lukavac Ljubuski Maglaj Modrica Mrkonjic N.G.Sarajevo
## 39 42 45 60 11 217
## N.Gorazde N.Sarajevo N.Travnik Nevesinje Novi Grad Orasje
## 16 67 80 12 61 33
## Pale Posusje Prijedor Prnjavor Prozor Ribnik
## 24 21 84 55 17 19
## Rogatica Rudo S.Brijeg S.G.Sarajevo Samac Sanski Most
## 22 22 59 51 85 90
## Sokolac Srbac Srebrenica Srebrenik Stolac Tesanj
## 12 19 11 11 17 74
## Teslic Tomislavgrad Travnik Trebinje Tuzla Ugljevik
## 47 21 136 41 104 25
## V.Kladusa Visegrad Visoko Vitez Vlasenica Vogosca
## 62 18 60 81 18 26
## Zavidovici Zenica Zepce Zivinice Zvornik
## 71 174 14 70 67#FAZA 7 BUILD UP A MODEL
Mincer earnings function log(wages)
#Isti model koji je winner kad uzmem sve opcije za Active
library(texreg)## Version: 1.39.3
## Date: 2023-11-09
## Author: Philip Leifeld (University of Essex)
##
## Consider submitting praise using the praise or praise_interactive functions.
## Please cite the JSS article in your publications -- see citation("texreg").##
## Attaching package: 'texreg'## The following object is masked from 'package:tidyr':
##
## extractWI.3 <- lm (log(wages)~ Age + as.numeric(Age^2)+ factor (Sex_) + Educ.1 + factor(Full.time) + factor (Munic) , data = WI.df) # No rur. urb THE WINNING MODEL
WI.3a <- lm (log(wages)~ Age + as.numeric(Age^2)+ factor (Sex_) + Educ.1 + factor(Full.time) + factor (Munic) , data = WI.dfa) # No rur. urb THE WINNING MODEL
texreg::screenreg(list(WI.3, WI.3a))##
## ====================================================
## Model 1 Model 2
## ----------------------------------------------------
## (Intercept) 4.37 *** 4.88 ***
## (0.16) (0.14)
## Age 0.06 *** 0.04 ***
## (0.01) (0.00)
## as.numeric(Age^2) -0.00 *** -0.00 ***
## (0.00) (0.00)
## factor(Sex_)2 -0.20 *** -0.15 ***
## (0.02) (0.02)
## Educ.1elementary 0.33 *** 0.17 *
## (0.08) (0.07)
## Educ.1high schl 0.61 *** 0.48 ***
## (0.08) (0.07)
## Educ.1faculty 1.08 *** 0.93 ***
## (0.08) (0.07)
## factor(Full.time)part time -0.85 *** -0.74 ***
## (0.03) (0.03)
## factor(Munic)B.Petrovac 0.07 -0.18
## (0.23) (0.19)
## factor(Munic)Banovici 0.33 0.30 *
## (0.18) (0.14)
## factor(Munic)Banja Luka 0.07 0.12
## (0.12) (0.10)
## factor(Munic)Bihac -0.04 -0.09
## (0.13) (0.11)
## factor(Munic)Bijeljina -0.00 -0.11
## (0.12) (0.10)
## factor(Munic)Bileca 0.15 0.05
## (0.22) (0.18)
## factor(Munic)Bratunac 0.39 0.35 *
## (0.22) (0.18)
## factor(Munic)Brcko 0.25 * 0.12
## (0.12) (0.10)
## factor(Munic)Breza 0.30 0.23
## (0.20) (0.16)
## factor(Munic)Brod -0.71 *** -0.26
## (0.16) (0.14)
## factor(Munic)Bugojno 0.06 0.04
## (0.14) (0.11)
## factor(Munic)Busovaca 0.19 0.06
## (0.17) (0.14)
## factor(Munic)Buzim -0.74 *** -0.67 ***
## (0.21) (0.18)
## factor(Munic)C.Sarajevo 0.44 ** 0.34 **
## (0.14) (0.12)
## factor(Munic)Cajnice -0.21 -0.18
## (0.21) (0.17)
## factor(Munic)Capljina 0.54 ** 0.48 ***
## (0.17) (0.14)
## factor(Munic)Cazin -0.11 -0.13
## (0.14) (0.11)
## factor(Munic)Celinac -0.02 -0.12
## (0.24) (0.19)
## factor(Munic)Citluk 0.32 0.29 *
## (0.17) (0.14)
## factor(Munic)D.Vakuf -0.01 -0.11
## (0.15) (0.13)
## factor(Munic)Derventa -0.47 ** -0.47 ***
## (0.15) (0.13)
## factor(Munic)Doboj 0.24 0.15
## (0.14) (0.11)
## factor(Munic)Doboj-Jug -0.06 -0.11
## (0.24) (0.19)
## factor(Munic)Drvar 0.15 0.03
## (0.19) (0.15)
## factor(Munic)Foca - RS 0.07 -0.03
## (0.17) (0.14)
## factor(Munic)Fojnica 0.38 0.25
## (0.22) (0.17)
## factor(Munic)G.Vakuf -0.01 -0.10
## (0.15) (0.13)
## factor(Munic)Gacko 0.31 0.25
## (0.17) (0.14)
## factor(Munic)Glamoc 0.34 0.24
## (0.22) (0.18)
## factor(Munic)Gorazde 0.31 0.24
## (0.16) (0.13)
## factor(Munic)Gracanica -0.03 -0.13
## (0.15) (0.13)
## factor(Munic)Grad Mostar 0.45 *** 0.33 **
## (0.13) (0.10)
## factor(Munic)Gradacac 0.08 0.07
## (0.15) (0.12)
## factor(Munic)Gradiska -0.15 -0.05
## (0.19) (0.15)
## factor(Munic)Grude 0.49 ** 0.40 **
## (0.17) (0.14)
## factor(Munic)I.Ilidza 0.09 -0.01
## (0.20) (0.16)
## factor(Munic)Ilidza 0.66 *** 0.51 ***
## (0.14) (0.12)
## factor(Munic)Ilijas 0.17 0.18
## (0.19) (0.16)
## factor(Munic)Jajce 0.46 * 0.28
## (0.19) (0.16)
## factor(Munic)K.Dubica -0.07 -0.14
## (0.21) (0.17)
## factor(Munic)Kakanj 0.13 0.08
## (0.16) (0.13)
## factor(Munic)Kalesija -0.23 -0.13
## (0.16) (0.14)
## factor(Munic)Kiseljak 0.36 * 0.25
## (0.16) (0.13)
## factor(Munic)Kladanj 0.16 0.13
## (0.21) (0.18)
## factor(Munic)Kljuc 0.08 -0.01
## (0.20) (0.16)
## factor(Munic)Knezevo 0.03 -0.03
## (0.23) (0.18)
## factor(Munic)Konjic 0.35 * 0.23 *
## (0.14) (0.12)
## factor(Munic)Kotor Varos 0.21 0.14
## (0.22) (0.18)
## factor(Munic)Kresevo 0.37 * 0.37 *
## (0.18) (0.15)
## factor(Munic)Kupres 0.29 0.20
## (0.21) (0.16)
## factor(Munic)Laktasi 0.10 0.02
## (0.15) (0.12)
## factor(Munic)Livno 0.53 *** 0.41 **
## (0.16) (0.13)
## factor(Munic)Lopare 0.10 -0.01
## (0.15) (0.13)
## factor(Munic)Lukavac 0.28 0.12
## (0.16) (0.13)
## factor(Munic)Ljubuski 0.66 *** 0.42 **
## (0.16) (0.13)
## factor(Munic)Maglaj 0.16 0.07
## (0.15) (0.14)
## factor(Munic)Modrica 0.06 -0.06
## (0.14) (0.12)
## factor(Munic)Mrkonjic 0.36 0.27
## (0.24) (0.20)
## factor(Munic)N.G.Sarajevo 0.44 *** 0.35 ***
## (0.12) (0.10)
## factor(Munic)N.Gorazde -0.28 -0.41 *
## (0.21) (0.19)
## factor(Munic)N.Sarajevo 0.57 *** 0.43 ***
## (0.14) (0.12)
## factor(Munic)N.Travnik -0.03 -0.17
## (0.14) (0.11)
## factor(Munic)Nevesinje 0.02 -0.03
## (0.23) (0.18)
## factor(Munic)Novi Grad 0.05 -0.05
## (0.14) (0.12)
## factor(Munic)Orasje 0.46 ** 0.34 *
## (0.17) (0.14)
## factor(Munic)Pale 0.14 0.07
## (0.18) (0.15)
## factor(Munic)Posusje 0.56 ** 0.41 **
## (0.19) (0.16)
## factor(Munic)Prijedor 0.21 0.14
## (0.14) (0.11)
## factor(Munic)Prnjavor 0.16 0.02
## (0.15) (0.12)
## factor(Munic)Prozor 0.35 0.22
## (0.20) (0.18)
## factor(Munic)Ribnik -0.15 -0.02
## (0.19) (0.16)
## factor(Munic)Rogatica 0.23 0.18
## (0.19) (0.15)
## factor(Munic)Rudo -0.07 -0.16
## (0.19) (0.16)
## factor(Munic)S.Brijeg 0.70 *** 0.61 ***
## (0.15) (0.12)
## factor(Munic)S.G.Sarajevo 0.46 ** 0.36 **
## (0.15) (0.12)
## factor(Munic)Samac 0.13 0.02
## (0.14) (0.11)
## factor(Munic)Sanski Most 0.27 * -0.01
## (0.14) (0.11)
## factor(Munic)Sokolac -0.11 -0.16
## (0.23) (0.18)
## factor(Munic)Srbac 0.13 0.06
## (0.19) (0.16)
## factor(Munic)Srebrenica 0.31 0.33
## (0.24) (0.19)
## factor(Munic)Srebrenik 0.21 0.16
## (0.24) (0.19)
## factor(Munic)Stolac 0.45 * 0.31
## (0.20) (0.16)
## factor(Munic)Tesanj 0.07 0.00
## (0.14) (0.12)
## factor(Munic)Teslic 0.06 -0.04
## (0.15) (0.13)
## factor(Munic)Tomislavgrad 0.72 *** 0.17
## (0.19) (0.18)
## factor(Munic)Travnik 0.40 ** 0.33 **
## (0.13) (0.11)
## factor(Munic)Trebinje 0.13 0.03
## (0.16) (0.13)
## factor(Munic)Tuzla 0.01 0.17
## (0.13) (0.11)
## factor(Munic)Ugljevik 0.31 0.17
## (0.18) (0.14)
## factor(Munic)V.Kladusa 0.03 -0.25 *
## (0.14) (0.12)
## factor(Munic)Visegrad 0.05 -0.15
## (0.20) (0.17)
## factor(Munic)Visoko 0.22 0.13
## (0.14) (0.12)
## factor(Munic)Vitez 0.04 -0.07
## (0.14) (0.11)
## factor(Munic)Vlasenica 0.18 0.12
## (0.20) (0.16)
## factor(Munic)Vogosca 0.20 0.14
## (0.18) (0.14)
## factor(Munic)Zavidovici 0.11 -0.04
## (0.14) (0.12)
## factor(Munic)Zenica 0.18 0.10
## (0.13) (0.10)
## factor(Munic)Zepce 0.17
## (0.22)
## factor(Munic)Zivinice -0.19 -0.05
## (0.14) (0.12)
## factor(Munic)Zvornik -0.24 -0.33 **
## (0.14) (0.12)
## ----------------------------------------------------
## R^2 0.36 0.40
## Adj. R^2 0.35 0.39
## Num. obs. 5095 4511
## ====================================================
## *** p < 0.001; ** p < 0.01; * p < 0.05screenreg(WI.3)##
## =======================================
## Model 1
## ---------------------------------------
## (Intercept) 4.37 ***
## (0.16)
## Age 0.06 ***
## (0.01)
## as.numeric(Age^2) -0.00 ***
## (0.00)
## factor(Sex_)2 -0.20 ***
## (0.02)
## Educ.1elementary 0.33 ***
## (0.08)
## Educ.1high schl 0.61 ***
## (0.08)
## Educ.1faculty 1.08 ***
## (0.08)
## factor(Full.time)part time -0.85 ***
## (0.03)
## factor(Munic)B.Petrovac 0.07
## (0.23)
## factor(Munic)Banovici 0.33
## (0.18)
## factor(Munic)Banja Luka 0.07
## (0.12)
## factor(Munic)Bihac -0.04
## (0.13)
## factor(Munic)Bijeljina -0.00
## (0.12)
## factor(Munic)Bileca 0.15
## (0.22)
## factor(Munic)Bratunac 0.39
## (0.22)
## factor(Munic)Brcko 0.25 *
## (0.12)
## factor(Munic)Breza 0.30
## (0.20)
## factor(Munic)Brod -0.71 ***
## (0.16)
## factor(Munic)Bugojno 0.06
## (0.14)
## factor(Munic)Busovaca 0.19
## (0.17)
## factor(Munic)Buzim -0.74 ***
## (0.21)
## factor(Munic)C.Sarajevo 0.44 **
## (0.14)
## factor(Munic)Cajnice -0.21
## (0.21)
## factor(Munic)Capljina 0.54 **
## (0.17)
## factor(Munic)Cazin -0.11
## (0.14)
## factor(Munic)Celinac -0.02
## (0.24)
## factor(Munic)Citluk 0.32
## (0.17)
## factor(Munic)D.Vakuf -0.01
## (0.15)
## factor(Munic)Derventa -0.47 **
## (0.15)
## factor(Munic)Doboj 0.24
## (0.14)
## factor(Munic)Doboj-Jug -0.06
## (0.24)
## factor(Munic)Drvar 0.15
## (0.19)
## factor(Munic)Foca - RS 0.07
## (0.17)
## factor(Munic)Fojnica 0.38
## (0.22)
## factor(Munic)G.Vakuf -0.01
## (0.15)
## factor(Munic)Gacko 0.31
## (0.17)
## factor(Munic)Glamoc 0.34
## (0.22)
## factor(Munic)Gorazde 0.31
## (0.16)
## factor(Munic)Gracanica -0.03
## (0.15)
## factor(Munic)Grad Mostar 0.45 ***
## (0.13)
## factor(Munic)Gradacac 0.08
## (0.15)
## factor(Munic)Gradiska -0.15
## (0.19)
## factor(Munic)Grude 0.49 **
## (0.17)
## factor(Munic)I.Ilidza 0.09
## (0.20)
## factor(Munic)Ilidza 0.66 ***
## (0.14)
## factor(Munic)Ilijas 0.17
## (0.19)
## factor(Munic)Jajce 0.46 *
## (0.19)
## factor(Munic)K.Dubica -0.07
## (0.21)
## factor(Munic)Kakanj 0.13
## (0.16)
## factor(Munic)Kalesija -0.23
## (0.16)
## factor(Munic)Kiseljak 0.36 *
## (0.16)
## factor(Munic)Kladanj 0.16
## (0.21)
## factor(Munic)Kljuc 0.08
## (0.20)
## factor(Munic)Knezevo 0.03
## (0.23)
## factor(Munic)Konjic 0.35 *
## (0.14)
## factor(Munic)Kotor Varos 0.21
## (0.22)
## factor(Munic)Kresevo 0.37 *
## (0.18)
## factor(Munic)Kupres 0.29
## (0.21)
## factor(Munic)Laktasi 0.10
## (0.15)
## factor(Munic)Livno 0.53 ***
## (0.16)
## factor(Munic)Lopare 0.10
## (0.15)
## factor(Munic)Lukavac 0.28
## (0.16)
## factor(Munic)Ljubuski 0.66 ***
## (0.16)
## factor(Munic)Maglaj 0.16
## (0.15)
## factor(Munic)Modrica 0.06
## (0.14)
## factor(Munic)Mrkonjic 0.36
## (0.24)
## factor(Munic)N.G.Sarajevo 0.44 ***
## (0.12)
## factor(Munic)N.Gorazde -0.28
## (0.21)
## factor(Munic)N.Sarajevo 0.57 ***
## (0.14)
## factor(Munic)N.Travnik -0.03
## (0.14)
## factor(Munic)Nevesinje 0.02
## (0.23)
## factor(Munic)Novi Grad 0.05
## (0.14)
## factor(Munic)Orasje 0.46 **
## (0.17)
## factor(Munic)Pale 0.14
## (0.18)
## factor(Munic)Posusje 0.56 **
## (0.19)
## factor(Munic)Prijedor 0.21
## (0.14)
## factor(Munic)Prnjavor 0.16
## (0.15)
## factor(Munic)Prozor 0.35
## (0.20)
## factor(Munic)Ribnik -0.15
## (0.19)
## factor(Munic)Rogatica 0.23
## (0.19)
## factor(Munic)Rudo -0.07
## (0.19)
## factor(Munic)S.Brijeg 0.70 ***
## (0.15)
## factor(Munic)S.G.Sarajevo 0.46 **
## (0.15)
## factor(Munic)Samac 0.13
## (0.14)
## factor(Munic)Sanski Most 0.27 *
## (0.14)
## factor(Munic)Sokolac -0.11
## (0.23)
## factor(Munic)Srbac 0.13
## (0.19)
## factor(Munic)Srebrenica 0.31
## (0.24)
## factor(Munic)Srebrenik 0.21
## (0.24)
## factor(Munic)Stolac 0.45 *
## (0.20)
## factor(Munic)Tesanj 0.07
## (0.14)
## factor(Munic)Teslic 0.06
## (0.15)
## factor(Munic)Tomislavgrad 0.72 ***
## (0.19)
## factor(Munic)Travnik 0.40 **
## (0.13)
## factor(Munic)Trebinje 0.13
## (0.16)
## factor(Munic)Tuzla 0.01
## (0.13)
## factor(Munic)Ugljevik 0.31
## (0.18)
## factor(Munic)V.Kladusa 0.03
## (0.14)
## factor(Munic)Visegrad 0.05
## (0.20)
## factor(Munic)Visoko 0.22
## (0.14)
## factor(Munic)Vitez 0.04
## (0.14)
## factor(Munic)Vlasenica 0.18
## (0.20)
## factor(Munic)Vogosca 0.20
## (0.18)
## factor(Munic)Zavidovici 0.11
## (0.14)
## factor(Munic)Zenica 0.18
## (0.13)
## factor(Munic)Zepce 0.17
## (0.22)
## factor(Munic)Zivinice -0.19
## (0.14)
## factor(Munic)Zvornik -0.24
## (0.14)
## ---------------------------------------
## R^2 0.36
## Adj. R^2 0.35
## Num. obs. 5095
## =======================================
## *** p < 0.001; ** p < 0.01; * p < 0.05summary(WI.3)##
## Call:
## lm(formula = log(wages) ~ Age + as.numeric(Age^2) + factor(Sex_) +
## Educ.1 + factor(Full.time) + factor(Munic), data = WI.df)
##
## Residuals:
## Min 1Q Median 3Q Max
## -8.9246 -0.2729 0.0715 0.3721 3.0283
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 4.366e+00 1.644e-01 26.562 < 2e-16 ***
## Age 5.626e-02 5.706e-03 9.859 < 2e-16 ***
## as.numeric(Age^2) -5.562e-04 6.711e-05 -8.288 < 2e-16 ***
## factor(Sex_)2 -1.977e-01 2.085e-02 -9.482 < 2e-16 ***
## Educ.1elementary 3.267e-01 7.768e-02 4.205 2.66e-05 ***
## Educ.1high schl 6.099e-01 7.564e-02 8.063 9.21e-16 ***
## Educ.1faculty 1.078e+00 7.843e-02 13.746 < 2e-16 ***
## factor(Full.time)part time -8.548e-01 3.331e-02 -25.659 < 2e-16 ***
## factor(Munic)B.Petrovac 7.057e-02 2.282e-01 0.309 0.757169
## factor(Munic)Banovici 3.309e-01 1.765e-01 1.875 0.060910 .
## factor(Munic)Banja Luka 7.302e-02 1.234e-01 0.592 0.554113
## factor(Munic)Bihac -3.746e-02 1.341e-01 -0.279 0.780044
## factor(Munic)Bijeljina -3.139e-03 1.222e-01 -0.026 0.979509
## factor(Munic)Bileca 1.464e-01 2.217e-01 0.661 0.508835
## factor(Munic)Bratunac 3.949e-01 2.213e-01 1.784 0.074437 .
## factor(Munic)Brcko 2.470e-01 1.209e-01 2.042 0.041156 *
## factor(Munic)Breza 3.049e-01 2.017e-01 1.511 0.130793
## factor(Munic)Brod -7.106e-01 1.592e-01 -4.463 8.26e-06 ***
## factor(Munic)Bugojno 6.004e-02 1.353e-01 0.444 0.657183
## factor(Munic)Busovaca 1.881e-01 1.652e-01 1.138 0.255098
## factor(Munic)Buzim -7.420e-01 2.060e-01 -3.602 0.000318 ***
## factor(Munic)C.Sarajevo 4.378e-01 1.424e-01 3.074 0.002120 **
## factor(Munic)Cajnice -2.111e-01 2.061e-01 -1.024 0.305726
## factor(Munic)Capljina 5.351e-01 1.723e-01 3.106 0.001908 **
## factor(Munic)Cazin -1.105e-01 1.395e-01 -0.792 0.428332
## factor(Munic)Celinac -1.803e-02 2.358e-01 -0.076 0.939052
## factor(Munic)Citluk 3.186e-01 1.696e-01 1.878 0.060397 .
## factor(Munic)D.Vakuf -1.033e-02 1.496e-01 -0.069 0.944943
## factor(Munic)Derventa -4.695e-01 1.541e-01 -3.047 0.002323 **
## factor(Munic)Doboj 2.389e-01 1.385e-01 1.725 0.084510 .
## factor(Munic)Doboj-Jug -6.272e-02 2.358e-01 -0.266 0.790286
## factor(Munic)Drvar 1.469e-01 1.914e-01 0.768 0.442688
## factor(Munic)Foca - RS 7.452e-02 1.682e-01 0.443 0.657720
## factor(Munic)Fojnica 3.810e-01 2.156e-01 1.767 0.077315 .
## factor(Munic)G.Vakuf -1.106e-02 1.547e-01 -0.071 0.943005
## factor(Munic)Gacko 3.146e-01 1.728e-01 1.821 0.068673 .
## factor(Munic)Glamoc 3.423e-01 2.219e-01 1.543 0.122984
## factor(Munic)Gorazde 3.090e-01 1.576e-01 1.960 0.050030 .
## factor(Munic)Gracanica -3.216e-02 1.500e-01 -0.214 0.830204
## factor(Munic)Grad Mostar 4.517e-01 1.252e-01 3.607 0.000313 ***
## factor(Munic)Gradacac 7.671e-02 1.489e-01 0.515 0.606402
## factor(Munic)Gradiska -1.547e-01 1.853e-01 -0.835 0.403786
## factor(Munic)Grude 4.921e-01 1.667e-01 2.953 0.003165 **
## factor(Munic)I.Ilidza 8.632e-02 1.978e-01 0.436 0.662619
## factor(Munic)Ilidza 6.618e-01 1.441e-01 4.591 4.51e-06 ***
## factor(Munic)Ilijas 1.746e-01 1.942e-01 0.899 0.368551
## factor(Munic)Jajce 4.614e-01 1.879e-01 2.455 0.014117 *
## factor(Munic)K.Dubica -7.096e-02 2.104e-01 -0.337 0.735935
## factor(Munic)Kakanj 1.296e-01 1.556e-01 0.833 0.404906
## factor(Munic)Kalesija -2.267e-01 1.562e-01 -1.451 0.146805
## factor(Munic)Kiseljak 3.583e-01 1.597e-01 2.244 0.024860 *
## factor(Munic)Kladanj 1.564e-01 2.101e-01 0.744 0.456721
## factor(Munic)Kljuc 7.723e-02 1.977e-01 0.391 0.696088
## factor(Munic)Knezevo 3.065e-02 2.281e-01 0.134 0.893147
## factor(Munic)Konjic 3.498e-01 1.411e-01 2.478 0.013242 *
## factor(Munic)Kotor Varos 2.126e-01 2.213e-01 0.961 0.336733
## factor(Munic)Kresevo 3.699e-01 1.808e-01 2.046 0.040790 *
## factor(Munic)Kupres 2.924e-01 2.054e-01 1.424 0.154633
## factor(Munic)Laktasi 1.004e-01 1.534e-01 0.655 0.512755
## factor(Munic)Livno 5.323e-01 1.615e-01 3.296 0.000987 ***
## factor(Munic)Lopare 1.024e-01 1.520e-01 0.674 0.500635
## factor(Munic)Lukavac 2.770e-01 1.585e-01 1.747 0.080700 .
## factor(Munic)Ljubuski 6.586e-01 1.561e-01 4.220 2.49e-05 ***
## factor(Munic)Maglaj 1.622e-01 1.534e-01 1.058 0.290327
## factor(Munic)Modrica 5.719e-02 1.448e-01 0.395 0.692976
## factor(Munic)Mrkonjic 3.583e-01 2.361e-01 1.518 0.129093
## factor(Munic)N.G.Sarajevo 4.364e-01 1.241e-01 3.517 0.000440 ***
## factor(Munic)N.Gorazde -2.803e-01 2.058e-01 -1.362 0.173256
## factor(Munic)N.Sarajevo 5.652e-01 1.425e-01 3.968 7.36e-05 ***
## factor(Munic)N.Travnik -2.957e-02 1.377e-01 -0.215 0.829979
## factor(Munic)Nevesinje 2.100e-02 2.282e-01 0.092 0.926689
## factor(Munic)Novi Grad 4.820e-02 1.443e-01 0.334 0.738276
## factor(Munic)Orasje 4.634e-01 1.653e-01 2.803 0.005087 **
## factor(Munic)Pale 1.427e-01 1.807e-01 0.789 0.429882
## factor(Munic)Posusje 5.560e-01 1.882e-01 2.955 0.003146 **
## factor(Munic)Prijedor 2.122e-01 1.368e-01 1.551 0.120889
## factor(Munic)Prnjavor 1.613e-01 1.472e-01 1.096 0.273280
## factor(Munic)Prozor 3.468e-01 2.018e-01 1.719 0.085722 .
## factor(Munic)Ribnik -1.538e-01 1.943e-01 -0.792 0.428510
## factor(Munic)Rogatica 2.268e-01 1.857e-01 1.221 0.222057
## factor(Munic)Rudo -6.754e-02 1.855e-01 -0.364 0.715744
## factor(Munic)S.Brijeg 6.994e-01 1.453e-01 4.815 1.52e-06 ***
## factor(Munic)S.G.Sarajevo 4.555e-01 1.497e-01 3.042 0.002361 **
## factor(Munic)Samac 1.256e-01 1.367e-01 0.919 0.358399
## factor(Munic)Sanski Most 2.718e-01 1.353e-01 2.010 0.044534 *
## factor(Munic)Sokolac -1.078e-01 2.285e-01 -0.472 0.637095
## factor(Munic)Srbac 1.300e-01 1.941e-01 0.670 0.503094
## factor(Munic)Srebrenica 3.142e-01 2.359e-01 1.332 0.182936
## factor(Munic)Srebrenik 2.136e-01 2.355e-01 0.907 0.364478
## factor(Munic)Stolac 4.472e-01 2.017e-01 2.218 0.026634 *
## factor(Munic)Tesanj 7.376e-02 1.393e-01 0.530 0.596471
## factor(Munic)Teslic 6.114e-02 1.519e-01 0.403 0.687276
## factor(Munic)Tomislavgrad 7.188e-01 1.882e-01 3.820 0.000135 ***
## factor(Munic)Travnik 4.027e-01 1.286e-01 3.132 0.001746 **
## factor(Munic)Trebinje 1.325e-01 1.568e-01 0.845 0.398210
## factor(Munic)Tuzla 1.185e-02 1.330e-01 0.089 0.928997
## factor(Munic)Ugljevik 3.098e-01 1.784e-01 1.736 0.082540 .
## factor(Munic)V.Kladusa 2.842e-02 1.434e-01 0.198 0.842858
## factor(Munic)Visegrad 5.109e-02 1.980e-01 0.258 0.796355
## factor(Munic)Visoko 2.156e-01 1.447e-01 1.490 0.136255
## factor(Munic)Vitez 4.461e-02 1.376e-01 0.324 0.745823
## factor(Munic)Vlasenica 1.797e-01 1.975e-01 0.910 0.362989
## factor(Munic)Vogosca 2.010e-01 1.765e-01 1.139 0.254865
## factor(Munic)Zavidovici 1.088e-01 1.406e-01 0.774 0.438869
## factor(Munic)Zenica 1.751e-01 1.261e-01 1.389 0.165043
## factor(Munic)Zepce 1.727e-01 2.156e-01 0.801 0.423149
## factor(Munic)Zivinice -1.905e-01 1.407e-01 -1.354 0.175862
## factor(Munic)Zvornik -2.376e-01 1.416e-01 -1.677 0.093534 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.6827 on 4987 degrees of freedom
## Multiple R-squared: 0.3618, Adjusted R-squared: 0.3481
## F-statistic: 26.42 on 107 and 4987 DF, p-value: < 2.2e-16summary(WI.3a)##
## Call:
## lm(formula = log(wages) ~ Age + as.numeric(Age^2) + factor(Sex_) +
## Educ.1 + factor(Full.time) + factor(Munic), data = WI.dfa)
##
## Residuals:
## Min 1Q Median 3Q Max
## -2.54025 -0.23473 0.06767 0.32865 2.32918
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 4.882e+00 1.385e-01 35.239 < 2e-16 ***
## Age 4.106e-02 4.939e-03 8.315 < 2e-16 ***
## as.numeric(Age^2) -3.956e-04 5.834e-05 -6.781 1.35e-11 ***
## factor(Sex_)2 -1.462e-01 1.734e-02 -8.431 < 2e-16 ***
## Educ.1elementary 1.655e-01 6.756e-02 2.450 0.014338 *
## Educ.1high schl 4.823e-01 6.571e-02 7.340 2.53e-13 ***
## Educ.1faculty 9.292e-01 6.791e-02 13.683 < 2e-16 ***
## factor(Full.time)part time -7.410e-01 3.014e-02 -24.587 < 2e-16 ***
## factor(Munic)B.Petrovac -1.828e-01 1.881e-01 -0.972 0.331014
## factor(Munic)Banovici 3.017e-01 1.435e-01 2.102 0.035622 *
## factor(Munic)Banja Luka 1.150e-01 1.018e-01 1.130 0.258490
## factor(Munic)Bihac -8.512e-02 1.110e-01 -0.767 0.443091
## factor(Munic)Bijeljina -1.147e-01 1.011e-01 -1.135 0.256639
## factor(Munic)Bileca 4.909e-02 1.771e-01 0.277 0.781594
## factor(Munic)Bratunac 3.501e-01 1.768e-01 1.981 0.047707 *
## factor(Munic)Brcko 1.161e-01 1.001e-01 1.160 0.246258
## factor(Munic)Breza 2.315e-01 1.615e-01 1.433 0.151933
## factor(Munic)Brod -2.568e-01 1.402e-01 -1.832 0.067015 .
## factor(Munic)Bugojno 3.926e-02 1.140e-01 0.344 0.730521
## factor(Munic)Busovaca 6.027e-02 1.354e-01 0.445 0.656223
## factor(Munic)Buzim -6.659e-01 1.767e-01 -3.768 0.000167 ***
## factor(Munic)C.Sarajevo 3.420e-01 1.193e-01 2.867 0.004157 **
## factor(Munic)Cajnice -1.785e-01 1.686e-01 -1.059 0.289878
## factor(Munic)Capljina 4.760e-01 1.400e-01 3.400 0.000679 ***
## factor(Munic)Cazin -1.264e-01 1.145e-01 -1.104 0.269596
## factor(Munic)Celinac -1.168e-01 1.880e-01 -0.621 0.534620
## factor(Munic)Citluk 2.928e-01 1.391e-01 2.105 0.035380 *
## factor(Munic)D.Vakuf -1.137e-01 1.259e-01 -0.904 0.366228
## factor(Munic)Derventa -4.717e-01 1.294e-01 -3.646 0.000270 ***
## factor(Munic)Doboj 1.505e-01 1.133e-01 1.328 0.184221
## factor(Munic)Doboj-Jug -1.087e-01 1.950e-01 -0.557 0.577387
## factor(Munic)Drvar 2.669e-02 1.535e-01 0.174 0.861975
## factor(Munic)Foca - RS -3.240e-02 1.355e-01 -0.239 0.811095
## factor(Munic)Fojnica 2.497e-01 1.723e-01 1.449 0.147484
## factor(Munic)G.Vakuf -1.004e-01 1.326e-01 -0.757 0.449049
## factor(Munic)Gacko 2.474e-01 1.404e-01 1.762 0.078135 .
## factor(Munic)Glamoc 2.371e-01 1.773e-01 1.337 0.181347
## factor(Munic)Gorazde 2.352e-01 1.289e-01 1.825 0.068039 .
## factor(Munic)Gracanica -1.288e-01 1.260e-01 -1.022 0.306996
## factor(Munic)Grad Mostar 3.323e-01 1.032e-01 3.219 0.001296 **
## factor(Munic)Gradacac 7.144e-02 1.237e-01 0.578 0.563598
## factor(Munic)Gradiska -4.956e-02 1.509e-01 -0.328 0.742638
## factor(Munic)Grude 4.016e-01 1.354e-01 2.965 0.003044 **
## factor(Munic)I.Ilidza -1.296e-02 1.615e-01 -0.080 0.936055
## factor(Munic)Ilidza 5.109e-01 1.214e-01 4.209 2.61e-05 ***
## factor(Munic)Ilijas 1.766e-01 1.646e-01 1.073 0.283470
## factor(Munic)Jajce 2.768e-01 1.556e-01 1.779 0.075271 .
## factor(Munic)K.Dubica -1.379e-01 1.723e-01 -0.801 0.423451
## factor(Munic)Kakanj 7.807e-02 1.280e-01 0.610 0.541813
## factor(Munic)Kalesija -1.259e-01 1.436e-01 -0.876 0.380820
## factor(Munic)Kiseljak 2.476e-01 1.314e-01 1.884 0.059658 .
## factor(Munic)Kladanj 1.331e-01 1.768e-01 0.753 0.451400
## factor(Munic)Kljuc -1.006e-02 1.646e-01 -0.061 0.951270
## factor(Munic)Knezevo -2.828e-02 1.821e-01 -0.155 0.876559
## factor(Munic)Konjic 2.337e-01 1.161e-01 2.012 0.044256 *
## factor(Munic)Kotor Varos 1.407e-01 1.768e-01 0.796 0.426000
## factor(Munic)Kresevo 3.693e-01 1.471e-01 2.511 0.012066 *
## factor(Munic)Kupres 2.037e-01 1.643e-01 1.240 0.215175
## factor(Munic)Laktasi 2.353e-02 1.241e-01 0.190 0.849623
## factor(Munic)Livno 4.060e-01 1.331e-01 3.050 0.002305 **
## factor(Munic)Lopare -1.237e-02 1.261e-01 -0.098 0.921844
## factor(Munic)Lukavac 1.165e-01 1.343e-01 0.867 0.385883
## factor(Munic)Ljubuski 4.160e-01 1.345e-01 3.093 0.001995 **
## factor(Munic)Maglaj 6.616e-02 1.378e-01 0.480 0.631132
## factor(Munic)Modrica -6.353e-02 1.195e-01 -0.532 0.595000
## factor(Munic)Mrkonjic 2.657e-01 1.952e-01 1.361 0.173453
## factor(Munic)N.G.Sarajevo 3.454e-01 1.022e-01 3.379 0.000735 ***
## factor(Munic)N.Gorazde -4.083e-01 1.950e-01 -2.094 0.036340 *
## factor(Munic)N.Sarajevo 4.299e-01 1.194e-01 3.600 0.000321 ***
## factor(Munic)N.Travnik -1.673e-01 1.125e-01 -1.487 0.137135
## factor(Munic)Nevesinje -3.494e-02 1.822e-01 -0.192 0.847908
## factor(Munic)Novi Grad -5.045e-02 1.193e-01 -0.423 0.672377
## factor(Munic)Orasje 3.446e-01 1.378e-01 2.501 0.012413 *
## factor(Munic)Pale 7.309e-02 1.452e-01 0.503 0.614781
## factor(Munic)Posusje 4.060e-01 1.558e-01 2.607 0.009176 **
## factor(Munic)Prijedor 1.439e-01 1.120e-01 1.285 0.198752
## factor(Munic)Prnjavor 1.964e-02 1.242e-01 0.158 0.874379
## factor(Munic)Prozor 2.234e-01 1.772e-01 1.261 0.207409
## factor(Munic)Ribnik -1.957e-02 1.614e-01 -0.121 0.903513
## factor(Munic)Rogatica 1.831e-01 1.536e-01 1.192 0.233289
## factor(Munic)Rudo -1.642e-01 1.615e-01 -1.017 0.309369
## factor(Munic)S.Brijeg 6.098e-01 1.190e-01 5.124 3.13e-07 ***
## factor(Munic)S.G.Sarajevo 3.636e-01 1.218e-01 2.986 0.002844 **
## factor(Munic)Samac 2.103e-02 1.124e-01 0.187 0.851528
## factor(Munic)Sanski Most -1.254e-02 1.129e-01 -0.111 0.911530
## factor(Munic)Sokolac -1.643e-01 1.824e-01 -0.901 0.367796
## factor(Munic)Srbac 5.756e-02 1.583e-01 0.363 0.716253
## factor(Munic)Srebrenica 3.278e-01 1.949e-01 1.681 0.092769 .
## factor(Munic)Srebrenik 1.607e-01 1.879e-01 0.856 0.392305
## factor(Munic)Stolac 3.132e-01 1.648e-01 1.901 0.057387 .
## factor(Munic)Tesanj 1.299e-03 1.188e-01 0.011 0.991277
## factor(Munic)Teslic -3.769e-02 1.288e-01 -0.293 0.769767
## factor(Munic)Tomislavgrad 1.667e-01 1.821e-01 0.916 0.359887
## factor(Munic)Travnik 3.328e-01 1.056e-01 3.152 0.001635 **
## factor(Munic)Trebinje 3.394e-02 1.274e-01 0.266 0.790013
## factor(Munic)Tuzla 1.673e-01 1.119e-01 1.495 0.134965
## factor(Munic)Ugljevik 1.734e-01 1.434e-01 1.209 0.226649
## factor(Munic)V.Kladusa -2.492e-01 1.243e-01 -2.005 0.045074 *
## factor(Munic)Visegrad -1.453e-01 1.685e-01 -0.862 0.388509
## factor(Munic)Visoko 1.286e-01 1.181e-01 1.089 0.276360
## factor(Munic)Vitez -6.815e-02 1.131e-01 -0.602 0.546978
## factor(Munic)Vlasenica 1.152e-01 1.583e-01 0.728 0.466599
## factor(Munic)Vogosca 1.412e-01 1.435e-01 0.983 0.325438
## factor(Munic)Zavidovici -4.154e-02 1.198e-01 -0.347 0.728695
## factor(Munic)Zenica 1.013e-01 1.038e-01 0.975 0.329371
## factor(Munic)Zivinice -5.412e-02 1.208e-01 -0.448 0.654192
## factor(Munic)Zvornik -3.265e-01 1.163e-01 -2.807 0.005018 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.5385 on 4404 degrees of freedom
## Multiple R-squared: 0.4012, Adjusted R-squared: 0.3868
## F-statistic: 27.84 on 106 and 4404 DF, p-value: < 2.2e-16FAZA 8: IZRADI MAPE ZA Winner model
Mean map
open libraries
Import map
library(readr)
BIH_adm3 <- read_csv("BIH_adm/BIH_adm3.csv")## Rows: 141 Columns: 14
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (8): ISO, NAME_0, NAME_1, NAME_2, NAME_3, TYPE_3, ENGTYPE_3, VARNAME_3
## dbl (5): OBJECTID, ID_0, ID_1, ID_2, ID_3
## lgl (1): NL_NAME_3
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.amraGIS <- st_read("BIH_adm/BIH_adm3.shp")## Reading layer `BIH_adm3' from data source
## `C:\Users\Amra\Documents\HBS database inspecting\BIH_adm\BIH_adm3.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 141 features and 13 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 15.72739 ymin: 42.56531 xmax: 19.61471 ymax: 45.27468
## Geodetic CRS: WGS 84#change interpunction Nema potrebe za ovim jer cu formirati varijablu opcina koju sam prilagodial hbs bazi
Mean wages
U skladu sa preporukom Ismirovom od dana 30.10.2023. koje su se odnosile an rents, je cu da uradim mean vrijednosti za ove opstine wages i to cu mapirati
tapply(WI.df$wages, WI.df$Munic, summary)## $B.Krupa
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.0 287.5 400.0 576.9 797.4 1858.3
##
## $B.Petrovac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 175.0 335.0 700.0 805.6 1143.4 2083.3
##
## $Banovici
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 530.2 920.8 857.2 1106.7 1556.7
##
## $`Banja Luka`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.083 500.000 690.000 790.656 1000.000 3353.333
##
## $Bihac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 10.0 360.0 600.0 676.6 900.0 3000.0
##
## $Bijeljina
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.0 300.0 500.0 744.7 866.7 8333.3
##
## $Bileca
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 200.0 370.0 500.0 598.5 800.0 1200.0
##
## $Bratunac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 200.0 520.0 700.0 800.8 1050.0 1500.0
##
## $Brcko
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 33.33 400.00 600.00 707.53 850.00 3500.00
##
## $Breza
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 100 860 1000 891 1100 1500
##
## $Brod
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 8.333 70.833 414.583 407.961 660.417 916.667
##
## $Bugojno
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 33.33 250.00 400.00 534.18 700.00 2030.00
##
## $Busovaca
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 75.0 208.3 800.0 703.7 1000.0 2000.0
##
## $Buzim
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 20.00 66.67 104.17 346.88 290.62 1391.67
##
## $C.Sarajevo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 53.33 615.00 1032.50 1403.79 1712.50 18350.00
##
## $Cajnice
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 62.5 218.8 550.0 501.6 725.0 950.0
##
## $Capljina
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 208.3 625.0 800.0 932.8 1248.8 2061.7
##
## $Cazin
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 33.33 237.50 500.00 554.41 700.00 3500.00
##
## $Celinac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 166.7 450.0 568.3 583.7 632.7 1200.0
##
## $Citluk
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 58.33 425.00 1037.50 1005.39 1350.00 3220.00
##
## $D.Vakuf
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 145.8 300.0 436.2 550.0 1777.2
##
## $Derventa
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 41.67 183.33 275.00 342.39 458.33 1400.00
##
## $Doboj
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 187.5 493.9 670.0 724.0 847.5 1700.0
##
## $`Doboj-Jug`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 291.7 400.0 410.0 491.7 550.0 920.0
##
## $Drvar
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 150.0 345.8 533.3 580.7 725.0 1300.0
##
## $`Foca - RS`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 58.33 385.00 600.00 719.68 935.00 2000.00
##
## $Fojnica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 266.7 420.0 500.0 623.6 875.0 1200.0
##
## $G.Vakuf
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.0 100.0 331.0 443.1 600.0 1875.0
##
## $Gacko
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 400.0 575.0 800.0 826.2 1062.5 1500.0
##
## $Glamoc
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 200.0 400.0 400.0 498.5 650.0 1100.0
##
## $Gorazde
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 250.0 522.2 660.0 716.8 800.0 2268.0
##
## $Gracanica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 29.17 293.75 450.00 499.11 637.50 1756.67
##
## $`Grad Mostar`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 58.33 600.00 900.00 1041.70 1281.67 3300.00
##
## $Gradacac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 33.33 300.00 500.00 551.46 775.00 1536.00
##
## $Gradiska
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 8.0 360.0 433.3 482.3 595.0 1000.0
##
## $Grude
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 58.33 800.00 900.00 1042.76 1250.00 2000.00
##
## $I.Ilidza
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 166.7 400.0 575.0 667.7 837.5 1500.0
##
## $Ilidza
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 83.33 650.00 885.00 1029.04 1250.00 3873.00
##
## $Ilijas
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 71.67 380.00 641.67 607.47 837.50 1000.00
##
## $Jajce
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 227.5 500.0 725.0 828.2 1010.0 2060.0
##
## $K.Dubica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 308.3 370.0 450.0 562.4 555.0 1733.3
##
## $Kakanj
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 12.5 350.0 700.0 900.6 1256.2 5000.0
##
## $Kalesija
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 33.33 87.50 346.67 474.34 700.00 2250.00
##
## $Kiseljak
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 100.0 600.0 760.0 734.8 900.0 1200.0
##
## $Kladanj
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 100.0 530.0 690.0 680.0 882.5 1096.7
##
## $Kljuc
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 80.0 350.0 762.5 682.3 937.5 1260.0
##
## $Knezevo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 137.5 487.5 620.8 679.2 949.0 1155.0
##
## $Konjic
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 100.0 400.0 700.0 791.8 1116.7 2700.0
##
## $`Kotor Varos`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 370.0 504.2 580.0 630.5 765.0 962.5
##
## $Kresevo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 41.67 687.50 800.00 740.00 877.50 1050.00
##
## $Kupres
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 133.3 208.3 970.0 870.7 1251.7 2000.0
##
## $Laktasi
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 70.0 458.3 550.0 636.4 800.0 1530.0
##
## $Livno
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 291.7 650.0 880.0 943.3 1200.0 2000.0
##
## $Lopare
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.0 262.5 400.0 486.1 600.0 1766.7
##
## $Lukavac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 86.67 430.00 600.00 750.91 1062.83 2500.00
##
## $Ljubuski
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 166.7 800.0 1200.0 1495.2 1587.5 12000.0
##
## $Maglaj
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 100.0 400.0 610.0 684.1 900.0 2000.0
##
## $Modrica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 83.33 243.75 408.33 491.21 666.67 2000.00
##
## $Mrkonjic
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 400.0 525.0 800.0 745.5 900.0 1200.0
##
## $N.G.Sarajevo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50 700 1000 1020 1253 3000
##
## $N.Gorazde
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 44.5 333.3 370.0 505.8 425.5 2000.0
##
## $N.Sarajevo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50 675 1000 1254 1539 5000
##
## $N.Travnik
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 66.67 200.00 442.08 481.61 655.25 3000.00
##
## $Nevesinje
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 370.0 400.0 580.0 715.0 962.5 1500.0
##
## $`Novi Grad`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 66.67 400.00 600.00 692.00 850.00 4000.00
##
## $Orasje
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.0 600.0 800.0 996.5 1200.0 2750.0
##
## $Pale
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 275.0 445.0 625.0 656.3 825.0 1250.0
##
## $Posusje
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 500 800 800 1012 1200 2000
##
## $Prijedor
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 452.1 666.7 672.0 936.7 1237.5
##
## $Prnjavor
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 400.0 600.0 729.2 750.0 3000.0
##
## $Prozor
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 125.0 400.0 500.0 751.9 1250.0 2000.0
##
## $Ribnik
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 8.333 425.000 600.000 589.035 750.000 1500.000
##
## $Rogatica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 512.5 775.0 993.9 1187.5 5390.0
##
## $Rudo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 185.0 370.0 500.0 633.9 746.2 2041.7
##
## $S.Brijeg
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 500 940 1000 1104 1300 2000
##
## $S.G.Sarajevo
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 300.0 652.1 1000.0 989.9 1241.5 2600.0
##
## $Samac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 333.3 533.3 630.9 800.0 2960.0
##
## $`Sanski Most`
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 83.33 458.33 695.75 908.12 990.62 4583.33
##
## $Sokolac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 70.0 505.0 685.0 785.8 1050.0 1620.0
##
## $Srbac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 166.7 400.0 480.0 585.8 740.0 1380.0
##
## $Srebrenica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 400 540 720 840 1100 1680
##
## $Srebrenik
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 300.0 533.3 600.0 571.2 633.3 850.0
##
## $Stolac
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 200.0 600.0 916.7 932.9 1183.3 1683.3
##
## $Tesanj
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 54.17 400.00 505.00 604.33 675.00 2200.00
##
## $Teslic
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 58.33 366.67 476.67 540.46 715.00 1458.33
##
## $Tomislavgrad
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 150 700 1000 1335 1400 4000
##
## $Travnik
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50.0 600.0 700.0 792.0 953.8 2000.0
##
## $Trebinje
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 217.0 400.0 600.0 684.1 800.0 2800.0
##
## $Tuzla
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 4.0 385.0 576.7 651.7 905.0 2100.0
##
## $Ugljevik
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 120 400 650 673 960 1470
##
## $V.Kladusa
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 52.5 214.2 525.0 578.3 787.5 1950.0
##
## $Visegrad
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 370.0 377.5 410.0 705.6 875.8 2800.0
##
## $Visoko
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 62.5 400.0 600.0 652.4 760.0 2018.0
##
## $Vitez
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 66.67 250.00 550.00 592.42 650.00 4000.00
##
## $Vlasenica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 133.3 407.5 620.0 701.9 815.0 1513.3
##
## $Vogosca
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 83.33 431.25 606.50 815.19 1073.50 2300.00
##
## $Zavidovici
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 104.2 400.0 600.0 638.9 800.0 2000.0
##
## $Zenica
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 37.5 450.0 700.0 736.6 993.3 2220.0
##
## $Zepce
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 33.33 400.00 600.00 859.52 1000.00 3000.00
##
## $Zivinice
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 12.5 150.0 397.5 544.4 877.1 2625.0
##
## $Zvornik
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 50 250 400 676 750 12833dim (WI.df)## [1] 5095 16mean.wages <- WI.df %>% dplyr::group_by(Munic) %>%summarise_at(dplyr::vars("wages"), list(~n(), ~mean(., na.rm = TRUE), ~sd(., na.rm = TRUE)))
mean.wages %>% arrange (desc(mean)) %>% print(n=103)## # A tibble: 101 × 4
## Munic n mean sd
## <chr> <int> <dbl> <dbl>
## 1 Ljubuski 42 1495. 1793.
## 2 C.Sarajevo 68 1404. 2193.
## 3 Tomislavgrad 21 1335. 1002.
## 4 N.Sarajevo 67 1254. 860.
## 5 S.Brijeg 59 1104. 317.
## 6 Grude 32 1043. 469.
## 7 Grad Mostar 190 1042. 612.
## 8 Ilidza 61 1029. 733.
## 9 N.G.Sarajevo 217 1020. 532.
## 10 Posusje 21 1012. 444.
## 11 Citluk 30 1005. 677.
## 12 Orasje 33 997. 628.
## 13 Rogatica 22 994. 1048.
## 14 S.G.Sarajevo 51 990. 470.
## 15 Livno 36 943. 399.
## 16 Stolac 17 933. 436.
## 17 Capljina 28 933. 431.
## 18 Sanski Most 90 908. 818.
## 19 Kakanj 42 901. 830.
## 20 Breza 17 891. 374.
## 21 Kupres 16 871. 575.
## 22 Zepce 14 860. 789.
## 23 Banovici 26 857. 350.
## 24 Srebrenica 11 840 392.
## 25 Jajce 21 828. 468.
## 26 Gacko 28 826. 307.
## 27 Vogosca 26 815. 545.
## 28 B.Petrovac 12 806. 592.
## 29 Bratunac 13 801. 379.
## 30 Travnik 136 792. 305.
## 31 Konjic 69 792. 538.
## 32 Banja Luka 236 791. 510.
## 33 Sokolac 12 786. 481.
## 34 Prozor 17 752. 526.
## 35 Lukavac 39 751. 502.
## 36 Mrkonjic 11 745. 246.
## 37 Bijeljina 245 745. 866.
## 38 Kresevo 24 740 236.
## 39 Zenica 174 737. 398.
## 40 Kiseljak 38 735. 249.
## 41 Prnjavor 55 729. 582.
## 42 Doboj 78 724. 321.
## 43 Foca - RS 31 720. 460.
## 44 Gorazde 40 717. 349.
## 45 Nevesinje 12 715 379.
## 46 Brcko 312 708. 502.
## 47 Visegrad 18 706. 594.
## 48 Busovaca 33 704. 486.
## 49 Vlasenica 18 702. 397.
## 50 Novi Grad 61 692. 560.
## 51 Maglaj 45 684. 397.
## 52 Trebinje 41 684. 444.
## 53 Kljuc 18 682. 377.
## 54 Kladanj 15 680. 288.
## 55 Knezevo 12 679. 311.
## 56 Bihac 97 677. 500.
## 57 Zvornik 67 676. 1542.
## 58 Ugljevik 25 673 408.
## 59 Prijedor 84 672. 303.
## 60 I.Ilidza 18 668. 361.
## 61 Pale 24 656. 282.
## 62 Visoko 60 652. 371.
## 63 Tuzla 104 652. 430.
## 64 Zavidovici 71 639. 330.
## 65 Laktasi 45 636. 320.
## 66 Rudo 22 634. 449.
## 67 Samac 85 631. 487.
## 68 Kotor Varos 13 631. 180.
## 69 Fojnica 14 624. 281.
## 70 Ilijas 19 607. 269.
## 71 Tesanj 74 604. 384.
## 72 Bileca 13 598. 319.
## 73 Vitez 81 592. 541.
## 74 Ribnik 19 589. 326.
## 75 Srbac 19 586. 285.
## 76 Celinac 11 584. 298.
## 77 Drvar 20 581. 301.
## 78 V.Kladusa 62 578. 455.
## 79 B.Krupa 36 577. 433.
## 80 Srebrenik 11 571. 148.
## 81 K.Dubica 15 562. 365.
## 82 Cazin 72 554. 496.
## 83 Gradacac 51 551. 348.
## 84 Zivinice 70 544. 517.
## 85 Teslic 47 540. 282.
## 86 Bugojno 93 534. 430.
## 87 N.Gorazde 16 506. 487.
## 88 Cajnice 16 502. 299.
## 89 Gracanica 50 499. 328.
## 90 Glamoc 13 498. 248.
## 91 Doboj-Jug 11 492. 182.
## 92 Modrica 60 491. 376.
## 93 Lopare 48 486. 383.
## 94 Gradiska 22 482. 246.
## 95 N.Travnik 80 482. 410.
## 96 Kalesija 42 474. 493.
## 97 G.Vakuf 44 443. 452.
## 98 D.Vakuf 52 436. 412.
## 99 Brod 38 408. 321.
## 100 Buzim 16 347. 471.
## 101 Derventa 44 342. 256.mean.wages.a <- WI.dfa %>% dplyr::group_by(Munic) %>% summarise_at(dplyr::vars("wages"), list(~n(), ~mean(., na.rm = TRUE), ~sd(., na.rm = TRUE)))
mean.wages.a %>% arrange (desc(mean)) %>% print(n=103)## # A tibble: 100 × 4
## Munic n mean sd
## <chr> <int> <dbl> <dbl>
## 1 C.Sarajevo 57 1089. 555.
## 2 S.Brijeg 56 1087. 313.
## 3 Ljubuski 32 1073. 523.
## 4 N.Sarajevo 56 1070. 518.
## 5 Grude 31 1008. 440.
## 6 N.G.Sarajevo 201 950. 398.
## 7 S.G.Sarajevo 50 947. 405.
## 8 Capljina 27 945 426.
## 9 Stolac 16 944. 431.
## 10 Grad Mostar 176 941. 443.
## 11 Citluk 28 940. 505.
## 12 Posusje 19 919. 349.
## 13 Livno 33 892. 361.
## 14 Srebrenica 10 884 383.
## 15 Breza 17 881. 370.
## 16 Ilidza 50 874. 441.
## 17 Kupres 16 869. 573.
## 18 Orasje 29 856. 401.
## 19 Rogatica 20 821. 343.
## 20 Vogosca 25 815. 515.
## 21 Gacko 27 812. 303.
## 22 Banovici 25 806. 243.
## 23 Bratunac 13 801. 379.
## 24 Sokolac 12 786. 481.
## 25 Kakanj 39 783. 468.
## 26 Travnik 130 778. 286.
## 27 Banja Luka 218 777. 416.
## 28 Mrkonjic 10 770 245.
## 29 Kresevo 23 757. 178.
## 30 Tomislavgrad 12 752. 278.
## 31 Kiseljak 35 742. 257.
## 32 Jajce 19 742. 438.
## 33 Tuzla 82 741. 399.
## 34 Zenica 163 738. 388.
## 35 Konjic 64 735. 446.
## 36 Kladanj 13 731. 256.
## 37 Gorazde 38 725. 352.
## 38 Doboj 75 718. 304.
## 39 Nevesinje 12 715 379.
## 40 Prozor 13 700. 451.
## 41 Vlasenica 18 698. 393.
## 42 Maglaj 29 690. 375.
## 43 Lukavac 32 689. 392.
## 44 I.Ilidza 17 686. 363.
## 45 Foca - RS 31 683. 390.
## 46 Knezevo 12 679. 311.
## 47 Kljuc 16 679. 342.
## 48 Busovaca 31 668. 472.
## 49 Prijedor 81 660. 282.
## 50 Pale 24 656. 282.
## 51 Visoko 58 652. 378.
## 52 Rudo 17 649. 489.
## 53 Bihac 87 647. 381.
## 54 Ugljevik 25 645. 397.
## 55 Novi Grad 55 640. 311.
## 56 Ilijas 16 637. 262.
## 57 Laktasi 45 636. 320.
## 58 Brcko 267 634. 352.
## 59 B.Petrovac 11 633. 417.
## 60 Trebinje 40 631. 291.
## 61 Kotor Varos 13 631. 180.
## 62 Sanski Most 76 619. 330.
## 63 Fojnica 14 612. 290.
## 64 Prnjavor 45 610. 255.
## 65 Zavidovici 54 610. 300.
## 66 Bijeljina 213 607. 425.
## 67 Ribnik 17 607. 306.
## 68 B.Krupa 33 602. 427.
## 69 Bileca 13 598. 319.
## 70 Srbac 18 596. 289.
## 71 Zivinice 51 595. 376.
## 72 Samac 80 593. 415.
## 73 Kalesija 25 591. 417.
## 74 Tesanj 56 587. 335.
## 75 Celinac 11 572. 305.
## 76 Srebrenik 11 571. 148.
## 77 Gradacac 45 571. 314.
## 78 Bugojno 73 568. 416.
## 79 Drvar 20 562. 293.
## 80 Visegrad 15 557. 277.
## 81 Teslic 38 545. 220.
## 82 Vitez 76 540. 378.
## 83 K.Dubica 14 538. 256.
## 84 Cazin 68 531. 340.
## 85 Cajnice 15 525 279.
## 86 Brod 27 514. 277.
## 87 Gradiska 21 505. 228.
## 88 Zvornik 63 504. 323.
## 89 Doboj-Jug 10 500. 190.
## 90 Glamoc 13 498. 248.
## 91 V.Kladusa 44 487. 407.
## 92 D.Vakuf 43 481. 421.
## 93 Lopare 43 471. 338.
## 94 Modrica 55 468. 348.
## 95 Gracanica 42 463. 251.
## 96 N.Travnik 78 436. 273.
## 97 Buzim 13 419. 492.
## 98 G.Vakuf 34 416. 394.
## 99 N.Gorazde 10 389. 145.
## 100 Derventa 37 345. 192.Spasi mean wages
#spasiti a u excel
library(openxlsx)
wb = createWorkbook()
addWorksheet(wb, "mean.wages.onlyempl") #ovo je naziv sheeta
addWorksheet(wb, "mean.notall_wages.onlyempl") #ovo je naziv sheeta
library(expss)## Loading required package: maditr##
## To select columns from data: columns(mtcars, mpg, vs:carb)##
## Attaching package: 'maditr'## The following object is masked from 'package:raster':
##
## shift## The following object is masked from 'package:purrr':
##
## transpose## The following object is masked from 'package:readr':
##
## cols## The following objects are masked from 'package:dplyr':
##
## between, coalesce, first, last##
## Attaching package: 'expss'## The following object is masked from 'package:raster':
##
## calc## The following objects are masked from 'package:stringr':
##
## fixed, regex## The following objects are masked from 'package:purrr':
##
## keep, modify, modify_if, when## The following object is masked from 'package:ggplot2':
##
## vars## The following objects are masked from 'package:tidyr':
##
## contains, nest## The following objects are masked from 'package:dplyr':
##
## compute, contains, na_if, recode, vars, where## The following objects are masked from 'package:haven':
##
## is.labelled, read_spsswriteData(wb,"mean.wages.onlyempl",mean.wages)
writeData(wb,"mean.notall_wages.onlyempl",mean.wages.a)
saveWorkbook(wb, "mean.wages.onlyempl.11.4.2024.xlsx", overwrite = T)map.wages <- merge (mean.wages, amraGIS, by.x='Munic', by.y = "NAME_4", all=T)
#142 muncicip
map.wages.a <- merge (mean.wages.a, amraGIS, by.x='Munic', by.y = "NAME_4", all=T)
#142 muncicipg = ggplot() + geom_sf(data = map.wages, aes(geometry=geometry, fill = mean)) + scale_fill_viridis_c() + scale_x_continuous(breaks = c(18, 48))
g + labs (title = "Mean wages")
g.a = ggplot() + geom_sf(data = map.wages.a, aes(geometry=geometry, fill = mean)) + scale_fill_viridis_c() + scale_x_continuous(breaks = c(18, 48))
g.a + labs (title = "Mean wages")
Maping coefficienst of models
Coeff and maping model WI.3
Extracting the coefficients from WINNING model
model.coef3 <- summary(WI.3)$coefficients
model.coef3 # vidi koliko redova su koef.za druge varijable (8redova) i B.Krupa je referentna opstina## Estimate Std. Error t value
## (Intercept) 4.3659883884 1.643692e-01 26.56208178
## Age 0.0562606318 5.706278e-03 9.85942729
## as.numeric(Age^2) -0.0005561833 6.710941e-05 -8.28770976
## factor(Sex_)2 -0.1977001270 2.085084e-02 -9.48163996
## Educ.1elementary 0.3266644379 7.768435e-02 4.20502247
## Educ.1high schl 0.6098989185 7.563748e-02 8.06344803
## Educ.1faculty 1.0781752145 7.843336e-02 13.74638495
## factor(Full.time)part time -0.8548068085 3.331397e-02 -25.65911204
## factor(Munic)B.Petrovac 0.0705671062 2.282118e-01 0.30921761
## factor(Munic)Banovici 0.3308523627 1.764948e-01 1.87457278
## factor(Munic)Banja Luka 0.0730207936 1.234192e-01 0.59164846
## factor(Munic)Bihac -0.0374563860 1.341190e-01 -0.27927731
## factor(Munic)Bijeljina -0.0031394567 1.222250e-01 -0.02568588
## factor(Munic)Bileca 0.1464476751 2.216551e-01 0.66070068
## factor(Munic)Bratunac 0.3948522744 2.212935e-01 1.78429218
## factor(Munic)Brcko 0.2469952639 1.209289e-01 2.04248285
## factor(Munic)Breza 0.3048600123 2.017294e-01 1.51123243
## factor(Munic)Brod -0.7105596047 1.592071e-01 -4.46311574
## factor(Munic)Bugojno 0.0600375326 1.352707e-01 0.44383246
## factor(Munic)Busovaca 0.1880624501 1.652301e-01 1.13818537
## factor(Munic)Buzim -0.7420081828 2.059813e-01 -3.60230845
## factor(Munic)C.Sarajevo 0.4377795698 1.423936e-01 3.07443219
## factor(Munic)Cajnice -0.2111280915 2.061121e-01 -1.02433631
## factor(Munic)Capljina 0.5350666802 1.722737e-01 3.10591117
## factor(Munic)Cazin -0.1105391303 1.395495e-01 -0.79211440
## factor(Munic)Celinac -0.0180275135 2.357614e-01 -0.07646507
## factor(Munic)Citluk 0.3185605215 1.695997e-01 1.87830837
## factor(Munic)D.Vakuf -0.0103312781 1.495939e-01 -0.06906217
## factor(Munic)Derventa -0.4695186451 1.540906e-01 -3.04703062
## factor(Munic)Doboj 0.2389127664 1.384646e-01 1.72544277
## factor(Munic)Doboj-Jug -0.0627222892 2.358392e-01 -0.26595359
## factor(Munic)Drvar 0.1469117239 1.913604e-01 0.76772256
## factor(Munic)Foca - RS 0.0745178915 1.681781e-01 0.44308913
## factor(Munic)Fojnica 0.3809556599 2.156132e-01 1.76684755
## factor(Munic)G.Vakuf -0.0110627051 1.547287e-01 -0.07149745
## factor(Munic)Gacko 0.3146141433 1.727738e-01 1.82095979
## factor(Munic)Glamoc 0.3423180660 2.219057e-01 1.54262847
## factor(Munic)Gorazde 0.3089617711 1.576185e-01 1.96018690
## factor(Munic)Gracanica -0.0321625146 1.499759e-01 -0.21445118
## factor(Munic)Grad Mostar 0.4517069832 1.252362e-01 3.60683962
## factor(Munic)Gradacac 0.0767090454 1.488777e-01 0.51524885
## factor(Munic)Gradiska -0.1547353027 1.853229e-01 -0.83494952
## factor(Munic)Grude 0.4921280157 1.666709e-01 2.95269333
## factor(Munic)I.Ilidza 0.0863229995 1.978406e-01 0.43632612
## factor(Munic)Ilidza 0.6617860766 1.441358e-01 4.59140592
## factor(Munic)Ilijas 0.1746241686 1.941839e-01 0.89927230
## factor(Munic)Jajce 0.4614114963 1.879375e-01 2.45513236
## factor(Munic)K.Dubica -0.0709565681 2.103909e-01 -0.33726069
## factor(Munic)Kakanj 0.1296385629 1.556357e-01 0.83296135
## factor(Munic)Kalesija -0.2266678449 1.562002e-01 -1.45113648
## factor(Munic)Kiseljak 0.3583194267 1.596614e-01 2.24424547
## factor(Munic)Kladanj 0.1563606899 2.100732e-01 0.74431534
## factor(Munic)Kljuc 0.0772279155 1.977016e-01 0.39062874
## factor(Munic)Knezevo 0.0306463233 2.281420e-01 0.13433005
## factor(Munic)Konjic 0.3497787816 1.411487e-01 2.47808685
## factor(Munic)Kotor Varos 0.2126302150 2.213206e-01 0.96073386
## factor(Munic)Kresevo 0.3698792651 1.807646e-01 2.04619335
## factor(Munic)Kupres 0.2924203493 2.054129e-01 1.42357342
## factor(Munic)Laktasi 0.1003922861 1.533642e-01 0.65460035
## factor(Munic)Livno 0.5323425604 1.615047e-01 3.29614253
## factor(Munic)Lopare 0.1023503584 1.519588e-01 0.67354013
## factor(Munic)Lukavac 0.2769705791 1.585412e-01 1.74699465
## factor(Munic)Ljubuski 0.6585852603 1.560632e-01 4.21998951
## factor(Munic)Maglaj 0.1622368012 1.534128e-01 1.05751771
## factor(Munic)Modrica 0.0571851011 1.448299e-01 0.39484318
## factor(Munic)Mrkonjic 0.3583132832 2.360522e-01 1.51794057
## factor(Munic)N.G.Sarajevo 0.4363612589 1.240567e-01 3.51743348
## factor(Munic)N.Gorazde -0.2802858760 2.057885e-01 -1.36200969
## factor(Munic)N.Sarajevo 0.5652466009 1.424654e-01 3.96760560
## factor(Munic)N.Travnik -0.0295701104 1.377018e-01 -0.21474012
## factor(Munic)Nevesinje 0.0210007713 2.282302e-01 0.09201575
## factor(Munic)Novi Grad 0.0482048374 1.442585e-01 0.33415587
## factor(Munic)Orasje 0.4633726383 1.653308e-01 2.80269920
## factor(Munic)Pale 0.1426644128 1.807119e-01 0.78945797
## factor(Munic)Posusje 0.5560415278 1.881981e-01 2.95455398
## factor(Munic)Prijedor 0.2122313116 1.368072e-01 1.55131731
## factor(Munic)Prnjavor 0.1612972537 1.472146e-01 1.09566044
## factor(Munic)Prozor 0.3467681162 2.017559e-01 1.71875113
## factor(Munic)Ribnik -0.1538460514 1.942969e-01 -0.79180919
## factor(Munic)Rogatica 0.2268020544 1.857162e-01 1.22122901
## factor(Munic)Rudo -0.0675427752 1.854685e-01 -0.36417392
## factor(Munic)S.Brijeg 0.6994496266 1.452742e-01 4.81468508
## factor(Munic)S.G.Sarajevo 0.4554604551 1.497142e-01 3.04219913
## factor(Munic)Samac 0.1255903694 1.367336e-01 0.91850422
## factor(Munic)Sanski Most 0.2718470039 1.352784e-01 2.00953682
## factor(Munic)Sokolac -0.1078236039 2.285403e-01 -0.47179259
## factor(Munic)Srbac 0.1299775432 1.940895e-01 0.66967856
## factor(Munic)Srebrenica 0.3141764735 2.358762e-01 1.33195493
## factor(Munic)Srebrenik 0.2136140785 2.355306e-01 0.90694816
## factor(Munic)Stolac 0.4472067254 2.016712e-01 2.21750440
## factor(Munic)Tesanj 0.0737589190 1.392948e-01 0.52951672
## factor(Munic)Teslic 0.0611353025 1.518596e-01 0.40257775
## factor(Munic)Tomislavgrad 0.7188404067 1.881819e-01 3.81992339
## factor(Munic)Travnik 0.4026913177 1.285736e-01 3.13198976
## factor(Munic)Trebinje 0.1324997382 1.568238e-01 0.84489565
## factor(Munic)Tuzla 0.0118487274 1.329657e-01 0.08911113
## factor(Munic)Ugljevik 0.3097543410 1.783804e-01 1.73648209
## factor(Munic)V.Kladusa 0.0284219244 1.433648e-01 0.19824902
## factor(Munic)Visegrad 0.0510917133 1.979678e-01 0.25808088
## factor(Munic)Visoko 0.2155745610 1.446688e-01 1.49012505
## factor(Munic)Vitez 0.0446100308 1.376133e-01 0.32416949
## factor(Munic)Vlasenica 0.1797038889 1.975273e-01 0.90976755
## factor(Munic)Vogosca 0.2009559193 1.764718e-01 1.13874238
## factor(Munic)Zavidovici 0.1088124342 1.405544e-01 0.77416621
## factor(Munic)Zenica 0.1751121675 1.261150e-01 1.38851208
## factor(Munic)Zepce 0.1726908697 2.155843e-01 0.80103630
## factor(Munic)Zivinice -0.1905094734 1.407222e-01 -1.35379865
## factor(Munic)Zvornik -0.2375947797 1.416479e-01 -1.67736142
## Pr(>|t|)
## (Intercept) 1.628195e-145
## Age 1.005724e-22
## as.numeric(Age^2) 1.469344e-16
## factor(Sex_)2 3.769003e-21
## Educ.1elementary 2.656310e-05
## Educ.1high schl 9.210427e-16
## Educ.1faculty 3.126198e-42
## factor(Full.time)part time 1.708215e-136
## factor(Munic)B.Petrovac 7.571689e-01
## factor(Munic)Banovici 6.090993e-02
## factor(Munic)Banja Luka 5.541128e-01
## factor(Munic)Bihac 7.800436e-01
## factor(Munic)Bijeljina 9.795089e-01
## factor(Munic)Bileca 5.088348e-01
## factor(Munic)Bratunac 7.443698e-02
## factor(Munic)Brcko 4.115614e-02
## factor(Munic)Breza 1.307926e-01
## factor(Munic)Brod 8.255629e-06
## factor(Munic)Bugojno 6.571830e-01
## factor(Munic)Busovaca 2.550978e-01
## factor(Munic)Buzim 3.184768e-04
## factor(Munic)C.Sarajevo 2.120443e-03
## factor(Munic)Cajnice 3.057261e-01
## factor(Munic)Capljina 1.907587e-03
## factor(Munic)Cazin 4.283317e-01
## factor(Munic)Celinac 9.390522e-01
## factor(Munic)Citluk 6.039729e-02
## factor(Munic)D.Vakuf 9.449429e-01
## factor(Munic)Derventa 2.323237e-03
## factor(Munic)Doboj 8.450966e-02
## factor(Munic)Doboj-Jug 7.902859e-01
## factor(Munic)Drvar 4.426884e-01
## factor(Munic)Foca - RS 6.577205e-01
## factor(Munic)Fojnica 7.731492e-02
## factor(Munic)G.Vakuf 9.430047e-01
## factor(Munic)Gacko 6.867286e-02
## factor(Munic)Glamoc 1.229844e-01
## factor(Munic)Gorazde 5.002955e-02
## factor(Munic)Gracanica 8.302040e-01
## factor(Munic)Grad Mostar 3.129832e-04
## factor(Munic)Gradacac 6.064020e-01
## factor(Munic)Gradiska 4.037861e-01
## factor(Munic)Grude 3.164844e-03
## factor(Munic)I.Ilidza 6.626190e-01
## factor(Munic)Ilidza 4.510947e-06
## factor(Munic)Ilijas 3.685511e-01
## factor(Munic)Jajce 1.411717e-02
## factor(Munic)K.Dubica 7.359346e-01
## factor(Munic)Kakanj 4.049064e-01
## factor(Munic)Kalesija 1.468048e-01
## factor(Munic)Kiseljak 2.486030e-02
## factor(Munic)Kladanj 4.567208e-01
## factor(Munic)Kljuc 6.960884e-01
## factor(Munic)Knezevo 8.931470e-01
## factor(Munic)Konjic 1.324176e-02
## factor(Munic)Kotor Varos 3.367326e-01
## factor(Munic)Kresevo 4.078968e-02
## factor(Munic)Kupres 1.546326e-01
## factor(Munic)Laktasi 5.127553e-01
## factor(Munic)Livno 9.870775e-04
## factor(Munic)Lopare 5.006349e-01
## factor(Munic)Lukavac 8.069982e-02
## factor(Munic)Ljubuski 2.486538e-05
## factor(Munic)Maglaj 2.903266e-01
## factor(Munic)Modrica 6.929755e-01
## factor(Munic)Mrkonjic 1.290928e-01
## factor(Munic)N.G.Sarajevo 4.396238e-04
## factor(Munic)N.Gorazde 1.732563e-01
## factor(Munic)N.Sarajevo 7.361742e-05
## factor(Munic)N.Travnik 8.299787e-01
## factor(Munic)Nevesinje 9.266892e-01
## factor(Munic)Novi Grad 7.382760e-01
## factor(Munic)Orasje 5.087257e-03
## factor(Munic)Pale 4.298819e-01
## factor(Munic)Posusje 3.145855e-03
## factor(Munic)Prijedor 1.208891e-01
## factor(Munic)Prnjavor 2.732803e-01
## factor(Munic)Prozor 8.572176e-02
## factor(Munic)Ribnik 4.285096e-01
## factor(Munic)Rogatica 2.220571e-01
## factor(Munic)Rudo 7.157436e-01
## factor(Munic)S.Brijeg 1.517911e-06
## factor(Munic)S.G.Sarajevo 2.360785e-03
## factor(Munic)Samac 3.583994e-01
## factor(Munic)Sanski Most 4.453401e-02
## factor(Munic)Sokolac 6.370955e-01
## factor(Munic)Srbac 5.030937e-01
## factor(Munic)Srebrenica 1.829359e-01
## factor(Munic)Srebrenik 3.644780e-01
## factor(Munic)Stolac 2.663355e-02
## factor(Munic)Tesanj 5.964706e-01
## factor(Munic)Teslic 6.872761e-01
## factor(Munic)Tomislavgrad 1.351158e-04
## factor(Munic)Travnik 1.746312e-03
## factor(Munic)Trebinje 3.982096e-01
## factor(Munic)Tuzla 9.289972e-01
## factor(Munic)Ugljevik 8.254038e-02
## factor(Munic)V.Kladusa 8.428583e-01
## factor(Munic)Visegrad 7.963551e-01
## factor(Munic)Visoko 1.362546e-01
## factor(Munic)Vitez 7.458233e-01
## factor(Munic)Vlasenica 3.629891e-01
## factor(Munic)Vogosca 2.548653e-01
## factor(Munic)Zavidovici 4.388692e-01
## factor(Munic)Zenica 1.650432e-01
## factor(Munic)Zepce 4.231488e-01
## factor(Munic)Zivinice 1.758620e-01
## factor(Munic)Zvornik 9.353450e-02###WI.3a model
model.coef3a <- summary(WI.3a)$coefficients
model.coef3a # vidi koliko redova su koef.za druge varijable (8redova) i B.Krupa je referentna opstina## Estimate Std. Error t value
## (Intercept) 4.8820174064 1.385390e-01 35.23929966
## Age 0.0410644316 4.938568e-03 8.31504816
## as.numeric(Age^2) -0.0003955788 5.833688e-05 -6.78093837
## factor(Sex_)2 -0.1462233294 1.734289e-02 -8.43131446
## Educ.1elementary 0.1654886912 6.755574e-02 2.44966130
## Educ.1high schl 0.4823290332 6.571441e-02 7.33977534
## Educ.1faculty 0.9292422135 6.791161e-02 13.68311171
## factor(Full.time)part time -0.7410491635 3.014004e-02 -24.58686674
## factor(Munic)B.Petrovac -0.1828358917 1.880676e-01 -0.97218162
## factor(Munic)Banovici 0.3016533953 1.435173e-01 2.10186063
## factor(Munic)Banja Luka 0.1150192556 1.017768e-01 1.13011324
## factor(Munic)Bihac -0.0851187048 1.109684e-01 -0.76705331
## factor(Munic)Bijeljina -0.1147259562 1.011230e-01 -1.13451837
## factor(Munic)Bileca 0.0490915703 1.770603e-01 0.27725901
## factor(Munic)Bratunac 0.3500611443 1.767523e-01 1.98051838
## factor(Munic)Brcko 0.1160503018 1.000743e-01 1.15964185
## factor(Munic)Breza 0.2314714318 1.615308e-01 1.43298605
## factor(Munic)Brod -0.2568216114 1.401847e-01 -1.83202342
## factor(Munic)Bugojno 0.0392603818 1.139786e-01 0.34445398
## factor(Munic)Busovaca 0.0602734344 1.353958e-01 0.44516461
## factor(Munic)Buzim -0.6658529584 1.767191e-01 -3.76786037
## factor(Munic)C.Sarajevo 0.3420448207 1.192847e-01 2.86746574
## factor(Munic)Cajnice -0.1784560884 1.685907e-01 -1.05851665
## factor(Munic)Capljina 0.4760086679 1.399878e-01 3.40035750
## factor(Munic)Cazin -0.1263838686 1.144643e-01 -1.10413341
## factor(Munic)Celinac -0.1167661273 1.880239e-01 -0.62101763
## factor(Munic)Citluk 0.2927611936 1.391035e-01 2.10462823
## factor(Munic)D.Vakuf -0.1137453229 1.258726e-01 -0.90365442
## factor(Munic)Derventa -0.4716639191 1.293800e-01 -3.64556949
## factor(Munic)Doboj 0.1505261578 1.133413e-01 1.32807822
## factor(Munic)Doboj-Jug -0.1086578800 1.949907e-01 -0.55724648
## factor(Munic)Drvar 0.0266899413 1.535047e-01 0.17387054
## factor(Munic)Foca - RS -0.0323988425 1.355439e-01 -0.23902837
## factor(Munic)Fojnica 0.2496832908 1.723464e-01 1.44872951
## factor(Munic)G.Vakuf -0.1004225275 1.326466e-01 -0.75706827
## factor(Munic)Gacko 0.2474449393 1.404322e-01 1.76202422
## factor(Munic)Glamoc 0.2370568654 1.773276e-01 1.33682969
## factor(Munic)Gorazde 0.2352043382 1.288652e-01 1.82519620
## factor(Munic)Gracanica -0.1287612023 1.260310e-01 -1.02166316
## factor(Munic)Grad Mostar 0.3323133995 1.032354e-01 3.21898710
## factor(Munic)Gradacac 0.0714432098 1.237006e-01 0.57754940
## factor(Munic)Gradiska -0.0495593918 1.509196e-01 -0.32838271
## factor(Munic)Grude 0.4015864622 1.354478e-01 2.96487917
## factor(Munic)I.Ilidza -0.0129577300 1.615006e-01 -0.08023330
## factor(Munic)Ilidza 0.5109311837 1.213849e-01 4.20918174
## factor(Munic)Ilijas 0.1765869018 1.646209e-01 1.07268834
## factor(Munic)Jajce 0.2767950564 1.555703e-01 1.77922800
## factor(Munic)K.Dubica -0.1379123331 1.722778e-01 -0.80052296
## factor(Munic)Kakanj 0.0780721742 1.279619e-01 0.61012063
## factor(Munic)Kalesija -0.1258889572 1.436310e-01 -0.87647505
## factor(Munic)Kiseljak 0.2476207118 1.314476e-01 1.88379772
## factor(Munic)Kladanj 0.1331337133 1.767690e-01 0.75315075
## factor(Munic)Kljuc -0.0100623402 1.646436e-01 -0.06111590
## factor(Munic)Knezevo -0.0282840147 1.820768e-01 -0.15534115
## factor(Munic)Konjic 0.2337023837 1.161409e-01 2.01223094
## factor(Munic)Kotor Varos 0.1407311704 1.767694e-01 0.79612880
## factor(Munic)Kresevo 0.3692902596 1.470536e-01 2.51126361
## factor(Munic)Kupres 0.2036694262 1.642973e-01 1.23963968
## factor(Munic)Laktasi 0.0235345614 1.241203e-01 0.18961095
## factor(Munic)Livno 0.4060325482 1.331447e-01 3.04955810
## factor(Munic)Lopare -0.0123681621 1.260554e-01 -0.09811685
## factor(Munic)Lukavac 0.1164707504 1.343077e-01 0.86719366
## factor(Munic)Ljubuski 0.4160130281 1.345058e-01 3.09290034
## factor(Munic)Maglaj 0.0661564022 1.377775e-01 0.48016838
## factor(Munic)Modrica -0.0635310175 1.194994e-01 -0.53164297
## factor(Munic)Mrkonjic 0.2657356976 1.951913e-01 1.36141191
## factor(Munic)N.G.Sarajevo 0.3454167452 1.022333e-01 3.37871201
## factor(Munic)N.Gorazde -0.4082614358 1.949913e-01 -2.09374232
## factor(Munic)N.Sarajevo 0.4298509631 1.193930e-01 3.60030410
## factor(Munic)N.Travnik -0.1673090938 1.125284e-01 -1.48681702
## factor(Munic)Nevesinje -0.0349373647 1.821560e-01 -0.19179918
## factor(Munic)Novi Grad -0.0504538219 1.192999e-01 -0.42291603
## factor(Munic)Orasje 0.3446396132 1.377899e-01 2.50119590
## factor(Munic)Pale 0.0730904827 1.452237e-01 0.50329590
## factor(Munic)Posusje 0.4059994225 1.557602e-01 2.60656645
## factor(Munic)Prijedor 0.1439496058 1.119961e-01 1.28530924
## factor(Munic)Prnjavor 0.0196426833 1.242357e-01 0.15810814
## factor(Munic)Prozor 0.2233714878 1.771515e-01 1.26090636
## factor(Munic)Ribnik -0.0195693403 1.614201e-01 -0.12123238
## factor(Munic)Rogatica 0.1830572232 1.535595e-01 1.19209305
## factor(Munic)Rudo -0.1641534327 1.614626e-01 -1.01666533
## factor(Munic)S.Brijeg 0.6098130173 1.190223e-01 5.12351911
## factor(Munic)S.G.Sarajevo 0.3636204138 1.217854e-01 2.98574672
## factor(Munic)Samac 0.0210317727 1.123612e-01 0.18718009
## factor(Munic)Sanski Most -0.0125448276 1.128992e-01 -0.11111529
## factor(Munic)Sokolac -0.1642980421 1.824111e-01 -0.90070221
## factor(Munic)Srbac 0.0575557836 1.583404e-01 0.36349408
## factor(Munic)Srebrenica 0.3277554536 1.949379e-01 1.68133236
## factor(Munic)Srebrenik 0.1607166511 1.878565e-01 0.85552860
## factor(Munic)Stolac 0.3131976234 1.647674e-01 1.90084684
## factor(Munic)Tesanj 0.0012990188 1.188165e-01 0.01093298
## factor(Munic)Teslic -0.0376860150 1.287540e-01 -0.29269772
## factor(Munic)Tomislavgrad 0.1667336034 1.820880e-01 0.91567581
## factor(Munic)Travnik 0.3328100303 1.056003e-01 3.15160040
## factor(Munic)Trebinje 0.0339378713 1.274372e-01 0.26631046
## factor(Munic)Tuzla 0.1672624447 1.118753e-01 1.49507904
## factor(Munic)Ugljevik 0.1734186954 1.434156e-01 1.20920414
## factor(Munic)V.Kladusa -0.2491723445 1.243046e-01 -2.00453001
## factor(Munic)Visegrad -0.1452941187 1.684744e-01 -0.86241053
## factor(Munic)Visoko 0.1285958545 1.181222e-01 1.08866806
## factor(Munic)Vitez -0.0681471761 1.131371e-01 -0.60234159
## factor(Munic)Vlasenica 0.1152210905 1.582518e-01 0.72808699
## factor(Munic)Vogosca 0.1411544974 1.435292e-01 0.98345500
## factor(Munic)Zavidovici -0.0415413631 1.197557e-01 -0.34688423
## factor(Munic)Zenica 0.1012596903 1.038040e-01 0.97548960
## factor(Munic)Zivinice -0.0541166033 1.208023e-01 -0.44797642
## factor(Munic)Zvornik -0.3264987028 1.163042e-01 -2.80728324
## Pr(>|t|)
## (Intercept) 7.233233e-240
## Age 1.209020e-16
## as.numeric(Age^2) 1.352056e-11
## factor(Sex_)2 4.576044e-17
## Educ.1elementary 1.433774e-02
## Educ.1high schl 2.534592e-13
## Educ.1faculty 9.057415e-42
## factor(Full.time)part time 3.397987e-125
## factor(Munic)B.Petrovac 3.310136e-01
## factor(Munic)Banovici 3.562214e-02
## factor(Munic)Banja Luka 2.584901e-01
## factor(Munic)Bihac 4.430909e-01
## factor(Munic)Bijeljina 2.566389e-01
## factor(Munic)Bileca 7.815942e-01
## factor(Munic)Bratunac 4.770744e-02
## factor(Munic)Brcko 2.462575e-01
## factor(Munic)Breza 1.519328e-01
## factor(Munic)Brod 6.701543e-02
## factor(Munic)Bugojno 7.305213e-01
## factor(Munic)Busovaca 6.562227e-01
## factor(Munic)Buzim 1.668056e-04
## factor(Munic)C.Sarajevo 4.157389e-03
## factor(Munic)Cajnice 2.898780e-01
## factor(Munic)Capljina 6.789631e-04
## factor(Munic)Cazin 2.695956e-01
## factor(Munic)Celinac 5.346202e-01
## factor(Munic)Citluk 3.538020e-02
## factor(Munic)D.Vakuf 3.662281e-01
## factor(Munic)Derventa 2.698804e-04
## factor(Munic)Doboj 1.842211e-01
## factor(Munic)Doboj-Jug 5.773874e-01
## factor(Munic)Drvar 8.619752e-01
## factor(Munic)Foca - RS 8.110947e-01
## factor(Munic)Fojnica 1.474843e-01
## factor(Munic)G.Vakuf 4.490495e-01
## factor(Munic)Gacko 7.813457e-02
## factor(Munic)Glamoc 1.813473e-01
## factor(Munic)Gorazde 6.803913e-02
## factor(Munic)Gracanica 3.069965e-01
## factor(Munic)Grad Mostar 1.295777e-03
## factor(Munic)Gradacac 5.635979e-01
## factor(Munic)Gradiska 7.426379e-01
## factor(Munic)Grude 3.044260e-03
## factor(Munic)I.Ilidza 9.360554e-01
## factor(Munic)Ilidza 2.614068e-05
## factor(Munic)Ilijas 2.834698e-01
## factor(Munic)Jajce 7.527134e-02
## factor(Munic)K.Dubica 4.234510e-01
## factor(Munic)Kakanj 5.418134e-01
## factor(Munic)Kalesija 3.808196e-01
## factor(Munic)Kiseljak 5.965817e-02
## factor(Munic)Kladanj 4.513996e-01
## factor(Munic)Kljuc 9.512697e-01
## factor(Munic)Knezevo 8.765595e-01
## factor(Munic)Konjic 4.425637e-02
## factor(Munic)Kotor Varos 4.260001e-01
## factor(Munic)Kresevo 1.206550e-02
## factor(Munic)Kupres 2.151747e-01
## factor(Munic)Laktasi 8.496227e-01
## factor(Munic)Livno 2.305409e-03
## factor(Munic)Lopare 9.218440e-01
## factor(Munic)Lukavac 3.858832e-01
## factor(Munic)Ljubuski 1.994518e-03
## factor(Munic)Maglaj 6.311315e-01
## factor(Munic)Modrica 5.950001e-01
## factor(Munic)Mrkonjic 1.734532e-01
## factor(Munic)N.G.Sarajevo 7.345870e-04
## factor(Munic)N.Gorazde 3.633998e-02
## factor(Munic)N.Sarajevo 3.213461e-04
## factor(Munic)N.Travnik 1.371347e-01
## factor(Munic)Nevesinje 8.479084e-01
## factor(Munic)Novi Grad 6.723772e-01
## factor(Munic)Orasje 1.241350e-02
## factor(Munic)Pale 6.147814e-01
## factor(Munic)Posusje 9.176328e-03
## factor(Munic)Prijedor 1.987518e-01
## factor(Munic)Prnjavor 8.743788e-01
## factor(Munic)Prozor 2.074094e-01
## factor(Munic)Ribnik 9.035125e-01
## factor(Munic)Rogatica 2.332890e-01
## factor(Munic)Rudo 3.093685e-01
## factor(Munic)S.Brijeg 3.127144e-07
## factor(Munic)S.G.Sarajevo 2.844435e-03
## factor(Munic)Samac 8.515281e-01
## factor(Munic)Sanski Most 9.115300e-01
## factor(Munic)Sokolac 3.677959e-01
## factor(Munic)Srbac 7.162533e-01
## factor(Munic)Srebrenica 9.276928e-02
## factor(Munic)Srebrenik 3.923051e-01
## factor(Munic)Stolac 5.738730e-02
## factor(Munic)Tesanj 9.912774e-01
## factor(Munic)Teslic 7.697670e-01
## factor(Munic)Tomislavgrad 3.598871e-01
## factor(Munic)Travnik 1.634680e-03
## factor(Munic)Trebinje 7.900126e-01
## factor(Munic)Tuzla 1.349655e-01
## factor(Munic)Ugljevik 2.266493e-01
## factor(Munic)V.Kladusa 4.507442e-02
## factor(Munic)Visegrad 3.885086e-01
## factor(Munic)Visoko 2.763599e-01
## factor(Munic)Vitez 5.469778e-01
## factor(Munic)Vlasenica 4.665990e-01
## factor(Munic)Vogosca 3.254376e-01
## factor(Munic)Zavidovici 7.286949e-01
## factor(Munic)Zenica 3.293710e-01
## factor(Munic)Zivinice 6.541922e-01
## factor(Munic)Zvornik 5.018104e-03model.coef3[,c("Estimate","Pr(>|t|)")]## Estimate Pr(>|t|)
## (Intercept) 4.3659883884 1.628195e-145
## Age 0.0562606318 1.005724e-22
## as.numeric(Age^2) -0.0005561833 1.469344e-16
## factor(Sex_)2 -0.1977001270 3.769003e-21
## Educ.1elementary 0.3266644379 2.656310e-05
## Educ.1high schl 0.6098989185 9.210427e-16
## Educ.1faculty 1.0781752145 3.126198e-42
## factor(Full.time)part time -0.8548068085 1.708215e-136
## factor(Munic)B.Petrovac 0.0705671062 7.571689e-01
## factor(Munic)Banovici 0.3308523627 6.090993e-02
## factor(Munic)Banja Luka 0.0730207936 5.541128e-01
## factor(Munic)Bihac -0.0374563860 7.800436e-01
## factor(Munic)Bijeljina -0.0031394567 9.795089e-01
## factor(Munic)Bileca 0.1464476751 5.088348e-01
## factor(Munic)Bratunac 0.3948522744 7.443698e-02
## factor(Munic)Brcko 0.2469952639 4.115614e-02
## factor(Munic)Breza 0.3048600123 1.307926e-01
## factor(Munic)Brod -0.7105596047 8.255629e-06
## factor(Munic)Bugojno 0.0600375326 6.571830e-01
## factor(Munic)Busovaca 0.1880624501 2.550978e-01
## factor(Munic)Buzim -0.7420081828 3.184768e-04
## factor(Munic)C.Sarajevo 0.4377795698 2.120443e-03
## factor(Munic)Cajnice -0.2111280915 3.057261e-01
## factor(Munic)Capljina 0.5350666802 1.907587e-03
## factor(Munic)Cazin -0.1105391303 4.283317e-01
## factor(Munic)Celinac -0.0180275135 9.390522e-01
## factor(Munic)Citluk 0.3185605215 6.039729e-02
## factor(Munic)D.Vakuf -0.0103312781 9.449429e-01
## factor(Munic)Derventa -0.4695186451 2.323237e-03
## factor(Munic)Doboj 0.2389127664 8.450966e-02
## factor(Munic)Doboj-Jug -0.0627222892 7.902859e-01
## factor(Munic)Drvar 0.1469117239 4.426884e-01
## factor(Munic)Foca - RS 0.0745178915 6.577205e-01
## factor(Munic)Fojnica 0.3809556599 7.731492e-02
## factor(Munic)G.Vakuf -0.0110627051 9.430047e-01
## factor(Munic)Gacko 0.3146141433 6.867286e-02
## factor(Munic)Glamoc 0.3423180660 1.229844e-01
## factor(Munic)Gorazde 0.3089617711 5.002955e-02
## factor(Munic)Gracanica -0.0321625146 8.302040e-01
## factor(Munic)Grad Mostar 0.4517069832 3.129832e-04
## factor(Munic)Gradacac 0.0767090454 6.064020e-01
## factor(Munic)Gradiska -0.1547353027 4.037861e-01
## factor(Munic)Grude 0.4921280157 3.164844e-03
## factor(Munic)I.Ilidza 0.0863229995 6.626190e-01
## factor(Munic)Ilidza 0.6617860766 4.510947e-06
## factor(Munic)Ilijas 0.1746241686 3.685511e-01
## factor(Munic)Jajce 0.4614114963 1.411717e-02
## factor(Munic)K.Dubica -0.0709565681 7.359346e-01
## factor(Munic)Kakanj 0.1296385629 4.049064e-01
## factor(Munic)Kalesija -0.2266678449 1.468048e-01
## factor(Munic)Kiseljak 0.3583194267 2.486030e-02
## factor(Munic)Kladanj 0.1563606899 4.567208e-01
## factor(Munic)Kljuc 0.0772279155 6.960884e-01
## factor(Munic)Knezevo 0.0306463233 8.931470e-01
## factor(Munic)Konjic 0.3497787816 1.324176e-02
## factor(Munic)Kotor Varos 0.2126302150 3.367326e-01
## factor(Munic)Kresevo 0.3698792651 4.078968e-02
## factor(Munic)Kupres 0.2924203493 1.546326e-01
## factor(Munic)Laktasi 0.1003922861 5.127553e-01
## factor(Munic)Livno 0.5323425604 9.870775e-04
## factor(Munic)Lopare 0.1023503584 5.006349e-01
## factor(Munic)Lukavac 0.2769705791 8.069982e-02
## factor(Munic)Ljubuski 0.6585852603 2.486538e-05
## factor(Munic)Maglaj 0.1622368012 2.903266e-01
## factor(Munic)Modrica 0.0571851011 6.929755e-01
## factor(Munic)Mrkonjic 0.3583132832 1.290928e-01
## factor(Munic)N.G.Sarajevo 0.4363612589 4.396238e-04
## factor(Munic)N.Gorazde -0.2802858760 1.732563e-01
## factor(Munic)N.Sarajevo 0.5652466009 7.361742e-05
## factor(Munic)N.Travnik -0.0295701104 8.299787e-01
## factor(Munic)Nevesinje 0.0210007713 9.266892e-01
## factor(Munic)Novi Grad 0.0482048374 7.382760e-01
## factor(Munic)Orasje 0.4633726383 5.087257e-03
## factor(Munic)Pale 0.1426644128 4.298819e-01
## factor(Munic)Posusje 0.5560415278 3.145855e-03
## factor(Munic)Prijedor 0.2122313116 1.208891e-01
## factor(Munic)Prnjavor 0.1612972537 2.732803e-01
## factor(Munic)Prozor 0.3467681162 8.572176e-02
## factor(Munic)Ribnik -0.1538460514 4.285096e-01
## factor(Munic)Rogatica 0.2268020544 2.220571e-01
## factor(Munic)Rudo -0.0675427752 7.157436e-01
## factor(Munic)S.Brijeg 0.6994496266 1.517911e-06
## factor(Munic)S.G.Sarajevo 0.4554604551 2.360785e-03
## factor(Munic)Samac 0.1255903694 3.583994e-01
## factor(Munic)Sanski Most 0.2718470039 4.453401e-02
## factor(Munic)Sokolac -0.1078236039 6.370955e-01
## factor(Munic)Srbac 0.1299775432 5.030937e-01
## factor(Munic)Srebrenica 0.3141764735 1.829359e-01
## factor(Munic)Srebrenik 0.2136140785 3.644780e-01
## factor(Munic)Stolac 0.4472067254 2.663355e-02
## factor(Munic)Tesanj 0.0737589190 5.964706e-01
## factor(Munic)Teslic 0.0611353025 6.872761e-01
## factor(Munic)Tomislavgrad 0.7188404067 1.351158e-04
## factor(Munic)Travnik 0.4026913177 1.746312e-03
## factor(Munic)Trebinje 0.1324997382 3.982096e-01
## factor(Munic)Tuzla 0.0118487274 9.289972e-01
## factor(Munic)Ugljevik 0.3097543410 8.254038e-02
## factor(Munic)V.Kladusa 0.0284219244 8.428583e-01
## factor(Munic)Visegrad 0.0510917133 7.963551e-01
## factor(Munic)Visoko 0.2155745610 1.362546e-01
## factor(Munic)Vitez 0.0446100308 7.458233e-01
## factor(Munic)Vlasenica 0.1797038889 3.629891e-01
## factor(Munic)Vogosca 0.2009559193 2.548653e-01
## factor(Munic)Zavidovici 0.1088124342 4.388692e-01
## factor(Munic)Zenica 0.1751121675 1.650432e-01
## factor(Munic)Zepce 0.1726908697 4.231488e-01
## factor(Munic)Zivinice -0.1905094734 1.758620e-01
## factor(Munic)Zvornik -0.2375947797 9.353450e-02model.coef3 <- model.coef3
model.coef3 <- as.data.frame.array(model.coef3[-c(1:8),]) # da bude dataframe a ne matrix i da mi brise prvih 8 redova koji su drugi koeficijenti
Munic <- rownames(model.coef3)
model.coef3$Munic <- Municdim(model.coef3) #da vidim koliko redova## [1] 100 5rownames(model.coef3) <- 1:100 #nedostaje B.Krupa koje je 0, B.Grahovo nije u modelu
head(model.coef3)## Estimate Std. Error t value Pr(>|t|) Munic
## 1 0.070567106 0.2282118 0.30921761 0.75716892 factor(Munic)B.Petrovac
## 2 0.330852363 0.1764948 1.87457278 0.06090993 factor(Munic)Banovici
## 3 0.073020794 0.1234192 0.59164846 0.55411283 factor(Munic)Banja Luka
## 4 -0.037456386 0.1341190 -0.27927731 0.78004360 factor(Munic)Bihac
## 5 -0.003139457 0.1222250 -0.02568588 0.97950892 factor(Munic)Bijeljina
## 6 0.146447675 0.2216551 0.66070068 0.50883480 factor(Munic)BilecaMake a new variable with municipal names without
factor(Munic)
Munic1 <- sapply(strsplit(Munic, split=')', fixed=TRUE), function(x) (x[2]))
Munic1## [1] "B.Petrovac" "Banovici" "Banja Luka" "Bihac" "Bijeljina"
## [6] "Bileca" "Bratunac" "Brcko" "Breza" "Brod"
## [11] "Bugojno" "Busovaca" "Buzim" "C.Sarajevo" "Cajnice"
## [16] "Capljina" "Cazin" "Celinac" "Citluk" "D.Vakuf"
## [21] "Derventa" "Doboj" "Doboj-Jug" "Drvar" "Foca - RS"
## [26] "Fojnica" "G.Vakuf" "Gacko" "Glamoc" "Gorazde"
## [31] "Gracanica" "Grad Mostar" "Gradacac" "Gradiska" "Grude"
## [36] "I.Ilidza" "Ilidza" "Ilijas" "Jajce" "K.Dubica"
## [41] "Kakanj" "Kalesija" "Kiseljak" "Kladanj" "Kljuc"
## [46] "Knezevo" "Konjic" "Kotor Varos" "Kresevo" "Kupres"
## [51] "Laktasi" "Livno" "Lopare" "Lukavac" "Ljubuski"
## [56] "Maglaj" "Modrica" "Mrkonjic" "N.G.Sarajevo" "N.Gorazde"
## [61] "N.Sarajevo" "N.Travnik" "Nevesinje" "Novi Grad" "Orasje"
## [66] "Pale" "Posusje" "Prijedor" "Prnjavor" "Prozor"
## [71] "Ribnik" "Rogatica" "Rudo" "S.Brijeg" "S.G.Sarajevo"
## [76] "Samac" "Sanski Most" "Sokolac" "Srbac" "Srebrenica"
## [81] "Srebrenik" "Stolac" "Tesanj" "Teslic" "Tomislavgrad"
## [86] "Travnik" "Trebinje" "Tuzla" "Ugljevik" "V.Kladusa"
## [91] "Visegrad" "Visoko" "Vitez" "Vlasenica" "Vogosca"
## [96] "Zavidovici" "Zenica" "Zepce" "Zivinice" "Zvornik"model.coef3$Munic1 <- Munic1sada za model 3a
model.coef3a[,c("Estimate","Pr(>|t|)")]## Estimate Pr(>|t|)
## (Intercept) 4.8820174064 7.233233e-240
## Age 0.0410644316 1.209020e-16
## as.numeric(Age^2) -0.0003955788 1.352056e-11
## factor(Sex_)2 -0.1462233294 4.576044e-17
## Educ.1elementary 0.1654886912 1.433774e-02
## Educ.1high schl 0.4823290332 2.534592e-13
## Educ.1faculty 0.9292422135 9.057415e-42
## factor(Full.time)part time -0.7410491635 3.397987e-125
## factor(Munic)B.Petrovac -0.1828358917 3.310136e-01
## factor(Munic)Banovici 0.3016533953 3.562214e-02
## factor(Munic)Banja Luka 0.1150192556 2.584901e-01
## factor(Munic)Bihac -0.0851187048 4.430909e-01
## factor(Munic)Bijeljina -0.1147259562 2.566389e-01
## factor(Munic)Bileca 0.0490915703 7.815942e-01
## factor(Munic)Bratunac 0.3500611443 4.770744e-02
## factor(Munic)Brcko 0.1160503018 2.462575e-01
## factor(Munic)Breza 0.2314714318 1.519328e-01
## factor(Munic)Brod -0.2568216114 6.701543e-02
## factor(Munic)Bugojno 0.0392603818 7.305213e-01
## factor(Munic)Busovaca 0.0602734344 6.562227e-01
## factor(Munic)Buzim -0.6658529584 1.668056e-04
## factor(Munic)C.Sarajevo 0.3420448207 4.157389e-03
## factor(Munic)Cajnice -0.1784560884 2.898780e-01
## factor(Munic)Capljina 0.4760086679 6.789631e-04
## factor(Munic)Cazin -0.1263838686 2.695956e-01
## factor(Munic)Celinac -0.1167661273 5.346202e-01
## factor(Munic)Citluk 0.2927611936 3.538020e-02
## factor(Munic)D.Vakuf -0.1137453229 3.662281e-01
## factor(Munic)Derventa -0.4716639191 2.698804e-04
## factor(Munic)Doboj 0.1505261578 1.842211e-01
## factor(Munic)Doboj-Jug -0.1086578800 5.773874e-01
## factor(Munic)Drvar 0.0266899413 8.619752e-01
## factor(Munic)Foca - RS -0.0323988425 8.110947e-01
## factor(Munic)Fojnica 0.2496832908 1.474843e-01
## factor(Munic)G.Vakuf -0.1004225275 4.490495e-01
## factor(Munic)Gacko 0.2474449393 7.813457e-02
## factor(Munic)Glamoc 0.2370568654 1.813473e-01
## factor(Munic)Gorazde 0.2352043382 6.803913e-02
## factor(Munic)Gracanica -0.1287612023 3.069965e-01
## factor(Munic)Grad Mostar 0.3323133995 1.295777e-03
## factor(Munic)Gradacac 0.0714432098 5.635979e-01
## factor(Munic)Gradiska -0.0495593918 7.426379e-01
## factor(Munic)Grude 0.4015864622 3.044260e-03
## factor(Munic)I.Ilidza -0.0129577300 9.360554e-01
## factor(Munic)Ilidza 0.5109311837 2.614068e-05
## factor(Munic)Ilijas 0.1765869018 2.834698e-01
## factor(Munic)Jajce 0.2767950564 7.527134e-02
## factor(Munic)K.Dubica -0.1379123331 4.234510e-01
## factor(Munic)Kakanj 0.0780721742 5.418134e-01
## factor(Munic)Kalesija -0.1258889572 3.808196e-01
## factor(Munic)Kiseljak 0.2476207118 5.965817e-02
## factor(Munic)Kladanj 0.1331337133 4.513996e-01
## factor(Munic)Kljuc -0.0100623402 9.512697e-01
## factor(Munic)Knezevo -0.0282840147 8.765595e-01
## factor(Munic)Konjic 0.2337023837 4.425637e-02
## factor(Munic)Kotor Varos 0.1407311704 4.260001e-01
## factor(Munic)Kresevo 0.3692902596 1.206550e-02
## factor(Munic)Kupres 0.2036694262 2.151747e-01
## factor(Munic)Laktasi 0.0235345614 8.496227e-01
## factor(Munic)Livno 0.4060325482 2.305409e-03
## factor(Munic)Lopare -0.0123681621 9.218440e-01
## factor(Munic)Lukavac 0.1164707504 3.858832e-01
## factor(Munic)Ljubuski 0.4160130281 1.994518e-03
## factor(Munic)Maglaj 0.0661564022 6.311315e-01
## factor(Munic)Modrica -0.0635310175 5.950001e-01
## factor(Munic)Mrkonjic 0.2657356976 1.734532e-01
## factor(Munic)N.G.Sarajevo 0.3454167452 7.345870e-04
## factor(Munic)N.Gorazde -0.4082614358 3.633998e-02
## factor(Munic)N.Sarajevo 0.4298509631 3.213461e-04
## factor(Munic)N.Travnik -0.1673090938 1.371347e-01
## factor(Munic)Nevesinje -0.0349373647 8.479084e-01
## factor(Munic)Novi Grad -0.0504538219 6.723772e-01
## factor(Munic)Orasje 0.3446396132 1.241350e-02
## factor(Munic)Pale 0.0730904827 6.147814e-01
## factor(Munic)Posusje 0.4059994225 9.176328e-03
## factor(Munic)Prijedor 0.1439496058 1.987518e-01
## factor(Munic)Prnjavor 0.0196426833 8.743788e-01
## factor(Munic)Prozor 0.2233714878 2.074094e-01
## factor(Munic)Ribnik -0.0195693403 9.035125e-01
## factor(Munic)Rogatica 0.1830572232 2.332890e-01
## factor(Munic)Rudo -0.1641534327 3.093685e-01
## factor(Munic)S.Brijeg 0.6098130173 3.127144e-07
## factor(Munic)S.G.Sarajevo 0.3636204138 2.844435e-03
## factor(Munic)Samac 0.0210317727 8.515281e-01
## factor(Munic)Sanski Most -0.0125448276 9.115300e-01
## factor(Munic)Sokolac -0.1642980421 3.677959e-01
## factor(Munic)Srbac 0.0575557836 7.162533e-01
## factor(Munic)Srebrenica 0.3277554536 9.276928e-02
## factor(Munic)Srebrenik 0.1607166511 3.923051e-01
## factor(Munic)Stolac 0.3131976234 5.738730e-02
## factor(Munic)Tesanj 0.0012990188 9.912774e-01
## factor(Munic)Teslic -0.0376860150 7.697670e-01
## factor(Munic)Tomislavgrad 0.1667336034 3.598871e-01
## factor(Munic)Travnik 0.3328100303 1.634680e-03
## factor(Munic)Trebinje 0.0339378713 7.900126e-01
## factor(Munic)Tuzla 0.1672624447 1.349655e-01
## factor(Munic)Ugljevik 0.1734186954 2.266493e-01
## factor(Munic)V.Kladusa -0.2491723445 4.507442e-02
## factor(Munic)Visegrad -0.1452941187 3.885086e-01
## factor(Munic)Visoko 0.1285958545 2.763599e-01
## factor(Munic)Vitez -0.0681471761 5.469778e-01
## factor(Munic)Vlasenica 0.1152210905 4.665990e-01
## factor(Munic)Vogosca 0.1411544974 3.254376e-01
## factor(Munic)Zavidovici -0.0415413631 7.286949e-01
## factor(Munic)Zenica 0.1012596903 3.293710e-01
## factor(Munic)Zivinice -0.0541166033 6.541922e-01
## factor(Munic)Zvornik -0.3264987028 5.018104e-03model.coef3a <- model.coef3a
model.coef3a <- as.data.frame.array(model.coef3a[-c(1:8),]) # da bude dataframe a ne matrix i da mi brise prvih 8 redova koji su drugi koeficijenti
Munic <- rownames(model.coef3a)
model.coef3a$Munic <- Municza 3a
dim(model.coef3a) #da vidim koliko redova## [1] 99 5rownames(model.coef3a) <- 1:99 #nedostaje B.Krupa koje je 0, B.Grahovo nije u modelu
head(model.coef3a)## Estimate Std. Error t value Pr(>|t|) Munic
## 1 -0.18283589 0.1880676 -0.9721816 0.33101360 factor(Munic)B.Petrovac
## 2 0.30165340 0.1435173 2.1018606 0.03562214 factor(Munic)Banovici
## 3 0.11501926 0.1017768 1.1301132 0.25849005 factor(Munic)Banja Luka
## 4 -0.08511870 0.1109684 -0.7670533 0.44309094 factor(Munic)Bihac
## 5 -0.11472596 0.1011230 -1.1345184 0.25663892 factor(Munic)Bijeljina
## 6 0.04909157 0.1770603 0.2772590 0.78159424 factor(Munic)Bilecamodela
Munic1 <- sapply(strsplit(Munic, split=')', fixed=TRUE), function(x) (x[2]))
Munic1## [1] "B.Petrovac" "Banovici" "Banja Luka" "Bihac" "Bijeljina"
## [6] "Bileca" "Bratunac" "Brcko" "Breza" "Brod"
## [11] "Bugojno" "Busovaca" "Buzim" "C.Sarajevo" "Cajnice"
## [16] "Capljina" "Cazin" "Celinac" "Citluk" "D.Vakuf"
## [21] "Derventa" "Doboj" "Doboj-Jug" "Drvar" "Foca - RS"
## [26] "Fojnica" "G.Vakuf" "Gacko" "Glamoc" "Gorazde"
## [31] "Gracanica" "Grad Mostar" "Gradacac" "Gradiska" "Grude"
## [36] "I.Ilidza" "Ilidza" "Ilijas" "Jajce" "K.Dubica"
## [41] "Kakanj" "Kalesija" "Kiseljak" "Kladanj" "Kljuc"
## [46] "Knezevo" "Konjic" "Kotor Varos" "Kresevo" "Kupres"
## [51] "Laktasi" "Livno" "Lopare" "Lukavac" "Ljubuski"
## [56] "Maglaj" "Modrica" "Mrkonjic" "N.G.Sarajevo" "N.Gorazde"
## [61] "N.Sarajevo" "N.Travnik" "Nevesinje" "Novi Grad" "Orasje"
## [66] "Pale" "Posusje" "Prijedor" "Prnjavor" "Prozor"
## [71] "Ribnik" "Rogatica" "Rudo" "S.Brijeg" "S.G.Sarajevo"
## [76] "Samac" "Sanski Most" "Sokolac" "Srbac" "Srebrenica"
## [81] "Srebrenik" "Stolac" "Tesanj" "Teslic" "Tomislavgrad"
## [86] "Travnik" "Trebinje" "Tuzla" "Ugljevik" "V.Kladusa"
## [91] "Visegrad" "Visoko" "Vitez" "Vlasenica" "Vogosca"
## [96] "Zavidovici" "Zenica" "Zivinice" "Zvornik"model.coef3a$Munic1 <- Munic1add B.Krupa to dataframe
model.coef3[nrow(model.coef3) + 1,] = c(0, NA, NA, NA,"B.Krupa","B.Krupa" )
model.coef3$Estimate <- as.numeric(model.coef3$Estimate)add B.Krupa to dataframe 3a
model.coef3a[nrow(model.coef3a) + 1,] = c(0, NA, NA, NA,"B.Krupa","B.Krupa" )
model.coef3a$Estimate <- as.numeric(model.coef3a$Estimate)top10.m3 <- model.coef3 %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
head(10) #top ten municipalities
top10.m3## Estimate Munic1
## 1 0.7188404 Tomislavgrad
## 2 0.6994496 S.Brijeg
## 3 0.6617861 Ilidza
## 4 0.6585853 Ljubuski
## 5 0.5652466 N.Sarajevo
## 6 0.5560415 Posusje
## 7 0.5350667 Capljina
## 8 0.5323426 Livno
## 9 0.4921280 Grude
## 10 0.4633726 Orasjetop10.m3a <- model.coef3a %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
head(10) #top ten municipalities
top10.m3a## Estimate Munic1
## 1 0.6098130 S.Brijeg
## 2 0.5109312 Ilidza
## 3 0.4760087 Capljina
## 4 0.4298510 N.Sarajevo
## 5 0.4160130 Ljubuski
## 6 0.4060325 Livno
## 7 0.4059994 Posusje
## 8 0.4015865 Grude
## 9 0.3692903 Kresevo
## 10 0.3636204 S.G.Sarajevo#spasiti a u excel
library(openxlsx)
wb = createWorkbook()
sh1 = addWorksheet(wb, "Estimates_allmj") #ovo je naziv sheeta
sh2 = addWorksheet(wb, "mean.wages.onlyempl")#naziv drugog sheeta
sh3 = addWorksheet(wb, "Estimates_notallmj") #ovo je naziv 3. sheeta
sh4 = addWorksheet(wb, "mean.wages.onlyempl_reducwage")#naziv 4. sheeta
library(expss)
xl_write(model.coef3, wb, sh1) #shares je naziv objekta koji hocu spasiti
xl_write(mean.wages, wb, sh2)
xl_write(model.coef3a, wb, sh3) #shares je naziv objekta koji hocu spasiti
xl_write(mean.wages.a, wb, sh4)
saveWorkbook(wb, "Mean_Estim_wages11.04.2024.xlsx", overwrite = T)##Make a graf
merge estimates data and map data
WI.3 <- lm (log(wages)~ Age + as.numeric(Age^2)+ Sex_ + Educ.1 + factor(Full.time) + factor (Munic) , data = WI.df)
map.coef3 <- merge(amraGIS, model.coef3, by.y='Munic1', by.x = "NAME_4",all = TRUE) #100 munic ako nema argumenta all=T, sa njim je 144- sad stoji 142 vidi da nije usora ili slicnoWI.3a <- lm (log(wages)~ Age + as.numeric(Age^2)+ Sex_ + Educ.1 + factor(Full.time) + factor (Munic) , data = WI.dfa)
map.coef3a <- merge(amraGIS, model.coef3a, by.y='Munic1', by.x = "NAME_4",all = TRUE) #100 munic ako nema argumenta all=T, sa njim je 144- sad stoji 142 vidi da nije usora ili slicnonow WI.3 (dependent log(wages))
library(gridExtra)# radi funkcije tableGrob##
## Attaching package: 'gridExtra'## The following object is masked from 'package:dplyr':
##
## combinelibrary(patchwork) #radi spajanja mape i tabele##
## Attaching package: 'patchwork'## The following object is masked from 'package:raster':
##
## arealibrary(ggplot2)
g3 = ggplot() + geom_sf(data = map.coef3, aes(geometry=geometry, fill = Estimate)) + scale_fill_viridis_c() + scale_x_continuous()
g3
#grafikon ispod spaja mi mapu i top 10 municipalities
TABLE <- tableGrob(top10.m3, theme = ttheme_default(base_size = 5))
g3 + ggtitle("lm (log(wages)~ Age + as.numeric(Age^2)+ Sex_ + Educ.1 \n + factor(Full.time) + factor (Munic)") + theme(plot.title = element_text(color="gray", size=12, face="italic")) + plot_spacer() + TABLE + plot_layout(widths = c(12, -1.3 ,5)) #plot_layout prvi broj oznacava velicinu mape, drugi velicinu razmaka izmedju mape i tabele i treci oznacava prostor za tabelu
g3 + geom_sf_text(data= map.coef3, aes(geometry=geometry,label=ifelse(as.numeric(Estimate) >= 0.6,as.character(NAME_4),"")),hjust=0.5,vjust=0,cex =3, col="black") + geom_sf_text(data= map.coef3, aes(geometry=geometry,label=ifelse(as.numeric(Estimate) < -0.4,as.character(NAME_4),"")),hjust=0.5,vjust=0,cex =3, col="red")## Warning in st_point_on_surface.sfc(sf::st_zm(x)): st_point_on_surface may not
## give correct results for longitude/latitude data
## Warning in st_point_on_surface.sfc(sf::st_zm(x)): st_point_on_surface may not
## give correct results for longitude/latitude data## Warning: Removed 41 rows containing missing values (`geom_text()`).
## Removed 41 rows containing missing values (`geom_text()`).
##import new shape file!!
amraGIS.nova <- st_read ("C:/Users/Amra/Documents/HBS database inspecting/geoBoundaries-BIH-ADM2-all/geoBoundaries-BIH-ADM2_simplified.shp")## Reading layer `geoBoundaries-BIH-ADM2_simplified' from data source
## `C:\Users\Amra\Documents\HBS database inspecting\geoBoundaries-BIH-ADM2-all\geoBoundaries-BIH-ADM2_simplified.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 142 features and 5 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 15.73111 ymin: 42.54831 xmax: 19.62407 ymax: 45.27963
## Geodetic CRS: WGS 84##Adjust the names of the municipalites
amraGIS.nova$shapeName <- c("Brcko", "Vukosavlje","Banja Luka", "Bihac", "Bijeljina", "Doboj", "Jajce", "Mostar", "Prijedor", "Trebinje", "Zenica", "Banovici", "Berkovici" , "Bileca","Gradiska", "Kostajnica", "B.Krupa" , "K.Dubica", "Brod", "Novi Grad", "B.Petrovac", "Samac", "B.Grahovo", "Breza" , "Bugojno" , "Busovaca", "Buzim" , "Cajnice", "Capljina", "Cazin", "Celic", "Celinac", "C.Sarajevo", "Citluk", "Derventa", "Doboj-Istok", "Doboj-Jug", "Dobretici", "D-Samac", "D.Vakuf", "Donji Zabar", "Drvar", "I.Drvar", "I.Ilidza", "I.Mostar", "I.N.Sarajevo", "I.S.Grad", "Foca - RS", "Fojnica", "Gacko", "Glamoc","Gorazde", "G.Vakuf", "Gracanica", "Gradacac", "Grude", "Hadzici", "Han Pijesak", "Ilidza", "Ilijas", "Jablanica", "Jezero", "Kakanj", "Kalesija", "Kalinovik", "Kiseljak", "Kladanj", "Kljuc", "Konjic", "Kotor Varos", "Kresevo", "Krupa na Uni", "Kupres-RS", "Kupres", "Laktasi", "Livno", "Ljubinje", "Ljubuski", "Lopare", "Lukavac", "Maglaj", "Milici", "Modrica", "Mrkonjic" ,"Neum", "Nevesinje", "N.Gorazde", "N.Sarajevo", "Novi Grad", "N.Travnik", "Odzak", "S.G.Sarajevo", "Olovo", "Orasje", "Osmaci", "Ostra Luka", "Pale","Pale-FBiH", "Pelagicevo", "Petrovac", "Petrovo", "Posusje", "Prnjavor", "Prozor", "Ravno", "Ribnik", "Rogatica", "Rudo", "Sanski Most", "Sapna","Sekovici", "Sipovo", "S.Brijeg", "Knezevo", "Sokolac", "Srbac", "Srebrenica", "Srebrenik", "Stolac", "Teocak", "Tesanj", "Teslic", "Tomislavgrad", "Travnik", "Trnovo", "Tuzla", "Ugljevik", "Usora", "Foca", "Vares","V.Kladusa", "Visegrad", "Visoko", "Vitez", "Vlasenica", "Vogosca", "Zavidovici", "Zepce", "Zivinice", "Zvornik", "Bratunac", "Trnovo-RS") Make a model with amraGIS.nova
map.coef3.nova <- merge(amraGIS.nova, model.coef3, by.y='Munic1', by.x = "shapeName",all = TRUE)map.coef3a.nova <- merge(amraGIS.nova, model.coef3a, by.y='Munic1', by.x = "shapeName",all = TRUE)Map with nova i exp.values
library(ggplot2)
g5 = ggplot() + geom_sf(data = map.coef3.nova, aes(geometry=geometry, fill = Estimate)) + scale_fill_viridis_c() + scale_x_continuous()
g5 + ggtitle("Wages Index") + theme(plot.title = element_text(color="gray", size=12, face="italic"))
library(ggplot2)
g5a = ggplot() + geom_sf(data = map.coef3a.nova, aes(geometry=geometry, fill = Estimate)) + scale_fill_viridis_c() + scale_x_continuous()
g5a + ggtitle("Wages(notall) Index") + theme(plot.title = element_text(color="gray", size=12, face="italic"))