A partir del conjunt de dades anorexia del paquet MASS, creeu un altre data frame que s’anomeni anorexia_treat_df format per Treat i per un vector nou calculat a partir de la diferència Prewt-Postwt. D’aquesta manera, ens quedarà un data frame que contingui el tipus de tractament i el valor del pes guanyat o perdut després d’haver fet el tractament.
library(MASS)
data("anorexia")
head(anorexia)## Treat Prewt Postwt
## 1 Cont 80.7 80.2
## 2 Cont 89.4 80.1
## 3 Cont 91.8 86.4
## 4 Cont 74.0 86.3
## 5 Cont 78.1 76.1
## 6 Cont 88.3 78.1anorexia_tret_df <-data.frame(anorexia$Treat,c(anorexia$Prewt-
anorexia$Postwt))
head(anorexia_tret_df)## anorexia.Treat c.anorexia.Prewt...anorexia.Postwt.
## 1 Cont 0.5
## 2 Cont 9.3
## 3 Cont 5.4
## 4 Cont -12.3
## 5 Cont 2.0
## 6 Cont 10.2Seleccioneu aquells individus que han guanyat pes després del
tractament i creeu un nou conjunt anomenat anorexia_tret_C_df que
contingui només les dades d’aquells que han seguit el tractament
<
colnames(anorexia_tret_df)## [1] "anorexia.Treat" "c.anorexia.Prewt...anorexia.Postwt."colnames(anorexia_tret_df)[colnames(anorexia_tret_df) == "c.anorexia.Prewt...anorexia.Postwt."] <- "Dif"
colnames(anorexia_tret_df)## [1] "anorexia.Treat" "Dif"anorexia_tret_C_df <- subset(anorexia_tret_df, Dif<0 & anorexia.Treat=="Cont")
head(anorexia_tret_C_df)## anorexia.Treat Dif
## 4 Cont -12.3
## 8 Cont -11.6
## 10 Cont -6.2
## 13 Cont -8.3
## 14 Cont -3.3
## 15 Cont -11.3