#UNIVRSIDAD CENTRAL DEL ECUADOR #Tema:Estadística Descriptiva-Variables discretas #AUTOR:Grupo 6
tle: “Estadistica Descriptiva-Variables Discretas” author: “Grupo No. 6” output: html_document: toc: true toc_depth: 2 toc_float: true number_sections: true theme: united # Puedes usar: journal highlight: tango code_folding: none fontsize: 12pt
#CARGA DE DATOS
setwd(“/cloud/project”) datos <- read.csv(“Variables_Proyectoo.csv”, header = TRUE, sep = “;” , dec = “.”) str(datos)
#VARIABLES DISCRETAS 1:
##DERRAME DE PETRÓLEO OCURRIDO POR AÑO
###Extraer la variable Year_edit_Fecha_del_derrame <- datos\(Year_edit_Fecha_del_derrame TDF_Year_edit_Fecha_del_derrame <- table(Year_edit_Fecha_del_derrame) ###Tabla de Distribución de Frecuencias Simples y Acumuladas tabla_Year_edit_Fecha_del_derrame <- as.data.frame(TDF_Year_edit_Fecha_del_derrame) hi_Year_edit_Fecha_del_derrame <- (tabla_Year_edit_Fecha_del_derrame\)Freq / sum(tabla_Year_edit_Fecha_del_derrame\(Freq)) * 100 sum(hi_Year_edit_Fecha_del_derrame) tabla_Year_edit_Fecha_del_derrame\)hi <- round(hi_Year_edit_Fecha_del_derrame, 2) Niasc_Year_edit_Fecha_del_derrame <- cumsum(tabla_Year_edit_Fecha_del_derrame\(Freq) Hiasc_Year_edit_Fecha_del_derrame <- cumsum(hi_Year_edit_Fecha_del_derrame) Nidsc_Year_edit_Fecha_del_derrame <- rev(cumsum(rev(tabla_Year_edit_Fecha_del_derrame\)Freq))) Hidsc_Year_edit_Fecha_del_derrame <- rev(cumsum(rev(hi_Year_edit_Fecha_del_derrame))) tabla_Year_edit_Fecha_del_derrame_Final <- data.frame( Year_edit_Fecha_del_derrame = tabla_Year_edit_Fecha_del_derrame\(Year_edit_Fecha_del_derrame, Frecuencia = tabla_Year_edit_Fecha_del_derrame\)Freq, Porcentaje = tabla_Year_edit_Fecha_del_derrame$hi, Niasc = Niasc_Year_edit_Fecha_del_derrame, Hiasc = round(Hiasc_Year_edit_Fecha_del_derrame, 2), Nidsc = Nidsc_Year_edit_Fecha_del_derrame, Hidsc = round(Hidsc_Year_edit_Fecha_del_derrame, 2) ) ####Ver tabla final print(tabla_Year_edit_Fecha_del_derrame_Final)
##Gráfica de barras
###Gráfica de barras versión local de derrames de petróleo por Año barplot(tabla_Year_edit_Fecha_del_derrame_Final\(Frecuencia, names.arg = tabla_Year_edit_Fecha_del_derrame_Final\)Year_edit_Fecha_del_derrame, main = “Frecuencia de Derrames de Petróleo por Año de Ocurrencia”, cex.main = 0.9, xlab = “Año”, ylab = “Cantidad de Derrames”, col = “cornsilk”)
###Gráfica de barras versión global de derrames de petróleo por Año nombres_Year_edit_Fecha_del_derrame <- as.character(tabla_Year_edit_Fecha_del_derrame_Final\(Year_edit_Fecha_del_derrame) TDFYear_edit_Fecha_del_derrame <- tabla_Year_edit_Fecha_del_derrame_Final\)Frecuencia barplot(TDFYear_edit_Fecha_del_derrame, main = “Gráfica N°2: Frecuencia de Derrames de Petróleo por Año de Ocurrencia”, xlab = “Año”, ylab = “Cantidad”, names.arg = nombres_Year_edit_Fecha_del_derrame, col = “cornsilk2”, cex.main = 0.9, cex.lab = 1, cex.axis = 0.7, cex.names = 0.7, ylim = c(0, sum(TDFYear_edit_Fecha_del_derrame)))
###Gráfica de barras versión local porcentual de derrames de petróleo por Año barplot(tabla_Year_edit_Fecha_del_derrame_Final\(Porcentaje, main = "Gráfica N°3: Porcentaje de Frecuencia de Derrames de Petróleo por Año de Ocurrencia", xlab = "Año", ylab = "Porcentaje", col = "cyan", names.arg = tabla_Year_edit_Fecha_del_derrame_Final\)Year_edit_Fecha_del_derrame, cex.main = 0.9, cex.lab = 1, cex.axis = 0.7, cex.names = 0.7, ylim = c(0, max(tabla_Year_edit_Fecha_del_derrame_Final$Porcentaje) * 1.1))
###Grafica de barras versión global porcentual de derrames de petróleo por Año barplot(tabla_Year_edit_Fecha_del_derrame_Final\(Porcentaje, main = "Gráfica N°4: Porcentaje de Frecuencia de Derrames de Petróleo por Año de Ocurrencia", xlab = "Año", ylab = "Porcentaje", col = "cyan3", cex.main = 0.9, names.arg = tabla_Year_edit_Fecha_del_derrame_Final\)Year_edit_Fecha_del_derrame, ylim = c(0, 100))
##Diagrama de cajas ###Diagrama de cajas de derrames de petróleo por Año boxplot(Year_edit_Fecha_del_derrame, horizontal = TRUE, col = “azure2”, xlab = “Año”, cex.main = 0.9, main = “Gráfica N°5: Variabilidad de Frecuencia de Derrames de Petróleo por Año de Ocurrencia”)
##Ojivas de Frecuencia Acumulada ###Ojivas de frecuencia Acumulada de derrames de petróleo por Año x_Year_edit_Fecha_del_derrame <- as.numeric(as.character(tabla_Year_edit_Fecha_del_derrame$Year_edit_Fecha_del_derrame)) y_ni_asc_Year_edit_Fecha_del_derrame <- Niasc_Year_edit_Fecha_del_derrame y_ni_dsc_Year_edit_Fecha_del_derrame <- Nidsc_Year_edit_Fecha_del_derrame
plot(x_Year_edit_Fecha_del_derrame, y_ni_asc_Year_edit_Fecha_del_derrame, type = “b”, main = “Gráfica N°6:Ojivas de Frecuencia Acumulada de Derrames de Petróleo por Año”, cex.main = 0.9, xlab = “Año”, ylab = “Frecuencia acumulada”, col = “black”) lines(x_Year_edit_Fecha_del_derrame, y_ni_dsc_Year_edit_Fecha_del_derrame, col = “blue”, type = “b”) grid() legend(“topright”, legend = c(“Ascendente”, “Descendente”), col = c(“black”, “blue”), lty = 1, cex = 0.8)
#Extraer la variable Mes_edit_Fecha_del_derrame <- datos\(Mes_edit_Fecha_del_derrame TDF_Mes_edit_Fecha_del_derrame <- table(Mes_edit_Fecha_del_derrame) #Tabla de Distribución de Frecuencias Simples y Acumuladas tabla_Mes_edit_Fecha_del_derrame <- as.data.frame(TDF_Mes_edit_Fecha_del_derrame) hi_Mes_edit_Fecha_del_derrame <- (tabla_Mes_edit_Fecha_del_derrame\)Freq / sum(tabla_Mes_edit_Fecha_del_derrame\(Freq)) * 100 sum(hi_Mes_edit_Fecha_del_derrame) tabla_Mes_edit_Fecha_del_derrame\)hi <- round(hi_Mes_edit_Fecha_del_derrame, 2) Niasc_Mes_edit_Fecha_del_derrame <- cumsum(tabla_Mes_edit_Fecha_del_derrame\(Freq) Hiasc_Mes_edit_Fecha_del_derrame <- cumsum(hi_Mes_edit_Fecha_del_derrame) Nidsc_Mes_edit_Fecha_del_derrame <- rev(cumsum(rev(tabla_Mes_edit_Fecha_del_derrame\)Freq))) Hidsc_Mes_edit_Fecha_del_derrame <- rev(cumsum(rev(hi_Mes_edit_Fecha_del_derrame))) tabla_Mes_edit_Fecha_del_derrame_Final <- data.frame( Mes_edit_Fecha_del_derrame = tabla_Mes_edit_Fecha_del_derrame\(Mes_edit_Fecha_del_derrame, Frecuencia = tabla_Mes_edit_Fecha_del_derrame\)Freq, Porcentaje = tabla_Mes_edit_Fecha_del_derrame$hi, Niasc = Niasc_Mes_edit_Fecha_del_derrame, Hiasc = round(Hiasc_Mes_edit_Fecha_del_derrame, 2), Nidsc = Nidsc_Mes_edit_Fecha_del_derrame, Hidsc = round(Hidsc_Mes_edit_Fecha_del_derrame, 2) ) #Ver tabla final print(tabla_Mes_edit_Fecha_del_derrame_Final)
barplot(tabla_Mes_edit_Fecha_del_derrame_Final\(Frecuencia, names.arg = tabla_Mes_edit_Fecha_del_derrame_Final\)Mes_edit_Fecha_del_derrame, main = “Gráfica N°1: Frecuencia de Derrames de Petróleo por Mes”, cex.main = 0.9, xlab = “Mes”, ylab = “Cantidad”, col = “cornsilk”)
nombres_Mes_edit_Fecha_del_derrame <- as.character(tabla_Mes_edit_Fecha_del_derrame_Final\(Mes_edit_Fecha_del_derrame) TDFMes_edit_Fecha_del_derrame <- tabla_Mes_edit_Fecha_del_derrame_Final\)Frecuencia
barplot(TDFMes_edit_Fecha_del_derrame, main = “Gráfica N°2: Frecuencia Total de Derrames de Petróleo por Mes”, xlab = “Mes”, ylab = “Cantidad”, names.arg = nombres_Mes_edit_Fecha_del_derrame, col = “cornsilk2”, cex.main = 0.9, cex.lab = 1, cex.axis = 0.7, cex.names = 0.7, ylim = c(0, sum(TDFMes_edit_Fecha_del_derrame)))
barplot(tabla_Mes_edit_Fecha_del_derrame_Final\(Porcentaje, main = "Gráfica N°3:Porcentaje de Derrames de Petróleo por Mes", xlab = "Mes", ylab = "Porcentaje", col = "cyan", names.arg = tabla_Mes_edit_Fecha_del_derrame_Final\)Mes_edit_Fecha_del_derrame, cex.main = 0.9, cex.lab = 1, cex.axis = 0.7, cex.names = 0.7, ylim = c(0, max(tabla_Mes_edit_Fecha_del_derrame_Final$Porcentaje) * 1.1))
barplot(tabla_Mes_edit_Fecha_del_derrame_Final\(Porcentaje, main = "Gráfica N°4: Distribución Porcentual de Derrames de Petróleo por Mes", xlab = "Mes", ylab = "Porcentaje", col = "cyan3", cex.main = 0.9, names.arg = tabla_Mes_edit_Fecha_del_derrame_Final\)Mes_edit_Fecha_del_derrame, ylim = c(0, 100))
boxplot(Mes_edit_Fecha_del_derrame, horizontal = TRUE, col = “azure2”, xlab = “Mes”, cex.main = 0.9, main = “Gráfica N°5: Variabilidad en los Derrames de Petróleo por Mes”)
x_Mes_edit_Fecha_del_derrame <- as.numeric(as.character(tabla_Mes_edit_Fecha_del_derrame$Mes_edit_Fecha_del_derrame)) y_ni_asc_Mes_edit_Fecha_del_derrame <- Niasc_Mes_edit_Fecha_del_derrame y_ni_dsc_Mes_edit_Fecha_del_derrame <- Nidsc_Mes_edit_Fecha_del_derrame
plot(x_Mes_edit_Fecha_del_derrame, y_ni_asc_Mes_edit_Fecha_del_derrame, type = “b”, main = “Gráfica N°6: Ojivas de Frecuencia Acumulada de Derrames de Petróleo por Mes”, cex.main = 0.9, xlab = “Mes”, ylab = “Frecuencia acumulada”, col = “black”) lines(x_Mes_edit_Fecha_del_derrame, y_ni_dsc_Mes_edit_Fecha_del_derrame, col = “blue”, type = “b”) grid() legend(“topright”, legend = c(“Ascendente”, “Descendente”), col = c(“black”, “blue”), lty = 1, cex = 0.8)