R-tools

> wiki_username <- "Jouni"
> library(OpasnetUtils)
> library(reshape2)
> library(ggplot2)
> orig.data <- opbase.data("Op_fi5697", subset="kyselydata")
> kyselydata <- as.data.frame(lapply(orig.data, FUN=function(x) {x[x == "En osaa sanoa"] <- NA; x})) # En osaa sanoa -> NA
> #kyselydata <- as.data.frame(lapply(kyselydata, FUN=as.numeric))
> #oprint(kyselydata)
> melted.kysely <- melt(kyselydata,
+ id=c("Sukupuoli", "Opiskeluala", "Opiskeluvuosi"),
+ variable.name = "argumentti",
+ value.name="vastaus"
+ )
> melted.kysely$hypoteesi <- ifelse(melted.kysely$argumentti %in% c("rystykavely_energiatehotonta",
+ "pystykavely_energiatehokasta",
+ "pienet_muutokset",
+ "v_simpanssi_pysty_tehoton",
+ "v_ei_energiatehokasta",
+ "v_bent_hip"), "energiatehokkuus",
+ ifelse(melted.kysely$argumentti %in% c("orangin_liike",
+ "jalan_ojennus",
+ "lapaluu",
+ "pystysuora_kiipeaminen",
+ "tukipisteet",
+ "peukalovarvas",
+ "ranteen_rakenteet",
+ "v_rysty_kantamuoto",
+ "v_orangit_ei_maassa",
+ "v_simpanssi_gorilla_merkki"), "puissakavely",
+ ifelse(melted.kysely$argumentti %in% c("apinat_vedessa",
+ "muut_vesiapinat",
+ "ranta_ruoka",
+ "noste",
+ "virtaviivaisuus",
+ "v_kadelliset_ei_veteen",
+ "v_yksinkertaisuus"), "kahlaus",
+ ifelse(melted.kysely$argumentti %in% c("iho_aurinko_tuuli",
+ "pidempaan_savannilla",
+ "v_savannielaimet_nelijalkaisia",
+ "v_helpommin_viileaksi",
+ "v_ei_savannilla"), "lammonsaately",
+ ifelse(melted.kysely$argumentti %in% c("simpanssi_kantaa",
+ "naaraat_ruokaa",
+ "vauvan_kanto",
+ "v_epatasapaino_vauva",
+ "v_tyokalujen_ajoitus"), "kadet",
+ ifelse(melted.kysely$argumentti %in% c("apinat_keraavat_seisten",
+ "pystyasento_rakenteet",
+ "pysty_korkeammalle",
+ "lyhyt_matka",
+ "v_miksi_erikoistua"), "ruokailuasento", "hypoteesi"
+ ))))))
> #oprint(melted.kysely)
> oprint(head(melted.kysely))

> ##HISTOGRAMMEJA
> for(i in unique(melted.kysely$hypoteesi)){
+ kuvat <- ggplot(melted.kysely[melted.kysely$hypoteesi == i,], aes(vastaus)) +
+ geom_histogram(fill="red") +
+ labs(x="", y="määrä", title= paste(i)) +
+ theme_grey(base_size=24) +
+ #theme(axis.text.x=element_text(angle=90,hjust=1,vjust=0.4)) +
+ facet_wrap(~argumentti)
+
+ print(kuvat)
+
+ }