#Olx
#15.2. Ismir's mail to put const.year and withouth and numb of obs.are robusness check
setwd("C:/Users/Amra/Documents/OLX R inspecting")
library(readxl)
Amra_Fetahović_doktorski_rad_OLX_data <- read_excel("Amra Fetahović - doktorski rad OLX data.xlsx", sheet = "No duplicates")
olx <- Amra_Fetahović_doktorski_rad_OLX_data #32508
head (olx)
## # A tibble: 6 × 20
## ID date_of_Sale prop_type sqm text Okucnica
## <dbl> <dttm> <chr> <chr> <chr> <chr>
## 1 29699 2013-10-18 00:00:00 Kuće 300 300 <NA>
## 2 68498 2014-11-16 00:00:00 Kuće 80 80 <NA>
## 3 118136 2016-02-10 00:00:00 Stanovi 85.1 85.1 <NA>
## 4 118137 2016-02-10 00:00:00 Kuće 285 285 300
## 5 118138 2016-02-10 00:00:00 Stanovi 82 82 <NA>
## 6 118140 2016-02-10 00:00:00 Kuće 200 200 m2 <NA>
## # ℹ 14 more variables: `Kvadratura okućnice` <chr>, Sprat <chr>, Adresa <chr>,
## # Grad <chr>, Grad1 <chr>, Cijena <dbl>, Cijena_original <dbl>,
## # prop.age <dbl>, `Godina izgradnje` <dbl>, `Godina izgradnje1` <chr>,
## # `Nedavno adaptiran` <chr>, Balkon <chr>, `Vrsta grijanja` <chr>, Lift <chr>
str(olx)
## tibble [32,508 × 20] (S3: tbl_df/tbl/data.frame)
## $ ID : num [1:32508] 29699 68498 118136 118137 118138 ...
## $ date_of_Sale : POSIXct[1:32508], format: "2013-10-18" "2014-11-16" ...
## $ prop_type : chr [1:32508] "Kuće" "Kuće" "Stanovi" "Kuće" ...
## $ sqm : chr [1:32508] "300" "80" "85.1" "285" ...
## $ text : chr [1:32508] "300" "80" "85.1" "285" ...
## $ Okucnica : chr [1:32508] NA NA NA "300" ...
## $ Kvadratura okućnice: chr [1:32508] NA NA NA "300" ...
## $ Sprat : chr [1:32508] NA NA "1" NA ...
## $ Adresa : chr [1:32508] NA "Stara Ilićka" "Trg Ilidžanske brigade" "Oteška" ...
## $ Grad : chr [1:32508] "Bihać" "Mostar" "Istočna Ilidža" "Ilidža" ...
## $ Grad1 : chr [1:32508] "Bihać" "Mostar" "Istočna Ilidža" "Ilidža" ...
## $ Cijena : num [1:32508] 400000 85000 102000 330000 240000 240000 90000 86000 165000 80000 ...
## $ Cijena_original : num [1:32508] 400000 85000 102000 330000 240000 240000 90000 86000 165000 80000 ...
## $ prop.age : num [1:32508] NA NA NA 32 NA NA NA NA 42 32 ...
## $ Godina izgradnje : num [1:32508] NA NA NA 1985 NA ...
## $ Godina izgradnje1 : chr [1:32508] NA NA NA "1980 do 1989" ...
## $ Nedavno adaptiran : chr [1:32508] NA NA "Da" "Da" ...
## $ Balkon : chr [1:32508] NA NA "Da" "Da" ...
## $ Vrsta grijanja : chr [1:32508] NA NA "Struja" "Centralno (Kotlovnica)" ...
## $ Lift : chr [1:32508] NA NA NA NA ...
# format variables
olx$sqm <- as.numeric(olx$sqm)
## Warning: NAs introduced by coercion
summary(olx$sqm) #1683NA
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.1 53.0 70.0 159.2 120.0 68685.0 1683
olx$Cijena <- as.numeric(olx$Cijena)
summary(olx$Cijena)#51 NA
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 10000 64000 91800 124616 140000 5900000 51
unique(sort(olx$Grad)) #136 opstina
## [1] "Banovići" "Banja Luka"
## [3] "Berkovići" "Bihać"
## [5] "Bijeljina" "Bileća"
## [7] "Bjelasnica" "Bosanska Krupa"
## [9] "Bosanski Petrovac" "Bosansko Grahovo"
## [11] "Bratunac" "Brčko"
## [13] "Breza" "Brod"
## [15] "Bugojno" "Busovača"
## [17] "Bužim" "Cazin"
## [19] "Čajniče" "Čapljina"
## [21] "Čelić" "Čelinac"
## [23] "Čitluk" "Derventa"
## [25] "Doboj" "Doboj Istok"
## [27] "Doboj Jug" "Dobretići"
## [29] "Domaljevac-Šamac" "Donji Vakuf"
## [31] "Donji Žabar" "Drvar"
## [33] "Foča" "Foča - Ustikolina"
## [35] "Fojnica" "Gacko"
## [37] "Glamoč" "Goražde"
## [39] "Gornji Vakuf-Uskoplje" "Gračanica"
## [41] "Gradačac" "Gradiška"
## [43] "Grude" "Hadžići"
## [45] "Han Pijesak" "Ilidža"
## [47] "Ilijaš" "Istočna Ilidža"
## [49] "Istočni Stari Grad" "Istočno Novo Sarajevo"
## [51] "Jablanica" "Jajce"
## [53] "Jezero" "Kakanj"
## [55] "Kalesija" "Kalinovik"
## [57] "Kiseljak" "Kladanj"
## [59] "Ključ" "Kneževo"
## [61] "Konjic" "Kostajnica"
## [63] "Kotor Varoš" "Kozarska Dubica"
## [65] "Kreševo" "Kupres"
## [67] "Laktaši" "Livno"
## [69] "Lopare" "Lukavac"
## [71] "Ljubinje" "Ljubuški"
## [73] "Maglaj" "Milići"
## [75] "Modriča" "Mostar"
## [77] "Mrkonjic Grad" "Neum"
## [79] "Nevesinje" "Novi Grad"
## [81] "Novi Travnik" "Odžak"
## [83] "Olovo" "Orašje"
## [85] "Osmaci" "Oštra Luka"
## [87] "Pale - FBiH" "Pale - RS"
## [89] "Pelagićevo" "Petrovo"
## [91] "Posušje" "Prijedor"
## [93] "Prnjavor" "Prozor"
## [95] "Ravno" "Rogatica"
## [97] "Rudo" "Sanski Most"
## [99] "Sapna" "Sarajevo - Centar"
## [101] "Sarajevo - Novi Grad" "Sarajevo - Novo Sarajevo"
## [103] "Sarajevo - Stari Grad" "Sokolac"
## [105] "Srbac" "Srebrenica"
## [107] "Srebrenik" "Stolac"
## [109] "Šamac" "Šekovići"
## [111] "Šipovo" "Široki Brijeg"
## [113] "Teočak" "Teslić"
## [115] "Tešanj" "Tomislavgrad"
## [117] "Travnik" "Trebinje"
## [119] "Trnovo - FBiH" "Trnovo - RS"
## [121] "Tuzla" "Ugljevik"
## [123] "Usora" "Vareš"
## [125] "Velika Kladuša" "Visoko"
## [127] "Višegrad" "Vitez"
## [129] "Vlasenica" "Vogošća"
## [131] "Vukosavlje" "Zavidovići"
## [133] "Zenica" "Zvornik"
## [135] "Žepče" "Živinice"
#check num of entry per municipality
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, union
olx %>%
group_by(Grad) %>%
summarise(n=n()) %>%
arrange(n) %>%
print(n=35) #33 less then 10
## # A tibble: 137 × 2
## Grad n
## <chr> <int>
## 1 Berkovići 1
## 2 Bjelasnica 1
## 3 Dobretići 1
## 4 Foča - Ustikolina 1
## 5 Sapna 1
## 6 Čajniče 1
## 7 Šekovići 1
## 8 Domaljevac-Šamac 2
## 9 Jezero 2
## 10 Ljubinje 2
## 11 Oštra Luka 2
## 12 Posušje 2
## 13 Rudo 2
## 14 Teočak 2
## 15 Donji Žabar 3
## 16 Istočni Stari Grad 3
## 17 Kalinovik 3
## 18 Mrkonjic Grad 3
## 19 Pale - FBiH 3
## 20 Glamoč 4
## 21 Han Pijesak 4
## 22 Trnovo - RS 4
## 23 Bosansko Grahovo 5
## 24 Vukosavlje 5
## 25 Doboj Istok 6
## 26 Grude 6
## 27 Osmaci 6
## 28 Čitluk 6
## 29 Kreševo 7
## 30 Pelagićevo 7
## 31 Šipovo 7
## 32 Kostajnica 8
## 33 Milići 9
## 34 Petrovo 10
## 35 Bileća 11
## # ℹ 102 more rows
# ommit NAs in Kvadratura and Grad and filter to Stanovi i Kuce
library(tidyr)
olx1 <- olx %>% drop_na(sqm)
olx1 <- olx1 %>% drop_na(Grad) #30110 obs
table (is.na(olx1$sqm))
##
## FALSE
## 30110
table (is.na(olx1$Grad))
##
## FALSE
## 30110
olx1 <- olx1 %>% dplyr::filter(prop_type=="Kuće"| prop_type=="Stanovi")
table (olx1$prop_type)
##
## Kuće Stanovi
## 11472 18628
# Analiza olx1
table(sort(olx1$Grad))
##
## Banovići Banja Luka Berkovići
## 145 2392 1
## Bihać Bijeljina Bileća
## 493 936 11
## Bjelasnica Bosanska Krupa Bosanski Petrovac
## 1 119 35
## Bosansko Grahovo Bratunac Brčko
## 4 42 1533
## Breza Brod Bugojno
## 80 82 163
## Busovača Bužim Cazin
## 69 13 124
## Čajniče Čapljina Čelić
## 1 77 26
## Čelinac Čitluk Derventa
## 26 6 113
## Doboj Doboj Istok Doboj Jug
## 473 6 19
## Dobretići Domaljevac-Šamac Donji Vakuf
## 1 2 49
## Donji Žabar Drvar Foča
## 3 12 76
## Foča - Ustikolina Fojnica Gacko
## 1 45 15
## Glamoč Goražde Gornji Vakuf-Uskoplje
## 4 155 30
## Gračanica Gradačac Gradiška
## 150 118 210
## Grude Hadžići Han Pijesak
## 5 438 4
## Ilidža Ilijaš Istočna Ilidža
## 2057 238 147
## Istočni Stari Grad Istočno Novo Sarajevo Jablanica
## 3 245 41
## Jajce Jezero Kakanj
## 87 2 147
## Kalesija Kalinovik Kiseljak
## 71 3 113
## Kladanj Ključ Kneževo
## 16 37 13
## Konjic Kostajnica Kotor Varoš
## 145 8 20
## Kozarska Dubica Kreševo Kupres
## 58 7 21
## Laktaši Livno Lopare
## 161 22 18
## Lukavac Ljubinje Ljubuški
## 396 2 37
## Maglaj Milići Modriča
## 70 9 91
## Mostar Mrkonjic Grad Neum
## 890 3 47
## Nevesinje Novi Grad Novi Travnik
## 11 82 156
## Odžak Olovo Orašje
## 35 33 23
## Osmaci Oštra Luka Pale - FBiH
## 6 2 3
## Pale - RS Pelagićevo Petrovo
## 152 7 10
## Posušje Prijedor Prnjavor
## 2 334 86
## Prozor Ravno Rogatica
## 11 20 43
## Rudo Sanski Most Sarajevo - Centar
## 2 114 2546
## Sarajevo - Novi Grad Sarajevo - Novo Sarajevo Sarajevo - Stari Grad
## 3290 2380 900
## Sokolac Srbac Srebrenica
## 21 64 31
## Srebrenik Stolac Šamac
## 154 31 18
## Šekovići Šipovo Široki Brijeg
## 1 4 15
## Teočak Teslić Tešanj
## 2 77 91
## Tomislavgrad Travnik Trebinje
## 11 254 121
## Trnovo - FBiH Trnovo - RS Tuzla
## 8 4 1877
## Ugljevik Usora Vareš
## 31 17 36
## Velika Kladuša Visoko Višegrad
## 86 181 52
## Vitez Vlasenica Vogošća
## 118 34 685
## Vukosavlje Zavidovići Zenica
## 5 173 1633
## Zvornik Žepče Živinice
## 127 53 304
table (olx1$prop_type) #11472 kuca i 18628 stanova
##
## Kuće Stanovi
## 11472 18628
unique(sort(olx1$Grad)) #135 opcina
## [1] "Banovići" "Banja Luka"
## [3] "Berkovići" "Bihać"
## [5] "Bijeljina" "Bileća"
## [7] "Bjelasnica" "Bosanska Krupa"
## [9] "Bosanski Petrovac" "Bosansko Grahovo"
## [11] "Bratunac" "Brčko"
## [13] "Breza" "Brod"
## [15] "Bugojno" "Busovača"
## [17] "Bužim" "Cazin"
## [19] "Čajniče" "Čapljina"
## [21] "Čelić" "Čelinac"
## [23] "Čitluk" "Derventa"
## [25] "Doboj" "Doboj Istok"
## [27] "Doboj Jug" "Dobretići"
## [29] "Domaljevac-Šamac" "Donji Vakuf"
## [31] "Donji Žabar" "Drvar"
## [33] "Foča" "Foča - Ustikolina"
## [35] "Fojnica" "Gacko"
## [37] "Glamoč" "Goražde"
## [39] "Gornji Vakuf-Uskoplje" "Gračanica"
## [41] "Gradačac" "Gradiška"
## [43] "Grude" "Hadžići"
## [45] "Han Pijesak" "Ilidža"
## [47] "Ilijaš" "Istočna Ilidža"
## [49] "Istočni Stari Grad" "Istočno Novo Sarajevo"
## [51] "Jablanica" "Jajce"
## [53] "Jezero" "Kakanj"
## [55] "Kalesija" "Kalinovik"
## [57] "Kiseljak" "Kladanj"
## [59] "Ključ" "Kneževo"
## [61] "Konjic" "Kostajnica"
## [63] "Kotor Varoš" "Kozarska Dubica"
## [65] "Kreševo" "Kupres"
## [67] "Laktaši" "Livno"
## [69] "Lopare" "Lukavac"
## [71] "Ljubinje" "Ljubuški"
## [73] "Maglaj" "Milići"
## [75] "Modriča" "Mostar"
## [77] "Mrkonjic Grad" "Neum"
## [79] "Nevesinje" "Novi Grad"
## [81] "Novi Travnik" "Odžak"
## [83] "Olovo" "Orašje"
## [85] "Osmaci" "Oštra Luka"
## [87] "Pale - FBiH" "Pale - RS"
## [89] "Pelagićevo" "Petrovo"
## [91] "Posušje" "Prijedor"
## [93] "Prnjavor" "Prozor"
## [95] "Ravno" "Rogatica"
## [97] "Rudo" "Sanski Most"
## [99] "Sarajevo - Centar" "Sarajevo - Novi Grad"
## [101] "Sarajevo - Novo Sarajevo" "Sarajevo - Stari Grad"
## [103] "Sokolac" "Srbac"
## [105] "Srebrenica" "Srebrenik"
## [107] "Stolac" "Šamac"
## [109] "Šekovići" "Šipovo"
## [111] "Široki Brijeg" "Teočak"
## [113] "Teslić" "Tešanj"
## [115] "Tomislavgrad" "Travnik"
## [117] "Trebinje" "Trnovo - FBiH"
## [119] "Trnovo - RS" "Tuzla"
## [121] "Ugljevik" "Usora"
## [123] "Vareš" "Velika Kladuša"
## [125] "Visoko" "Višegrad"
## [127] "Vitez" "Vlasenica"
## [129] "Vogošća" "Vukosavlje"
## [131] "Zavidovići" "Zenica"
## [133] "Zvornik" "Žepče"
## [135] "Živinice"
# pregledaj kvadraturu
summary(olx1$sqm)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.1 53.0 70.0 158.8 120.0 68685.0
# vrijednost kuce
summary(olx1$Cijena)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 10000 64000 91000 123387 140000 5900000 44
# izbaci percentile vrijednosti kuce
stan.kuce2 <- olx1 %>%
dplyr::filter(Cijena < quantile (olx1$Cijena, 0.99, na.rm = T)
& Cijena > quantile (olx1$Cijena, 0.01, na.rm = T))
summary(stan.kuce2$Cijena) #ima minimum 20400 to je ok i max 595000
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 20400 65000 91650 116093 139000 590000
# Izbaci percentile kvadrature
summary(stan.kuce2$sqm)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.1 53.0 70.0 150.8 120.0 68685.0
stan.kuce2 <- stan.kuce2 %>%
dplyr::filter(sqm < quantile (stan.kuce2$sqm, 0.99, na.rm = T)
& sqm > quantile (stan.kuce2$sqm, 0.01, na.rm = T))
summary(stan.kuce2$sqm) #this is more reasnable 25.2 min and 980 max
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.24 53.00 70.00 105.48 119.00 980.00
# made price per sqm
stan.kuce2$p.sqm <- round (stan.kuce2$Cijena/stan.kuce2$sqm,3)
summary(stan.kuce2$p.sqm) #26 km po kvm i 15161..maximum nije smislen
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 26.32 900.00 1347.72 1392.34 1800.00 15161.29
dim(stan.kuce2)
## [1] 28686 21
# izbaci percentile cijena per sqm ali samo upper
stan.kuce3 <- stan.kuce2 %>%
dplyr::filter(p.sqm < quantile (stan.kuce2$p.sqm, 0.99, na.rm = T)
)
summary (stan.kuce3$p.sqm) #niska cijena po kv je i moguca jer se odnosi na kuce u nedodjijama
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 26.32 900.00 1337.84 1362.53 1788.23 3470.59
unique(stan.kuce3$Grad) #sad je 133 dakle ne gubimo mnogo
## [1] "Bihać" "Mostar"
## [3] "Istočna Ilidža" "Ilidža"
## [5] "Sarajevo - Centar" "Gacko"
## [7] "Srbac" "Sarajevo - Novi Grad"
## [9] "Tuzla" "Zenica"
## [11] "Sarajevo - Novo Sarajevo" "Sarajevo - Stari Grad"
## [13] "Doboj" "Brčko"
## [15] "Banja Luka" "Kakanj"
## [17] "Travnik" "Bugojno"
## [19] "Bijeljina" "Banovići"
## [21] "Visoko" "Bosanska Krupa"
## [23] "Vogošća" "Busovača"
## [25] "Živinice" "Cazin"
## [27] "Vitez" "Konjic"
## [29] "Breza" "Brod"
## [31] "Istočno Novo Sarajevo" "Laktaši"
## [33] "Goražde" "Novi Travnik"
## [35] "Velika Kladuša" "Hadžići"
## [37] "Novi Grad" "Lukavac"
## [39] "Trebinje" "Ugljevik"
## [41] "Rogatica" "Gradiška"
## [43] "Foča" "Zavidovići"
## [45] "Derventa" "Višegrad"
## [47] "Sanski Most" "Zvornik"
## [49] "Ilijaš" "Maglaj"
## [51] "Čelić" "Gradačac"
## [53] "Usora" "Čitluk"
## [55] "Prnjavor" "Čapljina"
## [57] "Prijedor" "Jajce"
## [59] "Kalesija" "Jablanica"
## [61] "Grude" "Drvar"
## [63] "Kozarska Dubica" "Pale - RS"
## [65] "Olovo" "Tešanj"
## [67] "Gračanica" "Teslić"
## [69] "Livno" "Srebrenik"
## [71] "Kiseljak" "Kupres"
## [73] "Pale - FBiH" "Stolac"
## [75] "Vlasenica" "Bileća"
## [77] "Bosanski Petrovac" "Odžak"
## [79] "Donji Vakuf" "Žepče"
## [81] "Orašje" "Široki Brijeg"
## [83] "Kneževo" "Milići"
## [85] "Fojnica" "Modriča"
## [87] "Kotor Varoš" "Foča - Ustikolina"
## [89] "Bratunac" "Ravno"
## [91] "Doboj Jug" "Ljubuški"
## [93] "Kostajnica" "Neum"
## [95] "Kladanj" "Osmaci"
## [97] "Šamac" "Rudo"
## [99] "Trnovo - FBiH" "Petrovo"
## [101] "Pelagićevo" "Nevesinje"
## [103] "Vareš" "Ključ"
## [105] "Gornji Vakuf-Uskoplje" "Istočni Stari Grad"
## [107] "Lopare" "Čelinac"
## [109] "Sokolac" "Jezero"
## [111] "Tomislavgrad" "Domaljevac-Šamac"
## [113] "Han Pijesak" "Glamoč"
## [115] "Bjelasnica" "Šipovo"
## [117] "Vukosavlje" "Srebrenica"
## [119] "Kreševo" "Bužim"
## [121] "Prozor" "Posušje"
## [123] "Mrkonjic Grad" "Bosansko Grahovo"
## [125] "Trnovo - RS" "Dobretići"
## [127] "Donji Žabar" "Teočak"
## [129] "Berkovići" "Ljubinje"
## [131] "Doboj Istok" "Oštra Luka"
## [133] "Kalinovik"
dim(stan.kuce3) #izgubili smo oko 300 obs
## [1] 28397 21
stan.kuce3 %>%
group_by(Grad) %>%
summarise (n= n()) %>%
arrange(n) %>%
print(n=40) #35 less then 10
## # A tibble: 133 × 2
## Grad n
## <chr> <int>
## 1 Berkovići 1
## 2 Bjelasnica 1
## 3 Dobretići 1
## 4 Foča - Ustikolina 1
## 5 Kalinovik 1
## 6 Domaljevac-Šamac 2
## 7 Jezero 2
## 8 Ljubinje 2
## 9 Mrkonjic Grad 2
## 10 Oštra Luka 2
## 11 Posušje 2
## 12 Rudo 2
## 13 Teočak 2
## 14 Bosansko Grahovo 3
## 15 Donji Žabar 3
## 16 Istočni Stari Grad 3
## 17 Pale - FBiH 3
## 18 Glamoč 4
## 19 Grude 4
## 20 Han Pijesak 4
## 21 Trnovo - RS 4
## 22 Šipovo 4
## 23 Doboj Istok 5
## 24 Osmaci 5
## 25 Vukosavlje 5
## 26 Čitluk 5
## 27 Kreševo 6
## 28 Kostajnica 7
## 29 Pelagićevo 7
## 30 Prozor 7
## 31 Trnovo - FBiH 7
## 32 Bužim 8
## 33 Milići 8
## 34 Petrovo 8
## 35 Drvar 9
## 36 Nevesinje 10
## 37 Bileća 11
## 38 Tomislavgrad 11
## 39 Široki Brijeg 11
## 40 Kneževo 12
## # ℹ 93 more rows
stan.kuce2 %>% #35 less then 10
group_by(Grad) %>%
summarise (n= n()) %>%
arrange(n) %>%
print(n=40)
## # A tibble: 133 × 2
## Grad n
## <chr> <int>
## 1 Berkovići 1
## 2 Bjelasnica 1
## 3 Dobretići 1
## 4 Foča - Ustikolina 1
## 5 Kalinovik 1
## 6 Domaljevac-Šamac 2
## 7 Jezero 2
## 8 Ljubinje 2
## 9 Mrkonjic Grad 2
## 10 Oštra Luka 2
## 11 Posušje 2
## 12 Rudo 2
## 13 Teočak 2
## 14 Bosansko Grahovo 3
## 15 Donji Žabar 3
## 16 Istočni Stari Grad 3
## 17 Pale - FBiH 3
## 18 Glamoč 4
## 19 Grude 4
## 20 Han Pijesak 4
## 21 Trnovo - RS 4
## 22 Šipovo 4
## 23 Doboj Istok 5
## 24 Osmaci 5
## 25 Vukosavlje 5
## 26 Čitluk 5
## 27 Kreševo 6
## 28 Kostajnica 7
## 29 Pelagićevo 7
## 30 Prozor 7
## 31 Trnovo - FBiH 7
## 32 Bužim 8
## 33 Milići 8
## 34 Petrovo 8
## 35 Drvar 9
## 36 Nevesinje 10
## 37 Bileća 11
## 38 Tomislavgrad 11
## 39 Široki Brijeg 11
## 40 Kneževo 12
## # ℹ 93 more rows
# Delete low frequency obs (<10)
#make vector with low frequencies
low.freq <- stan.kuce3 %>% group_by(Grad) %>% summarise (n=n()) %>% dplyr::filter(n<10) %>% select (Grad)
stan.kuce4 <- stan.kuce3[!stan.kuce3$Grad %in% low.freq$Grad,] #ovo je bolji nacin da ne brisem ja pojedinacno
dim(stan.kuce4) #28257
## [1] 28257 21
stan.kuce4 %>%
group_by(Grad) %>%
summarise (n= n()) %>%
arrange(n) %>%
print(n=20) #no low freq# total 98
## # A tibble: 98 × 2
## Grad n
## <chr> <int>
## 1 Nevesinje 10
## 2 Bileća 11
## 3 Tomislavgrad 11
## 4 Široki Brijeg 11
## 5 Kneževo 12
## 6 Gacko 15
## 7 Kladanj 15
## 8 Usora 15
## 9 Lopare 16
## 10 Ravno 16
## 11 Šamac 16
## 12 Doboj Jug 17
## 13 Sokolac 18
## 14 Kotor Varoš 20
## 15 Vareš 20
## 16 Kupres 21
## 17 Livno 21
## 18 Orašje 21
## 19 Čelinac 23
## 20 Čelić 25
## # ℹ 78 more rows
# delete duplicates
#upper code on https://www.datanovia.com/en/lessons/identify-and-remove-duplicate-data-in-r/
#I did it in Excel so I dont do it now
#names(stan.kuce5)
#dist.s.k <- stan.kuce5 %>% distinct(Adresa, Cijena,Kvadratura,`Vrsta nekretnine`,Balkon, Lift, `Vrsta grijanja`, Sprat, Grad, `Godina izgradnje`, .keep_all = TRUE)
#View (dist.s.k)
#Prilagodi isti sadrzaj u Godini izgradnje: Ipak uzimam intervale po Ismirovoj instrukciji
table (stan.kuce4$`Godina izgradnje1`)
##
## 1950-1959 1955 1960-1969 1960 do 1969 1970-1979 1970 do 1979
## 169 42 858 2 1545 4
## 1980-1989 1980 do 1989 1990-1999 1990 do 1999 2000-2009 2000 do 2009
## 2715 16 897 7 1904 14
## 2010+ 2015+ 2020+ Novogradnja Prije 1950
## 1325 8 24 1929 325
table (stan.kuce4$`Godina izgradnje`)
##
## 1910 1940 1955 1965 1975 1985 1995 2005 2012 2015 2016 2017
## 263 62 211 860 1549 2731 904 1918 1325 5 879 1077
dist.s.k <- stan.kuce4 %>% mutate(constr.y = case_when(`Godina izgradnje1` == "Prije 1950" ~ "Until 1950",
`Godina izgradnje1` == "1950-1959" ~ "1950-1959",
`Godina izgradnje1` == "1960 do 1969" ~ "1960-1969",
`Godina izgradnje1` == "1960-1969" ~ "1960-1969",
`Godina izgradnje1` == "1970 do 1979" ~"1970-1979",
`Godina izgradnje1` == "1970-1979" ~"1970-1979",
`Godina izgradnje1` == "1980 do 1989" ~"1980-1989",
`Godina izgradnje1` == "1980-1989" ~"1980-1989",
`Godina izgradnje1` == "1990 do 1999" ~"1990-1999",
`Godina izgradnje1` == "1990-1999" ~"1990-1999",
`Godina izgradnje1` == "2000 do 2009" ~ "2000-2009",
`Godina izgradnje1` == "2000-2009" ~ "2000-2009",
`Godina izgradnje1` == "2010+" ~ "2010+",
`Godina izgradnje1` == "2010+" ~ "2015+",
`Godina izgradnje1` == "2010+" ~ "2020+",
`Godina izgradnje1` == "Novogradnja" ~ "New"))
table(dist.s.k$constr.y)
##
## 1950-1959 1960-1969 1970-1979 1980-1989 1990-1999 2000-2009 2010+
## 169 860 1549 2731 904 1918 1325
## New Until 1950
## 1929 325
#Sada treba relevel varijablu da "Until 1950" bude referentna vrijednost
table (is.na(dist.s.k$constr.y)) #16547 are NA
##
## FALSE TRUE
## 11710 16547
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats 1.0.0 ✔ readr 2.1.4
## ✔ ggplot2 3.4.4 ✔ stringr 1.5.1
## ✔ lubridate 1.9.3 ✔ tibble 3.2.1
## ✔ purrr 1.0.2
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library (conflicted)
dist.s.k$constr.y <- fct_relevel (dist.s.k$constr.y, "New")
table(dist.s.k$constr.y)
##
## New 1950-1959 1960-1969 1970-1979 1980-1989 1990-1999 2000-2009
## 1929 169 860 1549 2731 904 1918
## 2010+ Until 1950
## 1325 325
#Sprat staviti u intervale spratnost
#dist.s.k1 <- dist.s.k1 %>% mutate(Floor = case_when(Sprat == "Prizemlje" ~ "Prizemlje",
#Sprat == "Suteren" ~ "Suteren",
#Sprat == "Visoko prizemlje" ~ "Visoko prizemlje",
#Sprat %in% c(1,2,3) ~ "1-3",
#Sprat %in% c(4,5,6) ~ "4-6",
#Sprat %in% c(7,8,9) ~ "7-9",
#Sprat %in% c(10,11,12,13,14,15,16,17,18,19,"20+") ~ "10+"))
#table (dist.s.k1$Floor)
# Relevel Floor
#dist.s.k1$Floor <- fct_relevel (dist.s.k1$Floor, "Suteren","Prizemlje", "Visoko prizemlje","1-3","4-6", "7-9","10+")
#table (dist.s.k1$Floor)
#Extract year
dist.s.k$Y <- format(dist.s.k$date_of_Sale, format="%Y")
table (dist.s.k$Y)
##
## 2013 2014 2016 2017
## 1 1 13089 15166
#filter only 2016 and 2017
dist.s.k1 <- dist.s.k %>% dplyr::filter (Y >= 2016)
table (dist.s.k1$Y)
##
## 2016 2017
## 13089 15166
#Relevel B.Krupa is referent
library(tidyverse)
library (forcats)
dist.s.k1$Grad <- as.factor(dist.s.k1$Grad)
dist.s.k1$Grad <- fct_relevel(dist.s.k1$Grad, "Bosanska Krupa")
table (dist.s.k1$Grad)
##
## Bosanska Krupa Banovići Banja Luka
## 114 130 2334
## Bihać Bijeljina Bileća
## 469 896 11
## Bosanski Petrovac Bratunac Brčko
## 27 37 1460
## Breza Brod Bugojno
## 73 78 147
## Busovača Cazin Čapljina
## 61 117 73
## Čelić Čelinac Derventa
## 25 23 101
## Doboj Doboj Jug Donji Vakuf
## 452 17 44
## Foča Fojnica Gacko
## 74 41 15
## Goražde Gornji Vakuf-Uskoplje Gračanica
## 147 28 138
## Gradačac Gradiška Hadžići
## 105 187 403
## Ilidža Ilijaš Istočna Ilidža
## 1932 220 144
## Istočno Novo Sarajevo Jablanica Jajce
## 240 35 84
## Kakanj Kalesija Kiseljak
## 138 60 103
## Kladanj Ključ Kneževo
## 15 34 12
## Konjic Kotor Varoš Kozarska Dubica
## 123 20 50
## Kupres Laktaši Livno
## 21 152 21
## Lopare Lukavac Ljubuški
## 16 376 36
## Maglaj Modriča Mostar
## 64 85 862
## Neum Nevesinje Novi Grad
## 35 10 72
## Novi Travnik Odžak Olovo
## 148 31 29
## Orašje Pale - RS Prijedor
## 21 142 314
## Prnjavor Ravno Rogatica
## 82 16 37
## Sanski Most Sarajevo - Centar Sarajevo - Novi Grad
## 104 2329 3213
## Sarajevo - Novo Sarajevo Sarajevo - Stari Grad Sokolac
## 2326 806 18
## Srbac Srebrenica Srebrenik
## 58 30 132
## Stolac Šamac Široki Brijeg
## 27 16 11
## Teslić Tešanj Tomislavgrad
## 67 80 11
## Travnik Trebinje Tuzla
## 240 112 1811
## Ugljevik Usora Vareš
## 30 15 20
## Velika Kladuša Visoko Višegrad
## 83 165 48
## Vitez Vlasenica Vogošća
## 115 26 643
## Zavidovići Zenica Zvornik
## 167 1515 115
## Žepče Živinice
## 47 268
table (dist.s.k1$constr.y)
##
## New 1950-1959 1960-1969 1970-1979 1980-1989 1990-1999 2000-2009
## 1929 169 860 1549 2731 904 1918
## 2010+ Until 1950
## 1325 325
###########################
#MODEL 1: with age of prop#
###########################
names (dist.s.k1)
## [1] "ID" "date_of_Sale" "prop_type"
## [4] "sqm" "text" "Okucnica"
## [7] "Kvadratura okućnice" "Sprat" "Adresa"
## [10] "Grad" "Grad1" "Cijena"
## [13] "Cijena_original" "prop.age" "Godina izgradnje"
## [16] "Godina izgradnje1" "Nedavno adaptiran" "Balkon"
## [19] "Vrsta grijanja" "Lift" "p.sqm"
## [22] "constr.y" "Y"
model1 <- lm (log(p.sqm) ~ factor (Y) + as.factor(constr.y) +factor (prop_type) + factor (`Vrsta grijanja`) + factor (Grad), data = dist.s.k1)
summary(model1)$coefficients
## Estimate Std. Error
## (Intercept) 6.565634109 0.066072169
## factor(Y)2017 0.038164178 0.007567146
## as.factor(constr.y)1950-1959 -0.005747513 0.033212356
## as.factor(constr.y)1960-1969 -0.048890339 0.017410617
## as.factor(constr.y)1970-1979 -0.070126662 0.014747864
## as.factor(constr.y)1980-1989 -0.115088151 0.012828533
## as.factor(constr.y)1990-1999 -0.182284966 0.017051254
## as.factor(constr.y)2000-2009 -0.031534135 0.013573741
## as.factor(constr.y)2010+ -0.008948840 0.014674099
## as.factor(constr.y)Until 1950 0.037934485 0.026016440
## factor(prop_type)Stanovi 0.509133201 0.011227635
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.019944344 0.015104062
## factor(`Vrsta grijanja`)Centralno (Plin) 0.017525386 0.016864097
## factor(`Vrsta grijanja`)Drva -0.265391897 0.016162168
## factor(`Vrsta grijanja`)Ostalo -0.282873130 0.017220251
## factor(`Vrsta grijanja`)Plin -0.054937554 0.017199089
## factor(`Vrsta grijanja`)Struja -0.079544069 0.015001516
## factor(Grad)Banovići -0.051878148 0.088267642
## factor(Grad)Banja Luka 0.535288888 0.064701454
## factor(Grad)Bihać 0.425798791 0.069396882
## factor(Grad)Bijeljina 0.169299574 0.067390479
## factor(Grad)Bileća 0.277358440 0.166183676
## factor(Grad)Bosanski Petrovac -0.378999119 0.143144364
## factor(Grad)Bratunac -0.281497145 0.104289561
## factor(Grad)Brčko 0.094179820 0.067448364
## factor(Grad)Breza 0.125523320 0.094205301
## factor(Grad)Brod -0.201915039 0.096772010
## factor(Grad)Bugojno 0.077311437 0.079242654
## factor(Grad)Busovača 0.062251132 0.094145666
## factor(Grad)Cazin 0.236929143 0.087669096
## factor(Grad)Čapljina 0.247099360 0.096705962
## factor(Grad)Čelić -0.104948850 0.125672680
## factor(Grad)Čelinac 0.155445020 0.122513064
## factor(Grad)Derventa 0.008443602 0.086940530
## factor(Grad)Doboj -0.016383099 0.069419289
## factor(Grad)Doboj Jug 0.058037723 0.129339822
## factor(Grad)Donji Vakuf -0.150163441 0.119774205
## factor(Grad)Foča -0.048739889 0.101731613
## factor(Grad)Fojnica 0.291268073 0.108936165
## factor(Grad)Gacko 0.166326048 0.192236149
## factor(Grad)Goražde 0.154069404 0.084450944
## factor(Grad)Gornji Vakuf-Uskoplje 0.083413199 0.129212404
## factor(Grad)Gračanica 0.093656987 0.080485963
## factor(Grad)Gradačac 0.134379239 0.083438265
## factor(Grad)Gradiška 0.007134140 0.076144766
## factor(Grad)Hadžići 0.376020836 0.070943995
## factor(Grad)Ilidža 0.450822467 0.065152132
## factor(Grad)Ilijaš 0.399930153 0.076669567
## factor(Grad)Istočna Ilidža -0.001797591 0.079623412
## factor(Grad)Istočno Novo Sarajevo 0.107438565 0.072551319
## factor(Grad)Jablanica 0.360848449 0.122535079
## factor(Grad)Jajce 0.130553595 0.092889361
## factor(Grad)Kakanj 0.177639352 0.088669590
## factor(Grad)Kalesija -0.082943642 0.100581828
## factor(Grad)Kiseljak 0.223743182 0.089129372
## factor(Grad)Kladanj 0.217577458 0.177392949
## factor(Grad)Ključ -0.301365070 0.119634133
## factor(Grad)Kneževo -0.595962949 0.192329287
## factor(Grad)Konjic 0.438458032 0.085144502
## factor(Grad)Kotor Varoš -0.304575011 0.149501972
## factor(Grad)Kozarska Dubica -0.171734444 0.101847948
## factor(Grad)Kupres 0.037380089 0.133239474
## factor(Grad)Laktaši 0.335356551 0.080016665
## factor(Grad)Livno -0.161060995 0.149549243
## factor(Grad)Lopare 0.003805227 0.212795453
## factor(Grad)Lukavac 0.009331751 0.070889940
## factor(Grad)Ljubuški -0.127231668 0.138037035
## factor(Grad)Maglaj -0.079742300 0.104326665
## factor(Grad)Modriča -0.008969180 0.087682791
## factor(Grad)Mostar 0.500686993 0.068244915
## factor(Grad)Neum 0.797363775 0.129295161
## factor(Grad)Nevesinje -0.123455257 0.212629992
## factor(Grad)Novi Grad -0.006459067 0.093466586
## factor(Grad)Novi Travnik -0.105794583 0.088336686
## factor(Grad)Odžak 0.213380913 0.125782549
## factor(Grad)Olovo 0.293342583 0.137813944
## factor(Grad)Orašje 0.118540557 0.137966651
## factor(Grad)Pale - RS 0.024650278 0.079922971
## factor(Grad)Prijedor 0.024275700 0.072888476
## factor(Grad)Prnjavor 0.158572384 0.085522559
## factor(Grad)Ravno 0.228681165 0.149588200
## factor(Grad)Rogatica -0.014003421 0.133310271
## factor(Grad)Sanski Most 0.022412903 0.088111546
## factor(Grad)Sarajevo - Centar 0.696263136 0.065605251
## factor(Grad)Sarajevo - Novi Grad 0.400650048 0.064889793
## factor(Grad)Sarajevo - Novo Sarajevo 0.550293817 0.065492492
## factor(Grad)Sarajevo - Stari Grad 0.594268072 0.067477798
## factor(Grad)Sokolac -0.258936232 0.143143584
## factor(Grad)Srbac -0.049689971 0.105726404
## factor(Grad)Srebrenica -0.241660246 0.122615506
## factor(Grad)Srebrenik 0.177251827 0.078826656
## factor(Grad)Stolac -0.166173377 0.137795919
## factor(Grad)Šamac -0.485802359 0.166091272
## factor(Grad)Široki Brijeg 0.391146805 0.156905963
## factor(Grad)Teslić 0.216560491 0.104342440
## factor(Grad)Tešanj 0.251977175 0.095772704
## factor(Grad)Tomislavgrad 0.349045777 0.192306936
## factor(Grad)Travnik 0.206793698 0.077418787
## factor(Grad)Trebinje 0.436371443 0.082194455
## factor(Grad)Tuzla 0.231695105 0.065572088
## factor(Grad)Ugljevik 0.202382106 0.129334115
## factor(Grad)Usora -0.238131726 0.138264737
## factor(Grad)Vareš -0.709000127 0.192327707
## factor(Grad)Velika Kladuša 0.176858252 0.091428027
## factor(Grad)Visoko 0.287068647 0.080761730
## factor(Grad)Višegrad -0.254226117 0.110740564
## factor(Grad)Vitez -0.066817577 0.085932082
## factor(Grad)Vlasenica -0.304678570 0.133209788
## factor(Grad)Vogošća 0.362489801 0.067768599
## factor(Grad)Zavidovići -0.035170452 0.084431574
## factor(Grad)Zenica 0.211522423 0.066919767
## factor(Grad)Zvornik -0.209338981 0.080514165
## factor(Grad)Žepče -0.125598391 0.103063619
## factor(Grad)Živinice 0.103319816 0.072201158
## t value Pr(>|t|)
## (Intercept) 99.37064637 0.000000e+00
## factor(Y)2017 5.04340475 4.642354e-07
## as.factor(constr.y)1950-1959 -0.17305344 8.626124e-01
## as.factor(constr.y)1960-1969 -2.80807627 4.992171e-03
## as.factor(constr.y)1970-1979 -4.75503872 2.007967e-06
## as.factor(constr.y)1980-1989 -8.97126332 3.380092e-19
## as.factor(constr.y)1990-1999 -10.69041421 1.500294e-26
## as.factor(constr.y)2000-2009 -2.32317195 2.018712e-02
## as.factor(constr.y)2010+ -0.60983917 5.419803e-01
## as.factor(constr.y)Until 1950 1.45809666 1.448410e-01
## factor(prop_type)Stanovi 45.34643371 0.000000e+00
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -1.32046229 1.867068e-01
## factor(`Vrsta grijanja`)Centralno (Plin) 1.03921282 2.987274e-01
## factor(`Vrsta grijanja`)Drva -16.42056263 6.453818e-60
## factor(`Vrsta grijanja`)Ostalo -16.42677193 5.839493e-60
## factor(`Vrsta grijanja`)Plin -3.19421302 1.405876e-03
## factor(`Vrsta grijanja`)Struja -5.30240186 1.163890e-07
## factor(Grad)Banovići -0.58773688 5.567204e-01
## factor(Grad)Banja Luka 8.27321260 1.446203e-16
## factor(Grad)Bihać 6.13570498 8.754737e-10
## factor(Grad)Bijeljina 2.51221800 1.201096e-02
## factor(Grad)Bileća 1.66898727 9.514690e-02
## factor(Grad)Bosanski Petrovac -2.64767056 8.115810e-03
## factor(Grad)Bratunac -2.69918812 6.960962e-03
## factor(Grad)Brčko 1.39632476 1.626434e-01
## factor(Grad)Breza 1.33244434 1.827404e-01
## factor(Grad)Brod -2.08650247 3.695492e-02
## factor(Grad)Bugojno 0.97562907 3.292687e-01
## factor(Grad)Busovača 0.66122143 5.084835e-01
## factor(Grad)Cazin 2.70253889 6.891231e-03
## factor(Grad)Čapljina 2.55516159 1.062650e-02
## factor(Grad)Čelić -0.83509678 4.036804e-01
## factor(Grad)Čelinac 1.26880362 2.045366e-01
## factor(Grad)Derventa 0.09711928 9.226333e-01
## factor(Grad)Doboj -0.23600211 8.134352e-01
## factor(Grad)Doboj Jug 0.44872277 6.536400e-01
## factor(Grad)Donji Vakuf -1.25372104 2.099687e-01
## factor(Grad)Foča -0.47910269 6.318746e-01
## factor(Grad)Fojnica 2.67375001 7.511343e-03
## factor(Grad)Gacko 0.86521733 3.869374e-01
## factor(Grad)Goražde 1.82436568 6.812253e-02
## factor(Grad)Gornji Vakuf-Uskoplje 0.64555102 5.185829e-01
## factor(Grad)Gračanica 1.16364373 2.445924e-01
## factor(Grad)Gradačac 1.61052293 1.073110e-01
## factor(Grad)Gradiška 0.09369179 9.253556e-01
## factor(Grad)Hadžići 5.30024897 1.177665e-07
## factor(Grad)Ilidža 6.91953511 4.769947e-12
## factor(Grad)Ilijaš 5.21628295 1.856962e-07
## factor(Grad)Istočna Ilidža -0.02257616 9.819887e-01
## factor(Grad)Istočno Novo Sarajevo 1.48086300 1.386703e-01
## factor(Grad)Jablanica 2.94485833 3.237441e-03
## factor(Grad)Jajce 1.40547414 1.599071e-01
## factor(Grad)Kakanj 2.00338528 4.515918e-02
## factor(Grad)Kalesija -0.82463844 4.095939e-01
## factor(Grad)Kiseljak 2.51031928 1.207571e-02
## factor(Grad)Kladanj 1.22652822 2.200249e-01
## factor(Grad)Ključ -2.51905591 1.178033e-02
## factor(Grad)Kneževo -3.09865937 1.948634e-03
## factor(Grad)Konjic 5.14957541 2.653574e-07
## factor(Grad)Kotor Varoš -2.03726417 4.164623e-02
## factor(Grad)Kozarska Dubica -1.68618463 9.178715e-02
## factor(Grad)Kupres 0.28054816 7.790620e-01
## factor(Grad)Laktaši 4.19108382 2.796835e-05
## factor(Grad)Livno -1.07697633 2.815132e-01
## factor(Grad)Lopare 0.01788209 9.857332e-01
## factor(Grad)Lukavac 0.13163716 8.952736e-01
## factor(Grad)Ljubuški -0.92172124 3.566932e-01
## factor(Grad)Maglaj -0.76435205 4.446730e-01
## factor(Grad)Modriča -0.10229123 9.185273e-01
## factor(Grad)Mostar 7.33661973 2.336660e-13
## factor(Grad)Neum 6.16700400 7.191193e-10
## factor(Grad)Nevesinje -0.58061074 5.615141e-01
## factor(Grad)Novi Grad -0.06910563 9.449067e-01
## factor(Grad)Novi Travnik -1.19762907 2.310859e-01
## factor(Grad)Odžak 1.69642701 8.983188e-02
## factor(Grad)Olovo 2.12854065 3.331330e-02
## factor(Grad)Orašje 0.85919718 3.902495e-01
## factor(Grad)Pale - RS 0.30842545 7.577642e-01
## factor(Grad)Prijedor 0.33305265 7.391006e-01
## factor(Grad)Prnjavor 1.85415855 6.374185e-02
## factor(Grad)Ravno 1.52873799 1.263567e-01
## factor(Grad)Rogatica -0.10504382 9.163429e-01
## factor(Grad)Sanski Most 0.25436964 7.992145e-01
## factor(Grad)Sarajevo - Centar 10.61291773 3.421954e-26
## factor(Grad)Sarajevo - Novi Grad 6.17431538 6.867246e-10
## factor(Grad)Sarajevo - Novo Sarajevo 8.40239540 4.883491e-17
## factor(Grad)Sarajevo - Stari Grad 8.80686824 1.469136e-18
## factor(Grad)Sokolac -1.80892657 7.048830e-02
## factor(Grad)Srbac -0.46998639 6.383736e-01
## factor(Grad)Srebrenica -1.97087835 4.876154e-02
## factor(Grad)Srebrenik 2.24862800 2.455487e-02
## factor(Grad)Stolac -1.20593831 2.278659e-01
## factor(Grad)Šamac -2.92491202 3.452211e-03
## factor(Grad)Široki Brijeg 2.49287406 1.268522e-02
## factor(Grad)Teslić 2.07547850 3.796420e-02
## factor(Grad)Tešanj 2.63099155 8.524875e-03
## factor(Grad)Tomislavgrad 1.81504518 6.954279e-02
## factor(Grad)Travnik 2.67110487 7.570757e-03
## factor(Grad)Trebinje 5.30901319 1.122556e-07
## factor(Grad)Tuzla 3.53344101 4.117836e-04
## factor(Grad)Ugljevik 1.56480064 1.176569e-01
## factor(Grad)Usora -1.72228821 8.504399e-02
## factor(Grad)Vareš -3.68641700 2.284709e-04
## factor(Grad)Velika Kladuša 1.93439865 5.308843e-02
## factor(Grad)Visoko 3.55451336 3.801858e-04
## factor(Grad)Višegrad -2.29569100 2.171131e-02
## factor(Grad)Vitez -0.77756265 4.368427e-01
## factor(Grad)Vlasenica -2.28720858 2.220158e-02
## factor(Grad)Vogošća 5.34893454 9.015569e-08
## factor(Grad)Zavidovići -0.41655568 6.770111e-01
## factor(Grad)Zenica 3.16083619 1.577217e-03
## factor(Grad)Zvornik -2.60002673 9.333469e-03
## factor(Grad)Žepče -1.21864915 2.230022e-01
## factor(Grad)Živinice 1.43099943 1.524573e-01
summary(model1) #16545 deleted due to missingness
##
## Call:
## lm(formula = log(p.sqm) ~ factor(Y) + as.factor(constr.y) + factor(prop_type) +
## factor(`Vrsta grijanja`) + factor(Grad), data = dist.s.k1)
##
## Residuals:
## Min 1Q Median 3Q Max
## -3.10185 -0.15745 0.01065 0.18481 1.73099
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 6.565634 0.066072 99.371
## factor(Y)2017 0.038164 0.007567 5.043
## as.factor(constr.y)1950-1959 -0.005748 0.033212 -0.173
## as.factor(constr.y)1960-1969 -0.048890 0.017411 -2.808
## as.factor(constr.y)1970-1979 -0.070127 0.014748 -4.755
## as.factor(constr.y)1980-1989 -0.115088 0.012829 -8.971
## as.factor(constr.y)1990-1999 -0.182285 0.017051 -10.690
## as.factor(constr.y)2000-2009 -0.031534 0.013574 -2.323
## as.factor(constr.y)2010+ -0.008949 0.014674 -0.610
## as.factor(constr.y)Until 1950 0.037934 0.026016 1.458
## factor(prop_type)Stanovi 0.509133 0.011228 45.346
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.019944 0.015104 -1.320
## factor(`Vrsta grijanja`)Centralno (Plin) 0.017525 0.016864 1.039
## factor(`Vrsta grijanja`)Drva -0.265392 0.016162 -16.421
## factor(`Vrsta grijanja`)Ostalo -0.282873 0.017220 -16.427
## factor(`Vrsta grijanja`)Plin -0.054938 0.017199 -3.194
## factor(`Vrsta grijanja`)Struja -0.079544 0.015002 -5.302
## factor(Grad)Banovići -0.051878 0.088268 -0.588
## factor(Grad)Banja Luka 0.535289 0.064701 8.273
## factor(Grad)Bihać 0.425799 0.069397 6.136
## factor(Grad)Bijeljina 0.169300 0.067390 2.512
## factor(Grad)Bileća 0.277358 0.166184 1.669
## factor(Grad)Bosanski Petrovac -0.378999 0.143144 -2.648
## factor(Grad)Bratunac -0.281497 0.104290 -2.699
## factor(Grad)Brčko 0.094180 0.067448 1.396
## factor(Grad)Breza 0.125523 0.094205 1.332
## factor(Grad)Brod -0.201915 0.096772 -2.087
## factor(Grad)Bugojno 0.077311 0.079243 0.976
## factor(Grad)Busovača 0.062251 0.094146 0.661
## factor(Grad)Cazin 0.236929 0.087669 2.703
## factor(Grad)Čapljina 0.247099 0.096706 2.555
## factor(Grad)Čelić -0.104949 0.125673 -0.835
## factor(Grad)Čelinac 0.155445 0.122513 1.269
## factor(Grad)Derventa 0.008444 0.086941 0.097
## factor(Grad)Doboj -0.016383 0.069419 -0.236
## factor(Grad)Doboj Jug 0.058038 0.129340 0.449
## factor(Grad)Donji Vakuf -0.150163 0.119774 -1.254
## factor(Grad)Foča -0.048740 0.101732 -0.479
## factor(Grad)Fojnica 0.291268 0.108936 2.674
## factor(Grad)Gacko 0.166326 0.192236 0.865
## factor(Grad)Goražde 0.154069 0.084451 1.824
## factor(Grad)Gornji Vakuf-Uskoplje 0.083413 0.129212 0.646
## factor(Grad)Gračanica 0.093657 0.080486 1.164
## factor(Grad)Gradačac 0.134379 0.083438 1.611
## factor(Grad)Gradiška 0.007134 0.076145 0.094
## factor(Grad)Hadžići 0.376021 0.070944 5.300
## factor(Grad)Ilidža 0.450822 0.065152 6.920
## factor(Grad)Ilijaš 0.399930 0.076670 5.216
## factor(Grad)Istočna Ilidža -0.001798 0.079623 -0.023
## factor(Grad)Istočno Novo Sarajevo 0.107439 0.072551 1.481
## factor(Grad)Jablanica 0.360848 0.122535 2.945
## factor(Grad)Jajce 0.130554 0.092889 1.405
## factor(Grad)Kakanj 0.177639 0.088670 2.003
## factor(Grad)Kalesija -0.082944 0.100582 -0.825
## factor(Grad)Kiseljak 0.223743 0.089129 2.510
## factor(Grad)Kladanj 0.217577 0.177393 1.227
## factor(Grad)Ključ -0.301365 0.119634 -2.519
## factor(Grad)Kneževo -0.595963 0.192329 -3.099
## factor(Grad)Konjic 0.438458 0.085145 5.150
## factor(Grad)Kotor Varoš -0.304575 0.149502 -2.037
## factor(Grad)Kozarska Dubica -0.171734 0.101848 -1.686
## factor(Grad)Kupres 0.037380 0.133239 0.281
## factor(Grad)Laktaši 0.335357 0.080017 4.191
## factor(Grad)Livno -0.161061 0.149549 -1.077
## factor(Grad)Lopare 0.003805 0.212795 0.018
## factor(Grad)Lukavac 0.009332 0.070890 0.132
## factor(Grad)Ljubuški -0.127232 0.138037 -0.922
## factor(Grad)Maglaj -0.079742 0.104327 -0.764
## factor(Grad)Modriča -0.008969 0.087683 -0.102
## factor(Grad)Mostar 0.500687 0.068245 7.337
## factor(Grad)Neum 0.797364 0.129295 6.167
## factor(Grad)Nevesinje -0.123455 0.212630 -0.581
## factor(Grad)Novi Grad -0.006459 0.093467 -0.069
## factor(Grad)Novi Travnik -0.105795 0.088337 -1.198
## factor(Grad)Odžak 0.213381 0.125783 1.696
## factor(Grad)Olovo 0.293343 0.137814 2.129
## factor(Grad)Orašje 0.118541 0.137967 0.859
## factor(Grad)Pale - RS 0.024650 0.079923 0.308
## factor(Grad)Prijedor 0.024276 0.072888 0.333
## factor(Grad)Prnjavor 0.158572 0.085523 1.854
## factor(Grad)Ravno 0.228681 0.149588 1.529
## factor(Grad)Rogatica -0.014003 0.133310 -0.105
## factor(Grad)Sanski Most 0.022413 0.088112 0.254
## factor(Grad)Sarajevo - Centar 0.696263 0.065605 10.613
## factor(Grad)Sarajevo - Novi Grad 0.400650 0.064890 6.174
## factor(Grad)Sarajevo - Novo Sarajevo 0.550294 0.065492 8.402
## factor(Grad)Sarajevo - Stari Grad 0.594268 0.067478 8.807
## factor(Grad)Sokolac -0.258936 0.143144 -1.809
## factor(Grad)Srbac -0.049690 0.105726 -0.470
## factor(Grad)Srebrenica -0.241660 0.122616 -1.971
## factor(Grad)Srebrenik 0.177252 0.078827 2.249
## factor(Grad)Stolac -0.166173 0.137796 -1.206
## factor(Grad)Šamac -0.485802 0.166091 -2.925
## factor(Grad)Široki Brijeg 0.391147 0.156906 2.493
## factor(Grad)Teslić 0.216560 0.104342 2.075
## factor(Grad)Tešanj 0.251977 0.095773 2.631
## factor(Grad)Tomislavgrad 0.349046 0.192307 1.815
## factor(Grad)Travnik 0.206794 0.077419 2.671
## factor(Grad)Trebinje 0.436371 0.082194 5.309
## factor(Grad)Tuzla 0.231695 0.065572 3.533
## factor(Grad)Ugljevik 0.202382 0.129334 1.565
## factor(Grad)Usora -0.238132 0.138265 -1.722
## factor(Grad)Vareš -0.709000 0.192328 -3.686
## factor(Grad)Velika Kladuša 0.176858 0.091428 1.934
## factor(Grad)Visoko 0.287069 0.080762 3.555
## factor(Grad)Višegrad -0.254226 0.110741 -2.296
## factor(Grad)Vitez -0.066818 0.085932 -0.778
## factor(Grad)Vlasenica -0.304679 0.133210 -2.287
## factor(Grad)Vogošća 0.362490 0.067769 5.349
## factor(Grad)Zavidovići -0.035170 0.084432 -0.417
## factor(Grad)Zenica 0.211522 0.066920 3.161
## factor(Grad)Zvornik -0.209339 0.080514 -2.600
## factor(Grad)Žepče -0.125598 0.103064 -1.219
## factor(Grad)Živinice 0.103320 0.072201 1.431
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## factor(Y)2017 4.64e-07 ***
## as.factor(constr.y)1950-1959 0.862612
## as.factor(constr.y)1960-1969 0.004992 **
## as.factor(constr.y)1970-1979 2.01e-06 ***
## as.factor(constr.y)1980-1989 < 2e-16 ***
## as.factor(constr.y)1990-1999 < 2e-16 ***
## as.factor(constr.y)2000-2009 0.020187 *
## as.factor(constr.y)2010+ 0.541980
## as.factor(constr.y)Until 1950 0.144841
## factor(prop_type)Stanovi < 2e-16 ***
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 0.186707
## factor(`Vrsta grijanja`)Centralno (Plin) 0.298727
## factor(`Vrsta grijanja`)Drva < 2e-16 ***
## factor(`Vrsta grijanja`)Ostalo < 2e-16 ***
## factor(`Vrsta grijanja`)Plin 0.001406 **
## factor(`Vrsta grijanja`)Struja 1.16e-07 ***
## factor(Grad)Banovići 0.556720
## factor(Grad)Banja Luka < 2e-16 ***
## factor(Grad)Bihać 8.75e-10 ***
## factor(Grad)Bijeljina 0.012011 *
## factor(Grad)Bileća 0.095147 .
## factor(Grad)Bosanski Petrovac 0.008116 **
## factor(Grad)Bratunac 0.006961 **
## factor(Grad)Brčko 0.162643
## factor(Grad)Breza 0.182740
## factor(Grad)Brod 0.036955 *
## factor(Grad)Bugojno 0.329269
## factor(Grad)Busovača 0.508483
## factor(Grad)Cazin 0.006891 **
## factor(Grad)Čapljina 0.010627 *
## factor(Grad)Čelić 0.403680
## factor(Grad)Čelinac 0.204537
## factor(Grad)Derventa 0.922633
## factor(Grad)Doboj 0.813435
## factor(Grad)Doboj Jug 0.653640
## factor(Grad)Donji Vakuf 0.209969
## factor(Grad)Foča 0.631875
## factor(Grad)Fojnica 0.007511 **
## factor(Grad)Gacko 0.386937
## factor(Grad)Goražde 0.068123 .
## factor(Grad)Gornji Vakuf-Uskoplje 0.518583
## factor(Grad)Gračanica 0.244592
## factor(Grad)Gradačac 0.107311
## factor(Grad)Gradiška 0.925356
## factor(Grad)Hadžići 1.18e-07 ***
## factor(Grad)Ilidža 4.77e-12 ***
## factor(Grad)Ilijaš 1.86e-07 ***
## factor(Grad)Istočna Ilidža 0.981989
## factor(Grad)Istočno Novo Sarajevo 0.138670
## factor(Grad)Jablanica 0.003237 **
## factor(Grad)Jajce 0.159907
## factor(Grad)Kakanj 0.045159 *
## factor(Grad)Kalesija 0.409594
## factor(Grad)Kiseljak 0.012076 *
## factor(Grad)Kladanj 0.220025
## factor(Grad)Ključ 0.011780 *
## factor(Grad)Kneževo 0.001949 **
## factor(Grad)Konjic 2.65e-07 ***
## factor(Grad)Kotor Varoš 0.041646 *
## factor(Grad)Kozarska Dubica 0.091787 .
## factor(Grad)Kupres 0.779062
## factor(Grad)Laktaši 2.80e-05 ***
## factor(Grad)Livno 0.281513
## factor(Grad)Lopare 0.985733
## factor(Grad)Lukavac 0.895274
## factor(Grad)Ljubuški 0.356693
## factor(Grad)Maglaj 0.444673
## factor(Grad)Modriča 0.918527
## factor(Grad)Mostar 2.34e-13 ***
## factor(Grad)Neum 7.19e-10 ***
## factor(Grad)Nevesinje 0.561514
## factor(Grad)Novi Grad 0.944907
## factor(Grad)Novi Travnik 0.231086
## factor(Grad)Odžak 0.089832 .
## factor(Grad)Olovo 0.033313 *
## factor(Grad)Orašje 0.390250
## factor(Grad)Pale - RS 0.757764
## factor(Grad)Prijedor 0.739101
## factor(Grad)Prnjavor 0.063742 .
## factor(Grad)Ravno 0.126357
## factor(Grad)Rogatica 0.916343
## factor(Grad)Sanski Most 0.799215
## factor(Grad)Sarajevo - Centar < 2e-16 ***
## factor(Grad)Sarajevo - Novi Grad 6.87e-10 ***
## factor(Grad)Sarajevo - Novo Sarajevo < 2e-16 ***
## factor(Grad)Sarajevo - Stari Grad < 2e-16 ***
## factor(Grad)Sokolac 0.070488 .
## factor(Grad)Srbac 0.638374
## factor(Grad)Srebrenica 0.048762 *
## factor(Grad)Srebrenik 0.024555 *
## factor(Grad)Stolac 0.227866
## factor(Grad)Šamac 0.003452 **
## factor(Grad)Široki Brijeg 0.012685 *
## factor(Grad)Teslić 0.037964 *
## factor(Grad)Tešanj 0.008525 **
## factor(Grad)Tomislavgrad 0.069543 .
## factor(Grad)Travnik 0.007571 **
## factor(Grad)Trebinje 1.12e-07 ***
## factor(Grad)Tuzla 0.000412 ***
## factor(Grad)Ugljevik 0.117657
## factor(Grad)Usora 0.085044 .
## factor(Grad)Vareš 0.000228 ***
## factor(Grad)Velika Kladuša 0.053088 .
## factor(Grad)Visoko 0.000380 ***
## factor(Grad)Višegrad 0.021711 *
## factor(Grad)Vitez 0.436843
## factor(Grad)Vlasenica 0.022202 *
## factor(Grad)Vogošća 9.02e-08 ***
## factor(Grad)Zavidovići 0.677011
## factor(Grad)Zenica 0.001577 **
## factor(Grad)Zvornik 0.009333 **
## factor(Grad)Žepče 0.223002
## factor(Grad)Živinice 0.152457
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4057 on 11596 degrees of freedom
## (16545 observations deleted due to missingness)
## Multiple R-squared: 0.5754, Adjusted R-squared: 0.5713
## F-statistic: 139.1 on 113 and 11596 DF, p-value: < 2.2e-16
model.coef1 <- summary(model1)$coefficients #
model.coef1 # vidi koliko redova su koef.za druge varijable (17redova) i Opstina Bosanska Krupa = Ref
## Estimate Std. Error
## (Intercept) 6.565634109 0.066072169
## factor(Y)2017 0.038164178 0.007567146
## as.factor(constr.y)1950-1959 -0.005747513 0.033212356
## as.factor(constr.y)1960-1969 -0.048890339 0.017410617
## as.factor(constr.y)1970-1979 -0.070126662 0.014747864
## as.factor(constr.y)1980-1989 -0.115088151 0.012828533
## as.factor(constr.y)1990-1999 -0.182284966 0.017051254
## as.factor(constr.y)2000-2009 -0.031534135 0.013573741
## as.factor(constr.y)2010+ -0.008948840 0.014674099
## as.factor(constr.y)Until 1950 0.037934485 0.026016440
## factor(prop_type)Stanovi 0.509133201 0.011227635
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.019944344 0.015104062
## factor(`Vrsta grijanja`)Centralno (Plin) 0.017525386 0.016864097
## factor(`Vrsta grijanja`)Drva -0.265391897 0.016162168
## factor(`Vrsta grijanja`)Ostalo -0.282873130 0.017220251
## factor(`Vrsta grijanja`)Plin -0.054937554 0.017199089
## factor(`Vrsta grijanja`)Struja -0.079544069 0.015001516
## factor(Grad)Banovići -0.051878148 0.088267642
## factor(Grad)Banja Luka 0.535288888 0.064701454
## factor(Grad)Bihać 0.425798791 0.069396882
## factor(Grad)Bijeljina 0.169299574 0.067390479
## factor(Grad)Bileća 0.277358440 0.166183676
## factor(Grad)Bosanski Petrovac -0.378999119 0.143144364
## factor(Grad)Bratunac -0.281497145 0.104289561
## factor(Grad)Brčko 0.094179820 0.067448364
## factor(Grad)Breza 0.125523320 0.094205301
## factor(Grad)Brod -0.201915039 0.096772010
## factor(Grad)Bugojno 0.077311437 0.079242654
## factor(Grad)Busovača 0.062251132 0.094145666
## factor(Grad)Cazin 0.236929143 0.087669096
## factor(Grad)Čapljina 0.247099360 0.096705962
## factor(Grad)Čelić -0.104948850 0.125672680
## factor(Grad)Čelinac 0.155445020 0.122513064
## factor(Grad)Derventa 0.008443602 0.086940530
## factor(Grad)Doboj -0.016383099 0.069419289
## factor(Grad)Doboj Jug 0.058037723 0.129339822
## factor(Grad)Donji Vakuf -0.150163441 0.119774205
## factor(Grad)Foča -0.048739889 0.101731613
## factor(Grad)Fojnica 0.291268073 0.108936165
## factor(Grad)Gacko 0.166326048 0.192236149
## factor(Grad)Goražde 0.154069404 0.084450944
## factor(Grad)Gornji Vakuf-Uskoplje 0.083413199 0.129212404
## factor(Grad)Gračanica 0.093656987 0.080485963
## factor(Grad)Gradačac 0.134379239 0.083438265
## factor(Grad)Gradiška 0.007134140 0.076144766
## factor(Grad)Hadžići 0.376020836 0.070943995
## factor(Grad)Ilidža 0.450822467 0.065152132
## factor(Grad)Ilijaš 0.399930153 0.076669567
## factor(Grad)Istočna Ilidža -0.001797591 0.079623412
## factor(Grad)Istočno Novo Sarajevo 0.107438565 0.072551319
## factor(Grad)Jablanica 0.360848449 0.122535079
## factor(Grad)Jajce 0.130553595 0.092889361
## factor(Grad)Kakanj 0.177639352 0.088669590
## factor(Grad)Kalesija -0.082943642 0.100581828
## factor(Grad)Kiseljak 0.223743182 0.089129372
## factor(Grad)Kladanj 0.217577458 0.177392949
## factor(Grad)Ključ -0.301365070 0.119634133
## factor(Grad)Kneževo -0.595962949 0.192329287
## factor(Grad)Konjic 0.438458032 0.085144502
## factor(Grad)Kotor Varoš -0.304575011 0.149501972
## factor(Grad)Kozarska Dubica -0.171734444 0.101847948
## factor(Grad)Kupres 0.037380089 0.133239474
## factor(Grad)Laktaši 0.335356551 0.080016665
## factor(Grad)Livno -0.161060995 0.149549243
## factor(Grad)Lopare 0.003805227 0.212795453
## factor(Grad)Lukavac 0.009331751 0.070889940
## factor(Grad)Ljubuški -0.127231668 0.138037035
## factor(Grad)Maglaj -0.079742300 0.104326665
## factor(Grad)Modriča -0.008969180 0.087682791
## factor(Grad)Mostar 0.500686993 0.068244915
## factor(Grad)Neum 0.797363775 0.129295161
## factor(Grad)Nevesinje -0.123455257 0.212629992
## factor(Grad)Novi Grad -0.006459067 0.093466586
## factor(Grad)Novi Travnik -0.105794583 0.088336686
## factor(Grad)Odžak 0.213380913 0.125782549
## factor(Grad)Olovo 0.293342583 0.137813944
## factor(Grad)Orašje 0.118540557 0.137966651
## factor(Grad)Pale - RS 0.024650278 0.079922971
## factor(Grad)Prijedor 0.024275700 0.072888476
## factor(Grad)Prnjavor 0.158572384 0.085522559
## factor(Grad)Ravno 0.228681165 0.149588200
## factor(Grad)Rogatica -0.014003421 0.133310271
## factor(Grad)Sanski Most 0.022412903 0.088111546
## factor(Grad)Sarajevo - Centar 0.696263136 0.065605251
## factor(Grad)Sarajevo - Novi Grad 0.400650048 0.064889793
## factor(Grad)Sarajevo - Novo Sarajevo 0.550293817 0.065492492
## factor(Grad)Sarajevo - Stari Grad 0.594268072 0.067477798
## factor(Grad)Sokolac -0.258936232 0.143143584
## factor(Grad)Srbac -0.049689971 0.105726404
## factor(Grad)Srebrenica -0.241660246 0.122615506
## factor(Grad)Srebrenik 0.177251827 0.078826656
## factor(Grad)Stolac -0.166173377 0.137795919
## factor(Grad)Šamac -0.485802359 0.166091272
## factor(Grad)Široki Brijeg 0.391146805 0.156905963
## factor(Grad)Teslić 0.216560491 0.104342440
## factor(Grad)Tešanj 0.251977175 0.095772704
## factor(Grad)Tomislavgrad 0.349045777 0.192306936
## factor(Grad)Travnik 0.206793698 0.077418787
## factor(Grad)Trebinje 0.436371443 0.082194455
## factor(Grad)Tuzla 0.231695105 0.065572088
## factor(Grad)Ugljevik 0.202382106 0.129334115
## factor(Grad)Usora -0.238131726 0.138264737
## factor(Grad)Vareš -0.709000127 0.192327707
## factor(Grad)Velika Kladuša 0.176858252 0.091428027
## factor(Grad)Visoko 0.287068647 0.080761730
## factor(Grad)Višegrad -0.254226117 0.110740564
## factor(Grad)Vitez -0.066817577 0.085932082
## factor(Grad)Vlasenica -0.304678570 0.133209788
## factor(Grad)Vogošća 0.362489801 0.067768599
## factor(Grad)Zavidovići -0.035170452 0.084431574
## factor(Grad)Zenica 0.211522423 0.066919767
## factor(Grad)Zvornik -0.209338981 0.080514165
## factor(Grad)Žepče -0.125598391 0.103063619
## factor(Grad)Živinice 0.103319816 0.072201158
## t value Pr(>|t|)
## (Intercept) 99.37064637 0.000000e+00
## factor(Y)2017 5.04340475 4.642354e-07
## as.factor(constr.y)1950-1959 -0.17305344 8.626124e-01
## as.factor(constr.y)1960-1969 -2.80807627 4.992171e-03
## as.factor(constr.y)1970-1979 -4.75503872 2.007967e-06
## as.factor(constr.y)1980-1989 -8.97126332 3.380092e-19
## as.factor(constr.y)1990-1999 -10.69041421 1.500294e-26
## as.factor(constr.y)2000-2009 -2.32317195 2.018712e-02
## as.factor(constr.y)2010+ -0.60983917 5.419803e-01
## as.factor(constr.y)Until 1950 1.45809666 1.448410e-01
## factor(prop_type)Stanovi 45.34643371 0.000000e+00
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -1.32046229 1.867068e-01
## factor(`Vrsta grijanja`)Centralno (Plin) 1.03921282 2.987274e-01
## factor(`Vrsta grijanja`)Drva -16.42056263 6.453818e-60
## factor(`Vrsta grijanja`)Ostalo -16.42677193 5.839493e-60
## factor(`Vrsta grijanja`)Plin -3.19421302 1.405876e-03
## factor(`Vrsta grijanja`)Struja -5.30240186 1.163890e-07
## factor(Grad)Banovići -0.58773688 5.567204e-01
## factor(Grad)Banja Luka 8.27321260 1.446203e-16
## factor(Grad)Bihać 6.13570498 8.754737e-10
## factor(Grad)Bijeljina 2.51221800 1.201096e-02
## factor(Grad)Bileća 1.66898727 9.514690e-02
## factor(Grad)Bosanski Petrovac -2.64767056 8.115810e-03
## factor(Grad)Bratunac -2.69918812 6.960962e-03
## factor(Grad)Brčko 1.39632476 1.626434e-01
## factor(Grad)Breza 1.33244434 1.827404e-01
## factor(Grad)Brod -2.08650247 3.695492e-02
## factor(Grad)Bugojno 0.97562907 3.292687e-01
## factor(Grad)Busovača 0.66122143 5.084835e-01
## factor(Grad)Cazin 2.70253889 6.891231e-03
## factor(Grad)Čapljina 2.55516159 1.062650e-02
## factor(Grad)Čelić -0.83509678 4.036804e-01
## factor(Grad)Čelinac 1.26880362 2.045366e-01
## factor(Grad)Derventa 0.09711928 9.226333e-01
## factor(Grad)Doboj -0.23600211 8.134352e-01
## factor(Grad)Doboj Jug 0.44872277 6.536400e-01
## factor(Grad)Donji Vakuf -1.25372104 2.099687e-01
## factor(Grad)Foča -0.47910269 6.318746e-01
## factor(Grad)Fojnica 2.67375001 7.511343e-03
## factor(Grad)Gacko 0.86521733 3.869374e-01
## factor(Grad)Goražde 1.82436568 6.812253e-02
## factor(Grad)Gornji Vakuf-Uskoplje 0.64555102 5.185829e-01
## factor(Grad)Gračanica 1.16364373 2.445924e-01
## factor(Grad)Gradačac 1.61052293 1.073110e-01
## factor(Grad)Gradiška 0.09369179 9.253556e-01
## factor(Grad)Hadžići 5.30024897 1.177665e-07
## factor(Grad)Ilidža 6.91953511 4.769947e-12
## factor(Grad)Ilijaš 5.21628295 1.856962e-07
## factor(Grad)Istočna Ilidža -0.02257616 9.819887e-01
## factor(Grad)Istočno Novo Sarajevo 1.48086300 1.386703e-01
## factor(Grad)Jablanica 2.94485833 3.237441e-03
## factor(Grad)Jajce 1.40547414 1.599071e-01
## factor(Grad)Kakanj 2.00338528 4.515918e-02
## factor(Grad)Kalesija -0.82463844 4.095939e-01
## factor(Grad)Kiseljak 2.51031928 1.207571e-02
## factor(Grad)Kladanj 1.22652822 2.200249e-01
## factor(Grad)Ključ -2.51905591 1.178033e-02
## factor(Grad)Kneževo -3.09865937 1.948634e-03
## factor(Grad)Konjic 5.14957541 2.653574e-07
## factor(Grad)Kotor Varoš -2.03726417 4.164623e-02
## factor(Grad)Kozarska Dubica -1.68618463 9.178715e-02
## factor(Grad)Kupres 0.28054816 7.790620e-01
## factor(Grad)Laktaši 4.19108382 2.796835e-05
## factor(Grad)Livno -1.07697633 2.815132e-01
## factor(Grad)Lopare 0.01788209 9.857332e-01
## factor(Grad)Lukavac 0.13163716 8.952736e-01
## factor(Grad)Ljubuški -0.92172124 3.566932e-01
## factor(Grad)Maglaj -0.76435205 4.446730e-01
## factor(Grad)Modriča -0.10229123 9.185273e-01
## factor(Grad)Mostar 7.33661973 2.336660e-13
## factor(Grad)Neum 6.16700400 7.191193e-10
## factor(Grad)Nevesinje -0.58061074 5.615141e-01
## factor(Grad)Novi Grad -0.06910563 9.449067e-01
## factor(Grad)Novi Travnik -1.19762907 2.310859e-01
## factor(Grad)Odžak 1.69642701 8.983188e-02
## factor(Grad)Olovo 2.12854065 3.331330e-02
## factor(Grad)Orašje 0.85919718 3.902495e-01
## factor(Grad)Pale - RS 0.30842545 7.577642e-01
## factor(Grad)Prijedor 0.33305265 7.391006e-01
## factor(Grad)Prnjavor 1.85415855 6.374185e-02
## factor(Grad)Ravno 1.52873799 1.263567e-01
## factor(Grad)Rogatica -0.10504382 9.163429e-01
## factor(Grad)Sanski Most 0.25436964 7.992145e-01
## factor(Grad)Sarajevo - Centar 10.61291773 3.421954e-26
## factor(Grad)Sarajevo - Novi Grad 6.17431538 6.867246e-10
## factor(Grad)Sarajevo - Novo Sarajevo 8.40239540 4.883491e-17
## factor(Grad)Sarajevo - Stari Grad 8.80686824 1.469136e-18
## factor(Grad)Sokolac -1.80892657 7.048830e-02
## factor(Grad)Srbac -0.46998639 6.383736e-01
## factor(Grad)Srebrenica -1.97087835 4.876154e-02
## factor(Grad)Srebrenik 2.24862800 2.455487e-02
## factor(Grad)Stolac -1.20593831 2.278659e-01
## factor(Grad)Šamac -2.92491202 3.452211e-03
## factor(Grad)Široki Brijeg 2.49287406 1.268522e-02
## factor(Grad)Teslić 2.07547850 3.796420e-02
## factor(Grad)Tešanj 2.63099155 8.524875e-03
## factor(Grad)Tomislavgrad 1.81504518 6.954279e-02
## factor(Grad)Travnik 2.67110487 7.570757e-03
## factor(Grad)Trebinje 5.30901319 1.122556e-07
## factor(Grad)Tuzla 3.53344101 4.117836e-04
## factor(Grad)Ugljevik 1.56480064 1.176569e-01
## factor(Grad)Usora -1.72228821 8.504399e-02
## factor(Grad)Vareš -3.68641700 2.284709e-04
## factor(Grad)Velika Kladuša 1.93439865 5.308843e-02
## factor(Grad)Visoko 3.55451336 3.801858e-04
## factor(Grad)Višegrad -2.29569100 2.171131e-02
## factor(Grad)Vitez -0.77756265 4.368427e-01
## factor(Grad)Vlasenica -2.28720858 2.220158e-02
## factor(Grad)Vogošća 5.34893454 9.015569e-08
## factor(Grad)Zavidovići -0.41655568 6.770111e-01
## factor(Grad)Zenica 3.16083619 1.577217e-03
## factor(Grad)Zvornik -2.60002673 9.333469e-03
## factor(Grad)Žepče -1.21864915 2.230022e-01
## factor(Grad)Živinice 1.43099943 1.524573e-01
model.coef1[,c("Estimate","Pr(>|t|)")]
## Estimate Pr(>|t|)
## (Intercept) 6.565634109 0.000000e+00
## factor(Y)2017 0.038164178 4.642354e-07
## as.factor(constr.y)1950-1959 -0.005747513 8.626124e-01
## as.factor(constr.y)1960-1969 -0.048890339 4.992171e-03
## as.factor(constr.y)1970-1979 -0.070126662 2.007967e-06
## as.factor(constr.y)1980-1989 -0.115088151 3.380092e-19
## as.factor(constr.y)1990-1999 -0.182284966 1.500294e-26
## as.factor(constr.y)2000-2009 -0.031534135 2.018712e-02
## as.factor(constr.y)2010+ -0.008948840 5.419803e-01
## as.factor(constr.y)Until 1950 0.037934485 1.448410e-01
## factor(prop_type)Stanovi 0.509133201 0.000000e+00
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.019944344 1.867068e-01
## factor(`Vrsta grijanja`)Centralno (Plin) 0.017525386 2.987274e-01
## factor(`Vrsta grijanja`)Drva -0.265391897 6.453818e-60
## factor(`Vrsta grijanja`)Ostalo -0.282873130 5.839493e-60
## factor(`Vrsta grijanja`)Plin -0.054937554 1.405876e-03
## factor(`Vrsta grijanja`)Struja -0.079544069 1.163890e-07
## factor(Grad)Banovići -0.051878148 5.567204e-01
## factor(Grad)Banja Luka 0.535288888 1.446203e-16
## factor(Grad)Bihać 0.425798791 8.754737e-10
## factor(Grad)Bijeljina 0.169299574 1.201096e-02
## factor(Grad)Bileća 0.277358440 9.514690e-02
## factor(Grad)Bosanski Petrovac -0.378999119 8.115810e-03
## factor(Grad)Bratunac -0.281497145 6.960962e-03
## factor(Grad)Brčko 0.094179820 1.626434e-01
## factor(Grad)Breza 0.125523320 1.827404e-01
## factor(Grad)Brod -0.201915039 3.695492e-02
## factor(Grad)Bugojno 0.077311437 3.292687e-01
## factor(Grad)Busovača 0.062251132 5.084835e-01
## factor(Grad)Cazin 0.236929143 6.891231e-03
## factor(Grad)Čapljina 0.247099360 1.062650e-02
## factor(Grad)Čelić -0.104948850 4.036804e-01
## factor(Grad)Čelinac 0.155445020 2.045366e-01
## factor(Grad)Derventa 0.008443602 9.226333e-01
## factor(Grad)Doboj -0.016383099 8.134352e-01
## factor(Grad)Doboj Jug 0.058037723 6.536400e-01
## factor(Grad)Donji Vakuf -0.150163441 2.099687e-01
## factor(Grad)Foča -0.048739889 6.318746e-01
## factor(Grad)Fojnica 0.291268073 7.511343e-03
## factor(Grad)Gacko 0.166326048 3.869374e-01
## factor(Grad)Goražde 0.154069404 6.812253e-02
## factor(Grad)Gornji Vakuf-Uskoplje 0.083413199 5.185829e-01
## factor(Grad)Gračanica 0.093656987 2.445924e-01
## factor(Grad)Gradačac 0.134379239 1.073110e-01
## factor(Grad)Gradiška 0.007134140 9.253556e-01
## factor(Grad)Hadžići 0.376020836 1.177665e-07
## factor(Grad)Ilidža 0.450822467 4.769947e-12
## factor(Grad)Ilijaš 0.399930153 1.856962e-07
## factor(Grad)Istočna Ilidža -0.001797591 9.819887e-01
## factor(Grad)Istočno Novo Sarajevo 0.107438565 1.386703e-01
## factor(Grad)Jablanica 0.360848449 3.237441e-03
## factor(Grad)Jajce 0.130553595 1.599071e-01
## factor(Grad)Kakanj 0.177639352 4.515918e-02
## factor(Grad)Kalesija -0.082943642 4.095939e-01
## factor(Grad)Kiseljak 0.223743182 1.207571e-02
## factor(Grad)Kladanj 0.217577458 2.200249e-01
## factor(Grad)Ključ -0.301365070 1.178033e-02
## factor(Grad)Kneževo -0.595962949 1.948634e-03
## factor(Grad)Konjic 0.438458032 2.653574e-07
## factor(Grad)Kotor Varoš -0.304575011 4.164623e-02
## factor(Grad)Kozarska Dubica -0.171734444 9.178715e-02
## factor(Grad)Kupres 0.037380089 7.790620e-01
## factor(Grad)Laktaši 0.335356551 2.796835e-05
## factor(Grad)Livno -0.161060995 2.815132e-01
## factor(Grad)Lopare 0.003805227 9.857332e-01
## factor(Grad)Lukavac 0.009331751 8.952736e-01
## factor(Grad)Ljubuški -0.127231668 3.566932e-01
## factor(Grad)Maglaj -0.079742300 4.446730e-01
## factor(Grad)Modriča -0.008969180 9.185273e-01
## factor(Grad)Mostar 0.500686993 2.336660e-13
## factor(Grad)Neum 0.797363775 7.191193e-10
## factor(Grad)Nevesinje -0.123455257 5.615141e-01
## factor(Grad)Novi Grad -0.006459067 9.449067e-01
## factor(Grad)Novi Travnik -0.105794583 2.310859e-01
## factor(Grad)Odžak 0.213380913 8.983188e-02
## factor(Grad)Olovo 0.293342583 3.331330e-02
## factor(Grad)Orašje 0.118540557 3.902495e-01
## factor(Grad)Pale - RS 0.024650278 7.577642e-01
## factor(Grad)Prijedor 0.024275700 7.391006e-01
## factor(Grad)Prnjavor 0.158572384 6.374185e-02
## factor(Grad)Ravno 0.228681165 1.263567e-01
## factor(Grad)Rogatica -0.014003421 9.163429e-01
## factor(Grad)Sanski Most 0.022412903 7.992145e-01
## factor(Grad)Sarajevo - Centar 0.696263136 3.421954e-26
## factor(Grad)Sarajevo - Novi Grad 0.400650048 6.867246e-10
## factor(Grad)Sarajevo - Novo Sarajevo 0.550293817 4.883491e-17
## factor(Grad)Sarajevo - Stari Grad 0.594268072 1.469136e-18
## factor(Grad)Sokolac -0.258936232 7.048830e-02
## factor(Grad)Srbac -0.049689971 6.383736e-01
## factor(Grad)Srebrenica -0.241660246 4.876154e-02
## factor(Grad)Srebrenik 0.177251827 2.455487e-02
## factor(Grad)Stolac -0.166173377 2.278659e-01
## factor(Grad)Šamac -0.485802359 3.452211e-03
## factor(Grad)Široki Brijeg 0.391146805 1.268522e-02
## factor(Grad)Teslić 0.216560491 3.796420e-02
## factor(Grad)Tešanj 0.251977175 8.524875e-03
## factor(Grad)Tomislavgrad 0.349045777 6.954279e-02
## factor(Grad)Travnik 0.206793698 7.570757e-03
## factor(Grad)Trebinje 0.436371443 1.122556e-07
## factor(Grad)Tuzla 0.231695105 4.117836e-04
## factor(Grad)Ugljevik 0.202382106 1.176569e-01
## factor(Grad)Usora -0.238131726 8.504399e-02
## factor(Grad)Vareš -0.709000127 2.284709e-04
## factor(Grad)Velika Kladuša 0.176858252 5.308843e-02
## factor(Grad)Visoko 0.287068647 3.801858e-04
## factor(Grad)Višegrad -0.254226117 2.171131e-02
## factor(Grad)Vitez -0.066817577 4.368427e-01
## factor(Grad)Vlasenica -0.304678570 2.220158e-02
## factor(Grad)Vogošća 0.362489801 9.015569e-08
## factor(Grad)Zavidovići -0.035170452 6.770111e-01
## factor(Grad)Zenica 0.211522423 1.577217e-03
## factor(Grad)Zvornik -0.209338981 9.333469e-03
## factor(Grad)Žepče -0.125598391 2.230022e-01
## factor(Grad)Živinice 0.103319816 1.524573e-01
model.coef1 <- model.coef1
model.coef1 <- as.data.frame.array(model.coef1[-c(1:17),]) # da bude dataframe a ne matrix i da mi brise prvih 8 redova koji su drugi koeficijenti
Munic <- rownames(model.coef1)
model.coef1$Munic <- Munic
dim(model.coef1) #da vidim koliko redova
## [1] 97 5
rownames(model.coef1) <- 1:97 #nedostaje ref.opstina
head(model.coef1)
## Estimate Std. Error t value Pr(>|t|) Munic
## 1 -0.05187815 0.08826764 -0.5877369 5.567204e-01 factor(Grad)Banovići
## 2 0.53528889 0.06470145 8.2732126 1.446203e-16 factor(Grad)Banja Luka
## 3 0.42579879 0.06939688 6.1357050 8.754737e-10 factor(Grad)Bihać
## 4 0.16929957 0.06739048 2.5122180 1.201096e-02 factor(Grad)Bijeljina
## 5 0.27735844 0.16618368 1.6689873 9.514690e-02 factor(Grad)Bileća
## 6 -0.37899912 0.14314436 -2.6476706 8.115810e-03 factor(Grad)Bosanski Petrovac
### Make a new variable with municipal names without `factor(Munic)`
Munic1 <- sapply(strsplit(Munic, split=')', fixed=TRUE), function(x) (x[2]))
Munic1
## [1] "Banovići" "Banja Luka"
## [3] "Bihać" "Bijeljina"
## [5] "Bileća" "Bosanski Petrovac"
## [7] "Bratunac" "Brčko"
## [9] "Breza" "Brod"
## [11] "Bugojno" "Busovača"
## [13] "Cazin" "Čapljina"
## [15] "Čelić" "Čelinac"
## [17] "Derventa" "Doboj"
## [19] "Doboj Jug" "Donji Vakuf"
## [21] "Foča" "Fojnica"
## [23] "Gacko" "Goražde"
## [25] "Gornji Vakuf-Uskoplje" "Gračanica"
## [27] "Gradačac" "Gradiška"
## [29] "Hadžići" "Ilidža"
## [31] "Ilijaš" "Istočna Ilidža"
## [33] "Istočno Novo Sarajevo" "Jablanica"
## [35] "Jajce" "Kakanj"
## [37] "Kalesija" "Kiseljak"
## [39] "Kladanj" "Ključ"
## [41] "Kneževo" "Konjic"
## [43] "Kotor Varoš" "Kozarska Dubica"
## [45] "Kupres" "Laktaši"
## [47] "Livno" "Lopare"
## [49] "Lukavac" "Ljubuški"
## [51] "Maglaj" "Modriča"
## [53] "Mostar" "Neum"
## [55] "Nevesinje" "Novi Grad"
## [57] "Novi Travnik" "Odžak"
## [59] "Olovo" "Orašje"
## [61] "Pale - RS" "Prijedor"
## [63] "Prnjavor" "Ravno"
## [65] "Rogatica" "Sanski Most"
## [67] "Sarajevo - Centar" "Sarajevo - Novi Grad"
## [69] "Sarajevo - Novo Sarajevo" "Sarajevo - Stari Grad"
## [71] "Sokolac" "Srbac"
## [73] "Srebrenica" "Srebrenik"
## [75] "Stolac" "Šamac"
## [77] "Široki Brijeg" "Teslić"
## [79] "Tešanj" "Tomislavgrad"
## [81] "Travnik" "Trebinje"
## [83] "Tuzla" "Ugljevik"
## [85] "Usora" "Vareš"
## [87] "Velika Kladuša" "Visoko"
## [89] "Višegrad" "Vitez"
## [91] "Vlasenica" "Vogošća"
## [93] "Zavidovići" "Zenica"
## [95] "Zvornik" "Žepče"
## [97] "Živinice"
model.coef1$Munic1 <- Munic1
#TREBA DA JE B.KRUPA
model.coef1[nrow(model.coef1) + 1,] = c(0, NA, NA, NA,"Bosanska Krupa","Bosanska Krupa" )
model.coef1$Estimate <- as.numeric(model.coef1$Estimate)
top50.m1 <- model.coef1 %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
head(50) #top ten municipalities
top50.m1
## Estimate Munic1
## 1 0.79736378 Neum
## 2 0.69626314 Sarajevo - Centar
## 3 0.59426807 Sarajevo - Stari Grad
## 4 0.55029382 Sarajevo - Novo Sarajevo
## 5 0.53528889 Banja Luka
## 6 0.50068699 Mostar
## 7 0.45082247 Ilidža
## 8 0.43845803 Konjic
## 9 0.43637144 Trebinje
## 10 0.42579879 Bihać
## 11 0.40065005 Sarajevo - Novi Grad
## 12 0.39993015 Ilijaš
## 13 0.39114681 Široki Brijeg
## 14 0.37602084 Hadžići
## 15 0.36248980 Vogošća
## 16 0.36084845 Jablanica
## 17 0.34904578 Tomislavgrad
## 18 0.33535655 Laktaši
## 19 0.29334258 Olovo
## 20 0.29126807 Fojnica
## 21 0.28706865 Visoko
## 22 0.27735844 Bileća
## 23 0.25197717 Tešanj
## 24 0.24709936 Čapljina
## 25 0.23692914 Cazin
## 26 0.23169510 Tuzla
## 27 0.22868116 Ravno
## 28 0.22374318 Kiseljak
## 29 0.21757746 Kladanj
## 30 0.21656049 Teslić
## 31 0.21338091 Odžak
## 32 0.21152242 Zenica
## 33 0.20679370 Travnik
## 34 0.20238211 Ugljevik
## 35 0.17763935 Kakanj
## 36 0.17725183 Srebrenik
## 37 0.17685825 Velika Kladuša
## 38 0.16929957 Bijeljina
## 39 0.16632605 Gacko
## 40 0.15857238 Prnjavor
## 41 0.15544502 Čelinac
## 42 0.15406940 Goražde
## 43 0.13437924 Gradačac
## 44 0.13055359 Jajce
## 45 0.12552332 Breza
## 46 0.11854056 Orašje
## 47 0.10743856 Istočno Novo Sarajevo
## 48 0.10331982 Živinice
## 49 0.09417982 Brčko
## 50 0.09365699 Gračanica
bottom50.m1 <- model.coef1 %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
tail(50) #top ten municipalities
bottom50.m1
## Estimate Munic1
## 49 0.094179820 Brčko
## 50 0.093656987 Gračanica
## 51 0.083413199 Gornji Vakuf-Uskoplje
## 52 0.077311437 Bugojno
## 53 0.062251132 Busovača
## 54 0.058037723 Doboj Jug
## 55 0.037380089 Kupres
## 56 0.024650278 Pale - RS
## 57 0.024275700 Prijedor
## 58 0.022412903 Sanski Most
## 59 0.009331751 Lukavac
## 60 0.008443602 Derventa
## 61 0.007134140 Gradiška
## 62 0.003805227 Lopare
## 63 0.000000000 Bosanska Krupa
## 64 -0.001797591 Istočna Ilidža
## 65 -0.006459067 Novi Grad
## 66 -0.008969180 Modriča
## 67 -0.014003421 Rogatica
## 68 -0.016383099 Doboj
## 69 -0.035170452 Zavidovići
## 70 -0.048739889 Foča
## 71 -0.049689971 Srbac
## 72 -0.051878148 Banovići
## 73 -0.066817577 Vitez
## 74 -0.079742300 Maglaj
## 75 -0.082943642 Kalesija
## 76 -0.104948850 Čelić
## 77 -0.105794583 Novi Travnik
## 78 -0.123455257 Nevesinje
## 79 -0.125598391 Žepče
## 80 -0.127231668 Ljubuški
## 81 -0.150163441 Donji Vakuf
## 82 -0.161060995 Livno
## 83 -0.166173377 Stolac
## 84 -0.171734444 Kozarska Dubica
## 85 -0.201915039 Brod
## 86 -0.209338981 Zvornik
## 87 -0.238131726 Usora
## 88 -0.241660246 Srebrenica
## 89 -0.254226117 Višegrad
## 90 -0.258936232 Sokolac
## 91 -0.281497145 Bratunac
## 92 -0.301365070 Ključ
## 93 -0.304575011 Kotor Varoš
## 94 -0.304678570 Vlasenica
## 95 -0.378999119 Bosanski Petrovac
## 96 -0.485802359 Šamac
## 97 -0.595962949 Kneževo
## 98 -0.709000127 Vareš
#######################
#MODEL 1a SAMO STANOVI#
#######################
#2.4. 2024 sam sama ovaj model uradila koji ukljucuje samo stanove jer su mi malo kuce problematicne. Mozda je neko uracunavao i zemlju koja nije svaki put unesena.
#Ovo jos nisam provjerila sa Ismirom
model1a <- lm (log(p.sqm) ~ factor (Y) + as.factor(constr.y) + factor (`Vrsta grijanja`) + factor (Grad), data = dist.s.k1 %>% dplyr::filter(prop_type=="Stanovi"))
summary(model1a)$coefficients
## Estimate Std. Error
## (Intercept) 7.2060359341 0.065929883
## factor(Y)2017 0.0387259753 0.004715369
## as.factor(constr.y)1950-1959 -0.1160020113 0.018205932
## as.factor(constr.y)1960-1969 -0.1293911274 0.010224904
## as.factor(constr.y)1970-1979 -0.1751977725 0.008840236
## as.factor(constr.y)1980-1989 -0.1582435142 0.007786950
## as.factor(constr.y)1990-1999 -0.1718097199 0.012611046
## as.factor(constr.y)2000-2009 -0.1188278675 0.009026438
## as.factor(constr.y)2010+ -0.0435002002 0.008926697
## as.factor(constr.y)Until 1950 -0.0603167263 0.014458008
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.0300345690 0.009582765
## factor(`Vrsta grijanja`)Centralno (Plin) 0.0050348408 0.009180818
## factor(`Vrsta grijanja`)Drva -0.2169975562 0.015439336
## factor(`Vrsta grijanja`)Ostalo -0.1433315730 0.012595917
## factor(`Vrsta grijanja`)Plin -0.0492288248 0.009634789
## factor(`Vrsta grijanja`)Struja -0.0638718468 0.008425495
## factor(Grad)Banovići -0.1541401379 0.077080368
## factor(Grad)Banja Luka 0.4575410514 0.065788559
## factor(Grad)Bihać 0.3374652637 0.068701228
## factor(Grad)Bijeljina 0.0660492728 0.067128899
## factor(Grad)Bileća -0.0235962604 0.103023335
## factor(Grad)Bosanski Petrovac -0.5960722977 0.205743024
## factor(Grad)Bratunac -0.2245573582 0.109158644
## factor(Grad)Brčko -0.1022770522 0.066826192
## factor(Grad)Breza -0.1919830088 0.095102025
## factor(Grad)Brod -0.3679017079 0.078675297
## factor(Grad)Bugojno -0.0318434580 0.078441234
## factor(Grad)Busovača 0.0467363448 0.130142019
## factor(Grad)Cazin 0.2092211373 0.082532068
## factor(Grad)Čapljina 0.0579547653 0.089801218
## factor(Grad)Čelić -0.2328763561 0.117496721
## factor(Grad)Čelinac 0.0180991377 0.108896561
## factor(Grad)Derventa -0.2001379647 0.086451488
## factor(Grad)Doboj -0.0290847219 0.067911982
## factor(Grad)Donji Vakuf -0.2105721778 0.130296355
## factor(Grad)Foča -0.0116082937 0.098639570
## factor(Grad)Fojnica 0.0474992963 0.109157695
## factor(Grad)Gacko 0.0330174406 0.152496037
## factor(Grad)Goražde 0.0293289491 0.081397318
## factor(Grad)Gračanica -0.0631873072 0.086231493
## factor(Grad)Gradačac -0.0393205413 0.086082032
## factor(Grad)Gradiška -0.1196921913 0.074744248
## factor(Grad)Hadžići 0.0574511566 0.074592606
## factor(Grad)Ilidža 0.2980243698 0.065995670
## factor(Grad)Ilijaš 0.1302953810 0.080505460
## factor(Grad)Istočna Ilidža -0.1045076349 0.070931266
## factor(Grad)Istočno Novo Sarajevo -0.0188888226 0.068017395
## factor(Grad)Jablanica 0.1241168642 0.130171719
## factor(Grad)Jajce -0.0027484917 0.077790599
## factor(Grad)Kakanj 0.1082168392 0.081608922
## factor(Grad)Kalesija 0.0254937947 0.102909824
## factor(Grad)Kiseljak 0.1174945895 0.098645962
## factor(Grad)Kladanj -0.1720668044 0.152599181
## factor(Grad)Ključ 0.0845066595 0.206066869
## factor(Grad)Kneževo -0.4365122080 0.205900609
## factor(Grad)Konjic 0.0324874210 0.079805183
## factor(Grad)Kotor Varoš -0.0020597362 0.117370246
## factor(Grad)Kozarska Dubica -0.1946789787 0.098640807
## factor(Grad)Laktaši 0.2082458887 0.084688064
## factor(Grad)Livno -0.0006227113 0.130369156
## factor(Grad)Lopare -0.3151420072 0.130361643
## factor(Grad)Lukavac -0.0300477918 0.069342148
## factor(Grad)Ljubuški -0.2339008537 0.103118905
## factor(Grad)Maglaj -0.0218034126 0.130329917
## factor(Grad)Modriča -0.1585112160 0.077367952
## factor(Grad)Mostar 0.4331837905 0.066876410
## factor(Grad)Neum 0.7143208297 0.117440066
## factor(Grad)Nevesinje -0.0251903267 0.152509138
## factor(Grad)Novi Grad -0.0340491206 0.087704893
## factor(Grad)Novi Travnik -0.1459389154 0.074624952
## factor(Grad)Orašje -0.0462574363 0.098445364
## factor(Grad)Pale - RS -0.0552075378 0.072187012
## factor(Grad)Prijedor -0.0325032318 0.070359844
## factor(Grad)Prnjavor -0.0026287066 0.081509666
## factor(Grad)Rogatica -0.3475300957 0.109069151
## factor(Grad)Sanski Most -0.0100851070 0.083398565
## factor(Grad)Sarajevo - Centar 0.6302726762 0.065993510
## factor(Grad)Sarajevo - Novi Grad 0.3186763015 0.065786241
## factor(Grad)Sarajevo - Novo Sarajevo 0.5019179732 0.065966740
## factor(Grad)Sarajevo - Stari Grad 0.5579441600 0.067212941
## factor(Grad)Sokolac -0.3443857447 0.117313526
## factor(Grad)Srbac 0.0590790045 0.152647512
## factor(Grad)Srebrenica -0.3201266802 0.095054374
## factor(Grad)Srebrenik 0.0389474289 0.084882145
## factor(Grad)Stolac -0.2822916212 0.205986032
## factor(Grad)Široki Brijeg 0.2554750626 0.205818103
## factor(Grad)Teslić -0.0304758692 0.089642955
## factor(Grad)Tešanj 0.1376101068 0.130320582
## factor(Grad)Tomislavgrad 0.1513583333 0.152656512
## factor(Grad)Travnik 0.0279439939 0.073620642
## factor(Grad)Trebinje 0.2427015529 0.072691895
## factor(Grad)Tuzla 0.1874563102 0.066220620
## factor(Grad)Ugljevik -0.2088883978 0.109155608
## factor(Grad)Vareš -0.4414955011 0.152747169
## factor(Grad)Velika Kladuša 0.3263577248 0.080494199
## factor(Grad)Visoko 0.1051121399 0.086314469
## factor(Grad)Višegrad -0.3549196956 0.117376351
## factor(Grad)Vitez -0.1316447660 0.086141001
## factor(Grad)Vlasenica -0.3944335654 0.130080387
## factor(Grad)Vogošća 0.1624299656 0.067156399
## factor(Grad)Zavidovići -0.0750488280 0.095082061
## factor(Grad)Zenica 0.1022831784 0.066695101
## factor(Grad)Zvornik -0.1807071130 0.074183728
## factor(Grad)Žepče -0.3971418375 0.152701740
## factor(Grad)Živinice 0.0785354641 0.072866035
## t value Pr(>|t|)
## (Intercept) 109.298479823 0.000000e+00
## factor(Y)2017 8.212713855 2.553923e-16
## as.factor(constr.y)1950-1959 -6.371660302 1.988676e-10
## as.factor(constr.y)1960-1969 -12.654508054 2.635943e-36
## as.factor(constr.y)1970-1979 -19.818223191 4.342929e-85
## as.factor(constr.y)1980-1989 -20.321629458 3.011330e-89
## as.factor(constr.y)1990-1999 -13.623748562 9.833920e-42
## as.factor(constr.y)2000-2009 -13.164425859 4.093536e-39
## as.factor(constr.y)2010+ -4.873045438 1.123112e-06
## as.factor(constr.y)Until 1950 -4.171855964 3.057856e-05
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -3.134227743 1.730202e-03
## factor(`Vrsta grijanja`)Centralno (Plin) 0.548408765 5.834288e-01
## factor(`Vrsta grijanja`)Drva -14.054849956 2.869007e-44
## factor(`Vrsta grijanja`)Ostalo -11.379208830 9.684585e-30
## factor(`Vrsta grijanja`)Plin -5.109486423 3.315929e-07
## factor(`Vrsta grijanja`)Struja -7.580782983 3.878504e-14
## factor(Grad)Banovići -1.999732775 4.556775e-02
## factor(Grad)Banja Luka 6.954720711 3.852485e-12
## factor(Grad)Bihać 4.912070325 9.216707e-07
## factor(Grad)Bijeljina 0.983917124 3.251903e-01
## factor(Grad)Bileća -0.229038018 8.188461e-01
## factor(Grad)Bosanski Petrovac -2.897168937 3.777158e-03
## factor(Grad)Bratunac -2.057165153 3.970735e-02
## factor(Grad)Brčko -1.530493499 1.259399e-01
## factor(Grad)Breza -2.018705783 4.355591e-02
## factor(Grad)Brod -4.676203626 2.977198e-06
## factor(Grad)Bugojno -0.405953047 6.847895e-01
## factor(Grad)Busovača 0.359118026 7.195176e-01
## factor(Grad)Cazin 2.535028421 1.126543e-02
## factor(Grad)Čapljina 0.645367253 5.187104e-01
## factor(Grad)Čelić -1.981981739 4.752038e-02
## factor(Grad)Čelinac 0.166204860 8.680005e-01
## factor(Grad)Derventa -2.315032041 2.063999e-02
## factor(Grad)Doboj -0.428270846 6.684671e-01
## factor(Grad)Donji Vakuf -1.616101834 1.061174e-01
## factor(Grad)Foča -0.117683945 9.063215e-01
## factor(Grad)Fojnica 0.435143817 6.634715e-01
## factor(Grad)Gacko 0.216513434 8.285939e-01
## factor(Grad)Goražde 0.360318374 7.186199e-01
## factor(Grad)Gračanica -0.732763689 4.637271e-01
## factor(Grad)Gradačac -0.456779893 6.478434e-01
## factor(Grad)Gradiška -1.601356555 1.093431e-01
## factor(Grad)Hadžići 0.770199084 4.412078e-01
## factor(Grad)Ilidža 4.515817041 6.410775e-06
## factor(Grad)Ilijaš 1.618466388 1.056072e-01
## factor(Grad)Istočna Ilidža -1.473364866 1.406977e-01
## factor(Grad)Istočno Novo Sarajevo -0.277705763 7.812464e-01
## factor(Grad)Jablanica 0.953485635 3.403770e-01
## factor(Grad)Jajce -0.035331926 9.718161e-01
## factor(Grad)Kakanj 1.326041768 1.848692e-01
## factor(Grad)Kalesija 0.247729455 8.043510e-01
## factor(Grad)Kiseljak 1.191073486 2.336652e-01
## factor(Grad)Kladanj -1.127573580 2.595387e-01
## factor(Grad)Ključ 0.410093383 6.817500e-01
## factor(Grad)Kneževo -2.120014168 3.403998e-02
## factor(Grad)Konjic 0.407084100 6.839587e-01
## factor(Grad)Kotor Varoš -0.017549049 9.859991e-01
## factor(Grad)Kozarska Dubica -1.973615023 4.846484e-02
## factor(Grad)Laktaši 2.458975680 1.395746e-02
## factor(Grad)Livno -0.004776523 9.961890e-01
## factor(Grad)Lopare -2.417444281 1.565534e-02
## factor(Grad)Lukavac -0.433326521 6.647909e-01
## factor(Grad)Ljubuški -2.268263543 2.334351e-02
## factor(Grad)Maglaj -0.167293996 8.671436e-01
## factor(Grad)Modriča -2.048796832 4.051924e-02
## factor(Grad)Mostar 6.477378087 9.966304e-11
## factor(Grad)Neum 6.082428719 1.246160e-09
## factor(Grad)Nevesinje -0.165172573 8.688129e-01
## factor(Grad)Novi Grad -0.388223727 6.978623e-01
## factor(Grad)Novi Travnik -1.955631617 5.054834e-02
## factor(Grad)Orašje -0.469879277 6.384559e-01
## factor(Grad)Pale - RS -0.764784919 4.444254e-01
## factor(Grad)Prijedor -0.461957133 6.441264e-01
## factor(Grad)Prnjavor -0.032250244 9.742734e-01
## factor(Grad)Rogatica -3.186328041 1.447248e-03
## factor(Grad)Sanski Most -0.120926625 9.037526e-01
## factor(Grad)Sarajevo - Centar 9.550525092 1.744234e-21
## factor(Grad)Sarajevo - Novi Grad 4.844117831 1.299124e-06
## factor(Grad)Sarajevo - Novo Sarajevo 7.608652119 3.133031e-14
## factor(Grad)Sarajevo - Stari Grad 8.301141932 1.226573e-16
## factor(Grad)Sokolac -2.935601334 3.339860e-03
## factor(Grad)Srbac 0.387028939 6.987466e-01
## factor(Grad)Srebrenica -3.367826936 7.617221e-04
## factor(Grad)Srebrenik 0.458841244 6.463624e-01
## factor(Grad)Stolac -1.370440598 1.705934e-01
## factor(Grad)Široki Brijeg 1.241266237 2.145491e-01
## factor(Grad)Teslić -0.339969484 7.338897e-01
## factor(Grad)Tešanj 1.055935334 2.910342e-01
## factor(Grad)Tomislavgrad 0.991496078 3.214777e-01
## factor(Grad)Travnik 0.379567379 7.042781e-01
## factor(Grad)Trebinje 3.338770492 8.459002e-04
## factor(Grad)Tuzla 2.830784581 4.656651e-03
## factor(Grad)Ugljevik -1.913675354 5.570244e-02
## factor(Grad)Vareš -2.890367814 3.859757e-03
## factor(Grad)Velika Kladuša 4.054425394 5.080326e-05
## factor(Grad)Visoko 1.217781219 2.233482e-01
## factor(Grad)Višegrad -3.023775176 2.505476e-03
## factor(Grad)Vitez -1.528247464 1.264963e-01
## factor(Grad)Vlasenica -3.032229341 2.436445e-03
## factor(Grad)Vogošća 2.418681875 1.560223e-02
## factor(Grad)Zavidovići -0.789305860 4.299600e-01
## factor(Grad)Zenica 1.533593584 1.251749e-01
## factor(Grad)Zvornik -2.435940052 1.487795e-02
## factor(Grad)Žepče -2.600768250 9.321109e-03
## factor(Grad)Živinice 1.077806205 2.811574e-01
summary(model1a) #16545 deleted due to missingness
##
## Call:
## lm(formula = log(p.sqm) ~ factor(Y) + as.factor(constr.y) + factor(`Vrsta grijanja`) +
## factor(Grad), data = dist.s.k1 %>% dplyr::filter(prop_type ==
## "Stanovi"))
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.97334 -0.10439 0.00743 0.11661 0.98441
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 7.2060359 0.0659299 109.298
## factor(Y)2017 0.0387260 0.0047154 8.213
## as.factor(constr.y)1950-1959 -0.1160020 0.0182059 -6.372
## as.factor(constr.y)1960-1969 -0.1293911 0.0102249 -12.655
## as.factor(constr.y)1970-1979 -0.1751978 0.0088402 -19.818
## as.factor(constr.y)1980-1989 -0.1582435 0.0077870 -20.322
## as.factor(constr.y)1990-1999 -0.1718097 0.0126110 -13.624
## as.factor(constr.y)2000-2009 -0.1188279 0.0090264 -13.164
## as.factor(constr.y)2010+ -0.0435002 0.0089267 -4.873
## as.factor(constr.y)Until 1950 -0.0603167 0.0144580 -4.172
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.0300346 0.0095828 -3.134
## factor(`Vrsta grijanja`)Centralno (Plin) 0.0050348 0.0091808 0.548
## factor(`Vrsta grijanja`)Drva -0.2169976 0.0154393 -14.055
## factor(`Vrsta grijanja`)Ostalo -0.1433316 0.0125959 -11.379
## factor(`Vrsta grijanja`)Plin -0.0492288 0.0096348 -5.109
## factor(`Vrsta grijanja`)Struja -0.0638718 0.0084255 -7.581
## factor(Grad)Banovići -0.1541401 0.0770804 -2.000
## factor(Grad)Banja Luka 0.4575411 0.0657886 6.955
## factor(Grad)Bihać 0.3374653 0.0687012 4.912
## factor(Grad)Bijeljina 0.0660493 0.0671289 0.984
## factor(Grad)Bileća -0.0235963 0.1030233 -0.229
## factor(Grad)Bosanski Petrovac -0.5960723 0.2057430 -2.897
## factor(Grad)Bratunac -0.2245574 0.1091586 -2.057
## factor(Grad)Brčko -0.1022771 0.0668262 -1.530
## factor(Grad)Breza -0.1919830 0.0951020 -2.019
## factor(Grad)Brod -0.3679017 0.0786753 -4.676
## factor(Grad)Bugojno -0.0318435 0.0784412 -0.406
## factor(Grad)Busovača 0.0467363 0.1301420 0.359
## factor(Grad)Cazin 0.2092211 0.0825321 2.535
## factor(Grad)Čapljina 0.0579548 0.0898012 0.645
## factor(Grad)Čelić -0.2328764 0.1174967 -1.982
## factor(Grad)Čelinac 0.0180991 0.1088966 0.166
## factor(Grad)Derventa -0.2001380 0.0864515 -2.315
## factor(Grad)Doboj -0.0290847 0.0679120 -0.428
## factor(Grad)Donji Vakuf -0.2105722 0.1302964 -1.616
## factor(Grad)Foča -0.0116083 0.0986396 -0.118
## factor(Grad)Fojnica 0.0474993 0.1091577 0.435
## factor(Grad)Gacko 0.0330174 0.1524960 0.217
## factor(Grad)Goražde 0.0293289 0.0813973 0.360
## factor(Grad)Gračanica -0.0631873 0.0862315 -0.733
## factor(Grad)Gradačac -0.0393205 0.0860820 -0.457
## factor(Grad)Gradiška -0.1196922 0.0747442 -1.601
## factor(Grad)Hadžići 0.0574512 0.0745926 0.770
## factor(Grad)Ilidža 0.2980244 0.0659957 4.516
## factor(Grad)Ilijaš 0.1302954 0.0805055 1.618
## factor(Grad)Istočna Ilidža -0.1045076 0.0709313 -1.473
## factor(Grad)Istočno Novo Sarajevo -0.0188888 0.0680174 -0.278
## factor(Grad)Jablanica 0.1241169 0.1301717 0.953
## factor(Grad)Jajce -0.0027485 0.0777906 -0.035
## factor(Grad)Kakanj 0.1082168 0.0816089 1.326
## factor(Grad)Kalesija 0.0254938 0.1029098 0.248
## factor(Grad)Kiseljak 0.1174946 0.0986460 1.191
## factor(Grad)Kladanj -0.1720668 0.1525992 -1.128
## factor(Grad)Ključ 0.0845067 0.2060669 0.410
## factor(Grad)Kneževo -0.4365122 0.2059006 -2.120
## factor(Grad)Konjic 0.0324874 0.0798052 0.407
## factor(Grad)Kotor Varoš -0.0020597 0.1173702 -0.018
## factor(Grad)Kozarska Dubica -0.1946790 0.0986408 -1.974
## factor(Grad)Laktaši 0.2082459 0.0846881 2.459
## factor(Grad)Livno -0.0006227 0.1303692 -0.005
## factor(Grad)Lopare -0.3151420 0.1303616 -2.417
## factor(Grad)Lukavac -0.0300478 0.0693421 -0.433
## factor(Grad)Ljubuški -0.2339009 0.1031189 -2.268
## factor(Grad)Maglaj -0.0218034 0.1303299 -0.167
## factor(Grad)Modriča -0.1585112 0.0773680 -2.049
## factor(Grad)Mostar 0.4331838 0.0668764 6.477
## factor(Grad)Neum 0.7143208 0.1174401 6.082
## factor(Grad)Nevesinje -0.0251903 0.1525091 -0.165
## factor(Grad)Novi Grad -0.0340491 0.0877049 -0.388
## factor(Grad)Novi Travnik -0.1459389 0.0746250 -1.956
## factor(Grad)Orašje -0.0462574 0.0984454 -0.470
## factor(Grad)Pale - RS -0.0552075 0.0721870 -0.765
## factor(Grad)Prijedor -0.0325032 0.0703598 -0.462
## factor(Grad)Prnjavor -0.0026287 0.0815097 -0.032
## factor(Grad)Rogatica -0.3475301 0.1090692 -3.186
## factor(Grad)Sanski Most -0.0100851 0.0833986 -0.121
## factor(Grad)Sarajevo - Centar 0.6302727 0.0659935 9.551
## factor(Grad)Sarajevo - Novi Grad 0.3186763 0.0657862 4.844
## factor(Grad)Sarajevo - Novo Sarajevo 0.5019180 0.0659667 7.609
## factor(Grad)Sarajevo - Stari Grad 0.5579442 0.0672129 8.301
## factor(Grad)Sokolac -0.3443857 0.1173135 -2.936
## factor(Grad)Srbac 0.0590790 0.1526475 0.387
## factor(Grad)Srebrenica -0.3201267 0.0950544 -3.368
## factor(Grad)Srebrenik 0.0389474 0.0848821 0.459
## factor(Grad)Stolac -0.2822916 0.2059860 -1.370
## factor(Grad)Široki Brijeg 0.2554751 0.2058181 1.241
## factor(Grad)Teslić -0.0304759 0.0896430 -0.340
## factor(Grad)Tešanj 0.1376101 0.1303206 1.056
## factor(Grad)Tomislavgrad 0.1513583 0.1526565 0.991
## factor(Grad)Travnik 0.0279440 0.0736206 0.380
## factor(Grad)Trebinje 0.2427016 0.0726919 3.339
## factor(Grad)Tuzla 0.1874563 0.0662206 2.831
## factor(Grad)Ugljevik -0.2088884 0.1091556 -1.914
## factor(Grad)Vareš -0.4414955 0.1527472 -2.890
## factor(Grad)Velika Kladuša 0.3263577 0.0804942 4.054
## factor(Grad)Visoko 0.1051121 0.0863145 1.218
## factor(Grad)Višegrad -0.3549197 0.1173764 -3.024
## factor(Grad)Vitez -0.1316448 0.0861410 -1.528
## factor(Grad)Vlasenica -0.3944336 0.1300804 -3.032
## factor(Grad)Vogošća 0.1624300 0.0671564 2.419
## factor(Grad)Zavidovići -0.0750488 0.0950821 -0.789
## factor(Grad)Zenica 0.1022832 0.0666951 1.534
## factor(Grad)Zvornik -0.1807071 0.0741837 -2.436
## factor(Grad)Žepče -0.3971418 0.1527017 -2.601
## factor(Grad)Živinice 0.0785355 0.0728660 1.078
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## factor(Y)2017 2.55e-16 ***
## as.factor(constr.y)1950-1959 1.99e-10 ***
## as.factor(constr.y)1960-1969 < 2e-16 ***
## as.factor(constr.y)1970-1979 < 2e-16 ***
## as.factor(constr.y)1980-1989 < 2e-16 ***
## as.factor(constr.y)1990-1999 < 2e-16 ***
## as.factor(constr.y)2000-2009 < 2e-16 ***
## as.factor(constr.y)2010+ 1.12e-06 ***
## as.factor(constr.y)Until 1950 3.06e-05 ***
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 0.001730 **
## factor(`Vrsta grijanja`)Centralno (Plin) 0.583429
## factor(`Vrsta grijanja`)Drva < 2e-16 ***
## factor(`Vrsta grijanja`)Ostalo < 2e-16 ***
## factor(`Vrsta grijanja`)Plin 3.32e-07 ***
## factor(`Vrsta grijanja`)Struja 3.88e-14 ***
## factor(Grad)Banovići 0.045568 *
## factor(Grad)Banja Luka 3.85e-12 ***
## factor(Grad)Bihać 9.22e-07 ***
## factor(Grad)Bijeljina 0.325190
## factor(Grad)Bileća 0.818846
## factor(Grad)Bosanski Petrovac 0.003777 **
## factor(Grad)Bratunac 0.039707 *
## factor(Grad)Brčko 0.125940
## factor(Grad)Breza 0.043556 *
## factor(Grad)Brod 2.98e-06 ***
## factor(Grad)Bugojno 0.684790
## factor(Grad)Busovača 0.719518
## factor(Grad)Cazin 0.011265 *
## factor(Grad)Čapljina 0.518710
## factor(Grad)Čelić 0.047520 *
## factor(Grad)Čelinac 0.868001
## factor(Grad)Derventa 0.020640 *
## factor(Grad)Doboj 0.668467
## factor(Grad)Donji Vakuf 0.106117
## factor(Grad)Foča 0.906321
## factor(Grad)Fojnica 0.663471
## factor(Grad)Gacko 0.828594
## factor(Grad)Goražde 0.718620
## factor(Grad)Gračanica 0.463727
## factor(Grad)Gradačac 0.647843
## factor(Grad)Gradiška 0.109343
## factor(Grad)Hadžići 0.441208
## factor(Grad)Ilidža 6.41e-06 ***
## factor(Grad)Ilijaš 0.105607
## factor(Grad)Istočna Ilidža 0.140698
## factor(Grad)Istočno Novo Sarajevo 0.781246
## factor(Grad)Jablanica 0.340377
## factor(Grad)Jajce 0.971816
## factor(Grad)Kakanj 0.184869
## factor(Grad)Kalesija 0.804351
## factor(Grad)Kiseljak 0.233665
## factor(Grad)Kladanj 0.259539
## factor(Grad)Ključ 0.681750
## factor(Grad)Kneževo 0.034040 *
## factor(Grad)Konjic 0.683959
## factor(Grad)Kotor Varoš 0.985999
## factor(Grad)Kozarska Dubica 0.048465 *
## factor(Grad)Laktaši 0.013957 *
## factor(Grad)Livno 0.996189
## factor(Grad)Lopare 0.015655 *
## factor(Grad)Lukavac 0.664791
## factor(Grad)Ljubuški 0.023344 *
## factor(Grad)Maglaj 0.867144
## factor(Grad)Modriča 0.040519 *
## factor(Grad)Mostar 9.97e-11 ***
## factor(Grad)Neum 1.25e-09 ***
## factor(Grad)Nevesinje 0.868813
## factor(Grad)Novi Grad 0.697862
## factor(Grad)Novi Travnik 0.050548 .
## factor(Grad)Orašje 0.638456
## factor(Grad)Pale - RS 0.444425
## factor(Grad)Prijedor 0.644126
## factor(Grad)Prnjavor 0.974273
## factor(Grad)Rogatica 0.001447 **
## factor(Grad)Sanski Most 0.903753
## factor(Grad)Sarajevo - Centar < 2e-16 ***
## factor(Grad)Sarajevo - Novi Grad 1.30e-06 ***
## factor(Grad)Sarajevo - Novo Sarajevo 3.13e-14 ***
## factor(Grad)Sarajevo - Stari Grad < 2e-16 ***
## factor(Grad)Sokolac 0.003340 **
## factor(Grad)Srbac 0.698747
## factor(Grad)Srebrenica 0.000762 ***
## factor(Grad)Srebrenik 0.646362
## factor(Grad)Stolac 0.170593
## factor(Grad)Široki Brijeg 0.214549
## factor(Grad)Teslić 0.733890
## factor(Grad)Tešanj 0.291034
## factor(Grad)Tomislavgrad 0.321478
## factor(Grad)Travnik 0.704278
## factor(Grad)Trebinje 0.000846 ***
## factor(Grad)Tuzla 0.004657 **
## factor(Grad)Ugljevik 0.055702 .
## factor(Grad)Vareš 0.003860 **
## factor(Grad)Velika Kladuša 5.08e-05 ***
## factor(Grad)Visoko 0.223348
## factor(Grad)Višegrad 0.002505 **
## factor(Grad)Vitez 0.126496
## factor(Grad)Vlasenica 0.002436 **
## factor(Grad)Vogošća 0.015602 *
## factor(Grad)Zavidovići 0.429960
## factor(Grad)Zenica 0.125175
## factor(Grad)Zvornik 0.014878 *
## factor(Grad)Žepče 0.009321 **
## factor(Grad)Živinice 0.281157
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.195 on 6992 degrees of freedom
## (10819 observations deleted due to missingness)
## Multiple R-squared: 0.631, Adjusted R-squared: 0.6255
## F-statistic: 114.9 on 104 and 6992 DF, p-value: < 2.2e-16
model.coef1a <- summary(model1a)$coefficients #
model.coef1a # vidi koliko redova su koef.za druge varijable (16redova) i Opstina Bosanska Krupa = Ref
## Estimate Std. Error
## (Intercept) 7.2060359341 0.065929883
## factor(Y)2017 0.0387259753 0.004715369
## as.factor(constr.y)1950-1959 -0.1160020113 0.018205932
## as.factor(constr.y)1960-1969 -0.1293911274 0.010224904
## as.factor(constr.y)1970-1979 -0.1751977725 0.008840236
## as.factor(constr.y)1980-1989 -0.1582435142 0.007786950
## as.factor(constr.y)1990-1999 -0.1718097199 0.012611046
## as.factor(constr.y)2000-2009 -0.1188278675 0.009026438
## as.factor(constr.y)2010+ -0.0435002002 0.008926697
## as.factor(constr.y)Until 1950 -0.0603167263 0.014458008
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.0300345690 0.009582765
## factor(`Vrsta grijanja`)Centralno (Plin) 0.0050348408 0.009180818
## factor(`Vrsta grijanja`)Drva -0.2169975562 0.015439336
## factor(`Vrsta grijanja`)Ostalo -0.1433315730 0.012595917
## factor(`Vrsta grijanja`)Plin -0.0492288248 0.009634789
## factor(`Vrsta grijanja`)Struja -0.0638718468 0.008425495
## factor(Grad)Banovići -0.1541401379 0.077080368
## factor(Grad)Banja Luka 0.4575410514 0.065788559
## factor(Grad)Bihać 0.3374652637 0.068701228
## factor(Grad)Bijeljina 0.0660492728 0.067128899
## factor(Grad)Bileća -0.0235962604 0.103023335
## factor(Grad)Bosanski Petrovac -0.5960722977 0.205743024
## factor(Grad)Bratunac -0.2245573582 0.109158644
## factor(Grad)Brčko -0.1022770522 0.066826192
## factor(Grad)Breza -0.1919830088 0.095102025
## factor(Grad)Brod -0.3679017079 0.078675297
## factor(Grad)Bugojno -0.0318434580 0.078441234
## factor(Grad)Busovača 0.0467363448 0.130142019
## factor(Grad)Cazin 0.2092211373 0.082532068
## factor(Grad)Čapljina 0.0579547653 0.089801218
## factor(Grad)Čelić -0.2328763561 0.117496721
## factor(Grad)Čelinac 0.0180991377 0.108896561
## factor(Grad)Derventa -0.2001379647 0.086451488
## factor(Grad)Doboj -0.0290847219 0.067911982
## factor(Grad)Donji Vakuf -0.2105721778 0.130296355
## factor(Grad)Foča -0.0116082937 0.098639570
## factor(Grad)Fojnica 0.0474992963 0.109157695
## factor(Grad)Gacko 0.0330174406 0.152496037
## factor(Grad)Goražde 0.0293289491 0.081397318
## factor(Grad)Gračanica -0.0631873072 0.086231493
## factor(Grad)Gradačac -0.0393205413 0.086082032
## factor(Grad)Gradiška -0.1196921913 0.074744248
## factor(Grad)Hadžići 0.0574511566 0.074592606
## factor(Grad)Ilidža 0.2980243698 0.065995670
## factor(Grad)Ilijaš 0.1302953810 0.080505460
## factor(Grad)Istočna Ilidža -0.1045076349 0.070931266
## factor(Grad)Istočno Novo Sarajevo -0.0188888226 0.068017395
## factor(Grad)Jablanica 0.1241168642 0.130171719
## factor(Grad)Jajce -0.0027484917 0.077790599
## factor(Grad)Kakanj 0.1082168392 0.081608922
## factor(Grad)Kalesija 0.0254937947 0.102909824
## factor(Grad)Kiseljak 0.1174945895 0.098645962
## factor(Grad)Kladanj -0.1720668044 0.152599181
## factor(Grad)Ključ 0.0845066595 0.206066869
## factor(Grad)Kneževo -0.4365122080 0.205900609
## factor(Grad)Konjic 0.0324874210 0.079805183
## factor(Grad)Kotor Varoš -0.0020597362 0.117370246
## factor(Grad)Kozarska Dubica -0.1946789787 0.098640807
## factor(Grad)Laktaši 0.2082458887 0.084688064
## factor(Grad)Livno -0.0006227113 0.130369156
## factor(Grad)Lopare -0.3151420072 0.130361643
## factor(Grad)Lukavac -0.0300477918 0.069342148
## factor(Grad)Ljubuški -0.2339008537 0.103118905
## factor(Grad)Maglaj -0.0218034126 0.130329917
## factor(Grad)Modriča -0.1585112160 0.077367952
## factor(Grad)Mostar 0.4331837905 0.066876410
## factor(Grad)Neum 0.7143208297 0.117440066
## factor(Grad)Nevesinje -0.0251903267 0.152509138
## factor(Grad)Novi Grad -0.0340491206 0.087704893
## factor(Grad)Novi Travnik -0.1459389154 0.074624952
## factor(Grad)Orašje -0.0462574363 0.098445364
## factor(Grad)Pale - RS -0.0552075378 0.072187012
## factor(Grad)Prijedor -0.0325032318 0.070359844
## factor(Grad)Prnjavor -0.0026287066 0.081509666
## factor(Grad)Rogatica -0.3475300957 0.109069151
## factor(Grad)Sanski Most -0.0100851070 0.083398565
## factor(Grad)Sarajevo - Centar 0.6302726762 0.065993510
## factor(Grad)Sarajevo - Novi Grad 0.3186763015 0.065786241
## factor(Grad)Sarajevo - Novo Sarajevo 0.5019179732 0.065966740
## factor(Grad)Sarajevo - Stari Grad 0.5579441600 0.067212941
## factor(Grad)Sokolac -0.3443857447 0.117313526
## factor(Grad)Srbac 0.0590790045 0.152647512
## factor(Grad)Srebrenica -0.3201266802 0.095054374
## factor(Grad)Srebrenik 0.0389474289 0.084882145
## factor(Grad)Stolac -0.2822916212 0.205986032
## factor(Grad)Široki Brijeg 0.2554750626 0.205818103
## factor(Grad)Teslić -0.0304758692 0.089642955
## factor(Grad)Tešanj 0.1376101068 0.130320582
## factor(Grad)Tomislavgrad 0.1513583333 0.152656512
## factor(Grad)Travnik 0.0279439939 0.073620642
## factor(Grad)Trebinje 0.2427015529 0.072691895
## factor(Grad)Tuzla 0.1874563102 0.066220620
## factor(Grad)Ugljevik -0.2088883978 0.109155608
## factor(Grad)Vareš -0.4414955011 0.152747169
## factor(Grad)Velika Kladuša 0.3263577248 0.080494199
## factor(Grad)Visoko 0.1051121399 0.086314469
## factor(Grad)Višegrad -0.3549196956 0.117376351
## factor(Grad)Vitez -0.1316447660 0.086141001
## factor(Grad)Vlasenica -0.3944335654 0.130080387
## factor(Grad)Vogošća 0.1624299656 0.067156399
## factor(Grad)Zavidovići -0.0750488280 0.095082061
## factor(Grad)Zenica 0.1022831784 0.066695101
## factor(Grad)Zvornik -0.1807071130 0.074183728
## factor(Grad)Žepče -0.3971418375 0.152701740
## factor(Grad)Živinice 0.0785354641 0.072866035
## t value Pr(>|t|)
## (Intercept) 109.298479823 0.000000e+00
## factor(Y)2017 8.212713855 2.553923e-16
## as.factor(constr.y)1950-1959 -6.371660302 1.988676e-10
## as.factor(constr.y)1960-1969 -12.654508054 2.635943e-36
## as.factor(constr.y)1970-1979 -19.818223191 4.342929e-85
## as.factor(constr.y)1980-1989 -20.321629458 3.011330e-89
## as.factor(constr.y)1990-1999 -13.623748562 9.833920e-42
## as.factor(constr.y)2000-2009 -13.164425859 4.093536e-39
## as.factor(constr.y)2010+ -4.873045438 1.123112e-06
## as.factor(constr.y)Until 1950 -4.171855964 3.057856e-05
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -3.134227743 1.730202e-03
## factor(`Vrsta grijanja`)Centralno (Plin) 0.548408765 5.834288e-01
## factor(`Vrsta grijanja`)Drva -14.054849956 2.869007e-44
## factor(`Vrsta grijanja`)Ostalo -11.379208830 9.684585e-30
## factor(`Vrsta grijanja`)Plin -5.109486423 3.315929e-07
## factor(`Vrsta grijanja`)Struja -7.580782983 3.878504e-14
## factor(Grad)Banovići -1.999732775 4.556775e-02
## factor(Grad)Banja Luka 6.954720711 3.852485e-12
## factor(Grad)Bihać 4.912070325 9.216707e-07
## factor(Grad)Bijeljina 0.983917124 3.251903e-01
## factor(Grad)Bileća -0.229038018 8.188461e-01
## factor(Grad)Bosanski Petrovac -2.897168937 3.777158e-03
## factor(Grad)Bratunac -2.057165153 3.970735e-02
## factor(Grad)Brčko -1.530493499 1.259399e-01
## factor(Grad)Breza -2.018705783 4.355591e-02
## factor(Grad)Brod -4.676203626 2.977198e-06
## factor(Grad)Bugojno -0.405953047 6.847895e-01
## factor(Grad)Busovača 0.359118026 7.195176e-01
## factor(Grad)Cazin 2.535028421 1.126543e-02
## factor(Grad)Čapljina 0.645367253 5.187104e-01
## factor(Grad)Čelić -1.981981739 4.752038e-02
## factor(Grad)Čelinac 0.166204860 8.680005e-01
## factor(Grad)Derventa -2.315032041 2.063999e-02
## factor(Grad)Doboj -0.428270846 6.684671e-01
## factor(Grad)Donji Vakuf -1.616101834 1.061174e-01
## factor(Grad)Foča -0.117683945 9.063215e-01
## factor(Grad)Fojnica 0.435143817 6.634715e-01
## factor(Grad)Gacko 0.216513434 8.285939e-01
## factor(Grad)Goražde 0.360318374 7.186199e-01
## factor(Grad)Gračanica -0.732763689 4.637271e-01
## factor(Grad)Gradačac -0.456779893 6.478434e-01
## factor(Grad)Gradiška -1.601356555 1.093431e-01
## factor(Grad)Hadžići 0.770199084 4.412078e-01
## factor(Grad)Ilidža 4.515817041 6.410775e-06
## factor(Grad)Ilijaš 1.618466388 1.056072e-01
## factor(Grad)Istočna Ilidža -1.473364866 1.406977e-01
## factor(Grad)Istočno Novo Sarajevo -0.277705763 7.812464e-01
## factor(Grad)Jablanica 0.953485635 3.403770e-01
## factor(Grad)Jajce -0.035331926 9.718161e-01
## factor(Grad)Kakanj 1.326041768 1.848692e-01
## factor(Grad)Kalesija 0.247729455 8.043510e-01
## factor(Grad)Kiseljak 1.191073486 2.336652e-01
## factor(Grad)Kladanj -1.127573580 2.595387e-01
## factor(Grad)Ključ 0.410093383 6.817500e-01
## factor(Grad)Kneževo -2.120014168 3.403998e-02
## factor(Grad)Konjic 0.407084100 6.839587e-01
## factor(Grad)Kotor Varoš -0.017549049 9.859991e-01
## factor(Grad)Kozarska Dubica -1.973615023 4.846484e-02
## factor(Grad)Laktaši 2.458975680 1.395746e-02
## factor(Grad)Livno -0.004776523 9.961890e-01
## factor(Grad)Lopare -2.417444281 1.565534e-02
## factor(Grad)Lukavac -0.433326521 6.647909e-01
## factor(Grad)Ljubuški -2.268263543 2.334351e-02
## factor(Grad)Maglaj -0.167293996 8.671436e-01
## factor(Grad)Modriča -2.048796832 4.051924e-02
## factor(Grad)Mostar 6.477378087 9.966304e-11
## factor(Grad)Neum 6.082428719 1.246160e-09
## factor(Grad)Nevesinje -0.165172573 8.688129e-01
## factor(Grad)Novi Grad -0.388223727 6.978623e-01
## factor(Grad)Novi Travnik -1.955631617 5.054834e-02
## factor(Grad)Orašje -0.469879277 6.384559e-01
## factor(Grad)Pale - RS -0.764784919 4.444254e-01
## factor(Grad)Prijedor -0.461957133 6.441264e-01
## factor(Grad)Prnjavor -0.032250244 9.742734e-01
## factor(Grad)Rogatica -3.186328041 1.447248e-03
## factor(Grad)Sanski Most -0.120926625 9.037526e-01
## factor(Grad)Sarajevo - Centar 9.550525092 1.744234e-21
## factor(Grad)Sarajevo - Novi Grad 4.844117831 1.299124e-06
## factor(Grad)Sarajevo - Novo Sarajevo 7.608652119 3.133031e-14
## factor(Grad)Sarajevo - Stari Grad 8.301141932 1.226573e-16
## factor(Grad)Sokolac -2.935601334 3.339860e-03
## factor(Grad)Srbac 0.387028939 6.987466e-01
## factor(Grad)Srebrenica -3.367826936 7.617221e-04
## factor(Grad)Srebrenik 0.458841244 6.463624e-01
## factor(Grad)Stolac -1.370440598 1.705934e-01
## factor(Grad)Široki Brijeg 1.241266237 2.145491e-01
## factor(Grad)Teslić -0.339969484 7.338897e-01
## factor(Grad)Tešanj 1.055935334 2.910342e-01
## factor(Grad)Tomislavgrad 0.991496078 3.214777e-01
## factor(Grad)Travnik 0.379567379 7.042781e-01
## factor(Grad)Trebinje 3.338770492 8.459002e-04
## factor(Grad)Tuzla 2.830784581 4.656651e-03
## factor(Grad)Ugljevik -1.913675354 5.570244e-02
## factor(Grad)Vareš -2.890367814 3.859757e-03
## factor(Grad)Velika Kladuša 4.054425394 5.080326e-05
## factor(Grad)Visoko 1.217781219 2.233482e-01
## factor(Grad)Višegrad -3.023775176 2.505476e-03
## factor(Grad)Vitez -1.528247464 1.264963e-01
## factor(Grad)Vlasenica -3.032229341 2.436445e-03
## factor(Grad)Vogošća 2.418681875 1.560223e-02
## factor(Grad)Zavidovići -0.789305860 4.299600e-01
## factor(Grad)Zenica 1.533593584 1.251749e-01
## factor(Grad)Zvornik -2.435940052 1.487795e-02
## factor(Grad)Žepče -2.600768250 9.321109e-03
## factor(Grad)Živinice 1.077806205 2.811574e-01
model.coef1a[,c("Estimate","Pr(>|t|)")]
## Estimate Pr(>|t|)
## (Intercept) 7.2060359341 0.000000e+00
## factor(Y)2017 0.0387259753 2.553923e-16
## as.factor(constr.y)1950-1959 -0.1160020113 1.988676e-10
## as.factor(constr.y)1960-1969 -0.1293911274 2.635943e-36
## as.factor(constr.y)1970-1979 -0.1751977725 4.342929e-85
## as.factor(constr.y)1980-1989 -0.1582435142 3.011330e-89
## as.factor(constr.y)1990-1999 -0.1718097199 9.833920e-42
## as.factor(constr.y)2000-2009 -0.1188278675 4.093536e-39
## as.factor(constr.y)2010+ -0.0435002002 1.123112e-06
## as.factor(constr.y)Until 1950 -0.0603167263 3.057856e-05
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.0300345690 1.730202e-03
## factor(`Vrsta grijanja`)Centralno (Plin) 0.0050348408 5.834288e-01
## factor(`Vrsta grijanja`)Drva -0.2169975562 2.869007e-44
## factor(`Vrsta grijanja`)Ostalo -0.1433315730 9.684585e-30
## factor(`Vrsta grijanja`)Plin -0.0492288248 3.315929e-07
## factor(`Vrsta grijanja`)Struja -0.0638718468 3.878504e-14
## factor(Grad)Banovići -0.1541401379 4.556775e-02
## factor(Grad)Banja Luka 0.4575410514 3.852485e-12
## factor(Grad)Bihać 0.3374652637 9.216707e-07
## factor(Grad)Bijeljina 0.0660492728 3.251903e-01
## factor(Grad)Bileća -0.0235962604 8.188461e-01
## factor(Grad)Bosanski Petrovac -0.5960722977 3.777158e-03
## factor(Grad)Bratunac -0.2245573582 3.970735e-02
## factor(Grad)Brčko -0.1022770522 1.259399e-01
## factor(Grad)Breza -0.1919830088 4.355591e-02
## factor(Grad)Brod -0.3679017079 2.977198e-06
## factor(Grad)Bugojno -0.0318434580 6.847895e-01
## factor(Grad)Busovača 0.0467363448 7.195176e-01
## factor(Grad)Cazin 0.2092211373 1.126543e-02
## factor(Grad)Čapljina 0.0579547653 5.187104e-01
## factor(Grad)Čelić -0.2328763561 4.752038e-02
## factor(Grad)Čelinac 0.0180991377 8.680005e-01
## factor(Grad)Derventa -0.2001379647 2.063999e-02
## factor(Grad)Doboj -0.0290847219 6.684671e-01
## factor(Grad)Donji Vakuf -0.2105721778 1.061174e-01
## factor(Grad)Foča -0.0116082937 9.063215e-01
## factor(Grad)Fojnica 0.0474992963 6.634715e-01
## factor(Grad)Gacko 0.0330174406 8.285939e-01
## factor(Grad)Goražde 0.0293289491 7.186199e-01
## factor(Grad)Gračanica -0.0631873072 4.637271e-01
## factor(Grad)Gradačac -0.0393205413 6.478434e-01
## factor(Grad)Gradiška -0.1196921913 1.093431e-01
## factor(Grad)Hadžići 0.0574511566 4.412078e-01
## factor(Grad)Ilidža 0.2980243698 6.410775e-06
## factor(Grad)Ilijaš 0.1302953810 1.056072e-01
## factor(Grad)Istočna Ilidža -0.1045076349 1.406977e-01
## factor(Grad)Istočno Novo Sarajevo -0.0188888226 7.812464e-01
## factor(Grad)Jablanica 0.1241168642 3.403770e-01
## factor(Grad)Jajce -0.0027484917 9.718161e-01
## factor(Grad)Kakanj 0.1082168392 1.848692e-01
## factor(Grad)Kalesija 0.0254937947 8.043510e-01
## factor(Grad)Kiseljak 0.1174945895 2.336652e-01
## factor(Grad)Kladanj -0.1720668044 2.595387e-01
## factor(Grad)Ključ 0.0845066595 6.817500e-01
## factor(Grad)Kneževo -0.4365122080 3.403998e-02
## factor(Grad)Konjic 0.0324874210 6.839587e-01
## factor(Grad)Kotor Varoš -0.0020597362 9.859991e-01
## factor(Grad)Kozarska Dubica -0.1946789787 4.846484e-02
## factor(Grad)Laktaši 0.2082458887 1.395746e-02
## factor(Grad)Livno -0.0006227113 9.961890e-01
## factor(Grad)Lopare -0.3151420072 1.565534e-02
## factor(Grad)Lukavac -0.0300477918 6.647909e-01
## factor(Grad)Ljubuški -0.2339008537 2.334351e-02
## factor(Grad)Maglaj -0.0218034126 8.671436e-01
## factor(Grad)Modriča -0.1585112160 4.051924e-02
## factor(Grad)Mostar 0.4331837905 9.966304e-11
## factor(Grad)Neum 0.7143208297 1.246160e-09
## factor(Grad)Nevesinje -0.0251903267 8.688129e-01
## factor(Grad)Novi Grad -0.0340491206 6.978623e-01
## factor(Grad)Novi Travnik -0.1459389154 5.054834e-02
## factor(Grad)Orašje -0.0462574363 6.384559e-01
## factor(Grad)Pale - RS -0.0552075378 4.444254e-01
## factor(Grad)Prijedor -0.0325032318 6.441264e-01
## factor(Grad)Prnjavor -0.0026287066 9.742734e-01
## factor(Grad)Rogatica -0.3475300957 1.447248e-03
## factor(Grad)Sanski Most -0.0100851070 9.037526e-01
## factor(Grad)Sarajevo - Centar 0.6302726762 1.744234e-21
## factor(Grad)Sarajevo - Novi Grad 0.3186763015 1.299124e-06
## factor(Grad)Sarajevo - Novo Sarajevo 0.5019179732 3.133031e-14
## factor(Grad)Sarajevo - Stari Grad 0.5579441600 1.226573e-16
## factor(Grad)Sokolac -0.3443857447 3.339860e-03
## factor(Grad)Srbac 0.0590790045 6.987466e-01
## factor(Grad)Srebrenica -0.3201266802 7.617221e-04
## factor(Grad)Srebrenik 0.0389474289 6.463624e-01
## factor(Grad)Stolac -0.2822916212 1.705934e-01
## factor(Grad)Široki Brijeg 0.2554750626 2.145491e-01
## factor(Grad)Teslić -0.0304758692 7.338897e-01
## factor(Grad)Tešanj 0.1376101068 2.910342e-01
## factor(Grad)Tomislavgrad 0.1513583333 3.214777e-01
## factor(Grad)Travnik 0.0279439939 7.042781e-01
## factor(Grad)Trebinje 0.2427015529 8.459002e-04
## factor(Grad)Tuzla 0.1874563102 4.656651e-03
## factor(Grad)Ugljevik -0.2088883978 5.570244e-02
## factor(Grad)Vareš -0.4414955011 3.859757e-03
## factor(Grad)Velika Kladuša 0.3263577248 5.080326e-05
## factor(Grad)Visoko 0.1051121399 2.233482e-01
## factor(Grad)Višegrad -0.3549196956 2.505476e-03
## factor(Grad)Vitez -0.1316447660 1.264963e-01
## factor(Grad)Vlasenica -0.3944335654 2.436445e-03
## factor(Grad)Vogošća 0.1624299656 1.560223e-02
## factor(Grad)Zavidovići -0.0750488280 4.299600e-01
## factor(Grad)Zenica 0.1022831784 1.251749e-01
## factor(Grad)Zvornik -0.1807071130 1.487795e-02
## factor(Grad)Žepče -0.3971418375 9.321109e-03
## factor(Grad)Živinice 0.0785354641 2.811574e-01
model.coef1a <- model.coef1a
model.coef1a <- as.data.frame.array(model.coef1a[-c(1:16),]) # da bude dataframe a ne matrix i da mi brise prvih 8 redova koji su drugi koeficijenti
Munic <- rownames(model.coef1a)
model.coef1a$Munic <- Munic
dim(model.coef1a) #da vidim koliko redova
## [1] 89 5
rownames(model.coef1a) <- 1:89 #nedostaje B.Krupa koja je 0
head(model.coef1a)
## Estimate Std. Error t value Pr(>|t|) Munic
## 1 -0.15414014 0.07708037 -1.9997328 4.556775e-02 factor(Grad)Banovići
## 2 0.45754105 0.06578856 6.9547207 3.852485e-12 factor(Grad)Banja Luka
## 3 0.33746526 0.06870123 4.9120703 9.216707e-07 factor(Grad)Bihać
## 4 0.06604927 0.06712890 0.9839171 3.251903e-01 factor(Grad)Bijeljina
## 5 -0.02359626 0.10302333 -0.2290380 8.188461e-01 factor(Grad)Bileća
## 6 -0.59607230 0.20574302 -2.8971689 3.777158e-03 factor(Grad)Bosanski Petrovac
Munic1 <- sapply(strsplit(Munic, split=')', fixed=TRUE), function(x) (x[2]))
Munic1
## [1] "Banovići" "Banja Luka"
## [3] "Bihać" "Bijeljina"
## [5] "Bileća" "Bosanski Petrovac"
## [7] "Bratunac" "Brčko"
## [9] "Breza" "Brod"
## [11] "Bugojno" "Busovača"
## [13] "Cazin" "Čapljina"
## [15] "Čelić" "Čelinac"
## [17] "Derventa" "Doboj"
## [19] "Donji Vakuf" "Foča"
## [21] "Fojnica" "Gacko"
## [23] "Goražde" "Gračanica"
## [25] "Gradačac" "Gradiška"
## [27] "Hadžići" "Ilidža"
## [29] "Ilijaš" "Istočna Ilidža"
## [31] "Istočno Novo Sarajevo" "Jablanica"
## [33] "Jajce" "Kakanj"
## [35] "Kalesija" "Kiseljak"
## [37] "Kladanj" "Ključ"
## [39] "Kneževo" "Konjic"
## [41] "Kotor Varoš" "Kozarska Dubica"
## [43] "Laktaši" "Livno"
## [45] "Lopare" "Lukavac"
## [47] "Ljubuški" "Maglaj"
## [49] "Modriča" "Mostar"
## [51] "Neum" "Nevesinje"
## [53] "Novi Grad" "Novi Travnik"
## [55] "Orašje" "Pale - RS"
## [57] "Prijedor" "Prnjavor"
## [59] "Rogatica" "Sanski Most"
## [61] "Sarajevo - Centar" "Sarajevo - Novi Grad"
## [63] "Sarajevo - Novo Sarajevo" "Sarajevo - Stari Grad"
## [65] "Sokolac" "Srbac"
## [67] "Srebrenica" "Srebrenik"
## [69] "Stolac" "Široki Brijeg"
## [71] "Teslić" "Tešanj"
## [73] "Tomislavgrad" "Travnik"
## [75] "Trebinje" "Tuzla"
## [77] "Ugljevik" "Vareš"
## [79] "Velika Kladuša" "Visoko"
## [81] "Višegrad" "Vitez"
## [83] "Vlasenica" "Vogošća"
## [85] "Zavidovići" "Zenica"
## [87] "Zvornik" "Žepče"
## [89] "Živinice"
model.coef1a$Munic1 <- Munic1
model.coef1a[nrow(model.coef1a) + 1,] = c(0, NA, NA, NA,"Bosanska Krupa","Bosanska Krupa" )
model.coef1a$Estimate <- as.numeric(model.coef1a$Estimate)
top50.m1a <- model.coef1a %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
head(50) #top ten municipalities
top50.m1a
## Estimate Munic1
## 1 0.7143208297 Neum
## 2 0.6302726762 Sarajevo - Centar
## 3 0.5579441600 Sarajevo - Stari Grad
## 4 0.5019179732 Sarajevo - Novo Sarajevo
## 5 0.4575410514 Banja Luka
## 6 0.4331837905 Mostar
## 7 0.3374652637 Bihać
## 8 0.3263577248 Velika Kladuša
## 9 0.3186763015 Sarajevo - Novi Grad
## 10 0.2980243698 Ilidža
## 11 0.2554750626 Široki Brijeg
## 12 0.2427015529 Trebinje
## 13 0.2092211373 Cazin
## 14 0.2082458887 Laktaši
## 15 0.1874563102 Tuzla
## 16 0.1624299656 Vogošća
## 17 0.1513583333 Tomislavgrad
## 18 0.1376101068 Tešanj
## 19 0.1302953810 Ilijaš
## 20 0.1241168642 Jablanica
## 21 0.1174945895 Kiseljak
## 22 0.1082168392 Kakanj
## 23 0.1051121399 Visoko
## 24 0.1022831784 Zenica
## 25 0.0845066595 Ključ
## 26 0.0785354641 Živinice
## 27 0.0660492728 Bijeljina
## 28 0.0590790045 Srbac
## 29 0.0579547653 Čapljina
## 30 0.0574511566 Hadžići
## 31 0.0474992963 Fojnica
## 32 0.0467363448 Busovača
## 33 0.0389474289 Srebrenik
## 34 0.0330174406 Gacko
## 35 0.0324874210 Konjic
## 36 0.0293289491 Goražde
## 37 0.0279439939 Travnik
## 38 0.0254937947 Kalesija
## 39 0.0180991377 Čelinac
## 40 0.0000000000 Bosanska Krupa
## 41 -0.0006227113 Livno
## 42 -0.0020597362 Kotor Varoš
## 43 -0.0026287066 Prnjavor
## 44 -0.0027484917 Jajce
## 45 -0.0100851070 Sanski Most
## 46 -0.0116082937 Foča
## 47 -0.0188888226 Istočno Novo Sarajevo
## 48 -0.0218034126 Maglaj
## 49 -0.0235962604 Bileća
## 50 -0.0251903267 Nevesinje
#Analyse db with only flats
only_flats = dist.s.k1 %>% dplyr::filter(prop_type=="Stanovi")
prop.table(table(only_flats$constr.y))
##
## New 1950-1959 1960-1969 1970-1979 1980-1989 1990-1999 2000-2009
## 0.19289841 0.01902212 0.08355643 0.14527265 0.25052839 0.04494857 0.11596449
## 2010+ Until 1950
## 0.11075102 0.03705791
prop.table(table (only_flats$`Vrsta grijanja`))
##
## Centralno (gradsko) Centralno (Kotlovnica) Centralno (Plin)
## 0.46880063 0.06424222 0.08921289
## Drva Ostalo Plin
## 0.05496900 0.05742696 0.09027429
## Struja
## 0.17507402
########
#MODEL 2 no constr year#
########
#This alternative model gives me less deleted observations
model2 <- lm (log(p.sqm) ~ factor (Y) +factor (prop_type) + factor (`Vrsta grijanja`) + factor (Grad), data = dist.s.k1)
summary(model2) # only 35 deleted due to missing
##
## Call:
## lm(formula = log(p.sqm) ~ factor(Y) + factor(prop_type) + factor(`Vrsta grijanja`) +
## factor(Grad), data = dist.s.k1)
##
## Residuals:
## Min 1Q Median 3Q Max
## -3.2500 -0.1612 0.0135 0.1898 1.7771
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 6.449974 0.040535 159.119
## factor(Y)2017 0.031994 0.005001 6.398
## factor(prop_type)Stanovi 0.521728 0.007147 73.003
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 0.018217 0.010274 1.773
## factor(`Vrsta grijanja`)Centralno (Plin) 0.032052 0.011522 2.782
## factor(`Vrsta grijanja`)Drva -0.237388 0.010513 -22.580
## factor(`Vrsta grijanja`)Ostalo -0.241241 0.010744 -22.454
## factor(`Vrsta grijanja`)Plin -0.038921 0.011196 -3.476
## factor(`Vrsta grijanja`)Struja -0.056741 0.009690 -5.856
## factor(Grad)Banovići -0.075895 0.053725 -1.413
## factor(Grad)Banja Luka 0.564299 0.040302 14.002
## factor(Grad)Bihać 0.403619 0.043670 9.242
## factor(Grad)Bijeljina 0.227187 0.041637 5.456
## factor(Grad)Bileća -0.187479 0.131923 -1.421
## factor(Grad)Bosanski Petrovac -0.277320 0.089432 -3.101
## factor(Grad)Bratunac -0.254004 0.079059 -3.213
## factor(Grad)Brčko 0.175813 0.040973 4.291
## factor(Grad)Breza 0.083801 0.062640 1.338
## factor(Grad)Brod -0.177756 0.061504 -2.890
## factor(Grad)Bugojno 0.020978 0.052154 0.402
## factor(Grad)Busovača 0.056545 0.066287 0.853
## factor(Grad)Cazin 0.212583 0.054991 3.866
## factor(Grad)Čapljina 0.137585 0.062678 2.195
## factor(Grad)Čelić -0.022710 0.092270 -0.246
## factor(Grad)Čelinac 0.242522 0.095526 2.539
## factor(Grad)Derventa 0.025853 0.057140 0.452
## factor(Grad)Doboj 0.025161 0.044016 0.572
## factor(Grad)Doboj Jug -0.054230 0.108692 -0.499
## factor(Grad)Donji Vakuf -0.147944 0.074156 -1.995
## factor(Grad)Foča 0.041105 0.062369 0.659
## factor(Grad)Fojnica 0.184841 0.076075 2.430
## factor(Grad)Gacko 0.195654 0.114761 1.705
## factor(Grad)Goražde 0.126266 0.052232 2.417
## factor(Grad)Gornji Vakuf-Uskoplje 0.055524 0.088151 0.630
## factor(Grad)Gračanica 0.014744 0.052909 0.279
## factor(Grad)Gradačac 0.083752 0.056545 1.481
## factor(Grad)Gradiška 0.060020 0.049688 1.208
## factor(Grad)Hadžići 0.358236 0.044393 8.070
## factor(Grad)Ilidža 0.486482 0.040543 11.999
## factor(Grad)Ilijaš 0.430028 0.048255 8.912
## factor(Grad)Istočna Ilidža 0.088575 0.052690 1.681
## factor(Grad)Istočno Novo Sarajevo 0.137697 0.047932 2.873
## factor(Grad)Jablanica 0.381424 0.080724 4.725
## factor(Grad)Jajce 0.017081 0.060342 0.283
## factor(Grad)Kakanj 0.140737 0.052961 2.657
## factor(Grad)Kalesija -0.051183 0.066681 -0.768
## factor(Grad)Kiseljak 0.125291 0.056805 2.206
## factor(Grad)Kladanj 0.086279 0.114739 0.752
## factor(Grad)Ključ -0.186767 0.081631 -2.288
## factor(Grad)Kneževo -0.467874 0.126798 -3.690
## factor(Grad)Konjic 0.250828 0.054319 4.618
## factor(Grad)Kotor Varoš -0.120176 0.101296 -1.186
## factor(Grad)Kozarska Dubica -0.176476 0.070860 -2.490
## factor(Grad)Kupres -0.021235 0.099210 -0.214
## factor(Grad)Laktaši 0.335005 0.051866 6.459
## factor(Grad)Livno -0.210964 0.099221 -2.126
## factor(Grad)Lopare -0.195918 0.111564 -1.756
## factor(Grad)Lukavac 0.025548 0.044834 0.570
## factor(Grad)Ljubuški 0.105311 0.080044 1.316
## factor(Grad)Maglaj 0.043637 0.065248 0.669
## factor(Grad)Modriča -0.026455 0.059910 -0.442
## factor(Grad)Mostar 0.466726 0.041995 11.114
## factor(Grad)Neum 0.950721 0.080920 11.749
## factor(Grad)Nevesinje -0.113939 0.137770 -0.827
## factor(Grad)Novi Grad 0.049726 0.062898 0.791
## factor(Grad)Novi Travnik -0.080299 0.052209 -1.538
## factor(Grad)Odžak -0.227391 0.084620 -2.687
## factor(Grad)Olovo 0.301763 0.086887 3.473
## factor(Grad)Orašje 0.025865 0.099268 0.261
## factor(Grad)Pale - RS 0.047806 0.052634 0.908
## factor(Grad)Prijedor 0.051669 0.045763 1.129
## factor(Grad)Prnjavor 0.041012 0.060739 0.675
## factor(Grad)Ravno 0.279641 0.114868 2.434
## factor(Grad)Rogatica -0.088186 0.079039 -1.116
## factor(Grad)Sanski Most 0.080302 0.056655 1.417
## factor(Grad)Sarajevo - Centar 0.736783 0.040621 18.138
## factor(Grad)Sarajevo - Novi Grad 0.403095 0.040219 10.023
## factor(Grad)Sarajevo - Novo Sarajevo 0.577959 0.040503 14.270
## factor(Grad)Sarajevo - Stari Grad 0.641205 0.042271 15.169
## factor(Grad)Sokolac -0.185978 0.105959 -1.755
## factor(Grad)Srbac -0.014632 0.067392 -0.217
## factor(Grad)Srebrenica -0.321718 0.085775 -3.751
## factor(Grad)Srebrenik 0.148766 0.053473 2.782
## factor(Grad)Stolac -0.008138 0.089427 -0.091
## factor(Grad)Šamac -0.463776 0.111542 -4.158
## factor(Grad)Široki Brijeg 0.372232 0.131899 2.822
## factor(Grad)Teslić 0.153231 0.064613 2.372
## factor(Grad)Tešanj 0.303095 0.060959 4.972
## factor(Grad)Tomislavgrad 0.318566 0.131895 2.415
## factor(Grad)Travnik 0.155690 0.047578 3.272
## factor(Grad)Trebinje 0.455071 0.055749 8.163
## factor(Grad)Tuzla 0.254319 0.040689 6.250
## factor(Grad)Ugljevik 0.113970 0.085824 1.328
## factor(Grad)Usora -0.141301 0.114964 -1.229
## factor(Grad)Vareš -0.580171 0.101340 -5.725
## factor(Grad)Velika Kladuša 0.188368 0.060315 3.123
## factor(Grad)Visoko 0.279002 0.050975 5.473
## factor(Grad)Višegrad -0.225207 0.071885 -3.133
## factor(Grad)Vitez -0.016293 0.055334 -0.294
## factor(Grad)Vlasenica -0.535964 0.090787 -5.904
## factor(Grad)Vogošća 0.446645 0.042667 10.468
## factor(Grad)Zavidovići -0.013844 0.050753 -0.273
## factor(Grad)Zenica 0.215000 0.040854 5.263
## factor(Grad)Zvornik -0.186914 0.055510 -3.367
## factor(Grad)Žepče -0.087548 0.072435 -1.209
## factor(Grad)Živinice 0.078388 0.046791 1.675
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## factor(Y)2017 1.60e-10 ***
## factor(prop_type)Stanovi < 2e-16 ***
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 0.076226 .
## factor(`Vrsta grijanja`)Centralno (Plin) 0.005412 **
## factor(`Vrsta grijanja`)Drva < 2e-16 ***
## factor(`Vrsta grijanja`)Ostalo < 2e-16 ***
## factor(`Vrsta grijanja`)Plin 0.000509 ***
## factor(`Vrsta grijanja`)Struja 4.80e-09 ***
## factor(Grad)Banovići 0.157767
## factor(Grad)Banja Luka < 2e-16 ***
## factor(Grad)Bihać < 2e-16 ***
## factor(Grad)Bijeljina 4.90e-08 ***
## factor(Grad)Bileća 0.155292
## factor(Grad)Bosanski Petrovac 0.001931 **
## factor(Grad)Bratunac 0.001316 **
## factor(Grad)Brčko 1.78e-05 ***
## factor(Grad)Breza 0.180965
## factor(Grad)Brod 0.003853 **
## factor(Grad)Bugojno 0.687513
## factor(Grad)Busovača 0.393653
## factor(Grad)Cazin 0.000111 ***
## factor(Grad)Čapljina 0.028165 *
## factor(Grad)Čelić 0.805589
## factor(Grad)Čelinac 0.011129 *
## factor(Grad)Derventa 0.650947
## factor(Grad)Doboj 0.567573
## factor(Grad)Doboj Jug 0.617831
## factor(Grad)Donji Vakuf 0.046048 *
## factor(Grad)Foča 0.509863
## factor(Grad)Fojnica 0.015117 *
## factor(Grad)Gacko 0.088228 .
## factor(Grad)Goražde 0.015639 *
## factor(Grad)Gornji Vakuf-Uskoplje 0.528786
## factor(Grad)Gračanica 0.780508
## factor(Grad)Gradačac 0.138577
## factor(Grad)Gradiška 0.227084
## factor(Grad)Hadžići 7.33e-16 ***
## factor(Grad)Ilidža < 2e-16 ***
## factor(Grad)Ilijaš < 2e-16 ***
## factor(Grad)Istočna Ilidža 0.092761 .
## factor(Grad)Istočno Novo Sarajevo 0.004072 **
## factor(Grad)Jablanica 2.31e-06 ***
## factor(Grad)Jajce 0.777126
## factor(Grad)Kakanj 0.007880 **
## factor(Grad)Kalesija 0.442744
## factor(Grad)Kiseljak 0.027416 *
## factor(Grad)Kladanj 0.452081
## factor(Grad)Ključ 0.022147 *
## factor(Grad)Kneževo 0.000225 ***
## factor(Grad)Konjic 3.90e-06 ***
## factor(Grad)Kotor Varoš 0.235478
## factor(Grad)Kozarska Dubica 0.012762 *
## factor(Grad)Kupres 0.830513
## factor(Grad)Laktaši 1.07e-10 ***
## factor(Grad)Livno 0.033495 *
## factor(Grad)Lopare 0.079082 .
## factor(Grad)Lukavac 0.568793
## factor(Grad)Ljubuški 0.188296
## factor(Grad)Maglaj 0.503637
## factor(Grad)Modriča 0.658805
## factor(Grad)Mostar < 2e-16 ***
## factor(Grad)Neum < 2e-16 ***
## factor(Grad)Nevesinje 0.408233
## factor(Grad)Novi Grad 0.429200
## factor(Grad)Novi Travnik 0.124055
## factor(Grad)Odžak 0.007210 **
## factor(Grad)Olovo 0.000515 ***
## factor(Grad)Orašje 0.794437
## factor(Grad)Pale - RS 0.363742
## factor(Grad)Prijedor 0.258879
## factor(Grad)Prnjavor 0.499540
## factor(Grad)Ravno 0.014920 *
## factor(Grad)Rogatica 0.264543
## factor(Grad)Sanski Most 0.156385
## factor(Grad)Sarajevo - Centar < 2e-16 ***
## factor(Grad)Sarajevo - Novi Grad < 2e-16 ***
## factor(Grad)Sarajevo - Novo Sarajevo < 2e-16 ***
## factor(Grad)Sarajevo - Stari Grad < 2e-16 ***
## factor(Grad)Sokolac 0.079238 .
## factor(Grad)Srbac 0.828118
## factor(Grad)Srebrenica 0.000177 ***
## factor(Grad)Srebrenik 0.005405 **
## factor(Grad)Stolac 0.927491
## factor(Grad)Šamac 3.22e-05 ***
## factor(Grad)Široki Brijeg 0.004775 **
## factor(Grad)Teslić 0.017722 *
## factor(Grad)Tešanj 6.66e-07 ***
## factor(Grad)Tomislavgrad 0.015728 *
## factor(Grad)Travnik 0.001068 **
## factor(Grad)Trebinje 3.41e-16 ***
## factor(Grad)Tuzla 4.15e-10 ***
## factor(Grad)Ugljevik 0.184206
## factor(Grad)Usora 0.219048
## factor(Grad)Vareš 1.04e-08 ***
## factor(Grad)Velika Kladuša 0.001792 **
## factor(Grad)Visoko 4.45e-08 ***
## factor(Grad)Višegrad 0.001733 **
## factor(Grad)Vitez 0.768416
## factor(Grad)Vlasenica 3.60e-09 ***
## factor(Grad)Vogošća < 2e-16 ***
## factor(Grad)Zavidovići 0.785032
## factor(Grad)Zenica 1.43e-07 ***
## factor(Grad)Zvornik 0.000760 ***
## factor(Grad)Žepče 0.226812
## factor(Grad)Živinice 0.093893 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4177 on 28114 degrees of freedom
## (35 observations deleted due to missingness)
## Multiple R-squared: 0.535, Adjusted R-squared: 0.5333
## F-statistic: 308 on 105 and 28114 DF, p-value: < 2.2e-16
model.coef2 <- summary(model2)$coefficients #
model.coef2 # vidi koliko redova su koef.za druge varijable (9redova) i Opstina BAnovici = Ref
## Estimate Std. Error
## (Intercept) 6.449974137 0.040535438
## factor(Y)2017 0.031994084 0.005000769
## factor(prop_type)Stanovi 0.521727638 0.007146675
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 0.018217197 0.010274296
## factor(`Vrsta grijanja`)Centralno (Plin) 0.032051568 0.011522422
## factor(`Vrsta grijanja`)Drva -0.237388408 0.010513026
## factor(`Vrsta grijanja`)Ostalo -0.241240891 0.010743928
## factor(`Vrsta grijanja`)Plin -0.038921162 0.011195604
## factor(`Vrsta grijanja`)Struja -0.056741376 0.009689942
## factor(Grad)Banovići -0.075894939 0.053724811
## factor(Grad)Banja Luka 0.564299071 0.040301522
## factor(Grad)Bihać 0.403619096 0.043670457
## factor(Grad)Bijeljina 0.227186534 0.041637036
## factor(Grad)Bileća -0.187478597 0.131922957
## factor(Grad)Bosanski Petrovac -0.277319644 0.089432120
## factor(Grad)Bratunac -0.254003752 0.079059267
## factor(Grad)Brčko 0.175812940 0.040972778
## factor(Grad)Breza 0.083801206 0.062640012
## factor(Grad)Brod -0.177756334 0.061504087
## factor(Grad)Bugojno 0.020978243 0.052154023
## factor(Grad)Busovača 0.056544581 0.066287129
## factor(Grad)Cazin 0.212583096 0.054990923
## factor(Grad)Čapljina 0.137584785 0.062678205
## factor(Grad)Čelić -0.022709722 0.092269786
## factor(Grad)Čelinac 0.242521828 0.095526258
## factor(Grad)Derventa 0.025853341 0.057140464
## factor(Grad)Doboj 0.025161194 0.044016130
## factor(Grad)Doboj Jug -0.054229973 0.108692139
## factor(Grad)Donji Vakuf -0.147944465 0.074156117
## factor(Grad)Foča 0.041104946 0.062369048
## factor(Grad)Fojnica 0.184841044 0.076075261
## factor(Grad)Gacko 0.195653736 0.114761164
## factor(Grad)Goražde 0.126265581 0.052232461
## factor(Grad)Gornji Vakuf-Uskoplje 0.055523719 0.088151289
## factor(Grad)Gračanica 0.014743518 0.052908644
## factor(Grad)Gradačac 0.083752031 0.056545164
## factor(Grad)Gradiška 0.060019893 0.049688165
## factor(Grad)Hadžići 0.358236304 0.044393275
## factor(Grad)Ilidža 0.486482438 0.040542687
## factor(Grad)Ilijaš 0.430027869 0.048255186
## factor(Grad)Istočna Ilidža 0.088575244 0.052689815
## factor(Grad)Istočno Novo Sarajevo 0.137697268 0.047932147
## factor(Grad)Jablanica 0.381424104 0.080724106
## factor(Grad)Jajce 0.017080968 0.060341926
## factor(Grad)Kakanj 0.140736516 0.052961036
## factor(Grad)Kalesija -0.051182953 0.066680982
## factor(Grad)Kiseljak 0.125291122 0.056804505
## factor(Grad)Kladanj 0.086279365 0.114739150
## factor(Grad)Ključ -0.186767047 0.081630506
## factor(Grad)Kneževo -0.467874477 0.126797871
## factor(Grad)Konjic 0.250827735 0.054319136
## factor(Grad)Kotor Varoš -0.120176162 0.101295672
## factor(Grad)Kozarska Dubica -0.176476289 0.070859997
## factor(Grad)Kupres -0.021235408 0.099209572
## factor(Grad)Laktaši 0.335004655 0.051865874
## factor(Grad)Livno -0.210964125 0.099221156
## factor(Grad)Lopare -0.195917897 0.111564270
## factor(Grad)Lukavac 0.025548310 0.044834383
## factor(Grad)Ljubuški 0.105311017 0.080043746
## factor(Grad)Maglaj 0.043637330 0.065248437
## factor(Grad)Modriča -0.026454506 0.059910283
## factor(Grad)Mostar 0.466726406 0.041995370
## factor(Grad)Neum 0.950720616 0.080920496
## factor(Grad)Nevesinje -0.113938635 0.137770225
## factor(Grad)Novi Grad 0.049725826 0.062898486
## factor(Grad)Novi Travnik -0.080298842 0.052209173
## factor(Grad)Odžak -0.227391014 0.084620403
## factor(Grad)Olovo 0.301762800 0.086887452
## factor(Grad)Orašje 0.025864978 0.099268324
## factor(Grad)Pale - RS 0.047805957 0.052633985
## factor(Grad)Prijedor 0.051669054 0.045762614
## factor(Grad)Prnjavor 0.041012325 0.060738960
## factor(Grad)Ravno 0.279641090 0.114867619
## factor(Grad)Rogatica -0.088186462 0.079038573
## factor(Grad)Sanski Most 0.080301724 0.056655327
## factor(Grad)Sarajevo - Centar 0.736783007 0.040621394
## factor(Grad)Sarajevo - Novi Grad 0.403094745 0.040218729
## factor(Grad)Sarajevo - Novo Sarajevo 0.577959492 0.040502624
## factor(Grad)Sarajevo - Stari Grad 0.641204740 0.042270842
## factor(Grad)Sokolac -0.185978193 0.105959067
## factor(Grad)Srbac -0.014631936 0.067391784
## factor(Grad)Srebrenica -0.321718241 0.085774542
## factor(Grad)Srebrenik 0.148765767 0.053472662
## factor(Grad)Stolac -0.008138064 0.089426881
## factor(Grad)Šamac -0.463775632 0.111542256
## factor(Grad)Široki Brijeg 0.372231661 0.131899213
## factor(Grad)Teslić 0.153231353 0.064613212
## factor(Grad)Tešanj 0.303094973 0.060958803
## factor(Grad)Tomislavgrad 0.318566254 0.131894718
## factor(Grad)Travnik 0.155690199 0.047578341
## factor(Grad)Trebinje 0.455070728 0.055748997
## factor(Grad)Tuzla 0.254319110 0.040688719
## factor(Grad)Ugljevik 0.113969888 0.085823984
## factor(Grad)Usora -0.141300834 0.114963618
## factor(Grad)Vareš -0.580171419 0.101339641
## factor(Grad)Velika Kladuša 0.188367800 0.060315441
## factor(Grad)Visoko 0.279001807 0.050974560
## factor(Grad)Višegrad -0.225206519 0.071885322
## factor(Grad)Vitez -0.016293116 0.055334054
## factor(Grad)Vlasenica -0.535963856 0.090786862
## factor(Grad)Vogošća 0.446644739 0.042667431
## factor(Grad)Zavidovići -0.013843894 0.050752964
## factor(Grad)Zenica 0.214999732 0.040854101
## factor(Grad)Zvornik -0.186913531 0.055509664
## factor(Grad)Žepče -0.087547982 0.072435157
## factor(Grad)Živinice 0.078387883 0.046791241
## t value Pr(>|t|)
## (Intercept) 159.11938899 0.000000e+00
## factor(Y)2017 6.39783321 1.600758e-10
## factor(prop_type)Stanovi 73.00284566 0.000000e+00
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 1.77308461 7.622551e-02
## factor(`Vrsta grijanja`)Centralno (Plin) 2.78166945 5.411611e-03
## factor(`Vrsta grijanja`)Drva -22.58040703 6.689138e-112
## factor(`Vrsta grijanja`)Ostalo -22.45369647 1.108963e-110
## factor(`Vrsta grijanja`)Plin -3.47646836 5.088310e-04
## factor(`Vrsta grijanja`)Struja -5.85569836 4.802872e-09
## factor(Grad)Banovići -1.41266088 1.577665e-01
## factor(Grad)Banja Luka 14.00192960 2.139104e-44
## factor(Grad)Bihać 9.24238323 2.575821e-20
## factor(Grad)Bijeljina 5.45635700 4.901012e-08
## factor(Grad)Bileća -1.42112186 1.552924e-01
## factor(Grad)Bosanski Petrovac -3.10089533 1.931272e-03
## factor(Grad)Bratunac -3.21282706 1.315835e-03
## factor(Grad)Brčko 4.29096950 1.784889e-05
## factor(Grad)Breza 1.33782233 1.809652e-01
## factor(Grad)Brod -2.89015482 3.853466e-03
## factor(Grad)Bugojno 0.40223632 6.875132e-01
## factor(Grad)Busovača 0.85302505 3.936527e-01
## factor(Grad)Cazin 3.86578520 1.109813e-04
## factor(Grad)Čapljina 2.19509773 2.816473e-02
## factor(Grad)Čelić -0.24612306 8.055888e-01
## factor(Grad)Čelinac 2.53879752 1.112876e-02
## factor(Grad)Derventa 0.45245242 6.509466e-01
## factor(Grad)Doboj 0.57163577 5.675733e-01
## factor(Grad)Doboj Jug -0.49893188 6.178313e-01
## factor(Grad)Donji Vakuf -1.99504060 4.604809e-02
## factor(Grad)Foča 0.65906001 5.098626e-01
## factor(Grad)Fojnica 2.42971291 1.511701e-02
## factor(Grad)Gacko 1.70487759 8.822831e-02
## factor(Grad)Goražde 2.41737759 1.563911e-02
## factor(Grad)Gornji Vakuf-Uskoplje 0.62986848 5.287858e-01
## factor(Grad)Gračanica 0.27865990 7.805079e-01
## factor(Grad)Gradačac 1.48115286 1.385771e-01
## factor(Grad)Gradiška 1.20793135 2.270838e-01
## factor(Grad)Hadžići 8.06960750 7.331451e-16
## factor(Grad)Ilidža 11.99926477 4.318598e-33
## factor(Grad)Ilijaš 8.91153693 5.330231e-19
## factor(Grad)Istočna Ilidža 1.68106956 9.276051e-02
## factor(Grad)Istočno Novo Sarajevo 2.87275401 4.072153e-03
## factor(Grad)Jablanica 4.72503347 2.311892e-06
## factor(Grad)Jajce 0.28306966 7.771255e-01
## factor(Grad)Kakanj 2.65735956 7.879991e-03
## factor(Grad)Kalesija -0.76757947 4.427435e-01
## factor(Grad)Kiseljak 2.20565466 2.741625e-02
## factor(Grad)Kladanj 0.75196100 4.520808e-01
## factor(Grad)Ključ -2.28795649 2.214744e-02
## factor(Grad)Kneževo -3.68992375 2.247444e-04
## factor(Grad)Konjic 4.61766797 3.897931e-06
## factor(Grad)Kotor Varoš -1.18638989 2.354784e-01
## factor(Grad)Kozarska Dubica -2.49049246 1.276235e-02
## factor(Grad)Kupres -0.21404596 8.305128e-01
## factor(Grad)Laktaši 6.45905736 1.070774e-10
## factor(Grad)Livno -2.12620103 3.349521e-02
## factor(Grad)Lopare -1.75609894 7.908240e-02
## factor(Grad)Lukavac 0.56983743 5.687925e-01
## factor(Grad)Ljubuški 1.31566827 1.882961e-01
## factor(Grad)Maglaj 0.66878737 5.036366e-01
## factor(Grad)Modriča -0.44156871 6.588047e-01
## factor(Grad)Mostar 11.11375867 1.233806e-28
## factor(Grad)Neum 11.74882342 8.499052e-32
## factor(Grad)Nevesinje -0.82701930 4.082331e-01
## factor(Grad)Novi Grad 0.79057270 4.292001e-01
## factor(Grad)Novi Travnik -1.53802173 1.240546e-01
## factor(Grad)Odžak -2.68718898 7.209856e-03
## factor(Grad)Olovo 3.47303084 5.153901e-04
## factor(Grad)Orašje 0.26055621 7.944367e-01
## factor(Grad)Pale - RS 0.90827166 3.637425e-01
## factor(Grad)Prijedor 1.12906693 2.588792e-01
## factor(Grad)Prnjavor 0.67522271 4.995398e-01
## factor(Grad)Ravno 2.43446406 1.492004e-02
## factor(Grad)Rogatica -1.11573954 2.645432e-01
## factor(Grad)Sanski Most 1.41737288 1.563850e-01
## factor(Grad)Sarajevo - Centar 18.13780718 4.189779e-73
## factor(Grad)Sarajevo - Novi Grad 10.02256299 1.329117e-23
## factor(Grad)Sarajevo - Novo Sarajevo 14.26968008 4.897765e-46
## factor(Grad)Sarajevo - Stari Grad 15.16896077 9.104014e-52
## factor(Grad)Sokolac -1.75518905 7.923787e-02
## factor(Grad)Srbac -0.21711751 8.281183e-01
## factor(Grad)Srebrenica -3.75074273 1.766652e-04
## factor(Grad)Srebrenik 2.78209017 5.404602e-03
## factor(Grad)Stolac -0.09100244 9.274913e-01
## factor(Grad)Šamac -4.15784700 3.222127e-05
## factor(Grad)Široki Brijeg 2.82209159 4.774502e-03
## factor(Grad)Teslić 2.37151734 1.772191e-02
## factor(Grad)Tešanj 4.97212804 6.661144e-07
## factor(Grad)Tomislavgrad 2.41530713 1.572828e-02
## factor(Grad)Travnik 3.27229148 1.068081e-03
## factor(Grad)Trebinje 8.16285055 3.407653e-16
## factor(Grad)Tuzla 6.25035926 4.153853e-10
## factor(Grad)Ugljevik 1.32794917 1.842057e-01
## factor(Grad)Usora -1.22909175 2.190477e-01
## factor(Grad)Vareš -5.72501947 1.044739e-08
## factor(Grad)Velika Kladuša 3.12304442 1.791724e-03
## factor(Grad)Visoko 5.47335390 4.453667e-08
## factor(Grad)Višegrad -3.13285817 1.732908e-03
## factor(Grad)Vitez -0.29445007 7.684162e-01
## factor(Grad)Vlasenica -5.90353985 3.598644e-09
## factor(Grad)Vogošća 10.46804858 1.349858e-25
## factor(Grad)Zavidovići -0.27277016 7.850319e-01
## factor(Grad)Zenica 5.26262304 1.430553e-07
## factor(Grad)Zvornik -3.36722505 7.603046e-04
## factor(Grad)Žepče -1.20863936 2.268116e-01
## factor(Grad)Živinice 1.67526830 9.389272e-02
model.coef2[,c("Estimate","Pr(>|t|)")]
## Estimate Pr(>|t|)
## (Intercept) 6.449974137 0.000000e+00
## factor(Y)2017 0.031994084 1.600758e-10
## factor(prop_type)Stanovi 0.521727638 0.000000e+00
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) 0.018217197 7.622551e-02
## factor(`Vrsta grijanja`)Centralno (Plin) 0.032051568 5.411611e-03
## factor(`Vrsta grijanja`)Drva -0.237388408 6.689138e-112
## factor(`Vrsta grijanja`)Ostalo -0.241240891 1.108963e-110
## factor(`Vrsta grijanja`)Plin -0.038921162 5.088310e-04
## factor(`Vrsta grijanja`)Struja -0.056741376 4.802872e-09
## factor(Grad)Banovići -0.075894939 1.577665e-01
## factor(Grad)Banja Luka 0.564299071 2.139104e-44
## factor(Grad)Bihać 0.403619096 2.575821e-20
## factor(Grad)Bijeljina 0.227186534 4.901012e-08
## factor(Grad)Bileća -0.187478597 1.552924e-01
## factor(Grad)Bosanski Petrovac -0.277319644 1.931272e-03
## factor(Grad)Bratunac -0.254003752 1.315835e-03
## factor(Grad)Brčko 0.175812940 1.784889e-05
## factor(Grad)Breza 0.083801206 1.809652e-01
## factor(Grad)Brod -0.177756334 3.853466e-03
## factor(Grad)Bugojno 0.020978243 6.875132e-01
## factor(Grad)Busovača 0.056544581 3.936527e-01
## factor(Grad)Cazin 0.212583096 1.109813e-04
## factor(Grad)Čapljina 0.137584785 2.816473e-02
## factor(Grad)Čelić -0.022709722 8.055888e-01
## factor(Grad)Čelinac 0.242521828 1.112876e-02
## factor(Grad)Derventa 0.025853341 6.509466e-01
## factor(Grad)Doboj 0.025161194 5.675733e-01
## factor(Grad)Doboj Jug -0.054229973 6.178313e-01
## factor(Grad)Donji Vakuf -0.147944465 4.604809e-02
## factor(Grad)Foča 0.041104946 5.098626e-01
## factor(Grad)Fojnica 0.184841044 1.511701e-02
## factor(Grad)Gacko 0.195653736 8.822831e-02
## factor(Grad)Goražde 0.126265581 1.563911e-02
## factor(Grad)Gornji Vakuf-Uskoplje 0.055523719 5.287858e-01
## factor(Grad)Gračanica 0.014743518 7.805079e-01
## factor(Grad)Gradačac 0.083752031 1.385771e-01
## factor(Grad)Gradiška 0.060019893 2.270838e-01
## factor(Grad)Hadžići 0.358236304 7.331451e-16
## factor(Grad)Ilidža 0.486482438 4.318598e-33
## factor(Grad)Ilijaš 0.430027869 5.330231e-19
## factor(Grad)Istočna Ilidža 0.088575244 9.276051e-02
## factor(Grad)Istočno Novo Sarajevo 0.137697268 4.072153e-03
## factor(Grad)Jablanica 0.381424104 2.311892e-06
## factor(Grad)Jajce 0.017080968 7.771255e-01
## factor(Grad)Kakanj 0.140736516 7.879991e-03
## factor(Grad)Kalesija -0.051182953 4.427435e-01
## factor(Grad)Kiseljak 0.125291122 2.741625e-02
## factor(Grad)Kladanj 0.086279365 4.520808e-01
## factor(Grad)Ključ -0.186767047 2.214744e-02
## factor(Grad)Kneževo -0.467874477 2.247444e-04
## factor(Grad)Konjic 0.250827735 3.897931e-06
## factor(Grad)Kotor Varoš -0.120176162 2.354784e-01
## factor(Grad)Kozarska Dubica -0.176476289 1.276235e-02
## factor(Grad)Kupres -0.021235408 8.305128e-01
## factor(Grad)Laktaši 0.335004655 1.070774e-10
## factor(Grad)Livno -0.210964125 3.349521e-02
## factor(Grad)Lopare -0.195917897 7.908240e-02
## factor(Grad)Lukavac 0.025548310 5.687925e-01
## factor(Grad)Ljubuški 0.105311017 1.882961e-01
## factor(Grad)Maglaj 0.043637330 5.036366e-01
## factor(Grad)Modriča -0.026454506 6.588047e-01
## factor(Grad)Mostar 0.466726406 1.233806e-28
## factor(Grad)Neum 0.950720616 8.499052e-32
## factor(Grad)Nevesinje -0.113938635 4.082331e-01
## factor(Grad)Novi Grad 0.049725826 4.292001e-01
## factor(Grad)Novi Travnik -0.080298842 1.240546e-01
## factor(Grad)Odžak -0.227391014 7.209856e-03
## factor(Grad)Olovo 0.301762800 5.153901e-04
## factor(Grad)Orašje 0.025864978 7.944367e-01
## factor(Grad)Pale - RS 0.047805957 3.637425e-01
## factor(Grad)Prijedor 0.051669054 2.588792e-01
## factor(Grad)Prnjavor 0.041012325 4.995398e-01
## factor(Grad)Ravno 0.279641090 1.492004e-02
## factor(Grad)Rogatica -0.088186462 2.645432e-01
## factor(Grad)Sanski Most 0.080301724 1.563850e-01
## factor(Grad)Sarajevo - Centar 0.736783007 4.189779e-73
## factor(Grad)Sarajevo - Novi Grad 0.403094745 1.329117e-23
## factor(Grad)Sarajevo - Novo Sarajevo 0.577959492 4.897765e-46
## factor(Grad)Sarajevo - Stari Grad 0.641204740 9.104014e-52
## factor(Grad)Sokolac -0.185978193 7.923787e-02
## factor(Grad)Srbac -0.014631936 8.281183e-01
## factor(Grad)Srebrenica -0.321718241 1.766652e-04
## factor(Grad)Srebrenik 0.148765767 5.404602e-03
## factor(Grad)Stolac -0.008138064 9.274913e-01
## factor(Grad)Šamac -0.463775632 3.222127e-05
## factor(Grad)Široki Brijeg 0.372231661 4.774502e-03
## factor(Grad)Teslić 0.153231353 1.772191e-02
## factor(Grad)Tešanj 0.303094973 6.661144e-07
## factor(Grad)Tomislavgrad 0.318566254 1.572828e-02
## factor(Grad)Travnik 0.155690199 1.068081e-03
## factor(Grad)Trebinje 0.455070728 3.407653e-16
## factor(Grad)Tuzla 0.254319110 4.153853e-10
## factor(Grad)Ugljevik 0.113969888 1.842057e-01
## factor(Grad)Usora -0.141300834 2.190477e-01
## factor(Grad)Vareš -0.580171419 1.044739e-08
## factor(Grad)Velika Kladuša 0.188367800 1.791724e-03
## factor(Grad)Visoko 0.279001807 4.453667e-08
## factor(Grad)Višegrad -0.225206519 1.732908e-03
## factor(Grad)Vitez -0.016293116 7.684162e-01
## factor(Grad)Vlasenica -0.535963856 3.598644e-09
## factor(Grad)Vogošća 0.446644739 1.349858e-25
## factor(Grad)Zavidovići -0.013843894 7.850319e-01
## factor(Grad)Zenica 0.214999732 1.430553e-07
## factor(Grad)Zvornik -0.186913531 7.603046e-04
## factor(Grad)Žepče -0.087547982 2.268116e-01
## factor(Grad)Živinice 0.078387883 9.389272e-02
model.coef2 <- model.coef2
model.coef2 <- as.data.frame.array(model.coef2[-c(1:9),]) # da bude dataframe a ne matrix i da mi brise prvih 8 redova koji su drugi koeficijenti
Munic <- rownames(model.coef2)
model.coef2$Munic <- Munic
dim(model.coef2) #da vidim koliko redova
## [1] 97 5
rownames(model.coef2) <- 1:97 #nedostaje Banovici koje je 0
head(model.coef2)
## Estimate Std. Error t value Pr(>|t|) Munic
## 1 -0.07589494 0.05372481 -1.412661 1.577665e-01 factor(Grad)Banovići
## 2 0.56429907 0.04030152 14.001930 2.139104e-44 factor(Grad)Banja Luka
## 3 0.40361910 0.04367046 9.242383 2.575821e-20 factor(Grad)Bihać
## 4 0.22718653 0.04163704 5.456357 4.901012e-08 factor(Grad)Bijeljina
## 5 -0.18747860 0.13192296 -1.421122 1.552924e-01 factor(Grad)Bileća
## 6 -0.27731964 0.08943212 -3.100895 1.931272e-03 factor(Grad)Bosanski Petrovac
### Make a new variable with municipal names without `factor(Munic)`
Munic1 <- sapply(strsplit(Munic, split=')', fixed=TRUE), function(x) (x[2]))
Munic1
## [1] "Banovići" "Banja Luka"
## [3] "Bihać" "Bijeljina"
## [5] "Bileća" "Bosanski Petrovac"
## [7] "Bratunac" "Brčko"
## [9] "Breza" "Brod"
## [11] "Bugojno" "Busovača"
## [13] "Cazin" "Čapljina"
## [15] "Čelić" "Čelinac"
## [17] "Derventa" "Doboj"
## [19] "Doboj Jug" "Donji Vakuf"
## [21] "Foča" "Fojnica"
## [23] "Gacko" "Goražde"
## [25] "Gornji Vakuf-Uskoplje" "Gračanica"
## [27] "Gradačac" "Gradiška"
## [29] "Hadžići" "Ilidža"
## [31] "Ilijaš" "Istočna Ilidža"
## [33] "Istočno Novo Sarajevo" "Jablanica"
## [35] "Jajce" "Kakanj"
## [37] "Kalesija" "Kiseljak"
## [39] "Kladanj" "Ključ"
## [41] "Kneževo" "Konjic"
## [43] "Kotor Varoš" "Kozarska Dubica"
## [45] "Kupres" "Laktaši"
## [47] "Livno" "Lopare"
## [49] "Lukavac" "Ljubuški"
## [51] "Maglaj" "Modriča"
## [53] "Mostar" "Neum"
## [55] "Nevesinje" "Novi Grad"
## [57] "Novi Travnik" "Odžak"
## [59] "Olovo" "Orašje"
## [61] "Pale - RS" "Prijedor"
## [63] "Prnjavor" "Ravno"
## [65] "Rogatica" "Sanski Most"
## [67] "Sarajevo - Centar" "Sarajevo - Novi Grad"
## [69] "Sarajevo - Novo Sarajevo" "Sarajevo - Stari Grad"
## [71] "Sokolac" "Srbac"
## [73] "Srebrenica" "Srebrenik"
## [75] "Stolac" "Šamac"
## [77] "Široki Brijeg" "Teslić"
## [79] "Tešanj" "Tomislavgrad"
## [81] "Travnik" "Trebinje"
## [83] "Tuzla" "Ugljevik"
## [85] "Usora" "Vareš"
## [87] "Velika Kladuša" "Visoko"
## [89] "Višegrad" "Vitez"
## [91] "Vlasenica" "Vogošća"
## [93] "Zavidovići" "Zenica"
## [95] "Zvornik" "Žepče"
## [97] "Živinice"
model.coef2$Munic1 <- Munic1
#add Bosanska Krupa to dataframe
model.coef2[nrow(model.coef2) + 1,] = c(0, NA, NA, NA,"Bosanska Krupa","Bosanska Krupa" )
model.coef2$Estimate <- as.numeric(model.coef2$Estimate)
top50.m2 <- model.coef2 %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
head(50) #top ten municipalities
top50.m2
## Estimate Munic1
## 1 0.95072062 Neum
## 2 0.73678301 Sarajevo - Centar
## 3 0.64120474 Sarajevo - Stari Grad
## 4 0.57795949 Sarajevo - Novo Sarajevo
## 5 0.56429907 Banja Luka
## 6 0.48648244 Ilidža
## 7 0.46672641 Mostar
## 8 0.45507073 Trebinje
## 9 0.44664474 Vogošća
## 10 0.43002787 Ilijaš
## 11 0.40361910 Bihać
## 12 0.40309474 Sarajevo - Novi Grad
## 13 0.38142410 Jablanica
## 14 0.37223166 Široki Brijeg
## 15 0.35823630 Hadžići
## 16 0.33500465 Laktaši
## 17 0.31856625 Tomislavgrad
## 18 0.30309497 Tešanj
## 19 0.30176280 Olovo
## 20 0.27964109 Ravno
## 21 0.27900181 Visoko
## 22 0.25431911 Tuzla
## 23 0.25082774 Konjic
## 24 0.24252183 Čelinac
## 25 0.22718653 Bijeljina
## 26 0.21499973 Zenica
## 27 0.21258310 Cazin
## 28 0.19565374 Gacko
## 29 0.18836780 Velika Kladuša
## 30 0.18484104 Fojnica
## 31 0.17581294 Brčko
## 32 0.15569020 Travnik
## 33 0.15323135 Teslić
## 34 0.14876577 Srebrenik
## 35 0.14073652 Kakanj
## 36 0.13769727 Istočno Novo Sarajevo
## 37 0.13758479 Čapljina
## 38 0.12626558 Goražde
## 39 0.12529112 Kiseljak
## 40 0.11396989 Ugljevik
## 41 0.10531102 Ljubuški
## 42 0.08857524 Istočna Ilidža
## 43 0.08627937 Kladanj
## 44 0.08380121 Breza
## 45 0.08375203 Gradačac
## 46 0.08030172 Sanski Most
## 47 0.07838788 Živinice
## 48 0.06001989 Gradiška
## 49 0.05654458 Busovača
## 50 0.05552372 Gornji Vakuf-Uskoplje
bottom50.m2 <- model.coef2 %>%
dplyr::select(Estimate, Munic1) %>%
arrange (desc(Estimate))%>%
tail(50) #top ten municipalities
bottom50.m2
## Estimate Munic1
## 49 0.056544581 Busovača
## 50 0.055523719 Gornji Vakuf-Uskoplje
## 51 0.051669054 Prijedor
## 52 0.049725826 Novi Grad
## 53 0.047805957 Pale - RS
## 54 0.043637330 Maglaj
## 55 0.041104946 Foča
## 56 0.041012325 Prnjavor
## 57 0.025864978 Orašje
## 58 0.025853341 Derventa
## 59 0.025548310 Lukavac
## 60 0.025161194 Doboj
## 61 0.020978243 Bugojno
## 62 0.017080968 Jajce
## 63 0.014743518 Gračanica
## 64 0.000000000 Bosanska Krupa
## 65 -0.008138064 Stolac
## 66 -0.013843894 Zavidovići
## 67 -0.014631936 Srbac
## 68 -0.016293116 Vitez
## 69 -0.021235408 Kupres
## 70 -0.022709722 Čelić
## 71 -0.026454506 Modriča
## 72 -0.051182953 Kalesija
## 73 -0.054229973 Doboj Jug
## 74 -0.075894939 Banovići
## 75 -0.080298842 Novi Travnik
## 76 -0.087547982 Žepče
## 77 -0.088186462 Rogatica
## 78 -0.113938635 Nevesinje
## 79 -0.120176162 Kotor Varoš
## 80 -0.141300834 Usora
## 81 -0.147944465 Donji Vakuf
## 82 -0.176476289 Kozarska Dubica
## 83 -0.177756334 Brod
## 84 -0.185978193 Sokolac
## 85 -0.186767047 Ključ
## 86 -0.186913531 Zvornik
## 87 -0.187478597 Bileća
## 88 -0.195917897 Lopare
## 89 -0.210964125 Livno
## 90 -0.225206519 Višegrad
## 91 -0.227391014 Odžak
## 92 -0.254003752 Bratunac
## 93 -0.277319644 Bosanski Petrovac
## 94 -0.321718241 Srebrenica
## 95 -0.463775632 Šamac
## 96 -0.467874477 Kneževo
## 97 -0.535963856 Vlasenica
## 98 -0.580171419 Vareš
###to do next : UPOREDI ESTIMATORS ZA OVE DVIJE OPSTINE
#spasiti a u excel
library(openxlsx)
wb = createWorkbook()
sh1 = addWorksheet(wb, "Estimates model1") #ovo je naziv sheeta
sh2 = addWorksheet(wb, "Estimates model2 no cons.y")#naziv drugog sheeta
sh3 = addWorksheet(wb, "Estimates model 1 no houses")
library(expss)
## Loading required package: maditr
##
## Use magrittr pipe '%>%' to chain several operations:
## mtcars %>%
## let(mpg_hp = mpg/hp) %>%
## take(mean(mpg_hp), by = am)
##
##
##
## Attaching package: 'maditr'
##
## 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
xl_write(model.coef1, wb, sh1) #model.coef1 je naziv objekta koji hocu spasiti
xl_write(model.coef2, wb, sh2)
xl_write(model.coef1a, wb, sh3)
saveWorkbook(wb, "HPI comp Olx 02042024.xlsx", overwrite = T)
## comparison of models
texreg::screenreg(list(model1, model2,model1a))
##
## =======================================================================================
## Model 1 Model 2 Model 3
## ---------------------------------------------------------------------------------------
## (Intercept) 6.57 *** 6.45 *** 7.21 ***
## (0.07) (0.04) (0.07)
## factor(Y)2017 0.04 *** 0.03 *** 0.04 ***
## (0.01) (0.01) (0.00)
## as.factor(constr.y)1950-1959 -0.01 -0.12 ***
## (0.03) (0.02)
## as.factor(constr.y)1960-1969 -0.05 ** -0.13 ***
## (0.02) (0.01)
## as.factor(constr.y)1970-1979 -0.07 *** -0.18 ***
## (0.01) (0.01)
## as.factor(constr.y)1980-1989 -0.12 *** -0.16 ***
## (0.01) (0.01)
## as.factor(constr.y)1990-1999 -0.18 *** -0.17 ***
## (0.02) (0.01)
## as.factor(constr.y)2000-2009 -0.03 * -0.12 ***
## (0.01) (0.01)
## as.factor(constr.y)2010+ -0.01 -0.04 ***
## (0.01) (0.01)
## as.factor(constr.y)Until 1950 0.04 -0.06 ***
## (0.03) (0.01)
## factor(prop_type)Stanovi 0.51 *** 0.52 ***
## (0.01) (0.01)
## factor(`Vrsta grijanja`)Centralno (Kotlovnica) -0.02 0.02 -0.03 **
## (0.02) (0.01) (0.01)
## factor(`Vrsta grijanja`)Centralno (Plin) 0.02 0.03 ** 0.01
## (0.02) (0.01) (0.01)
## factor(`Vrsta grijanja`)Drva -0.27 *** -0.24 *** -0.22 ***
## (0.02) (0.01) (0.02)
## factor(`Vrsta grijanja`)Ostalo -0.28 *** -0.24 *** -0.14 ***
## (0.02) (0.01) (0.01)
## factor(`Vrsta grijanja`)Plin -0.05 ** -0.04 *** -0.05 ***
## (0.02) (0.01) (0.01)
## factor(`Vrsta grijanja`)Struja -0.08 *** -0.06 *** -0.06 ***
## (0.02) (0.01) (0.01)
## factor(Grad)Banovići -0.05 -0.08 -0.15 *
## (0.09) (0.05) (0.08)
## factor(Grad)Banja Luka 0.54 *** 0.56 *** 0.46 ***
## (0.06) (0.04) (0.07)
## factor(Grad)Bihać 0.43 *** 0.40 *** 0.34 ***
## (0.07) (0.04) (0.07)
## factor(Grad)Bijeljina 0.17 * 0.23 *** 0.07
## (0.07) (0.04) (0.07)
## factor(Grad)Bileća 0.28 -0.19 -0.02
## (0.17) (0.13) (0.10)
## factor(Grad)Bosanski Petrovac -0.38 ** -0.28 ** -0.60 **
## (0.14) (0.09) (0.21)
## factor(Grad)Bratunac -0.28 ** -0.25 ** -0.22 *
## (0.10) (0.08) (0.11)
## factor(Grad)Brčko 0.09 0.18 *** -0.10
## (0.07) (0.04) (0.07)
## factor(Grad)Breza 0.13 0.08 -0.19 *
## (0.09) (0.06) (0.10)
## factor(Grad)Brod -0.20 * -0.18 ** -0.37 ***
## (0.10) (0.06) (0.08)
## factor(Grad)Bugojno 0.08 0.02 -0.03
## (0.08) (0.05) (0.08)
## factor(Grad)Busovača 0.06 0.06 0.05
## (0.09) (0.07) (0.13)
## factor(Grad)Cazin 0.24 ** 0.21 *** 0.21 *
## (0.09) (0.05) (0.08)
## factor(Grad)Čapljina 0.25 * 0.14 * 0.06
## (0.10) (0.06) (0.09)
## factor(Grad)Čelić -0.10 -0.02 -0.23 *
## (0.13) (0.09) (0.12)
## factor(Grad)Čelinac 0.16 0.24 * 0.02
## (0.12) (0.10) (0.11)
## factor(Grad)Derventa 0.01 0.03 -0.20 *
## (0.09) (0.06) (0.09)
## factor(Grad)Doboj -0.02 0.03 -0.03
## (0.07) (0.04) (0.07)
## factor(Grad)Doboj Jug 0.06 -0.05
## (0.13) (0.11)
## factor(Grad)Donji Vakuf -0.15 -0.15 * -0.21
## (0.12) (0.07) (0.13)
## factor(Grad)Foča -0.05 0.04 -0.01
## (0.10) (0.06) (0.10)
## factor(Grad)Fojnica 0.29 ** 0.18 * 0.05
## (0.11) (0.08) (0.11)
## factor(Grad)Gacko 0.17 0.20 0.03
## (0.19) (0.11) (0.15)
## factor(Grad)Goražde 0.15 0.13 * 0.03
## (0.08) (0.05) (0.08)
## factor(Grad)Gornji Vakuf-Uskoplje 0.08 0.06
## (0.13) (0.09)
## factor(Grad)Gračanica 0.09 0.01 -0.06
## (0.08) (0.05) (0.09)
## factor(Grad)Gradačac 0.13 0.08 -0.04
## (0.08) (0.06) (0.09)
## factor(Grad)Gradiška 0.01 0.06 -0.12
## (0.08) (0.05) (0.07)
## factor(Grad)Hadžići 0.38 *** 0.36 *** 0.06
## (0.07) (0.04) (0.07)
## factor(Grad)Ilidža 0.45 *** 0.49 *** 0.30 ***
## (0.07) (0.04) (0.07)
## factor(Grad)Ilijaš 0.40 *** 0.43 *** 0.13
## (0.08) (0.05) (0.08)
## factor(Grad)Istočna Ilidža -0.00 0.09 -0.10
## (0.08) (0.05) (0.07)
## factor(Grad)Istočno Novo Sarajevo 0.11 0.14 ** -0.02
## (0.07) (0.05) (0.07)
## factor(Grad)Jablanica 0.36 ** 0.38 *** 0.12
## (0.12) (0.08) (0.13)
## factor(Grad)Jajce 0.13 0.02 -0.00
## (0.09) (0.06) (0.08)
## factor(Grad)Kakanj 0.18 * 0.14 ** 0.11
## (0.09) (0.05) (0.08)
## factor(Grad)Kalesija -0.08 -0.05 0.03
## (0.10) (0.07) (0.10)
## factor(Grad)Kiseljak 0.22 * 0.13 * 0.12
## (0.09) (0.06) (0.10)
## factor(Grad)Kladanj 0.22 0.09 -0.17
## (0.18) (0.11) (0.15)
## factor(Grad)Ključ -0.30 * -0.19 * 0.08
## (0.12) (0.08) (0.21)
## factor(Grad)Kneževo -0.60 ** -0.47 *** -0.44 *
## (0.19) (0.13) (0.21)
## factor(Grad)Konjic 0.44 *** 0.25 *** 0.03
## (0.09) (0.05) (0.08)
## factor(Grad)Kotor Varoš -0.30 * -0.12 -0.00
## (0.15) (0.10) (0.12)
## factor(Grad)Kozarska Dubica -0.17 -0.18 * -0.19 *
## (0.10) (0.07) (0.10)
## factor(Grad)Kupres 0.04 -0.02
## (0.13) (0.10)
## factor(Grad)Laktaši 0.34 *** 0.34 *** 0.21 *
## (0.08) (0.05) (0.08)
## factor(Grad)Livno -0.16 -0.21 * -0.00
## (0.15) (0.10) (0.13)
## factor(Grad)Lopare 0.00 -0.20 -0.32 *
## (0.21) (0.11) (0.13)
## factor(Grad)Lukavac 0.01 0.03 -0.03
## (0.07) (0.04) (0.07)
## factor(Grad)Ljubuški -0.13 0.11 -0.23 *
## (0.14) (0.08) (0.10)
## factor(Grad)Maglaj -0.08 0.04 -0.02
## (0.10) (0.07) (0.13)
## factor(Grad)Modriča -0.01 -0.03 -0.16 *
## (0.09) (0.06) (0.08)
## factor(Grad)Mostar 0.50 *** 0.47 *** 0.43 ***
## (0.07) (0.04) (0.07)
## factor(Grad)Neum 0.80 *** 0.95 *** 0.71 ***
## (0.13) (0.08) (0.12)
## factor(Grad)Nevesinje -0.12 -0.11 -0.03
## (0.21) (0.14) (0.15)
## factor(Grad)Novi Grad -0.01 0.05 -0.03
## (0.09) (0.06) (0.09)
## factor(Grad)Novi Travnik -0.11 -0.08 -0.15
## (0.09) (0.05) (0.07)
## factor(Grad)Odžak 0.21 -0.23 **
## (0.13) (0.08)
## factor(Grad)Olovo 0.29 * 0.30 ***
## (0.14) (0.09)
## factor(Grad)Orašje 0.12 0.03 -0.05
## (0.14) (0.10) (0.10)
## factor(Grad)Pale - RS 0.02 0.05 -0.06
## (0.08) (0.05) (0.07)
## factor(Grad)Prijedor 0.02 0.05 -0.03
## (0.07) (0.05) (0.07)
## factor(Grad)Prnjavor 0.16 0.04 -0.00
## (0.09) (0.06) (0.08)
## factor(Grad)Ravno 0.23 0.28 *
## (0.15) (0.11)
## factor(Grad)Rogatica -0.01 -0.09 -0.35 **
## (0.13) (0.08) (0.11)
## factor(Grad)Sanski Most 0.02 0.08 -0.01
## (0.09) (0.06) (0.08)
## factor(Grad)Sarajevo - Centar 0.70 *** 0.74 *** 0.63 ***
## (0.07) (0.04) (0.07)
## factor(Grad)Sarajevo - Novi Grad 0.40 *** 0.40 *** 0.32 ***
## (0.06) (0.04) (0.07)
## factor(Grad)Sarajevo - Novo Sarajevo 0.55 *** 0.58 *** 0.50 ***
## (0.07) (0.04) (0.07)
## factor(Grad)Sarajevo - Stari Grad 0.59 *** 0.64 *** 0.56 ***
## (0.07) (0.04) (0.07)
## factor(Grad)Sokolac -0.26 -0.19 -0.34 **
## (0.14) (0.11) (0.12)
## factor(Grad)Srbac -0.05 -0.01 0.06
## (0.11) (0.07) (0.15)
## factor(Grad)Srebrenica -0.24 * -0.32 *** -0.32 ***
## (0.12) (0.09) (0.10)
## factor(Grad)Srebrenik 0.18 * 0.15 ** 0.04
## (0.08) (0.05) (0.08)
## factor(Grad)Stolac -0.17 -0.01 -0.28
## (0.14) (0.09) (0.21)
## factor(Grad)Šamac -0.49 ** -0.46 ***
## (0.17) (0.11)
## factor(Grad)Široki Brijeg 0.39 * 0.37 ** 0.26
## (0.16) (0.13) (0.21)
## factor(Grad)Teslić 0.22 * 0.15 * -0.03
## (0.10) (0.06) (0.09)
## factor(Grad)Tešanj 0.25 ** 0.30 *** 0.14
## (0.10) (0.06) (0.13)
## factor(Grad)Tomislavgrad 0.35 0.32 * 0.15
## (0.19) (0.13) (0.15)
## factor(Grad)Travnik 0.21 ** 0.16 ** 0.03
## (0.08) (0.05) (0.07)
## factor(Grad)Trebinje 0.44 *** 0.46 *** 0.24 ***
## (0.08) (0.06) (0.07)
## factor(Grad)Tuzla 0.23 *** 0.25 *** 0.19 **
## (0.07) (0.04) (0.07)
## factor(Grad)Ugljevik 0.20 0.11 -0.21
## (0.13) (0.09) (0.11)
## factor(Grad)Usora -0.24 -0.14
## (0.14) (0.11)
## factor(Grad)Vareš -0.71 *** -0.58 *** -0.44 **
## (0.19) (0.10) (0.15)
## factor(Grad)Velika Kladuša 0.18 0.19 ** 0.33 ***
## (0.09) (0.06) (0.08)
## factor(Grad)Visoko 0.29 *** 0.28 *** 0.11
## (0.08) (0.05) (0.09)
## factor(Grad)Višegrad -0.25 * -0.23 ** -0.35 **
## (0.11) (0.07) (0.12)
## factor(Grad)Vitez -0.07 -0.02 -0.13
## (0.09) (0.06) (0.09)
## factor(Grad)Vlasenica -0.30 * -0.54 *** -0.39 **
## (0.13) (0.09) (0.13)
## factor(Grad)Vogošća 0.36 *** 0.45 *** 0.16 *
## (0.07) (0.04) (0.07)
## factor(Grad)Zavidovići -0.04 -0.01 -0.08
## (0.08) (0.05) (0.10)
## factor(Grad)Zenica 0.21 ** 0.21 *** 0.10
## (0.07) (0.04) (0.07)
## factor(Grad)Zvornik -0.21 ** -0.19 *** -0.18 *
## (0.08) (0.06) (0.07)
## factor(Grad)Žepče -0.13 -0.09 -0.40 **
## (0.10) (0.07) (0.15)
## factor(Grad)Živinice 0.10 0.08 0.08
## (0.07) (0.05) (0.07)
## ---------------------------------------------------------------------------------------
## R^2 0.58 0.53 0.63
## Adj. R^2 0.57 0.53 0.63
## Num. obs. 11710 28220 7097
## =======================================================================================
## *** p < 0.001; ** p < 0.01; * p < 0.05
############################
#FOR DESCRIPTIVE STATISTICS#
############################
dim(dist.s.k1)
## [1] 28255 23
summary(dist.s.k1$Cijena)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 20400 65000 92000 114979 137000 590000
table(is.na(dist.s.k1$Cijena))
##
## FALSE
## 28255
sd (dist.s.k1$Cijena)
## [1] 78830.01
summary(dist.s.k1$sqm)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 25.24 53.00 70.00 105.62 120.00 980.00
table(is.na(dist.s.k1$sqm))
##
## FALSE
## 28255
sd (dist.s.k1$sqm)
## [1] 100.9103
summary(dist.s.k1$p.sqm)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 26.32 900.00 1340.43 1366.15 1791.05 3470.59
table(is.na(dist.s.k1$p.sqm))
##
## FALSE
## 28255
sd (dist.s.k1$p.sqm)
## [1] 643.5168
table(dist.s.k1$prop_type)
##
## Kuće Stanovi
## 10339 17916
table (is.na(dist.s.k1$prop_type))
##
## FALSE
## 28255
prop.table(table(dist.s.k1$prop_type))
##
## Kuće Stanovi
## 0.3659175 0.6340825
table(dist.s.k1$Y)
##
## 2016 2017
## 13089 15166
table(is.na(dist.s.k1$Y)) #FALSE
##
## FALSE
## 28255
prop.table(table(dist.s.k1$Y))
##
## 2016 2017
## 0.4632454 0.5367546
table (dist.s.k1$prop_type)
##
## Kuće Stanovi
## 10339 17916
table (is.na (dist.s.k1$prop_type)) #FALSE
##
## FALSE
## 28255
prop.table(table(dist.s.k1$prop_type))
##
## Kuće Stanovi
## 0.3659175 0.6340825
table (dist.s.k1$constr.y)
##
## New 1950-1959 1960-1969 1970-1979 1980-1989 1990-1999 2000-2009
## 1929 169 860 1549 2731 904 1918
## 2010+ Until 1950
## 1325 325
table(is.na(dist.s.k1$constr.y)) #T: 16545 F:11710
##
## FALSE TRUE
## 11710 16545
prop.table(table(dist.s.k1$constr.y))
##
## New 1950-1959 1960-1969 1970-1979 1980-1989 1990-1999 2000-2009
## 0.16473100 0.01443211 0.07344150 0.13228010 0.23321947 0.07719898 0.16379163
## 2010+ Until 1950
## 0.11315115 0.02775406
table (dist.s.k1$`Vrsta grijanja`)
##
## Centralno (gradsko) Centralno (Kotlovnica) Centralno (Plin)
## 8392 3324 2024
## Drva Ostalo Plin
## 5228 3043 2250
## Struja
## 3959
table(is.na(dist.s.k1$`Vrsta grijanja`)) #F: 28220 T: 35
##
## FALSE TRUE
## 28220 35
prop.table(table(dist.s.k1$`Vrsta grijanja`))
##
## Centralno (gradsko) Centralno (Kotlovnica) Centralno (Plin)
## 0.29737775 0.11778880 0.07172218
## Drva Ostalo Plin
## 0.18525868 0.10783133 0.07973069
## Struja
## 0.14029057
#########
#MAPPING#
#########
### open libraries
library(ggplot2)
library(ggmap)
## ℹ Google's Terms of Service: <https://mapsplatform.google.com>
## ℹ Please cite ggmap if you use it! Use `citation("ggmap")` for details.
# maping with sf and ggplot
library(sf)
## Linking to GEOS 3.11.2, GDAL 3.7.2, PROJ 9.3.0; sf_use_s2() is TRUE
library(tmap)
## Breaking News: tmap 3.x is retiring. Please test v4, e.g. with
## remotes::install_github('r-tmap/tmap')
library(maps)
library(raster)
## Loading required package: sp
library(dplyr)
library(spData)
#library(spDataLarge) samo mi ova nije instalirana, probacu bez nje
## Import map
library(readr)
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
amraGIS.nova$shapeName
## [1] "Brčko distrikt Bosne i Hercegovine" "Vukosavlje"
## [3] "Banja Luka" "Bihac"
## [5] "Bijeljina" "Doboj"
## [7] "Jajce" "Mostar"
## [9] "Prijedor" "Trebinje"
## [11] "Zenica" "Banovići"
## [13] "Berkovići" "Bileća"
## [15] "Bosanska Gradiška/Gradiška" "Bosanska Kostajnica/Kostajnica"
## [17] "Bosanska Krupa" "Bosanska/Kozarska Dubica"
## [19] "Bosanski Brod" "Bosanski Novi/Novi"
## [21] "Bosanski Petrovac" "Bosanski Šamac/Šamac"
## [23] "Bosansko Grahovo" "Breza"
## [25] "Bugojno" "Busovača"
## [27] "Bužim" "Čajniče"
## [29] "Čapljina" "Cazin"
## [31] "Čelić" "Čelinac"
## [33] "Centar" "Čitluk"
## [35] "Derventa" "Doboj Istok"
## [37] "Doboj Jug" "Dobretići"
## [39] "Domaljevac-Šamac" "Donji Vakuf"
## [41] "Donji Žabar" "Drvar"
## [43] "Istočni Drvar" "Istočna Ilidža"
## [45] "Istočni Mostar" "Istočno Novo Sarajevo"
## [47] "Istočni Stari Grad" "Foča"
## [49] "Fojnica" "Gacko"
## [51] "Glamoč" "Goražde"
## [53] "Gornji Vakuf/Uskoplje" "Gračanica"
## [55] "Gradačac" "Grude"
## [57] "Hadžići" "Han Pijesak"
## [59] "Ilidža" "Ilijaš"
## [61] "Jablanica" "Jezero"
## [63] "Kakanj" "Kalesija"
## [65] "Kalinovik" "Kiseljak"
## [67] "Kladanj" "Ključ"
## [69] "Konjic" "Kotor Varoš"
## [71] "Kreševo" "Krupa na Uni"
## [73] "Kupres (RS)" "Kupres"
## [75] "Laktaši" "Livno"
## [77] "Ljubinje" "Ljubuški"
## [79] "Lopare" "Lukavac"
## [81] "Maglaj" "Milići"
## [83] "Modriča" "Mrkonjic Grad"
## [85] "Neum" "Nevesinje"
## [87] "Novo Goražde" "Novo Sarajevo"
## [89] "Novi Grad" "Novi Travnik"
## [91] "Odžak" "Stari Grad"
## [93] "Olovo" "Orašje"
## [95] "Osmaci" "Oštra Luka"
## [97] "Pale" "Pale-Prača"
## [99] "Pelagićevo" "Petrovac"
## [101] "Petrovo" "Posušje"
## [103] "Prnjavor" "Prozor"
## [105] "Ravno" "Ribnik"
## [107] "Rogatica" "Rudo"
## [109] "Sanski Most" "Sapna"
## [111] "Šekovići" "Šipovo"
## [113] "Široki Brijeg" "Skender Vakuf/Kneževo"
## [115] "Sokolac" "Srbac"
## [117] "Srebrenica" "Srebrenik"
## [119] "Stolac" "Teočak"
## [121] "Tešanj" "Teslić"
## [123] "Tomislavgrad" "Travnik"
## [125] "Trnovo" "Tuzla"
## [127] "Ugljevik" "Usora"
## [129] "Ustikolina" "Vareš"
## [131] "Velika Kladuša" "Višegrad"
## [133] "Visoko" "Vitez"
## [135] "Vlasenica" "Vogošća"
## [137] "Zavidovići" "Žepče"
## [139] "Živinice" "Zvornik"
## [141] "Bratunac" "Trnovo (RS)"
amraGIS.nova$shapeName1 <- c("Brcko", "Vukosavlje","Banja Luka", "Bihac", "Bijeljina", "Doboj", "Jajce", "Mostar", "Prijedor", "Trebinje",
"Zenica", "Banovici", "Berkovici" , "Bileca","Gradiska", "Kostajnica",
"Bosanska Krupa" , "Kozarska Dubica", "Brod", "Novi Grad", "Bosanski Petrovac",
"Samac", "Bosansko Grahovo", "Breza" , "Bugojno" , "Busovaca", "Buzim" , "Cajnice",
"Capljina", "Cazin", "Celic", "Celinac", "Sarajevo - Centar", "Citluk",
"Derventa", "Doboj Istok", "Doboj Jug", "Dobretici", "Domaljevac-Samac",
"Donji Vakuf", "Donji Zabar", "Drvar", "Istocni Drvar", "Istocna Ilidza",
"Istocni Mostar", "Istocno Novo Sarajevo", "Istocni Stari Grad", "Foca (RS)",
"Fojnica", "Gacko", "Glamoc","Gorazde", "Gornji Vakuf-Uskoplje", "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 (FBiH)", "Laktasi", "Livno", "Ljubinje", "Ljubuski", "Lopare", "Lukavac",
"Maglaj", "Milici", "Modrica", "Mrkonjic Grad" ,"Neum", "Nevesinje", "Novo Gorazde",
"Sarajevo - Novo Sarajevo", "Sarajevo - Novi Grad", "Novi Travnik", "Odzak",
"Sarajevo - Stari Grad", "Olovo", "Orasje", "Osmaci", "Ostra Luka", "Pale (RS)",
"Pale (FBiH)", "Pelagicevo", "Petrovac", "Petrovo", "Posusje", "Prnjavor",
"Prozor-Rama", "Ravno", "Ribnik", "Rogatica", "Rudo", "Sanski Most", "Sapna",
"Sekovici", "Sipovo", "Siroki Brijeg", "Knezevo", "Sokolac", "Srbac", "Srebrenica",
"Srebrenik", "Stolac", "Teocak", "Tesanj", "Teslic", "Tomislavgrad", "Travnik",
"Trnovo (FBiH)", "Tuzla", "Ugljevik", "Usora", "Foca (FBiH)", "Vares",
"Velika Kladusa", "Visegrad", "Visoko", "Vitez", "Vlasenica", "Vogosca",
"Zavidovici", "Zepce", "Zivinice", "Zvornik", "Bratunac", "Trnovo (RS)")
## import datasets
library(readxl)
HPI_monetar_all <- read_excel("C:/Users/Amra/OneDrive - Direkcija za ekonomsko planiranje/Prijava teme doktorske disertacije/BHAAS Conference/HPI monetar all.xlsx",
col_types = c("numeric", "text", "text",
"text", "text", "text", "numeric",
"text", "numeric", "numeric", "numeric",
"numeric", "numeric", "numeric"))
## Warning: Expecting numeric in N7 / R7C14: got 'Top five cites: Neum,
## Sarajevo-Center, Sarajevo - Stari Grad, Sarajevo - Novo Sarajevo, Banja Luka,
## Mostar, Ilidza'
## Warning: Expecting numeric in N8 / R8C14: got 'Bottom five municipalities:
## Vares, Knezevo, Samac, Donji Zabar, Bosanski Petrovac, Vlasenica'
## New names:
## • `Estimate` -> `Estimate...1`
## • `Estimate` -> `Estimate...9`
## • `` -> `...13`
library(tidyverse)
library (conflicted)
HPI_monetar_all$Entity<- fct_relevel(HPI_monetar_all$Entity,"FBiH","RS")
table (HPI_monetar_all$Entity)
##
## FBiH RS BDBiH
## 79 62 1
HPI.map <- merge ( HPI_monetar_all, amraGIS.nova, by.y='shapeName1', by.x = "Munic1", all=T)
dim(HPI_monetar_all)
## [1] 142 14
dim(amraGIS.nova)
## [1] 142 7
dim (HPI.map)
## [1] 142 20
ggplot() + geom_sf(data = HPI.map, aes(geometry=geometry, fill = `Pricepersqm.no formula`),colour="white") +
scale_fill_viridis_c(name = "Price per sqm") +
scale_x_continuous(breaks = c(18, 48)) +
geom_sf_text(data=HPI.map, aes(geometry=geometry,label=ifelse(as.numeric(`Pricepersqm.no formula`) < 600,as.character(Munic1),"")),hjust=0.5,vjust=0,cex =3, col="white") +
geom_sf_text(data=HPI.map, aes(geometry=geometry,label=ifelse(as.numeric(`Pricepersqm.no formula`) > 2200,as.character(Munic1),"")),hjust=0.5,vjust=0,cex =3, col = "darkred")
## 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
##map entities and BD
ggplot() + geom_sf(data = HPI.map, aes(geometry=geometry, fill = Entity)) +
scale_fill_viridis_d(name = "ID_2")
#******************************************************
#ISMIROV INPUT ZA KORELATION BETWEEN DIVERSITY AND HPI*
#******************************************************
Diversity_index_correlation_BHAS <- read_excel("C:/Users/Amra/OneDrive - Direkcija za ekonomsko planiranje/Prijava teme doktorske disertacije/BHAAS Conference/Diversity index correlation.xlsx",
sheet = "BHAS")
#Omit BD because it is not relevant - only one obs for entity impact
diversity.Hpi.noBD <- Diversity_index_correlation_BHAS %>%
dplyr::filter(Entity != "BDBiH")
table(diversity.Hpi.noBD$Entity)
##
## FBiH RS
## 79 62
model.d.1 <- lm (log(Price.sqm) ~
`Diversity Index-ethnicity`,
data = diversity.Hpi.noBD)
summary(model.d.1)
##
## Call:
## lm(formula = log(Price.sqm) ~ `Diversity Index-ethnicity`, data = diversity.Hpi.noBD)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.71280 -0.15173 -0.01177 0.14066 0.68132
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.82991 0.08951 76.300 < 2e-16 ***
## `Diversity Index-ethnicity` 0.31829 0.11486 2.771 0.00635 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2411 on 139 degrees of freedom
## Multiple R-squared: 0.05235, Adjusted R-squared: 0.04554
## F-statistic: 7.679 on 1 and 139 DF, p-value: 0.006352
model.d.2 <- lm (log(Price.sqm) ~
`Diversity Index-mother language`,
data = diversity.Hpi.noBD)
summary(model.d.2)
##
## Call:
## lm(formula = log(Price.sqm) ~ `Diversity Index-mother language`,
## data = diversity.Hpi.noBD)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.71734 -0.15070 -0.01216 0.14454 0.65500
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.74389 0.08877 75.974 < 2e-16 ***
## `Diversity Index-mother language` 0.41700 0.11012 3.787 0.000226 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2358 on 139 degrees of freedom
## Multiple R-squared: 0.09351, Adjusted R-squared: 0.08699
## F-statistic: 14.34 on 1 and 139 DF, p-value: 0.0002264
model.d.3 <- lm (log(Price.sqm) ~
`Diversity Index-religion`,
data = diversity.Hpi.noBD)
summary(model.d.3)
##
## Call:
## lm(formula = log(Price.sqm) ~ `Diversity Index-religion`, data = diversity.Hpi.noBD)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.71087 -0.14750 -0.00561 0.13671 0.67682
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.82458 0.08826 77.325 < 2e-16 ***
## `Diversity Index-religion` 0.32361 0.11258 2.874 0.00468 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2406 on 139 degrees of freedom
## Multiple R-squared: 0.05611, Adjusted R-squared: 0.04932
## F-statistic: 8.262 on 1 and 139 DF, p-value: 0.004684
model.all3 <- lm(log(Price.sqm)~
`Diversity Index-ethnicity` + `Diversity Index-mother language`+ `Diversity Index-religion`,
data = diversity.Hpi.noBD)
summary(model.all3)
##
## Call:
## lm(formula = log(Price.sqm) ~ `Diversity Index-ethnicity` + `Diversity Index-mother language` +
## `Diversity Index-religion`, data = diversity.Hpi.noBD)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.73246 -0.13228 0.00112 0.16187 0.71511
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.70503 0.08514 78.753 < 2e-16 ***
## `Diversity Index-ethnicity` 0.73170 0.73449 0.996 0.32091
## `Diversity Index-mother language` 3.11198 0.58295 5.338 3.79e-07 ***
## `Diversity Index-religion` -3.45191 1.03606 -3.332 0.00111 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2203 on 137 degrees of freedom
## Multiple R-squared: 0.2197, Adjusted R-squared: 0.2027
## F-statistic: 12.86 on 3 and 137 DF, p-value: 1.873e-07
model.5.inter <- lm(log(Price.sqm)~
`Diversity Index-mother language`+ `Diversity Index-mother language`*factor(Entity),
data = diversity.Hpi.noBD)
summary(model.5.inter)
##
## Call:
## lm(formula = log(Price.sqm) ~ `Diversity Index-mother language` +
## `Diversity Index-mother language` * factor(Entity), data = diversity.Hpi.noBD)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.78760 -0.13176 -0.00702 0.12205 0.60582
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 6.9419 0.1148 60.476
## `Diversity Index-mother language` 0.2612 0.1392 1.877
## factor(Entity)RS -0.3446 0.1680 -2.052
## `Diversity Index-mother language`:factor(Entity)RS 0.2267 0.2097 1.081
## Pr(>|t|)
## (Intercept) <2e-16 ***
## `Diversity Index-mother language` 0.0627 .
## factor(Entity)RS 0.0421 *
## `Diversity Index-mother language`:factor(Entity)RS 0.2815
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2213 on 137 degrees of freedom
## Multiple R-squared: 0.2131, Adjusted R-squared: 0.1959
## F-statistic: 12.37 on 3 and 137 DF, p-value: 3.302e-07
model.6.no.int <- lm(log(Price.sqm) ~
`Diversity Index-mother language` + factor(Entity),
data = diversity.Hpi.noBD)
summary(model.6.no.int)
##
## Call:
## lm(formula = log(Price.sqm) ~ `Diversity Index-mother language` +
## factor(Entity), data = diversity.Hpi.noBD)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.75946 -0.12325 -0.00013 0.14170 0.59077
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 6.86150 0.08748 78.433 < 2e-16 ***
## `Diversity Index-mother language` 0.36113 0.10418 3.467 0.000703 ***
## factor(Entity)RS -0.16765 0.03784 -4.430 1.9e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2214 on 138 degrees of freedom
## Multiple R-squared: 0.2064, Adjusted R-squared: 0.1949
## F-statistic: 17.94 on 2 and 138 DF, p-value: 1.184e-07
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").
screenreg(list(model.d.1,model.d.2, model.d.3, model.all3,model.5.inter, model.6.no.int))
##
## ==========================================================================================================================
## Model 1 Model 2 Model 3 Model 4 Model 5 Model 6
## --------------------------------------------------------------------------------------------------------------------------
## (Intercept) 6.83 *** 6.74 *** 6.82 *** 6.71 *** 6.94 *** 6.86 ***
## (0.09) (0.09) (0.09) (0.09) (0.11) (0.09)
## `Diversity Index-ethnicity` 0.32 ** 0.73
## (0.11) (0.73)
## `Diversity Index-mother language` 0.42 *** 3.11 *** 0.26 0.36 ***
## (0.11) (0.58) (0.14) (0.10)
## `Diversity Index-religion` 0.32 ** -3.45 **
## (0.11) (1.04)
## factor(Entity)RS -0.34 * -0.17 ***
## (0.17) (0.04)
## `Diversity Index-mother language`:factor(Entity)RS 0.23
## (0.21)
## --------------------------------------------------------------------------------------------------------------------------
## R^2 0.05 0.09 0.06 0.22 0.21 0.21
## Adj. R^2 0.05 0.09 0.05 0.20 0.20 0.19
## Num. obs. 141 141 141 141 141 141
## ==========================================================================================================================
## *** p < 0.001; ** p < 0.01; * p < 0.05
##################################################
#correlation matrix - Ismirov prijedlog sastanak 8.3.
###################################################
##BEZ BDBIH
##select only those that are numerical
names(diversity.Hpi.noBD)
## [1] "Munic1" "Entity"
## [3] "Diversity Index-ethnicity" "Diversity Index-religion"
## [5] "Diversity Index-mother language" "Price.sqm"
library (Hmisc)
## Registered S3 methods overwritten by 'Hmisc':
## method from
## [.labelled expss
## print.labelled expss
## as.data.frame.labelled expss
rcorr(as.matrix(diversity.Hpi.noBD[,3:5]))
## Diversity Index-ethnicity
## Diversity Index-ethnicity 1.00
## Diversity Index-religion 0.99
## Diversity Index-mother language 0.96
## Diversity Index-religion
## Diversity Index-ethnicity 0.99
## Diversity Index-religion 1.00
## Diversity Index-mother language 0.98
## Diversity Index-mother language
## Diversity Index-ethnicity 0.96
## Diversity Index-religion 0.98
## Diversity Index-mother language 1.00
##
## n= 141
##
##
## P
## Diversity Index-ethnicity
## Diversity Index-ethnicity
## Diversity Index-religion 0
## Diversity Index-mother language 0
## Diversity Index-religion
## Diversity Index-ethnicity 0
## Diversity Index-religion
## Diversity Index-mother language 0
## Diversity Index-mother language
## Diversity Index-ethnicity 0
## Diversity Index-religion 0
## Diversity Index-mother language
##SA BDBIH
names(Diversity_index_correlation_BHAS)
## [1] "Munic1" "Entity"
## [3] "Diversity Index-ethnicity" "Diversity Index-religion"
## [5] "Diversity Index-mother language" "Price.sqm"
rcorr(as.matrix(Diversity_index_correlation_BHAS[,3:5]))
## Diversity Index-ethnicity
## Diversity Index-ethnicity 1.00
## Diversity Index-religion 0.99
## Diversity Index-mother language 0.97
## Diversity Index-religion
## Diversity Index-ethnicity 0.99
## Diversity Index-religion 1.00
## Diversity Index-mother language 0.98
## Diversity Index-mother language
## Diversity Index-ethnicity 0.97
## Diversity Index-religion 0.98
## Diversity Index-mother language 1.00
##
## n= 142
##
##
## P
## Diversity Index-ethnicity
## Diversity Index-ethnicity
## Diversity Index-religion 0
## Diversity Index-mother language 0
## Diversity Index-religion
## Diversity Index-ethnicity 0
## Diversity Index-religion
## Diversity Index-mother language 0
## Diversity Index-mother language
## Diversity Index-ethnicity 0
## Diversity Index-religion 0
## Diversity Index-mother language
####################
#POPULATION DENSITY#
####################
library(readxl)
povrsina_munic_density <- read_excel("C:/Users/Amra/OneDrive - Direkcija za ekonomsko planiranje/Prijava teme doktorske disertacije/data collection for BiH/povrsina munic-density.xlsx")
View(povrsina_munic_density)
names(povrsina_munic_density)
## [1] "Num." "Munic" "Povrsina_km2" "Pop.density" "Entity"
density.map <- merge ( povrsina_munic_density, amraGIS.nova, by.y='shapeName', by.x = "Munic", all=T)
dim (density.map)
## [1] 212 11
# min/max values
#vmin <- min(density.map$Pop.density, na.rm = T)
#vmax <- max(density.map$Pop.density, na.rm = T)
# bins
library (classInt)
#brk <- round(classIntervals(density.map$Pop.density,
# n = 12, style = "equal" )$brks, 0) |>
#head(-1) |>
#tail(-1) |>
#append(vmax)
# breaks
#breaks <- c(vmin, brk)
# colors
#cols <- rev(c(
# "#140e26","#451e15", "#451a40","#b32927",
# "#7d1d53", "#9d0d23","#b32957","#b32987","#b95999" ,"#ccaead",
#"#ccabcc","#eb804e", "#ffdc58"
#))
ggplot() + geom_sf(data = density.map, aes(geometry=geometry, fill = Pop.density),colour = "white") +
scale_fill_viridis_b(breaks = c(2,4,10,50,100,200,1000,10000),trans="log10" )
#####################
#MAP DIVERSITY INDEX#
#####################
diversity.map <- merge (Diversity_index_correlation_BHAS, amraGIS.nova, by.y='shapeName', by.x = "Munic1", all=T)
dim (diversity.map)
## [1] 212 12
ggplot() + geom_sf(data = diversity.map, aes(geometry=geometry, fill = `Diversity Index-mother language`),colour = "darkgray") +
scale_fill_viridis_b(limits = c(0,1), breaks = c(0,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1))
#***************
#LINIJA SA MAPE jos nisam uspjela
#***************
## Prebaci MULTIPOLYGON u long i lat data
#amraGIS.nova1 <- cbind(amraGIS.nova,st_coordinates(st_centroid(amraGIS.nova)))
library("tools")
#amraGIS.nova1$shapeName <- toTitleCase(amraGIS.nova1$shapeName)
#head(amraGIS.nova1)
#HPI.map$segment <- ifelse(HPI.map$Munic1 == "Neum" | HPI.map$Munic1 == "Sarajevo - Centar"| HPI.map$Munic1 == "Vares","yes", NA)
#Finally, we can plot the segments using geom_segment().
#ggplot(HPI.map) +
# geom_sf(data = HPI.map, aes(geometry=geometry, fill = `Pricepersqm.no formula`)) +
# scale_fill_viridis_c(name = "Price per sqm") +
# geom_segment(data = subset(HPI.map, segment == "yes"), aes(geometry[[1]]),xend=geometry[[1]])