histo_clinica = data.frame(
  
  Cedula = c(1005944792, 1007358169, 1002890765, 1006543287, 1004721983, 1008365921, 1001653478, 1009278106, 1007132546, 1005817392, 1006490273, 1003726159, 1009581624, 1004268701, 1007819352, 1002136479, 1005942861, 1008651392, 1007264815, 1003497168),

  Nombre = c("Wilson", "Estela", "Luis", "Ana", "Maria", "Carlos", "Laura", "Jose", "Isabel", "Juan", "Carmen", "Fernando", "Elena", "Andres", "Sofia", "Pedro", "Monica", "Ricardo", "Valeria", "Diego"),
  
  Apellido = c("Gonzalez", "Rodriguez", "Lopez", "Perez", "Martinez", "Sanchez", "Romero", "Fernandez", "Torres", "Gomez", "Diaz", "Vargas", "Jimenez", "Ruiz", "Ramirez", "Hernandez", "Garcia", "Silva", "Ortega", "Castro"),
  
  Edad = c(7, 32, 56, 19, 42, 68, 11, 25, 38, 73, 61, 15, 49, 30, 77, 22, 64, 9, 53, 35),
  
  Sexo =  c("Masculino", "Femenino", "Masculino", "Femenino", "Femenino", "Masculino", "Femenino", "Masculino", "Femenino", "Masculino", "Femenino", "Masculino", "Femenino", "Masculino", "Femenino", "Masculino", "Femenino", "Masculino", "Femenino", "Masculino"),
  
  Domicilio = c("Casa", "Apartamento", "Finca", "Apartamento", "Casa", "Finca", "Casa", "Finca", "Casa", "Apartamento", "Casa", "Casa", "Apartamento", "Finca", "Casa", "Casa", "Apartamento", "Casa", "Casa", "Casa"),
  
  Ciudad = c("Cali", "Cali", "Bogota", "Medellin", "Cali", "Cali", "Bogota", "Medellin", "Cali", "Neiva", "Bogota", "Cali", "Tunja", "Medellin", "Cali", "Bogota", "Cali", "Medellin", "Cali", "Medellin"),
  
  Departamento = c("Valle del Cauca", "Valle del Cauca", "Bogota, D.C.", "Antioquia", "Valle del Cauca", "Valle del Cauca", "Bogota, D.C.", "Antioquia", "Valle del Cauca", "Huila", "Bogota, D.C.", "Valle del Cauca", "Boyaca", "Antioquia", "Valle del Cauca", "Bogota, D.C.", "Valle del Cauca", "Antioquia", "Valle del Cauca", "Antioquia"),
  
  Telefono = c(3187465290, 3021984756, 3912345678, 3210987654, 3309876541, 3245678901, 3123456789, 3045678912, 3934567890, 3812345670, 3278901456, 3345678901, 3009876542, 3198765430, 3456789012, 3389012345, 3034567891, 3098765432, 3290876541, 3145678902),
  
  Num_historial = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20),
  
  Fuma = 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","Reubicacion","Domicilio","Accidente","Accidente","Domicilio","Domicilio","Accidente","Accidente","Accidente","Traslado","Traslado","Accidente","Domicilio","Reubicacion","Reubicacion","Accidente","Domicilio","Reubicacion"),

  Fecha_ingreso = 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"),
  
  Piso = c(5,4,3,2,1,1,2,3,4,5,5,2,1,4,3,3,1,5,2,4),
  
  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   Nombre  Apellido Edad      Sexo   Domicilio   Ciudad
