# Kod 4.1. Bezrobocie rejestrowane w roku 2019 w powiatach województw śląskiego i podkarpackiego -pozyskanie danych i weryfikacja hipotezy
library(bdl)
library(dplyr)
library(coin)
library(ggplot2)
rok=2019
no_uL=5 #
z<-as.data.frame(get_units(level=2))
woj1=z %>% dplyr::filter(name=="ŚLĄSKIE")
woj2=z %>% dplyr::filter(name=="PODKARPACKIE")
zmienna=60270 # stopa bezrobocia rejestrowanego
d1<-get_data_by_variable(zmienna,unitParentId=woj1[,'id'], unitLevel=no_uL,year=rok)
d2<-get_data_by_variable(zmienna,unitParentId=woj2[,'id'], unitLevel=no_uL,year=rok)
dane1=d1$val
dane2=d2$val
n1=length(d1$val)
n2=length(d2$val)
woj=factor(rep(c("ŚLĄSKIE", "PODKARPACKIE"), c(n1, n2)))
dane=c(dane1,dane2)
df=data.frame(dane,woj)
## Permutacyjny test dokładny
t_ex=oneway_test(dane ~ woj, data = df,distribution = "exact")
pvalue(t_ex)
statistic(t_ex)
województwo=woj
qplot(woj, dane, geom=c("boxplot"),
fill=województwo, main="Bezrobocie rejestrowane",
xlab="", ylab="%")
