Parcial02

library(MASS)

## Warning: package 'MASS' was built under R version 4.0.5

data(epil)
library(ggplot2)

## Warning: package 'ggplot2' was built under R version 4.0.5

#Punto 1: Numero total de convulsiones por paciente
vector = c()
for (i in 1:59) {
  vector[i] = sum(epil$y[epil$subject==i])
}
vector

##  [1]  14  14  11  13  55  22  12  93  22  33  66  30  16  42  59  16   6 123  15
## [20]  16  14  14  13  30 143   6  10  53  42  28   7  13  19  11  74  20  10  24
## [39]  29   4   6  12  65  26  39   7  32   3 302  13  26  10  70  13  15  51   6
## [58]   0  10

#Punto 2: Relación #Convulsiones vs Convulsiones de referencia
cor(epil$y,epil$base)

## [1] 0.7620134

ggplot(data= epil,aes(epil$y,epil$base))+geom_point()

## Warning: Use of `epil$y` is discouraged. Use `y` instead.

## Warning: Use of `epil$base` is discouraged. Use `base` instead.

#Punto 3: Pacientes con menos de 3 convulsiones en el último periodo
print("Pacientes con menos de 3 convulsiones y en el último periodo: ")

## [1] "Pacientes con menos de 3 convulsiones y en el último periodo: "

pacientes= nrow(epil[epil$y<3 & epil$period=='4',])
pacientes

## [1] 13

#punto 4: Grafique en 4 rangos de edad el número de convulsiones promedio
promedio1 = mean(epil$y[epil$age<=23])
promedio1

## [1] 10.85938

promedio2 = mean(epil$y[epil$age >=24 & epil$age<=29])
promedio2

## [1] 5.367647

promedio3 = mean(epil$y[epil$age >=30 & epil$age<=35])
promedio3

## [1] 9.546875

promedio4 = mean(epil$y[epil$age >=36 & epil$age<=42])
promedio4

## [1] 6.925

promedios=c(promedio1,promedio2,promedio3,promedio4)
colores = c("dimgrey","coral","pink","skyblue")
pie(c(promedios),labels = c(promedios),main = "Promedio de convulsiones por edades", col = colores)
legend("topright",c(paste(promedio1,": de 18 a 23 años"),paste(promedio2,": de 24 a 29 años"),paste(promedio3,": de 30 a 35 años"),paste(promedio3,": de 36 a 42 años")),fill = colores, cex = 0.5)