Control De Dieta

knitr::opts_chunk$set(echo = TRUE)



dieta1 <- subset(ChickWeight, Diet == 1) 
dieta2 <- subset(ChickWeight, Diet == 2) 
dieta3 <- subset(ChickWeight, Diet == 3) 
dieta4 <- subset(ChickWeight, Diet == 4)


factor1 <- as.factor(dieta1$Time)
factor2 <- as.factor(dieta2$Time)
factor3 <- as.factor(dieta3$Time)
factor4 <- as.factor(dieta4$Time)


promedios_dieta1 <- as.numeric(tapply(dieta1$weight, dieta1$Time, mean))
promedios_dieta2 <- as.numeric(tapply(dieta2$weight, dieta2$Time, mean))
promedios_dieta3 <- as.numeric(tapply(dieta3$weight, dieta3$Time, mean))
promedios_dieta4 <- as.numeric(tapply(dieta4$weight, dieta4$Time, mean))


df_promedios <- data.frame(
  Time = unique(dieta1$Time),
  Dieta1 = promedios_dieta1,
  Dieta2 = promedios_dieta2,
  Dieta3 = promedios_dieta3,
  Dieta4 = promedios_dieta4
)

df_promedios

##    Time    Dieta1 Dieta2 Dieta3   Dieta4
## 1     0  41.40000   40.7   40.8  41.0000
## 2     2  47.25000   49.4   50.4  51.8000
## 3     4  56.47368   59.8   62.2  64.5000
## 4     6  66.78947   75.4   77.9  83.9000
## 5     8  79.68421   91.7   98.4 105.6000
## 6    10  93.05263  108.5  117.1 126.0000
## 7    12 108.52632  131.3  144.4 151.4000
## 8    14 123.38889  141.9  164.5 161.8000
## 9    16 144.64706  164.7  197.4 182.0000
## 10   18 158.94118  187.7  233.1 202.9000
## 11   20 170.41176  205.6  258.9 233.8889
## 12   21 177.75000  214.7  270.3 238.5556

plot(df_promedios$Time, df_promedios$Dieta1, type = "b", col = "red", pch = 19,
     xlab = "Dias", ylab = "Gramos", main = "Crecimiento promedio en funcion de la dieta")


lines(df_promedios$Time, df_promedios$Dieta2, type = "b", col = "blue", pch = 19)
lines(df_promedios$Time, df_promedios$Dieta3, type = "b", col = "green", pch = 19)
lines(df_promedios$Time, df_promedios$Dieta4, type = "b", col = "purple", pch = 19)


legend("topright", legend = c("Dieta 1", "Dieta 2", "Dieta 3", "Dieta 4"),
       col = c("red", "blue", "green", "purple"), pch = 19, title = "Dietas")