| avm99963 | 50871ca | 2021-02-22 19:18:49 +0100 | [diff] [blame^] | 1 | # Classe 16/02/2021
|
| 2 | # Estadística, FME, presencial, curs 2020-21, 2n semestre
|
| 3 |
|
| 4 | library(tables)
|
| 5 | library(RcmdrMisc)
|
| 6 |
|
| 7 | # Conversió de dades
|
| 8 | dd = read.csv2("Enquesta.csv")
|
| 9 | dd$MO = factor(dd$MO, labels=c("No", "Sí"))
|
| 10 | dd$Envas = factor(dd$Envas, labels=c("No", "Sí"))
|
| 11 | dd$Vidre = factor(dd$Vidre, labels=c("No", "Sí"))
|
| 12 | dd$Sexe = factor(dd$Sexe)
|
| 13 | dd$Aspecte = factor(dd$Aspecte)
|
| 14 | dd$IM = dd$D2/dd$D4
|
| 15 |
|
| 16 |
|
| 17 | # Apartats 1 i 2
|
| 18 | dd$BMI = dd$Pes/(dd$H/100)^2
|
| 19 | dd$CBMI = cut(dd$BMI, c(0, 19, 27, 100), labels=c("1-baix", "2-mig", "3-alt"))
|
| 20 | head(dd)
|
| 21 |
|
| 22 | # Apartat 3
|
| 23 | summary(dd$Edat)
|
| 24 | mean(dd$Edat, trim=.1)
|
| 25 | boxplot(dd$Edat)
|
| 26 |
|
| 27 | boxplot(dd$Pes)
|
| 28 | boxplot(dd$H)
|
| 29 | boxplot(dd$D2)
|
| 30 |
|
| 31 | tsexe = table(dd$Sexe)
|
| 32 | tsexe
|
| 33 | barplot(tsexe)
|
| 34 |
|
| 35 | tmo = table(dd$MO)
|
| 36 | tmo
|
| 37 | barplot(tmo)
|
| 38 |
|
| 39 | tvidre = table(dd$Vidre)
|
| 40 | tvidre
|
| 41 | barplot(tvidre)
|
| 42 |
|
| 43 | tenvas = table(dd$Envas)
|
| 44 | tenvas
|
| 45 | barplot(tenvas)
|
| 46 |
|
| 47 | tpc = table(dd$PC)
|
| 48 | tpc
|
| 49 | barplot(tpc)
|
| 50 |
|
| 51 | taspecte = table(dd$Aspecte)
|
| 52 | taspecte
|
| 53 | barplot(taspecte)
|
| 54 |
|
| 55 | # Apartat 4
|
| 56 | boxplot(dd$BMI)
|
| 57 | abline(h=c(19, 27))
|
| 58 |
|
| 59 | # Apartat 5
|
| 60 | hist(dd$Pes, freq=F, ylim=c(0, 0.04), xlim=c(30, 130))
|
| 61 | curve(dnorm(x, mean=mean(dd$Pes), sd=sd(dd$Pes)), add=T)
|
| 62 |
|
| 63 | boxplot(dd$IM)
|
| 64 |
|
| 65 | # == BLOC 2 - ANÀLISI DESCRIPTIVA ==
|
| 66 | # Apartat 1
|
| 67 | # Tabular:
|
| 68 | # ~: Separa files i columnes
|
| 69 | # +: Concatena
|
| 70 | # *: Nesting
|
| 71 | tabular(mean*(Sexe+1)~H*(CBMI+1), dd)
|
| 72 | tabular(var*Sexe~H*CBMI, dd)
|
| 73 | interaction.plot(dd$CBMI, dd$Sexe, dd$H)
|
| 74 | plotMeans(dd$H, dd$CBMI, dd$Sexe, error.bars="sd")
|