In martinmodrak/revize-rs: Code for analysis of revize-rs project data
set.seed(181126)
knitr::opts_chunk$set(echo=FALSE)
library(tidyverse)
source(here::here("R/tools_sondy.R"))
sondy <- readRDS(here::here("private_data/sondy_180420.rds"))
Základní deskripce
Nejprve se podíváme na základní deskriptivní data o středisku.
Počty dat
Zda vubec nektere promenne ma smysl davat do analyz
table(sondy$je_rs_kmen_ve_stredisku) %>% as.tibble() %>% transmute(`Je kmen ve stredisku`=Var1,n) %>%
knitr::kable(caption = "Pocet odpovedi, kde existuje kmen ve stredisku")
table(sondy$je_rs_kmen_registrovany) %>% as.tibble() %>% transmute(`Je rs kmen registrovany`=Var1,n) %>% knitr::kable( caption = "Pocet odpovedi, zda je kmen registovany")
vlastni zvyklosti
rbind(table(sondy$vlastni_zvyklosti_obnova_slibu),
table(sondy$vlastni_zvyklosti_ritual_odchodu),
table(sondy$vlastni_zvyklosti_ritual_vstupu),
table(sondy$vlastni_zvyklosti_ritual_predavani_vudcovstvi),
table(sondy$vlastni_zvyklosti_nic)) %>% as_tibble() %>%
select(-`N/A`) %>% transmute(promenna = c("Obnova slibu", "Ritual odchodu", "Ritual vstupu", "Predavani vudcovstvi", "Nic"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro ritualy (803 lidi nevyplnilo)")
Ok, to znamena, ze ma smysl pouzivat jen ritual vstupu, ci rozdelovat lidi, zda maji/nemaji ritual.
Vychovne nastroje
rbind(table(sondy$vychovne_nastroje_vyzvy),
table(sondy$vychovne_nastroje_projekty),
table(sondy$vychovne_nastroje_rovernet),
table(sondy$vychovne_nastroje_kmen),
table(sondy$vychovne_nastroje_roversky_zacatek),
table(sondy$vychovne_nastroje_diar),
table(sondy$vychovne_nastroje_vlastni_stezka)) %>% as_tibble() %>%
select(-`N/A`) %>% transmute(promenna = c("Vyzvy", "Projekty", "Rovernet", "Kmen", "Roversky zacatek", "Diar", "Vlastni stezka"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro vychovne nastroje (803 lidi nevyplnilo)")
Tady to docela jde, pouze je teba davat pozor na konfidencni intervaly
Pomoc u dětských oddílů
TODO: jiný počet odpovědí než u členů RS.
rbind(table(sondy$pomoc_u_detskych_oddilu_jsem_stale_clenem),
table(sondy$pomoc_u_detskych_oddilu_nic),
table(sondy$pomoc_u_detskych_oddilu_vudce_oddilu),
table(sondy$pomoc_u_detskych_oddilu_vedu_oddil_papirove_nekdo_jiny),
table(sondy$pomoc_u_detskych_oddilu_zastupce_oddilu_nebo_oddilovy_radce),
table(sondy$pomoc_u_detskych_oddilu_obecna_pomoc_u_oddilu_bez_funkce),
table(sondy$pomoc_u_detskych_oddilu_radce_druziny),
table(sondy$pomoc_u_detskych_oddilu_pomoc_na_schuzkach),
table(sondy$pomoc_u_detskych_oddilu_pomoc_na_tabore),
table(sondy$pomoc_u_detskych_oddilu_pomoc_na_vypravach),
table(sondy$pomoc_u_detskych_oddilu_pomoc_na_klubovne_zakladne_brigadach),
table(sondy$pomoc_u_detskych_oddilu_jine)) %>% as_tibble() %>%
select(-`N/A`) %>%
transmute(promenna = c("Clennem u oddilu", "Nic", "Vudce oddilu", "Vedu oddil pouze papirove", "Zastupce oddilu nebo oddilovy radce", "Obecna pomoc u oddilu", "Radce druziny", "Pomoc na schuzkach", "Pomoc na tabore", "Pomoc na vypravach", "Pomoc na klubovne/zakladne/brigadach", "Jine"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro pomoc u detskych oddilu (424 nevyplnilo)")
Prekvapuje me, kolik lidi je ve vedeni oddilu! Nejak pomaha hodne lidi, jinak je asi taky ok
Probehle akce
rbind(table(sondy$roverske_akce_probehle_2roky_zavody),
table(sondy$roverske_akce_diskuze),
table(sondy$roverske_akce_probehle_expedice),
table(sondy$roverske_akce_probehle_tabor),
table(sondy$roverske_akce_duchovni_programy),
table(sondy$roverske_akce_pobyt_v_prirode),
table(sondy$roverske_akce_zahranicni_akce),
table(sondy$roverske_akce_kulturni_akce),
table(sondy$roverske_akce_tradicni_akce),
table(sondy$roverske_akce_setkani_bez_programu_hospoda),
table(sondy$roverske_akce_sluzba),
table(sondy$roverske_akce_velke_rs_akce),
table(sondy$roverske_akce_kurzy_vzdelavaci_akce)) %>%
as_tibble() %>%
select(-`N/A`) %>%
transmute(promenna = c("Probehle zavody behem 2 let", "Diskuze", "Expedice", "Tabor", "Duchovni programy", "Pobyt v prirode", "Zahranicni akce", "Kulturni akce", "Tradicni akce", "Setkani bez programu/hospoda", "Sluzba", "Velke RS akce", "Kurzy/Vzdelavaci akce"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro probehle akce (803 nevyplnilo)")
Taky asi mame dost programu
Proc neni registrovany
Nasledujici graf ukazje procenta uzivatelu, kteri vybrali dany duvod neregistrace kmene. Lide mohli zaskrtavat vice odpovedi, proto nesedi soucet do sta procent. Svisla cara oznacuje standardni chybu prumeru
sondy %>%
select(rs_kmen_registace_neni_registrovany_neni_vudce:rs_kmen_registrace_jine) %>%
gather(key = "duvod") %>%
filter(value!="N/A") %>%
mutate(value = if_else(value == "Ano", 1, 0)) %>%
ggplot(aes(x=duvod, y = value)) + stat_summary(fun.data = "mean_cl_boot") +
theme(axis.text.x = element_text(angle = 45, hjust = 1)) + ylim(0,1) +
scale_x_discrete(labels = c("Neni dost clenu", "Nema vudce", "Clenove se neregistruji", "Clenove se se registruji pod strediskem", "clenove se registruji u oddilu", "Jine", "Registace cinnosti nepomuze"))
Nejvice lidi se neregistruje, protoze pracuji u oddilu. Zajimava je i cast, ktera povazuje registraci za zbytecnou.
Ucast na rs - kmeans
V datech mame popsany duvody, proc se lide ucastni akci. Kazdy mohl procenta rozdelit dle sebe. Je zde tedy moznost, ze tu jsou urcite "skupiny" lidi, ktere chodi na akce z podobnych duvodu. Proto jsme zkusili udelat k-means, ktera nam data rozdeli do do k skupin, aby si byly co nejblize. Protoze nevime, kolik by melo byt takovychto skupin, budeme zkouset rozdeleni do 2 az 5 skupin
sondy_ucast_tmp <- sondy %>% select(id, ucast_rs_kmen_procenta_lide:ucast_rs_kmen_procenta_neco_jineho) %>% filter(!is.na(ucast_rs_kmen_procenta_lide))
sondy_ucast <- sondy_ucast_tmp %>%
mutate(skupina_k2 = kmeans(sondy_ucast_tmp,2)$cluster,
skupina_k3 = kmeans(sondy_ucast_tmp,3)$cluster,
skupina_k4 = kmeans(sondy_ucast_tmp,4)$cluster,
skupina_k5 = kmeans(sondy_ucast_tmp,5)$cluster)
lab_ucast <- c("Lide", "Zajimavy program", "Odpocinek od oddilu", "Moznost realizace", "Setrvacnost", "Jine")
summarize_group(sondy_ucast) %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100)
summarize_group(sondy_ucast, "skupina_k2") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100)
summarize_group(sondy_ucast, "skupina_k3") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100)
summarize_group(sondy_ucast, "skupina_k4") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100)
summarize_group(sondy_ucast, "skupina_k5") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100)
Nevychazi to nijak zajimave. Obecne je nejvic lidi registrovanych kvuli lidem.
Neucast na rs - kmeans
To stejne muzeme zjistit pro neucast
sondy_neucast_tmp <- sondy %>% select(neucast_rs_kmen_procenta_prace_na_stredisku:neucast_rs_kmen_procenta_jine) %>%
filter(!is.na(neucast_rs_kmen_procenta_prace_na_stredisku))
sondy_neucast <- sondy_neucast_tmp %>%
mutate(skupina_k2 = kmeans(sondy_neucast_tmp,2)$cluster,
skupina_k3 = kmeans(sondy_neucast_tmp,3)$cluster,
skupina_k4 = kmeans(sondy_neucast_tmp,4)$cluster,
skupina_k5 = kmeans(sondy_neucast_tmp,5)$cluster)
lab_neucast <- c("Cinovnicka funkce", "Jina parta", "Jine", "Nemam cas", "Nezajimavy program", "Odchod na SS/VS", "Prace na stredisku", "Rodinny zivot", "Skolni povinnosti", "Vztahy ve spolecenstvi")
summarize_group_neucast(sondy_neucast, "skupina_k2") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100)
summarize_group_neucast(sondy_neucast, "skupina_k3") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100)
summarize_group_neucast(sondy_neucast, "skupina_k4") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100)
summarize_group_neucast(sondy_neucast, "skupina_k5") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100)
Tady to vychazi lepe, jsou zde ruzne skupiny lidi, ktere se lid dle duvodu. Napr. pro pet clusteru je skupina 1 "cinovnicka", skupina 2 "obecny nedostatek casu", skupina 3 "zmena skoly", skupina 5 "jina parta". Zajimava je skupina 4, ktera se vyskytuje u vsech clusteru, u ktere nezname duvod.
TODO: tady by možná stálo zjistit, kolik lidí vůbec dalo výrazný podíl více věcem naráz...
Spokojenost s programem
v <- vars(spokojenost_s_frekvenci_akci2, celkova_spokojenost_s_programem_rs_kmenu, frekvence_rs_akci_celodenni_akce, frekvence_rs_akci_schuzky, frekvence_rs_akci_vicedenni_akce, frekvence_rs_akci_tabory_expedice,vychovne_nastroje_vyzvy:vychovne_nastroje_vlastni_stezka, vlastni_zvyklosti_nic,vlastni_zvyklosti_ritual_vstupu,rs_kmen_prvky_vudce:rs_kmen_prvky_ani_jedno)
aa <- sondy %>% mutate(celkova_spokojenost_s_programem_rs_kmenu = recode(celkova_spokojenost_s_programem_rs_kmenu, `Velmi nespokojen/a`= 1, `Spíše nespokojen/a` = 2, `Spíše spokojen/a` = 3, `Spíše spokojen/a` = 4, `Velmi spokojen/a` = 5, .default = NA_real_),
spokojenost_s_frekvenci_akci2 = recode(spokojenost_s_frekvenci_akci, `Jsem spokojen/a`= 1, `Nejsem spokojen/a - přál/a bych si akce častěji` = 0, `Nejsem spokojen/a - přál/a bych si akce méně často` = 0, .default = NA_real_),
frekvence_rs_akci_celodenni_akce = recode(frekvence_rs_akci_celodenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_),
frekvence_rs_akci_schuzky = recode(frekvence_rs_akci_schuzky, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_),
frekvence_rs_akci_vicedenni_akce = recode(frekvence_rs_akci_vicedenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, .default = NA_real_),
frekvence_rs_akci_tabory_expedice = recode(frekvence_rs_akci_tabory_expedice, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, .default = NA_real_),
vychovne_nastroje_vyzvy = recode(vychovne_nastroje_vyzvy, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_projekty = recode(vychovne_nastroje_projekty, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_rovernet = recode(vychovne_nastroje_rovernet, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_kmen = recode(vychovne_nastroje_kmen, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_roversky_zacatek = recode(vychovne_nastroje_roversky_zacatek, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_diar = recode(vychovne_nastroje_diar, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_vlastni_stezka = recode(vychovne_nastroje_vlastni_stezka, `Ano` = 1, `Ne` = 0, .default = NA_real_),
ma_kmen_vychovny_plan = recode(ma_kmen_vychovny_plan, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_nic = recode(vlastni_zvyklosti_nic, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_vstupu = recode(vlastni_zvyklosti_ritual_vstupu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_obnova_slibu = recode(vlastni_zvyklosti_obnova_slibu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_odchodu = recode(vlastni_zvyklosti_ritual_odchodu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_predavani_vudcovstvi = recode(vlastni_zvyklosti_ritual_predavani_vudcovstvi, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_vudce = recode(rs_kmen_prvky_vudce, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_rada = recode(rs_kmen_prvky_rada, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_ani_jedno = recode(rs_kmen_prvky_ani_jedno, `Ano` = 1, `Ne` = 0, .default = NA_real_)
)
c <- aa %>% select_at(v) %>%
filter(!is.na(spokojenost_s_frekvenci_akci2)) %>% as.matrix() %>% Hmisc::rcorr(type = "spearman")
format_corrtable(c) %>% knitr::kable()
Podil prace ve stredisku
TODO: ke clusterům počty lidí v nich
sondy_podil_casu_tmp <- sondy %>% select(id, podil_casu_ve_skautu_u_detskych_oddilu:podil_casu_ve_skautu_jine) %>%
filter(!is.na(podil_casu_ve_skautu_u_detskych_oddilu))
sondy_podil_casu <- sondy_podil_casu_tmp %>%
mutate(skupina_k2 = kmeans(sondy_podil_casu_tmp %>% select(-id),2)$cluster,
skupina_k3 = kmeans(sondy_podil_casu_tmp %>% select(-id),3)$cluster,
skupina_k4 = kmeans(sondy_podil_casu_tmp %>% select(-id),4)$cluster,
skupina_k5 = kmeans(sondy_podil_casu_tmp %>% select(-id),5)$cluster)
sondy_podil_casu %>% summarize_group_podil_prace() %>% plot_sondy_clusters(lab = c("Jine", "Roverska cinnost", "U detskych oddilu", "Ucast/pomoc celostatni a vzdelavaci"), ylim_min = 0, ylim_max = 100)
sondy_podil_casu %>% summarize_group_podil_prace("skupina_k2") %>% plot_sondy_podil_casu()
sondy_podil_casu %>% summarize_group_podil_prace("skupina_k3") %>% plot_sondy_podil_casu()
sondy_podil_casu %>% summarize_group_podil_prace("skupina_k4") %>% plot_sondy_podil_casu()
sondy_podil_casu %>% summarize_group_podil_prace("skupina_k5") %>% plot_sondy_podil_casu()
Jak se lisi korelacni matice u jednotlivych subskupin
Vezmeme podskupiny podle kmeans a replikujeme na nich korelacni matice
sondy2 <- sondy %>% left_join(sondy_podil_casu %>% select(-(podil_casu_ve_skautu_u_detskych_oddilu:podil_casu_ve_skautu_jine)), by = "id")
aa <- sondy2 %>% mutate(celkova_spokojenost_s_programem_rs_kmenu = recode(celkova_spokojenost_s_programem_rs_kmenu, `Velmi nespokojen/a`= 1, `Spíše nespokojen/a` = 2, `Spíše spokojen/a` = 3, `Spíše spokojen/a` = 4, `Velmi spokojen/a` = 5, .default = NA_real_),spokojenost_s_frekvenci_akci2 = recode(spokojenost_s_frekvenci_akci, `Jsem spokojen/a`= 1, `Nejsem spokojen/a - přál/a bych si akce častěji` = 0, `Nejsem spokojen/a - přál/a bych si akce méně často` = 0, .default = NA_real_),
frekvence_rs_akci_celodenni_akce = recode(frekvence_rs_akci_celodenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_),
frekvence_rs_akci_schuzky = recode(frekvence_rs_akci_schuzky, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_),
frekvence_rs_akci_vicedenni_akce = recode(frekvence_rs_akci_vicedenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, .default = NA_real_),
frekvence_rs_akci_tabory_expedice = recode(frekvence_rs_akci_tabory_expedice, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, .default = NA_real_),
vychovne_nastroje_vyzvy = recode(vychovne_nastroje_vyzvy, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_projekty = recode(vychovne_nastroje_projekty, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_rovernet = recode(vychovne_nastroje_rovernet, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_kmen = recode(vychovne_nastroje_kmen, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_roversky_zacatek = recode(vychovne_nastroje_roversky_zacatek, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_diar = recode(vychovne_nastroje_diar, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_vlastni_stezka = recode(vychovne_nastroje_vlastni_stezka, `Ano` = 1, `Ne` = 0, .default = NA_real_),
ma_kmen_vychovny_plan = recode(ma_kmen_vychovny_plan, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_nic = recode(vlastni_zvyklosti_nic, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_vstupu = recode(vlastni_zvyklosti_ritual_vstupu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_obnova_slibu = recode(vlastni_zvyklosti_obnova_slibu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_odchodu = recode(vlastni_zvyklosti_ritual_odchodu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_predavani_vudcovstvi = recode(vlastni_zvyklosti_ritual_predavani_vudcovstvi, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_vudce = recode(rs_kmen_prvky_vudce, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_rada = recode(rs_kmen_prvky_rada, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_ani_jedno = recode(rs_kmen_prvky_ani_jedno, `Ano` = 1, `Ne` = 0, .default = NA_real_)
)
4 skupiny - rozdily
aa2 <- aa %>% filter(!is.na(skupina_k5))
aa2$skupina_k5 <- recode(aa2$skupina_k5, `1`="detska-vzdelavaci",`2`="prevazne-detska",`3`="jina cinnost",`4`="roverska", `5`="detska")
vs <- aa %>% select_at(v) %>% colnames()
for (i in 1) {
p <- aa2 %>% ggplot(aes_string(x = "skupina_k5", y = vs[i])) + stat_summary(fun.data = "mean_cl_boot") + ylab("Procenta") + ylim(0,1) + ggtitle(vs[i])
print(p)
}
for (i in 2:6) {
p <- aa2 %>% ggplot(aes_string(x = "skupina_k5", y = vs[i])) + stat_summary(fun.data = "mean_cl_boot") + ylab("Procenta") + ylim(0,6) + ggtitle(vs[i])
print(p)
}
for (i in 7:length(vs)) {
p <- aa2 %>% ggplot(aes_string(x = "skupina_k5", y = vs[i])) + stat_summary(fun.data = "mean_cl_boot") + ylab("Procenta") + ylim(0,1) + ggtitle(vs[i])
print(p)
}
2 skupiny
subgroup_corr(aa,v,vars(skupina_k2), 1) %>% knitr::kable(caption="dve skupiny, skupina 1")
subgroup_corr(aa,v,vars(skupina_k2), 2) %>% knitr::kable(caption="dve skupiny, skupina 2")
3 skupiny
subgroup_corr(aa,v,vars(skupina_k3), 1) %>% knitr::kable(caption="tri skupiny, skupina 1")
subgroup_corr(aa,v,vars(skupina_k3), 2) %>% knitr::kable(caption="tri skupiny, skupina 2")
subgroup_corr(aa,v,vars(skupina_k3), 2) %>% knitr::kable(caption="tri skupiny, skupina 3")
4 skupiny
subgroup_corr(aa,v,vars(skupina_k4), 1) %>% knitr::kable(caption="ctyri skupiny, skupina 1")
subgroup_corr(aa,v,vars(skupina_k4), 2) %>% knitr::kable(caption="ctyri skupiny, skupina 2")
subgroup_corr(aa,v,vars(skupina_k4), 3) %>% knitr::kable(caption="ctyri skupiny, skupina 3")
subgroup_corr(aa,v,vars(skupina_k4), 4) %>% knitr::kable(caption="ctyri skupiny, skupina 4")
5 skupin
subgroup_corr(aa,v,vars(skupina_k5), 1) %>% knitr::kable(caption="pet skupin, skupina 1")
subgroup_corr(aa,v,vars(skupina_k5), 2) %>% knitr::kable(caption="pet skupin, skupina 2")
subgroup_corr(aa,v,vars(skupina_k5), 3) %>% knitr::kable(caption="pet skupin, skupina 3")
subgroup_corr(aa,v,vars(skupina_k5), 4) %>% knitr::kable(caption="pet skupin, skupina 4")
subgroup_corr(aa,v,vars(skupina_k5), 5) %>% knitr::kable(caption="pet skupin, skupina 5")
Clusterováni na základě typu proběhlých programů
cluster_var <- vars(roverske_akce_probehle_2roky_zavody:roverske_akce_kurzy_vzdelavaci_akce)
sondy_program_tmp <- sondy %>%
select_at(cluster_var) %>%
mutate_all(funs(recode(.,"Ano"=T,"Ne"=F,"N/A"=NA))) %>%
mutate(id = sondy$id) %>%
filter(!is.na(roverske_akce_probehle_2roky_zavody))
sondy_program <- sondy_program_tmp %>%
mutate(skupina_k2 = kmeans(sondy_program_tmp %>% select(-id),2)$cluster,
skupina_k3 = kmeans(sondy_program_tmp %>% select(-id),3)$cluster,
skupina_k4 = kmeans(sondy_program_tmp %>% select(-id),4)$cluster,
skupina_k5 = kmeans(sondy_program_tmp %>% select(-id),5)$cluster,
skupina_k6 = kmeans(sondy_program_tmp %>% select(-id),6)$cluster,
skupina_k7 = kmeans(sondy_program_tmp %>% select(-id),7)$cluster,
skupina_k8 = kmeans(sondy_program_tmp %>% select(-id),8)$cluster,
skupina_k9 = kmeans(sondy_program_tmp %>% select(-id),9)$cluster,
skupina_k10 = kmeans(sondy_program_tmp %>% select(-id),10)$cluster)
l <- c("zavody", "expedice", "tabor", "programy", "pobyt v prirode", "zahranicni akce", "kulturni akce", "diskuze", "tradicni akce", "setkani bez programu/hospoda", "sluzba", "velke rs akce", "vzdelavaci akce")
sondy_program %>% summarize_group_typ_programu() %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k2") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k3") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k4") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k5") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k6") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k7") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k8") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k9") %>% plot_sondy_clusters(l)
sondy_program %>% summarize_group_typ_programu("skupina_k10") %>% plot_sondy_clusters(l)
Korelace per akcni kmeny
ab <- sondy %>% left_join(sondy_program %>% select(id, skupina_k2:skupina_k10), by = "id") %>% mutate(celkova_spokojenost_s_programem_rs_kmenu = recode(celkova_spokojenost_s_programem_rs_kmenu, `Velmi nespokojen/a`= 1, `Spíše nespokojen/a` = 2, `Spíše spokojen/a` = 3, `Spíše spokojen/a` = 4, `Velmi spokojen/a` = 5, .default = NA_real_),
spokojenost_s_frekvenci_akci2 = recode(spokojenost_s_frekvenci_akci, `Jsem spokojen/a`= 1, `Nejsem spokojen/a - přál/a bych si akce častěji` = 0, `Nejsem spokojen/a - přál/a bych si akce méně často` = 0, .default = NA_real_),
frekvence_rs_akci_celodenni_akce = recode(frekvence_rs_akci_celodenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_),
frekvence_rs_akci_schuzky = recode(frekvence_rs_akci_schuzky, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_),
frekvence_rs_akci_vicedenni_akce = recode(frekvence_rs_akci_vicedenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, .default = NA_real_),
frekvence_rs_akci_tabory_expedice = recode(frekvence_rs_akci_tabory_expedice, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, .default = NA_real_),
vychovne_nastroje_vyzvy = recode(vychovne_nastroje_vyzvy, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_projekty = recode(vychovne_nastroje_projekty, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_rovernet = recode(vychovne_nastroje_rovernet, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_kmen = recode(vychovne_nastroje_kmen, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_roversky_zacatek = recode(vychovne_nastroje_roversky_zacatek, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_diar = recode(vychovne_nastroje_diar, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vychovne_nastroje_vlastni_stezka = recode(vychovne_nastroje_vlastni_stezka, `Ano` = 1, `Ne` = 0, .default = NA_real_),
ma_kmen_vychovny_plan = recode(ma_kmen_vychovny_plan, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_nic = recode(vlastni_zvyklosti_nic, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_vstupu = recode(vlastni_zvyklosti_ritual_vstupu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_obnova_slibu = recode(vlastni_zvyklosti_obnova_slibu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_odchodu = recode(vlastni_zvyklosti_ritual_odchodu, `Ano` = 1, `Ne` = 0, .default = NA_real_),
vlastni_zvyklosti_ritual_predavani_vudcovstvi = recode(vlastni_zvyklosti_ritual_predavani_vudcovstvi, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_vudce = recode(rs_kmen_prvky_vudce, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_rada = recode(rs_kmen_prvky_rada, `Ano` = 1, `Ne` = 0, .default = NA_real_),
rs_kmen_prvky_ani_jedno = recode(rs_kmen_prvky_ani_jedno, `Ano` = 1, `Ne` = 0, .default = NA_real_)
)
v <- vars(spokojenost_s_frekvenci_akci2, celkova_spokojenost_s_programem_rs_kmenu, frekvence_rs_akci_celodenni_akce, frekvence_rs_akci_schuzky, frekvence_rs_akci_vicedenni_akce, frekvence_rs_akci_tabory_expedice,vychovne_nastroje_vyzvy:vychovne_nastroje_vlastni_stezka, vlastni_zvyklosti_nic:vlastni_zvyklosti_ritual_predavani_vudcovstvi,rs_kmen_prvky_vudce:rs_kmen_prvky_ani_jedno)
subgroup_corr(ab,v,vars(skupina_k2), 1) %>% knitr::kable(caption="dve skupiny, skupina 1")
subgroup_corr(ab,v,vars(skupina_k2), 2) %>% knitr::kable(caption="dve skupiny, skupina 2")
subgroup_corr(ab, v, vars(skupina_k3), 1) %>% knitr::kable(caption="tri skupiny, skupina 1")
subgroup_corr(ab, v, vars(skupina_k3), 2) %>% knitr::kable(caption="tri skupiny, skupina 2")
subgroup_corr(ab, v, vars(skupina_k3), 3) %>% knitr::kable(caption="tri skupiny, skupina 3")
subgroup_corr(ab, v, vars(skupina_k4), 1) %>% knitr::kable(caption="ctyri skupiny, skupina 1")
subgroup_corr(ab, v, vars(skupina_k4), 2) %>% knitr::kable(caption="ctyri skupiny, skupina 2")
subgroup_corr(ab, v, vars(skupina_k4), 3) %>% knitr::kable(caption="ctyri skupiny, skupina 3")
subgroup_corr(ab, v, vars(skupina_k4), 4) %>% knitr::kable(caption="ctyri skupiny, skupina 4")
Korelace mezi podilem cinnosti a typem programu
podilcasu_program <- sondy_podil_casu %>% rename(podilcasu2 = skupina_k2,podilcasu3 = skupina_k3,
podilcasu4 = skupina_k4,podilcasu5 = skupina_k5) %>% left_join(sondy_program, by = "id")
podilcasu_program %>% select(podilcasu2:podilcasu5,skupina_k2:skupina_k10) %>% as.matrix() %>% Hmisc::rcorr() %>% format_corrtable() %>% knitr::kable()
#xtabs(~podilcasu5+skupina_k2, podilcasu_program) %>% knitr::kable()
#xtabs(~podilcasu3+skupina_k2, podilcasu_program)%>% knitr::kable()
Nvychazi tam zadna zajimava souvislost
martinmodrak/revize-rs documentation built on March 9, 2021, 5:30 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
set.seed(181126) knitr::opts_chunk$set(echo=FALSE) library(tidyverse) source(here::here("R/tools_sondy.R"))
sondy <- readRDS(here::here("private_data/sondy_180420.rds"))
Základní deskripce
Nejprve se podíváme na základní deskriptivní data o středisku.
Počty dat
Zda vubec nektere promenne ma smysl davat do analyz
table(sondy$je_rs_kmen_ve_stredisku) %>% as.tibble() %>% transmute(`Je kmen ve stredisku`=Var1,n) %>% knitr::kable(caption = "Pocet odpovedi, kde existuje kmen ve stredisku")
table(sondy$je_rs_kmen_registrovany) %>% as.tibble() %>% transmute(`Je rs kmen registrovany`=Var1,n) %>% knitr::kable( caption = "Pocet odpovedi, zda je kmen registovany")
vlastni zvyklosti
rbind(table(sondy$vlastni_zvyklosti_obnova_slibu), table(sondy$vlastni_zvyklosti_ritual_odchodu), table(sondy$vlastni_zvyklosti_ritual_vstupu), table(sondy$vlastni_zvyklosti_ritual_predavani_vudcovstvi), table(sondy$vlastni_zvyklosti_nic)) %>% as_tibble() %>% select(-`N/A`) %>% transmute(promenna = c("Obnova slibu", "Ritual odchodu", "Ritual vstupu", "Predavani vudcovstvi", "Nic"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro ritualy (803 lidi nevyplnilo)")
Ok, to znamena, ze ma smysl pouzivat jen ritual vstupu, ci rozdelovat lidi, zda maji/nemaji ritual.
Vychovne nastroje
rbind(table(sondy$vychovne_nastroje_vyzvy), table(sondy$vychovne_nastroje_projekty), table(sondy$vychovne_nastroje_rovernet), table(sondy$vychovne_nastroje_kmen), table(sondy$vychovne_nastroje_roversky_zacatek), table(sondy$vychovne_nastroje_diar), table(sondy$vychovne_nastroje_vlastni_stezka)) %>% as_tibble() %>% select(-`N/A`) %>% transmute(promenna = c("Vyzvy", "Projekty", "Rovernet", "Kmen", "Roversky zacatek", "Diar", "Vlastni stezka"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro vychovne nastroje (803 lidi nevyplnilo)")
Tady to docela jde, pouze je teba davat pozor na konfidencni intervaly
Pomoc u dětských oddílů
TODO: jiný počet odpovědí než u členů RS.
rbind(table(sondy$pomoc_u_detskych_oddilu_jsem_stale_clenem), table(sondy$pomoc_u_detskych_oddilu_nic), table(sondy$pomoc_u_detskych_oddilu_vudce_oddilu), table(sondy$pomoc_u_detskych_oddilu_vedu_oddil_papirove_nekdo_jiny), table(sondy$pomoc_u_detskych_oddilu_zastupce_oddilu_nebo_oddilovy_radce), table(sondy$pomoc_u_detskych_oddilu_obecna_pomoc_u_oddilu_bez_funkce), table(sondy$pomoc_u_detskych_oddilu_radce_druziny), table(sondy$pomoc_u_detskych_oddilu_pomoc_na_schuzkach), table(sondy$pomoc_u_detskych_oddilu_pomoc_na_tabore), table(sondy$pomoc_u_detskych_oddilu_pomoc_na_vypravach), table(sondy$pomoc_u_detskych_oddilu_pomoc_na_klubovne_zakladne_brigadach), table(sondy$pomoc_u_detskych_oddilu_jine)) %>% as_tibble() %>% select(-`N/A`) %>% transmute(promenna = c("Clennem u oddilu", "Nic", "Vudce oddilu", "Vedu oddil pouze papirove", "Zastupce oddilu nebo oddilovy radce", "Obecna pomoc u oddilu", "Radce druziny", "Pomoc na schuzkach", "Pomoc na tabore", "Pomoc na vypravach", "Pomoc na klubovne/zakladne/brigadach", "Jine"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro pomoc u detskych oddilu (424 nevyplnilo)")
Prekvapuje me, kolik lidi je ve vedeni oddilu! Nejak pomaha hodne lidi, jinak je asi taky ok
Probehle akce
rbind(table(sondy$roverske_akce_probehle_2roky_zavody), table(sondy$roverske_akce_diskuze), table(sondy$roverske_akce_probehle_expedice), table(sondy$roverske_akce_probehle_tabor), table(sondy$roverske_akce_duchovni_programy), table(sondy$roverske_akce_pobyt_v_prirode), table(sondy$roverske_akce_zahranicni_akce), table(sondy$roverske_akce_kulturni_akce), table(sondy$roverske_akce_tradicni_akce), table(sondy$roverske_akce_setkani_bez_programu_hospoda), table(sondy$roverske_akce_sluzba), table(sondy$roverske_akce_velke_rs_akce), table(sondy$roverske_akce_kurzy_vzdelavaci_akce)) %>% as_tibble() %>% select(-`N/A`) %>% transmute(promenna = c("Probehle zavody behem 2 let", "Diskuze", "Expedice", "Tabor", "Duchovni programy", "Pobyt v prirode", "Zahranicni akce", "Kulturni akce", "Tradicni akce", "Setkani bez programu/hospoda", "Sluzba", "Velke RS akce", "Kurzy/Vzdelavaci akce"), Ano, Ne) %>% knitr::kable(caption = "Pocty dat pro probehle akce (803 nevyplnilo)")
Taky asi mame dost programu
Proc neni registrovany
Nasledujici graf ukazje procenta uzivatelu, kteri vybrali dany duvod neregistrace kmene. Lide mohli zaskrtavat vice odpovedi, proto nesedi soucet do sta procent. Svisla cara oznacuje standardni chybu prumeru
sondy %>% select(rs_kmen_registace_neni_registrovany_neni_vudce:rs_kmen_registrace_jine) %>% gather(key = "duvod") %>% filter(value!="N/A") %>% mutate(value = if_else(value == "Ano", 1, 0)) %>% ggplot(aes(x=duvod, y = value)) + stat_summary(fun.data = "mean_cl_boot") + theme(axis.text.x = element_text(angle = 45, hjust = 1)) + ylim(0,1) + scale_x_discrete(labels = c("Neni dost clenu", "Nema vudce", "Clenove se neregistruji", "Clenove se se registruji pod strediskem", "clenove se registruji u oddilu", "Jine", "Registace cinnosti nepomuze"))
Nejvice lidi se neregistruje, protoze pracuji u oddilu. Zajimava je i cast, ktera povazuje registraci za zbytecnou.
Ucast na rs - kmeans
V datech mame popsany duvody, proc se lide ucastni akci. Kazdy mohl procenta rozdelit dle sebe. Je zde tedy moznost, ze tu jsou urcite "skupiny" lidi, ktere chodi na akce z podobnych duvodu. Proto jsme zkusili udelat k-means, ktera nam data rozdeli do do k skupin, aby si byly co nejblize. Protoze nevime, kolik by melo byt takovychto skupin, budeme zkouset rozdeleni do 2 az 5 skupin
sondy_ucast_tmp <- sondy %>% select(id, ucast_rs_kmen_procenta_lide:ucast_rs_kmen_procenta_neco_jineho) %>% filter(!is.na(ucast_rs_kmen_procenta_lide)) sondy_ucast <- sondy_ucast_tmp %>% mutate(skupina_k2 = kmeans(sondy_ucast_tmp,2)$cluster, skupina_k3 = kmeans(sondy_ucast_tmp,3)$cluster, skupina_k4 = kmeans(sondy_ucast_tmp,4)$cluster, skupina_k5 = kmeans(sondy_ucast_tmp,5)$cluster) lab_ucast <- c("Lide", "Zajimavy program", "Odpocinek od oddilu", "Moznost realizace", "Setrvacnost", "Jine") summarize_group(sondy_ucast) %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100) summarize_group(sondy_ucast, "skupina_k2") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100) summarize_group(sondy_ucast, "skupina_k3") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100) summarize_group(sondy_ucast, "skupina_k4") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100) summarize_group(sondy_ucast, "skupina_k5") %>% plot_sondy_clusters(lab_ucast, ylim_min = 0, ylim_max = 100)
Nevychazi to nijak zajimave. Obecne je nejvic lidi registrovanych kvuli lidem.
Neucast na rs - kmeans
To stejne muzeme zjistit pro neucast
sondy_neucast_tmp <- sondy %>% select(neucast_rs_kmen_procenta_prace_na_stredisku:neucast_rs_kmen_procenta_jine) %>% filter(!is.na(neucast_rs_kmen_procenta_prace_na_stredisku)) sondy_neucast <- sondy_neucast_tmp %>% mutate(skupina_k2 = kmeans(sondy_neucast_tmp,2)$cluster, skupina_k3 = kmeans(sondy_neucast_tmp,3)$cluster, skupina_k4 = kmeans(sondy_neucast_tmp,4)$cluster, skupina_k5 = kmeans(sondy_neucast_tmp,5)$cluster) lab_neucast <- c("Cinovnicka funkce", "Jina parta", "Jine", "Nemam cas", "Nezajimavy program", "Odchod na SS/VS", "Prace na stredisku", "Rodinny zivot", "Skolni povinnosti", "Vztahy ve spolecenstvi") summarize_group_neucast(sondy_neucast, "skupina_k2") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100) summarize_group_neucast(sondy_neucast, "skupina_k3") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100) summarize_group_neucast(sondy_neucast, "skupina_k4") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100) summarize_group_neucast(sondy_neucast, "skupina_k5") %>% plot_sondy_clusters(lab_neucast, ylim_min = 0, ylim_max = 100)
Tady to vychazi lepe, jsou zde ruzne skupiny lidi, ktere se lid dle duvodu. Napr. pro pet clusteru je skupina 1 "cinovnicka", skupina 2 "obecny nedostatek casu", skupina 3 "zmena skoly", skupina 5 "jina parta". Zajimava je skupina 4, ktera se vyskytuje u vsech clusteru, u ktere nezname duvod.
TODO: tady by možná stálo zjistit, kolik lidí vůbec dalo výrazný podíl více věcem naráz...
Spokojenost s programem
v <- vars(spokojenost_s_frekvenci_akci2, celkova_spokojenost_s_programem_rs_kmenu, frekvence_rs_akci_celodenni_akce, frekvence_rs_akci_schuzky, frekvence_rs_akci_vicedenni_akce, frekvence_rs_akci_tabory_expedice,vychovne_nastroje_vyzvy:vychovne_nastroje_vlastni_stezka, vlastni_zvyklosti_nic,vlastni_zvyklosti_ritual_vstupu,rs_kmen_prvky_vudce:rs_kmen_prvky_ani_jedno) aa <- sondy %>% mutate(celkova_spokojenost_s_programem_rs_kmenu = recode(celkova_spokojenost_s_programem_rs_kmenu, `Velmi nespokojen/a`= 1, `Spíše nespokojen/a` = 2, `Spíše spokojen/a` = 3, `Spíše spokojen/a` = 4, `Velmi spokojen/a` = 5, .default = NA_real_), spokojenost_s_frekvenci_akci2 = recode(spokojenost_s_frekvenci_akci, `Jsem spokojen/a`= 1, `Nejsem spokojen/a - přál/a bych si akce častěji` = 0, `Nejsem spokojen/a - přál/a bych si akce méně často` = 0, .default = NA_real_), frekvence_rs_akci_celodenni_akce = recode(frekvence_rs_akci_celodenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_), frekvence_rs_akci_schuzky = recode(frekvence_rs_akci_schuzky, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_), frekvence_rs_akci_vicedenni_akce = recode(frekvence_rs_akci_vicedenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, .default = NA_real_), frekvence_rs_akci_tabory_expedice = recode(frekvence_rs_akci_tabory_expedice, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, .default = NA_real_), vychovne_nastroje_vyzvy = recode(vychovne_nastroje_vyzvy, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_projekty = recode(vychovne_nastroje_projekty, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_rovernet = recode(vychovne_nastroje_rovernet, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_kmen = recode(vychovne_nastroje_kmen, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_roversky_zacatek = recode(vychovne_nastroje_roversky_zacatek, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_diar = recode(vychovne_nastroje_diar, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_vlastni_stezka = recode(vychovne_nastroje_vlastni_stezka, `Ano` = 1, `Ne` = 0, .default = NA_real_), ma_kmen_vychovny_plan = recode(ma_kmen_vychovny_plan, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_nic = recode(vlastni_zvyklosti_nic, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_vstupu = recode(vlastni_zvyklosti_ritual_vstupu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_obnova_slibu = recode(vlastni_zvyklosti_obnova_slibu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_odchodu = recode(vlastni_zvyklosti_ritual_odchodu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_predavani_vudcovstvi = recode(vlastni_zvyklosti_ritual_predavani_vudcovstvi, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_vudce = recode(rs_kmen_prvky_vudce, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_rada = recode(rs_kmen_prvky_rada, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_ani_jedno = recode(rs_kmen_prvky_ani_jedno, `Ano` = 1, `Ne` = 0, .default = NA_real_) ) c <- aa %>% select_at(v) %>% filter(!is.na(spokojenost_s_frekvenci_akci2)) %>% as.matrix() %>% Hmisc::rcorr(type = "spearman") format_corrtable(c) %>% knitr::kable()
Podil prace ve stredisku
TODO: ke clusterům počty lidí v nich
sondy_podil_casu_tmp <- sondy %>% select(id, podil_casu_ve_skautu_u_detskych_oddilu:podil_casu_ve_skautu_jine) %>% filter(!is.na(podil_casu_ve_skautu_u_detskych_oddilu)) sondy_podil_casu <- sondy_podil_casu_tmp %>% mutate(skupina_k2 = kmeans(sondy_podil_casu_tmp %>% select(-id),2)$cluster, skupina_k3 = kmeans(sondy_podil_casu_tmp %>% select(-id),3)$cluster, skupina_k4 = kmeans(sondy_podil_casu_tmp %>% select(-id),4)$cluster, skupina_k5 = kmeans(sondy_podil_casu_tmp %>% select(-id),5)$cluster) sondy_podil_casu %>% summarize_group_podil_prace() %>% plot_sondy_clusters(lab = c("Jine", "Roverska cinnost", "U detskych oddilu", "Ucast/pomoc celostatni a vzdelavaci"), ylim_min = 0, ylim_max = 100) sondy_podil_casu %>% summarize_group_podil_prace("skupina_k2") %>% plot_sondy_podil_casu() sondy_podil_casu %>% summarize_group_podil_prace("skupina_k3") %>% plot_sondy_podil_casu() sondy_podil_casu %>% summarize_group_podil_prace("skupina_k4") %>% plot_sondy_podil_casu() sondy_podil_casu %>% summarize_group_podil_prace("skupina_k5") %>% plot_sondy_podil_casu()
Jak se lisi korelacni matice u jednotlivych subskupin
Vezmeme podskupiny podle kmeans a replikujeme na nich korelacni matice
sondy2 <- sondy %>% left_join(sondy_podil_casu %>% select(-(podil_casu_ve_skautu_u_detskych_oddilu:podil_casu_ve_skautu_jine)), by = "id") aa <- sondy2 %>% mutate(celkova_spokojenost_s_programem_rs_kmenu = recode(celkova_spokojenost_s_programem_rs_kmenu, `Velmi nespokojen/a`= 1, `Spíše nespokojen/a` = 2, `Spíše spokojen/a` = 3, `Spíše spokojen/a` = 4, `Velmi spokojen/a` = 5, .default = NA_real_),spokojenost_s_frekvenci_akci2 = recode(spokojenost_s_frekvenci_akci, `Jsem spokojen/a`= 1, `Nejsem spokojen/a - přál/a bych si akce častěji` = 0, `Nejsem spokojen/a - přál/a bych si akce méně často` = 0, .default = NA_real_), frekvence_rs_akci_celodenni_akce = recode(frekvence_rs_akci_celodenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_), frekvence_rs_akci_schuzky = recode(frekvence_rs_akci_schuzky, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_), frekvence_rs_akci_vicedenni_akce = recode(frekvence_rs_akci_vicedenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, .default = NA_real_), frekvence_rs_akci_tabory_expedice = recode(frekvence_rs_akci_tabory_expedice, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, .default = NA_real_), vychovne_nastroje_vyzvy = recode(vychovne_nastroje_vyzvy, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_projekty = recode(vychovne_nastroje_projekty, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_rovernet = recode(vychovne_nastroje_rovernet, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_kmen = recode(vychovne_nastroje_kmen, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_roversky_zacatek = recode(vychovne_nastroje_roversky_zacatek, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_diar = recode(vychovne_nastroje_diar, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_vlastni_stezka = recode(vychovne_nastroje_vlastni_stezka, `Ano` = 1, `Ne` = 0, .default = NA_real_), ma_kmen_vychovny_plan = recode(ma_kmen_vychovny_plan, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_nic = recode(vlastni_zvyklosti_nic, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_vstupu = recode(vlastni_zvyklosti_ritual_vstupu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_obnova_slibu = recode(vlastni_zvyklosti_obnova_slibu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_odchodu = recode(vlastni_zvyklosti_ritual_odchodu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_predavani_vudcovstvi = recode(vlastni_zvyklosti_ritual_predavani_vudcovstvi, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_vudce = recode(rs_kmen_prvky_vudce, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_rada = recode(rs_kmen_prvky_rada, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_ani_jedno = recode(rs_kmen_prvky_ani_jedno, `Ano` = 1, `Ne` = 0, .default = NA_real_) )
4 skupiny - rozdily
aa2 <- aa %>% filter(!is.na(skupina_k5)) aa2$skupina_k5 <- recode(aa2$skupina_k5, `1`="detska-vzdelavaci",`2`="prevazne-detska",`3`="jina cinnost",`4`="roverska", `5`="detska") vs <- aa %>% select_at(v) %>% colnames() for (i in 1) { p <- aa2 %>% ggplot(aes_string(x = "skupina_k5", y = vs[i])) + stat_summary(fun.data = "mean_cl_boot") + ylab("Procenta") + ylim(0,1) + ggtitle(vs[i]) print(p) } for (i in 2:6) { p <- aa2 %>% ggplot(aes_string(x = "skupina_k5", y = vs[i])) + stat_summary(fun.data = "mean_cl_boot") + ylab("Procenta") + ylim(0,6) + ggtitle(vs[i]) print(p) } for (i in 7:length(vs)) { p <- aa2 %>% ggplot(aes_string(x = "skupina_k5", y = vs[i])) + stat_summary(fun.data = "mean_cl_boot") + ylab("Procenta") + ylim(0,1) + ggtitle(vs[i]) print(p) }
2 skupiny
subgroup_corr(aa,v,vars(skupina_k2), 1) %>% knitr::kable(caption="dve skupiny, skupina 1") subgroup_corr(aa,v,vars(skupina_k2), 2) %>% knitr::kable(caption="dve skupiny, skupina 2")
3 skupiny
subgroup_corr(aa,v,vars(skupina_k3), 1) %>% knitr::kable(caption="tri skupiny, skupina 1") subgroup_corr(aa,v,vars(skupina_k3), 2) %>% knitr::kable(caption="tri skupiny, skupina 2") subgroup_corr(aa,v,vars(skupina_k3), 2) %>% knitr::kable(caption="tri skupiny, skupina 3")
4 skupiny
subgroup_corr(aa,v,vars(skupina_k4), 1) %>% knitr::kable(caption="ctyri skupiny, skupina 1") subgroup_corr(aa,v,vars(skupina_k4), 2) %>% knitr::kable(caption="ctyri skupiny, skupina 2") subgroup_corr(aa,v,vars(skupina_k4), 3) %>% knitr::kable(caption="ctyri skupiny, skupina 3") subgroup_corr(aa,v,vars(skupina_k4), 4) %>% knitr::kable(caption="ctyri skupiny, skupina 4")
5 skupin
subgroup_corr(aa,v,vars(skupina_k5), 1) %>% knitr::kable(caption="pet skupin, skupina 1") subgroup_corr(aa,v,vars(skupina_k5), 2) %>% knitr::kable(caption="pet skupin, skupina 2") subgroup_corr(aa,v,vars(skupina_k5), 3) %>% knitr::kable(caption="pet skupin, skupina 3") subgroup_corr(aa,v,vars(skupina_k5), 4) %>% knitr::kable(caption="pet skupin, skupina 4") subgroup_corr(aa,v,vars(skupina_k5), 5) %>% knitr::kable(caption="pet skupin, skupina 5")
Clusterováni na základě typu proběhlých programů
cluster_var <- vars(roverske_akce_probehle_2roky_zavody:roverske_akce_kurzy_vzdelavaci_akce) sondy_program_tmp <- sondy %>% select_at(cluster_var) %>% mutate_all(funs(recode(.,"Ano"=T,"Ne"=F,"N/A"=NA))) %>% mutate(id = sondy$id) %>% filter(!is.na(roverske_akce_probehle_2roky_zavody)) sondy_program <- sondy_program_tmp %>% mutate(skupina_k2 = kmeans(sondy_program_tmp %>% select(-id),2)$cluster, skupina_k3 = kmeans(sondy_program_tmp %>% select(-id),3)$cluster, skupina_k4 = kmeans(sondy_program_tmp %>% select(-id),4)$cluster, skupina_k5 = kmeans(sondy_program_tmp %>% select(-id),5)$cluster, skupina_k6 = kmeans(sondy_program_tmp %>% select(-id),6)$cluster, skupina_k7 = kmeans(sondy_program_tmp %>% select(-id),7)$cluster, skupina_k8 = kmeans(sondy_program_tmp %>% select(-id),8)$cluster, skupina_k9 = kmeans(sondy_program_tmp %>% select(-id),9)$cluster, skupina_k10 = kmeans(sondy_program_tmp %>% select(-id),10)$cluster) l <- c("zavody", "expedice", "tabor", "programy", "pobyt v prirode", "zahranicni akce", "kulturni akce", "diskuze", "tradicni akce", "setkani bez programu/hospoda", "sluzba", "velke rs akce", "vzdelavaci akce") sondy_program %>% summarize_group_typ_programu() %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k2") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k3") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k4") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k5") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k6") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k7") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k8") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k9") %>% plot_sondy_clusters(l) sondy_program %>% summarize_group_typ_programu("skupina_k10") %>% plot_sondy_clusters(l)
Korelace per akcni kmeny
ab <- sondy %>% left_join(sondy_program %>% select(id, skupina_k2:skupina_k10), by = "id") %>% mutate(celkova_spokojenost_s_programem_rs_kmenu = recode(celkova_spokojenost_s_programem_rs_kmenu, `Velmi nespokojen/a`= 1, `Spíše nespokojen/a` = 2, `Spíše spokojen/a` = 3, `Spíše spokojen/a` = 4, `Velmi spokojen/a` = 5, .default = NA_real_), spokojenost_s_frekvenci_akci2 = recode(spokojenost_s_frekvenci_akci, `Jsem spokojen/a`= 1, `Nejsem spokojen/a - přál/a bych si akce častěji` = 0, `Nejsem spokojen/a - přál/a bych si akce méně často` = 0, .default = NA_real_), frekvence_rs_akci_celodenni_akce = recode(frekvence_rs_akci_celodenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_), frekvence_rs_akci_schuzky = recode(frekvence_rs_akci_schuzky, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, `Každý týden nebo častěji` = 6, .default = NA_real_), frekvence_rs_akci_vicedenni_akce = recode(frekvence_rs_akci_vicedenni_akce, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, `Asi jednou měsíčně` = 4, `Asi jednou za 14 dní` = 5, .default = NA_real_), frekvence_rs_akci_tabory_expedice = recode(frekvence_rs_akci_tabory_expedice, `Nikdy` = 0, `Méně často než 1 za rok` = 1, `Asi jednou za půl roku` = 2, `Asi jednou za čtvrt roku` = 3, .default = NA_real_), vychovne_nastroje_vyzvy = recode(vychovne_nastroje_vyzvy, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_projekty = recode(vychovne_nastroje_projekty, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_rovernet = recode(vychovne_nastroje_rovernet, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_kmen = recode(vychovne_nastroje_kmen, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_roversky_zacatek = recode(vychovne_nastroje_roversky_zacatek, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_diar = recode(vychovne_nastroje_diar, `Ano` = 1, `Ne` = 0, .default = NA_real_), vychovne_nastroje_vlastni_stezka = recode(vychovne_nastroje_vlastni_stezka, `Ano` = 1, `Ne` = 0, .default = NA_real_), ma_kmen_vychovny_plan = recode(ma_kmen_vychovny_plan, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_nic = recode(vlastni_zvyklosti_nic, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_vstupu = recode(vlastni_zvyklosti_ritual_vstupu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_obnova_slibu = recode(vlastni_zvyklosti_obnova_slibu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_odchodu = recode(vlastni_zvyklosti_ritual_odchodu, `Ano` = 1, `Ne` = 0, .default = NA_real_), vlastni_zvyklosti_ritual_predavani_vudcovstvi = recode(vlastni_zvyklosti_ritual_predavani_vudcovstvi, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_vudce = recode(rs_kmen_prvky_vudce, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_rada = recode(rs_kmen_prvky_rada, `Ano` = 1, `Ne` = 0, .default = NA_real_), rs_kmen_prvky_ani_jedno = recode(rs_kmen_prvky_ani_jedno, `Ano` = 1, `Ne` = 0, .default = NA_real_) ) v <- vars(spokojenost_s_frekvenci_akci2, celkova_spokojenost_s_programem_rs_kmenu, frekvence_rs_akci_celodenni_akce, frekvence_rs_akci_schuzky, frekvence_rs_akci_vicedenni_akce, frekvence_rs_akci_tabory_expedice,vychovne_nastroje_vyzvy:vychovne_nastroje_vlastni_stezka, vlastni_zvyklosti_nic:vlastni_zvyklosti_ritual_predavani_vudcovstvi,rs_kmen_prvky_vudce:rs_kmen_prvky_ani_jedno) subgroup_corr(ab,v,vars(skupina_k2), 1) %>% knitr::kable(caption="dve skupiny, skupina 1") subgroup_corr(ab,v,vars(skupina_k2), 2) %>% knitr::kable(caption="dve skupiny, skupina 2") subgroup_corr(ab, v, vars(skupina_k3), 1) %>% knitr::kable(caption="tri skupiny, skupina 1") subgroup_corr(ab, v, vars(skupina_k3), 2) %>% knitr::kable(caption="tri skupiny, skupina 2") subgroup_corr(ab, v, vars(skupina_k3), 3) %>% knitr::kable(caption="tri skupiny, skupina 3") subgroup_corr(ab, v, vars(skupina_k4), 1) %>% knitr::kable(caption="ctyri skupiny, skupina 1") subgroup_corr(ab, v, vars(skupina_k4), 2) %>% knitr::kable(caption="ctyri skupiny, skupina 2") subgroup_corr(ab, v, vars(skupina_k4), 3) %>% knitr::kable(caption="ctyri skupiny, skupina 3") subgroup_corr(ab, v, vars(skupina_k4), 4) %>% knitr::kable(caption="ctyri skupiny, skupina 4")
Korelace mezi podilem cinnosti a typem programu
podilcasu_program <- sondy_podil_casu %>% rename(podilcasu2 = skupina_k2,podilcasu3 = skupina_k3, podilcasu4 = skupina_k4,podilcasu5 = skupina_k5) %>% left_join(sondy_program, by = "id") podilcasu_program %>% select(podilcasu2:podilcasu5,skupina_k2:skupina_k10) %>% as.matrix() %>% Hmisc::rcorr() %>% format_corrtable() %>% knitr::kable() #xtabs(~podilcasu5+skupina_k2, podilcasu_program) %>% knitr::kable() #xtabs(~podilcasu3+skupina_k2, podilcasu_program)%>% knitr::kable()
Nvychazi tam zadna zajimava souvislost
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.