## 1  1005944792   Wilson  Gonzalez    7 Masculino        Casa     Cali
## 2  1007358169   Estela Rodriguez   32  Femenino Apartamento     Cali
## 3  1002890765     Luis     Lopez   56 Masculino       Finca   Bogota
## 4  1006543287      Ana     Perez   19  Femenino Apartamento Medellin
## 5  1004721983    Maria  Martinez   42  Femenino        Casa     Cali
## 6  1008365921   Carlos   Sanchez   68 Masculino       Finca     Cali
## 7  1001653478    Laura    Romero   11  Femenino        Casa   Bogota
## 8  1009278106     Jose Fernandez   25 Masculino       Finca Medellin
## 9  1007132546   Isabel    Torres   38  Femenino        Casa     Cali
## 10 1005817392     Juan     Gomez   73 Masculino Apartamento    Neiva
## 11 1006490273   Carmen      Diaz   61  Femenino        Casa   Bogota
## 12 1003726159 Fernando    Vargas   15 Masculino        Casa     Cali
## 13 1009581624    Elena   Jimenez   49  Femenino Apartamento    Tunja
## 14 1004268701   Andres      Ruiz   30 Masculino       Finca Medellin
## 15 1007819352    Sofia   Ramirez   77  Femenino        Casa     Cali
## 16 1002136479    Pedro Hernandez   22 Masculino        Casa   Bogota
## 17 1005942861   Monica    Garcia   64  Femenino Apartamento     Cali
## 18 1008651392  Ricardo     Silva    9 Masculino        Casa Medellin
## 19 1007264815  Valeria    Ortega   53  Femenino        Casa     Cali
## 20 1003497168    Diego    Castro   35 Masculino        Casa Medellin
##       Departamento   Telefono Num_historial Fuma Procedencia Fecha_ingreso Piso
## 1  Valle del Cauca 3187465290             1    N    Traslado    2023-12-01    5
## 2  Valle del Cauca 3021984756             2    S   Domicilio    2023-12-03    4
## 3     Bogota, D.C. 3912345678             3    S Reubicacion    2023-12-05    3
## 4        Antioquia 3210987654             4    N   Domicilio    2023-12-10    2
## 5  Valle del Cauca 3309876541             5    S   Accidente    2023-12-12    1
## 6  Valle del Cauca 3245678901             6    S   Accidente    2023-12-15    1
## 7     Bogota, D.C. 3123456789             7    N   Domicilio    2023-12-18    2
## 8        Antioquia 3045678912             8    S   Domicilio    2023-12-20    3
## 9  Valle del Cauca 3934567890             9    S   Accidente    2023-12-21    4
## 10           Huila 3812345670            10    S   Accidente    2023-12-22    5
## 11    Bogota, D.C. 3278901456            11    S   Accidente    2023-12-23    5
## 12 Valle del Cauca 3345678901            12    N    Traslado    2023-12-24    2
## 13          Boyaca 3009876542            13    S    Traslado    2023-12-26    1
## 14       Antioquia 3198765430            14    N   Accidente    2023-12-28    4
## 15 Valle del Cauca 3456789012            15    S   Domicilio    2023-12-29    3
## 16    Bogota, D.C. 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 Valle del Cauca 3290876541            19    S   Domicilio    2023-10-20    2
## 20       Antioquia 3145678902            20    S Reubicacion    2023-09-05    4
##    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
#Grafico hombres/mujeres

HM = table(histo_clinica$Sexo)
barplot(HM,main = "Sexo del Paciente",  names.arg=c("Femenino", "Masculino"), ylab="No. Personas", xlab="Sexo", col = c("yellow", "green"))

edad_fuma = table(histo_clinica$Fuma, histo_clinica$Edad)
barplot(edad_fuma, main = "Fumante vs Edad", ylab = "No. Personas", xlab = "Edad", las = 2, col = c("brown", "pink"))
legend(x = "topright", legend = c("No fumador", "Fumador"), fill = c("brown", "pink"), title = "Tipo de paciente")

# Grafica Px por piso y sexo
sex_piso = table(histo_clinica$Sexo, histo_clinica$Piso)
barplot(sex_piso, main = "Pacientes por piso", ylab = "No. Persona", xlab = "# piso", col = c("purple", "blue"))
legend(x = "topright", legend = c("Femenino", "Masculino"), fill = c("purple", "blue"), title = "Sexo")

# Grafica Edad y causa de ingreso
edad_ingreso = table(histo_clinica$Causa_ingreso, histo_clinica$Edad)
barplot(edad_ingreso, 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")

#Grafico promedio de edad por genero
prome_man = mean(histo_clinica$Edad[histo_clinica$Sexo == "Masculino"])
prome_women = mean(histo_clinica$Edad[histo_clinica$Sexo == "Femenino"])

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

#legend(x = "topright", legend = c(prom_man, prom_women), fill = c("#CAFF70","#FF34B3"), title = "Promedio de edades")
# Grafico datos independientes de 3 pacientes
barplot(c(histo_clinica$Num_historial[histo_clinica$Nombre == "Juan"], histo_clinica$Cama[histo_clinica$Nombre == "Sofia"], histo_clinica$Edad[histo_clinica$Nombre == "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")

Px1= c(histo_clinica$Edad[histo_clinica$Cedula == "1005944792"],histo_clinica$Edad[histo_clinica$Cedula == "1007132546"],histo_clinica$Edad[histo_clinica$Cedula == "1008651392"])
Px2= c(histo_clinica$Cama[histo_clinica$Cedula == "1005944792"],histo_clinica$Cama[histo_clinica$Cedula == "1007132546"],histo_clinica$Cama[histo_clinica$Cedula == "1008651392"])
Px3= c(histo_clinica$Piso[histo_clinica$Cedula =="1005944792"],histo_clinica$Piso[histo_clinica$Cedula =="1007132546"],histo_clinica$Piso[histo_clinica$Cedula =="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("#8B4513", "#FF7F24","#CD661D"), ylab = "Valor", xlab = "ID Px", ylim =c(0,40))
legend(x = "topright", legend = c("Edad","No. de cama","No. de piso"), fill = c("#8B4513", "#FF7F24","#CD661D"), title = "Datos") 

# Lista: 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"))

import = write.table(histo_clinica, file = "Practica_03_IM_2023B.csv", sep = ",",qmethod = "double")