histo_clinica = data.frame(
  
  cedula_p = c("1005944792", "1007358169", "1002890765", "1006543287", "1004721983",
                       "1008365921", "1001653478", "1009278106", "1007132546", "1005817392",
                       "1006490273", "1003726159", "1009581624", "1004268701", "1007819352",
                       "1002136479", "1005942861", "1008651392", "1007264815", "1003497168"),
  
  nombre_p = c("Juan", "Maria", "Pedro", "Ana", "Luis",
                           "Laura", "Carlos", "Sofia", "Diego", "Elena",
                           "Roberto", "Carolina", "Andres", "Valeria", "Martin",
                           "Lucia", "David", "Paula", "Felipe", "Catalina"),
  
  apellidos_p = c("Gomez", "Perez", "Lopez", "Martinez", "Rodriguez",
                             "Sanchez", "Gonzalez", "Ramirez", "Torres", "Flores",
                             "Diaz", "Garcia", "Reyes", "Vargas", "Alvarez",
                             "Jimenez", "Rojas", "Molina", "Ortiz", "Fernandez"),
  
  edad_p = c(7, 32, 56, 19, 42, 68, 11, 25, 38, 73, 61, 15, 49, 30, 77, 22, 64, 9, 53, 35),
  
  sexo_p_ = c("M", "F", "M", "F", "M",
           "F", "M", "F", "M", "F",
           "M", "F", "M", "F", "M",
           "F", "M", "F", "M", "F"),
  
  domicilio_p = c("Casa", "Apartamento", "Finca", "Apartamento", "Casa", "Finca", "Casa", "Finca", "Casa", "Apartamento", "Casa", "Casa", "Apartamento", "Finca", "Casa", "Casa", "Apartamento", "Casa", "Casa", "Casa"),
  
  ciudad = c("Cali", "Manizales", "Medellin", "Bogota", "Bogota", "Cali", "Bogota", "Manizales", "Medellin", "Cali", "Cali", "Bogota", "Bogota", "Bogota", "Manizales", "Medellin", "Cali", "Manizales", "Cali", "Medellin"),

  departamento = c("Antioquia", "Caldas", "Valle del Cauca", "Bogota D.C", "Valle del Cauca", "Antioquia", "Bogota D.C", "Caldas", "Antioquia", "Valle del Cauca", "Bogota D.C", "Bogota D.C", "Valle del Cauca", "Caldas", "Bogota D.C", "Caldas", "Valle del Cauca", "Antioquia", "Bogota D.C", "Antioquia"),
  
  telefono_p = c("3187465290", "3021984756", "3912345678", "3210987654", "3309876541",
                         "3245678901", "3123456789", "3045678912", "3934567890", "3812345670",
                         "3278901456", "3345678901", "3009876542", "3198765430", "3456789012",
                         "3389012345", "3034567891", "3098765432", "3290876541", "3145678902"),
  
  num_hist = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20),
  
  fumante = c("N", "S", "S", "N", "S", "S", "N", "S", "S", "S", "S", "N", "S", "N", "S", "N", "S", "N", "S", "S"),
  
  procedencia = c("Traslado","domicilio_p","Reubicacion","domicilio_p","Accidente","Accidente","domicilio_p","domicilio_p","Accidente","Accidente","Accidente","Traslado","Traslado","Accidente","domicilio_p","Reubicacion","Reubicacion","Accidente","domicilio_p","Reubicacion"),
  
  fecha_ing = c("2023-12-01", "2023-12-03", "2023-12-05", "2023-12-10", "2023-12-12","2023-12-15", "2023-12-18", "2023-12-20", "2023-12-21", "2023-12-22","2023-12-23", "2023-12-24", "2023-12-26", "2023-12-28", "2023-12-29","2023-11-05", "2023-11-10", "2023-11-15", "2023-10-20","2023-09-05"),
  
  numpiso_ = c(5,4,3,2,1,1,2,3,4,5,5,2,1,4,3,3,1,5,2,4),
  
  num_cama = c(34,23,45,36,11,25,46,58,33,31,30,28,42,56,39,24,25,18,37,20),
  
  causa_ingreso = c("Consulta", "Terapia", "Cita medica", "Cirugia", "Terapia", "Consulta", "Cita medica", "Terapia", "Cirugia", "Cita medica", "Consulta", "Cirugia", "Terapia", "Cirugia", "Cita medica", "Terapia", "Consulta", "Cita medica", "Cirugia", "Terapia"))



