| # Classe 16/02/2021 | |
| # Estadística, FME, presencial, curs 2020-21, 2n semestre | |
| library(tables) | |
| library(RcmdrMisc) | |
| # Conversió de dades | |
| dd = read.csv2("Enquesta.csv") | |
| dd$MO = factor(dd$MO, labels=c("No", "Sí")) | |
| dd$Envas = factor(dd$Envas, labels=c("No", "Sí")) | |
| dd$Vidre = factor(dd$Vidre, labels=c("No", "Sí")) | |
| dd$Sexe = factor(dd$Sexe) | |
| dd$Aspecte = factor(dd$Aspecte) | |
| dd$IM = dd$D2/dd$D4 | |
| # Apartats 1 i 2 | |
| dd$BMI = dd$Pes/(dd$H/100)^2 | |
| dd$CBMI = cut(dd$BMI, c(0, 19, 27, 100), labels=c("1-baix", "2-mig", "3-alt")) | |
| head(dd) | |
| # Apartat 3 | |
| summary(dd$Edat) | |
| mean(dd$Edat, trim=.1) | |
| boxplot(dd$Edat) | |
| boxplot(dd$Pes) | |
| boxplot(dd$H) | |
| boxplot(dd$D2) | |
| tsexe = table(dd$Sexe) | |
| tsexe | |
| barplot(tsexe) | |
| tmo = table(dd$MO) | |
| tmo | |
| barplot(tmo) | |
| tvidre = table(dd$Vidre) | |
| tvidre | |
| barplot(tvidre) | |
| tenvas = table(dd$Envas) | |
| tenvas | |
| barplot(tenvas) | |
| tpc = table(dd$PC) | |
| tpc | |
| barplot(tpc) | |
| taspecte = table(dd$Aspecte) | |
| taspecte | |
| barplot(taspecte) | |
| # Apartat 4 | |
| boxplot(dd$BMI) | |
| abline(h=c(19, 27)) | |
| # Apartat 5 | |
| hist(dd$Pes, freq=F, ylim=c(0, 0.04), xlim=c(30, 130)) | |
| curve(dnorm(x, mean=mean(dd$Pes), sd=sd(dd$Pes)), add=T) | |
| boxplot(dd$IM) | |
| # == BLOC 2 - ANÀLISI DESCRIPTIVA == | |
| # Apartat 1 | |
| # Tabular: | |
| # ~: Separa files i columnes | |
| # +: Concatena | |
| # *: Nesting | |
| tabular(mean*(Sexe+1)~H*(CBMI+1), dd) | |
| tabular(var*Sexe~H*CBMI, dd) | |
| interaction.plot(dd$CBMI, dd$Sexe, dd$H) | |
| plotMeans(dd$H, dd$CBMI, dd$Sexe, error.bars="sd") |