datosalumnos <- read.csv("https://raw.githubusercontent.com/rpizarrog/probabilidad-y-estad-stica/master/practicas%20R/unidad%204/Promediode%20alumnos%20inscritos%20en%20ISC.%20los%20que%20tienen%20promedio%20de%20todas%20las%20carreras.csv")
# datosalumnos
mediaADMINISTRACION <-mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ADMINISTRACION')])
mediaARQUITECTURA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ARQUITECTURA')])
mediaBIOQUIMICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'BIOQUIMICA')])
mediaCIVIL <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'CIVIL')])
mediaELECTRICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ELECTRICA')])
mediaELECTRONICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ELECTRONICA')])
mediaGESTIONEMPRESARIAL <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'GESTION EMPRESARIAL')])
mediaINFORMATICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'INFORMATICA')])
mediaMECANICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'MECANICA')])
mediaMECATRONICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'MECATRONICA')])
mediaQUIMICA <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'QUIMICA')])
mediaSISTEMAS <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'SISTEMAS')])
mediaTIC <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'TECNOLOGIAS DE INFORMACION Y COMUNICACIONES')])
mediaINDUSTRIAL <- mean(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'INDUSTRIAL')])
# mediaADMINISTRACION
# mediaARQUITECTURA
# mediaBIOQUIMICA
# mediaCIVIL
# mediaELECTRICA
# mediaELECTRONICA
# mediaGESTIONEMPRESARIAL
# mediaINFORMATICA
# mediaMECANICA
# mediaMECATRONICA
# mediaQUIMICA
# mediaSISTEMAS
# mediaTIC
# mediaINDUSTRIAL
desvADMINISTRACION <-sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ADMINISTRACION')])
desvARQUITECTURA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ARQUITECTURA')])
desvBIOQUIMICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'BIOQUIMICA')])
desvCIVIL <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'CIVIL')])
desvELECTRICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ELECTRICA')])
desvELECTRONICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'ELECTRONICA')])
desvGESTIONEMPRESARIAL <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'GESTION EMPRESARIAL')])
desvINFORMATICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'INFORMATICA')])
desvMECANICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'MECANICA')])
desvMECATRONICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'MECATRONICA')])
desvQUIMICA <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'QUIMICA')])
desvSISTEMAS <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'SISTEMAS')])
desvTIC <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'TECNOLOGIAS DE INFORMACION Y COMUNICACIONES')])
desvINDUSTRIAL <- sd(datosalumnos$Promedio[which(datosalumnos$CARRERA == 'INDUSTRIAL')])
# desvADMINISTRACION
# desvARQUITECTURA
# desvBIOQUIMICA
# desvCIVIL
# desvELECTRICA
# desvELECTRONICA
# desvGESTIONEMPRESARIAL
# desvINFORMATICA
# desvMECANICA
# desvMECATRONICA
# desvQUIMICA
# desvSISTEMAS
# desvTIC
# desvINDUSTRIAL
tablaPROMS <- data.frame(carrera = unique(datosalumnos$CARRERA), promedios = c(mediaADMINISTRACION, mediaARQUITECTURA, mediaBIOQUIMICA, mediaCIVIL, mediaELECTRICA, mediaELECTRONICA, mediaGESTIONEMPRESARIAL, mediaINFORMATICA, mediaMECANICA, mediaMECATRONICA, mediaQUIMICA, mediaSISTEMAS, mediaTIC, mediaINDUSTRIAL), desvstd = c(desvADMINISTRACION, desvARQUITECTURA, desvBIOQUIMICA, desvCIVIL, desvELECTRICA, desvELECTRONICA, desvGESTIONEMPRESARIAL, desvINFORMATICA, desvMECANICA, desvMECATRONICA, desvQUIMICA, desvSISTEMAS, desvTIC, desvINDUSTRIAL))
tablaPROMS
## carrera promedios desvstd
## 1 ADMINISTRACION 87.18570 4.330804
## 2 ARQUITECTURA 85.04016 3.753101
## 3 BIOQUIMICA 84.35411 4.332458
## 4 CIVIL 83.95414 4.263399
## 5 ELECTRICA 84.56220 3.536467
## 6 ELECTRONICA 83.45317 4.242269
## 7 GESTION EMPRESARIAL 86.81308 4.331956
## 8 INDUSTRIAL 84.35158 4.566525
## 9 INFORMATICA 82.43986 4.040545
## 10 MECANICA 83.39957 4.034543
## 11 MECATRONICA 84.68841 3.638157
## 12 QUIMICA 84.71179 4.704340
## 13 SISTEMAS 84.18871 4.163146
## 14 TECNOLOGIAS DE INFORMACION Y COMUNICACIONES 84.50989 4.573360
probnormal <- NULL
carreras <- NULL
for (i in 1:length(tablaPROMS$carrera)) {
probnormal[i] <- pnorm(88, tablaPROMS$promedios[i], tablaPROMS$desvstd[i], lower.tail = F) * 100
carreras[i] <- as.character(tablaPROMS$carrera[i])
# print(as.character(tablaPROMS$carrera[i]))
}
tablaprob <- data.frame(carrera = carreras, prob = probnormal)
tablaprob
## carrera prob
## 1 ADMINISTRACION 42.542849
## 2 ARQUITECTURA 21.516173
## 3 BIOQUIMICA 20.002604
## 4 CIVIL 17.131662
## 5 ELECTRICA 16.550061
## 6 ELECTRONICA 14.190678
## 7 GESTION EMPRESARIAL 39.204580
## 8 INDUSTRIAL 21.215999
## 9 INFORMATICA 8.439755
## 10 MECANICA 12.708912
## 11 MECATRONICA 18.134845
## 12 QUIMICA 24.228404
## 13 SISTEMAS 17.996932
## 14 TECNOLOGIAS DE INFORMACION Y COMUNICACIONES 22.268990