histo_clinica
##      cedula_p nombre_p apellidos_p edad_p sexo_p_ domicilio_p    ciudad
## 1  1005944792     Juan       Gomez      7       M        Casa      Cali
## 2  1007358169    Maria       Perez     32       F Apartamento Manizales
## 3  1002890765    Pedro       Lopez     56       M       Finca  Medellin
## 4  1006543287      Ana    Martinez     19       F Apartamento    Bogota
## 5  1004721983     Luis   Rodriguez     42       M        Casa    Bogota
## 6  1008365921    Laura     Sanchez     68       F       Finca      Cali
## 7  1001653478   Carlos    Gonzalez     11       M        Casa    Bogota
## 8  1009278106    Sofia     Ramirez     25       F       Finca Manizales
## 9  1007132546    Diego      Torres     38       M        Casa  Medellin
## 10 1005817392    Elena      Flores     73       F Apartamento      Cali
## 11 1006490273  Roberto        Diaz     61       M        Casa      Cali
## 12 1003726159 Carolina      Garcia     15       F        Casa    Bogota
## 13 1009581624   Andres       Reyes     49       M Apartamento    Bogota
## 14 1004268701  Valeria      Vargas     30       F       Finca    Bogota
## 15 1007819352   Martin     Alvarez     77       M        Casa Manizales
## 16 1002136479    Lucia     Jimenez     22       F        Casa  Medellin
## 17 1005942861    David       Rojas     64       M Apartamento      Cali
## 18 1008651392    Paula      Molina      9       F        Casa Manizales
## 19 1007264815   Felipe       Ortiz     53       M        Casa      Cali
## 20 1003497168 Catalina   Fernandez     35       F        Casa  Medellin
##       departamento telefono_p num_hist fumante procedencia  fecha_ing numpiso_
## 1        Antioquia 3187465290        1       N    Traslado 2023-12-01        5
## 2           Caldas 3021984756        2       S domicilio_p 2023-12-03        4
## 3  Valle del Cauca 3912345678        3       S Reubicacion 2023-12-05        3
## 4       Bogota D.C 3210987654        4       N domicilio_p 2023-12-10        2
## 5  Valle del Cauca 3309876541        5       S   Accidente 2023-12-12        1
## 6        Antioquia 3245678901        6       S   Accidente 2023-12-15        1
## 7       Bogota D.C 3123456789        7       N domicilio_p 2023-12-18        2
## 8           Caldas 3045678912        8       S domicilio_p 2023-12-20        3
## 9        Antioquia 3934567890        9       S   Accidente 2023-12-21        4
## 10 Valle del Cauca 3812345670       10       S   Accidente 2023-12-22        5
## 11      Bogota D.C 3278901456       11       S   Accidente 2023-12-23        5
## 12      Bogota D.C 3345678901       12       N    Traslado 2023-12-24        2
## 13 Valle del Cauca 3009876542       13       S    Traslado 2023-12-26        1
## 14          Caldas 3198765430       14       N   Accidente 2023-12-28        4
## 15      Bogota D.C 3456789012       15       S domicilio_p 2023-12-29        3
## 16          Caldas 3389012345       16       N Reubicacion 2023-11-05        3
## 17 Valle del Cauca 3034567891       17       S Reubicacion 2023-11-10        1
## 18       Antioquia 3098765432       18       N   Accidente 2023-11-15        5
## 19      Bogota D.C 3290876541       19       S domicilio_p 2023-10-20        2
## 20       Antioquia 3145678902       20       S Reubicacion 2023-09-05        4
##    num_cama causa_ingreso
## 1        34      Consulta
## 2        23       Terapia
## 3        45   Cita medica
## 4        36       Cirugia
## 5        11       Terapia
## 6        25      Consulta
## 7        46   Cita medica
## 8        58       Terapia
## 9        33       Cirugia
## 10       31   Cita medica
## 11       30      Consulta
## 12       28       Cirugia
## 13       42       Terapia
## 14       56       Cirugia
## 15       39   Cita medica
## 16       24       Terapia
## 17       25      Consulta
## 18       18   Cita medica
## 19       37       Cirugia
## 20       20       Terapia
#Pacientes mujeres y Hombres
HM = table(histo_clinica$sexo_p)
barplot(HM,main = "sexo del Paciente",  names.arg=c("Femenino", "Masculino"), ylab="No. Personas", xlab="sexo_p", col = c("pink", "blue"))

