# Kod 4.2. Realizacja testu permutacyjnego dokładnego dla dwóch wartości oczekiwanych, zmienna: mediana ceny m2 mieszkania w roku 2018 w województwach małopolskim, mazowieckim i śląskim
library(bdl)
library(dplyr)
library(ggplot2)
library(perm)
rok=2018
no_uL=5 #
z<-as.data.frame(get_units(level=2))
woj1=z %>% dplyr::filter(name=="MAŁOPOLSKIE")
woj2=z %>% dplyr::filter(name=="MAZOWIECKIE")
woj3=z %>% dplyr::filter(name=="ŚLĄSKIE")
zmienna=633687 # mediana ceny za 1m2 lokali mieszkalnych na rynku wtórnym ogółem
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)
d3<-get_data_by_variable(zmienna,unitParentId=woj3[,'id'], unitLevel=no_uL,year=rok)
dane1=d1$val
dane2=d2$val
dane3=d3$val
n1=length(d1$val)
n2=length(d2$val)
n3=length(d3$val)
województwo=factor(rep(c("MAŁOPOLSKIE", "MAZOWIECKIE","ŚLĄSKIE"), c(n1,n2,n3)))
dane=c(dane1,dane2,dane3)
df=data.frame(dane,województwo)
qplot(województwo, dane, geom=c("boxplot"),
fill=województwo, xlab="", ylab="Cena w zł.")
permKS(dane, województwo)
