A partir del conjunto de datos anorexia del paquete MASS:
Crea un data frame que se llame anorexia_treat_df formado por Treat y por un vector nuevo calculado a partir de la diferencia Prewt-Postwt.
Selecciona aquellos individuos que han ganado peso después del tratamiento y crea un nuevo conjunto llamado anorexia_treat_C_df que contenga solo los datos de aquellos que han seguido el tratamiento Cont y que han ganado peso después del tratamiento.
library(MASS)
anorexia_treat<- subset(anorexia, select= -c(Postwt, Prewt))
difwt<-c(anorexia$Postwt-anorexia$Prewt)
anorexia_treat_df<-data.frame( anorexia_treat, difwt)
names(anorexia_treat_df)
## [1] "Treat" "difwt"
anorexia_treat_C_d<-subset(anorexia_treat_df, difwt > 0 )
anorexia_treat_C_df<-subset(anorexia_treat_C_d, anorexia_treat_C_d$Treat=="Cont")
anorexia_treat_C_df
## Treat difwt
## 4 Cont 12.3
## 8 Cont 11.6
## 10 Cont 6.2
## 13 Cont 8.3
## 14 Cont 3.3
## 15 Cont 11.3
## 21 Cont 2.8
## 22 Cont 0.3
## 23 Cont 1.8
## 24 Cont 3.7
## 25 Cont 15.9
boxplot(difwt~Treat, anorexia_treat_C_d, main= "Pesos después del tratamiento")