#edad y si es fumante o no
edad_fuma = table(histo_clinica$fumante, histo_clinica$edad_p)
barplot(edad_fuma, main = "Fumante vs Edad", ylab = "No. Personas", xlab = "Edad", las = 2, col = c("white", "red"))
legend(x = "topright", legend = c("No fumador", "Fumador"), fill = c("white", "red"), title = "Tipo de paciente")

#Pacientes por piso divididos por sexo_p
sex_piso = table(histo_clinica$sexo_p, histo_clinica$numpiso)
barplot(sex_piso, main = "Pacientes por piso", ylab = "#. Persona", xlab = "# piso", col = c("purple", "green"))
legend(x = "topright", legend = c("Femenino", "Masculino"), fill = c("purple", "green"), title = "sexo_p")

#Edad y causa de ingreso
edad_pingreso = table(histo_clinica$causa_ingreso, histo_clinica$edad_p)
barplot(edad_pingreso, main = "Relacion Causa de Ingresos por edad", ylab = "No. Personas", xlab = "Edad", las = 2, col = c("#528B8B","#8DEEEE","#97FFFF","#2F4F4F"))

legend(x = "right", legend = c("Consulta", "Terapia", "Cita medica", "Cirugia"), fill = c("#528B8B","#8DEEEE","#97FFFF","#2F4F4F"), title = "Causa de ingreso")

#Promedio de edad (M,H) 

prome_man = mean(histo_clinica$edad_p[histo_clinica$sexo_p_ == "Masculino"])
prome_women = mean(histo_clinica$edad_p[histo_clinica$sexo_p_ == "Femenino"])

barplot(c(mean(histo_clinica$edad_p[histo_clinica$sexo_p_ == "Masculino"]), mean(histo_clinica$edad_p[histo_clinica$sexo_p_=="Femenino"])), main = "Promedio de edad Hombres y Mujeres", names.arg=c("M", "F"), col=c("yellow","purple"), ylab="Edad promedio", xlab="Sexo", las=1, ylim=c(0,50))

barplot(c(histo_clinica$num_hist[histo_clinica$nombre_p == "Juan"], histo_clinica$num_cama[histo_clinica$nombre_p == "Sofia"], histo_clinica$edad_p[histo_clinica$nombre_p == "Pedro"]), main ="Datos independientes de 3 pacientes",names.arg=c("Juan","Sofia","Pedro"),col=c("blue","darkgreen","red"), ylab="Valor", xlab="Nombre", las=1, ylim=c(0,40))
legend(x = "topright", legend = c("No. Historia Clinica", "Ni. cama", "Edad del Px"), fill = c("blue","darkgreen","red"), title = "Datos")

#Datos de 3 pacientes por identidad
Px1= c(histo_clinica$edad_p[histo_clinica$cedula_p == "1005944792"],histo_clinica$edad_p[histo_clinica$cedula_p == "1007132546"],histo_clinica$edad_p[histo_clinica$cedula_p == "1008651392"])

Px2= c(histo_clinica$num_cama[histo_clinica$cedula_p == "1005944792"],histo_clinica$num_cama[histo_clinica$cedula_p == "1007132546"],histo_clinica$num_cama[histo_clinica$cedula_p == "1008651392"])

Px3= c(histo_clinica$numpiso_[histo_clinica$cedula_p =="1005944792"],histo_clinica$numpiso_[histo_clinica$cedula_p =="1007132546"],histo_clinica$numpiso_[histo_clinica$cedula_p =="1008651392"])

ID = matrix(c(Px1,Px2,Px3),nrow=3,byrow=3)
dimnames(ID)=list(c("Edad","No. cama","No. piso"),c("1005944792","1007132546","1008651392"))
barplot(ID,beside= TRUE, col =c("grey", "blue","green"), ylab = "Valor", xlab = "ID Px", ylim =c(0,40))
legend(x = "topright", legend = c("Edad","No. de cama","No. de piso"), fill = c("grey", "blue","green"), title = "Datos") 

 #Lista de causas de ingreso
ingreso = table(histo_clinica$causa_ingreso)
ingreso
## 
##     Cirugia Cita medica    Consulta     Terapia 
##           5           5           4           6
barplot(ingreso, ylab = "No. Personsas", xlab = "Causa de Ingreso", col = c("#C0FF3E","#FFDEAD","#AB82FF","#FFE4E1"))

# Exportar Base de Datos a excel
import = write.table(histo_clinica, file = "practica_3_IM.xlsx", sep = ",",qmethod = "double")