Humedad relativa

Variable Cuantitativa Continua

Cargamos las librería

library(PASWR)
## Loading required package: lattice
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(readr)
library(e1071)

Carga los datos (Conjunto de datos)

setwd("/cloud/project")
read_csv("weatherdataANTISANA.csv")
## Rows: 366 Columns: 10
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (1): Date
## dbl (9): Longitude, Latitude, Elevation, Max Temperature, Min Temperature, P...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
## # A tibble: 366 × 10
##    Date       Longitude Latitude Elevation `Max Temperature` `Min Temperature`
##    <chr>          <dbl>    <dbl>     <dbl>             <dbl>             <dbl>
##  1 01/01/2012     -78.1   -0.468      4048              16.1              6.91
##  2 01/02/2012     -78.1   -0.468      4048              15.5              9.23
##  3 01/03/2012     -78.1   -0.468      4048              11.5              8.69
##  4 01/04/2012     -78.1   -0.468      4048              12.0              9.53
##  5 01/05/2012     -78.1   -0.468      4048              11.7              7.90
##  6 01/06/2012     -78.1   -0.468      4048              12.1              7.84
##  7 01/07/2012     -78.1   -0.468      4048              13.1              6.39
##  8 01/08/2012     -78.1   -0.468      4048              11.5              9.76
##  9 01/09/2012     -78.1   -0.468      4048              12.9             10.1 
## 10 01/10/2012     -78.1   -0.468      4048              13.4              8.41
## # ℹ 356 more rows
## # ℹ 4 more variables: Precipitation <dbl>, Wind <dbl>,
## #   `Relative Humidity` <dbl>, Solar <dbl>
datos<- read.csv("weatherdataANTISANA.csv", header = T, sep = ",", dec = ".")

EXTRAER LA VARIABLE CONTINUA

# 1) Humedad relativa
humedad<-datos$Relative.Humidity

Para una mejor análisis de la variable, utilizaremos los intervalos que nos proporciona R.

# DETERMINACIÓN DE INTERVALOS CON R

histo_humedad<-hist(humedad,main= "Grafica No. 5.1:Distribución de la cantidad de humedad relativa \nde cada uno de los registros de clima en el volcán Antisana",
                   xlab= "%", ylab= "Cantidad", col="green", las=2)

TABLA DE FRECUENCIAS

# CREAR LOS LIMITES
Limites<-histo_humedad$breaks
Limites
##  [1] 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00
LimiteSup<-Limites[2:10]
LimiteSup
## [1] 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00
LimiteInf<-Limites[1:9]
LimiteInf
## [1] 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95
# MARCA DE CLASE
MC<-histo_humedad$mids
MC
## [1] 0.575 0.625 0.675 0.725 0.775 0.825 0.875 0.925 0.975
#Tabla
ni<-histo_humedad$counts
ni
## [1]   2  14  18  20  18  24  35  61 174
sum(ni)
## [1] 366
hi<-ni/sum(ni)*100
hi
## [1]  0.5464481  3.8251366  4.9180328  5.4644809  4.9180328  6.5573770  9.5628415
## [8] 16.6666667 47.5409836
sum(hi)
## [1] 100
Ni_asc<-cumsum(ni)
Ni_asc
## [1]   2  16  34  54  72  96 131 192 366
Hi_asc<- cumsum(hi)
Hi_asc
## [1]   0.5464481   4.3715847   9.2896175  14.7540984  19.6721311  26.2295082
## [7]  35.7923497  52.4590164 100.0000000
Ni_dsc<- rev(cumsum(rev(ni)))  
Ni_dsc  
## [1] 366 364 350 332 312 294 270 235 174
Hi_dsc<- rev(cumsum(rev(hi)))  
Hi_dsc 
## [1] 100.00000  99.45355  95.62842  90.71038  85.24590  80.32787  73.77049
## [8]  64.20765  47.54098
Tabla_Humedad <- data.frame(LimiteInf,LimiteSup,MC,ni,round(hi,2),Ni_asc,round(Hi_asc,2),
                           Ni_dsc,round(Hi_dsc,2))  
colnames(Tabla_Humedad)<- c("LimiteInf","LimiteSup","MC","ni","hi (%)","Ni_asc","Hi_asc(%)","Ni_dsc","Hi_dsc(%)")
Tabla_Humedad
##   LimiteInf LimiteSup    MC  ni hi (%) Ni_asc Hi_asc(%) Ni_dsc Hi_dsc(%)
## 1      0.55      0.60 0.575   2   0.55      2      0.55    366    100.00
## 2      0.60      0.65 0.625  14   3.83     16      4.37    364     99.45
## 3      0.65      0.70 0.675  18   4.92     34      9.29    350     95.63
## 4      0.70      0.75 0.725  20   5.46     54     14.75    332     90.71
## 5      0.75      0.80 0.775  18   4.92     72     19.67    312     85.25
## 6      0.80      0.85 0.825  24   6.56     96     26.23    294     80.33
## 7      0.85      0.90 0.875  35   9.56    131     35.79    270     73.77
## 8      0.90      0.95 0.925  61  16.67    192     52.46    235     64.21
## 9      0.95      1.00 0.975 174  47.54    366    100.00    174     47.54
library(knitr)
kable(Tabla_Humedad, format = "markdown", caption = "Tabla 5.1:Tabla de frecuencias agrupadas de la humedad de los registros de clima")
Tabla 5.1:Tabla de frecuencias agrupadas de la humedad de los registros de clima
LimiteInf LimiteSup MC ni hi (%) Ni_asc Hi_asc(%) Ni_dsc Hi_dsc(%)
0.55 0.60 0.575 2 0.55 2 0.55 366 100.00
0.60 0.65 0.625 14 3.83 16 4.37 364 99.45
0.65 0.70 0.675 18 4.92 34 9.29 350 95.63
0.70 0.75 0.725 20 5.46 54 14.75 332 90.71
0.75 0.80 0.775 18 4.92 72 19.67 312 85.25
0.80 0.85 0.825 24 6.56 96 26.23 294 80.33
0.85 0.90 0.875 35 9.56 131 35.79 270 73.77
0.90 0.95 0.925 61 16.67 192 52.46 235 64.21
0.95 1.00 0.975 174 47.54 366 100.00 174 47.54

GRÁFICAS

#Gráfico No. 2
# DIAGRAMA DE BARRAS LOCAL
histograma_Humedad<-hist(humedad,main= "Grafica No.5.2: Distribución de la cantidad de humedad \nde cada uno de los registros de clima en el volcán Antisana",
                        xlab= "%", ylab= "Cantidad", col="green", las=2)

#Gráfico No. 3
# DIAGRAMA DE BARRAS LOCAL
histograma_Humedad<-hist(humedad,main= "Grafica No.5.3: Distribución de la cantidad de humedad \nde cada uno de los registros de clima en el volcán Antisana",
                        xlab= "%", ylab= "Cantidad", col="green", las=2,ylim = c(0,length(humedad)))

#Gráfico No. 4
# DIAGRAMA DE BARRAS relativo
barplot(hi,main= "Grafica No.5.4: Distribución de la cantidad relativa de la viento \nde cada uno de los registros de clima en el volcán Antisana",
        xlab= "%", ylab= "Cantidad (%)", col="green", las=2,names.arg =Tabla_Humedad$MC)

#Gráfico No. 5
# DIAGRMA DE CAJA Y BIGOTE
Caja<-boxplot(humedad,horizontal = T, col = "green",main="Gráfica No. 5.5:Distribución de la cantidad de humedad \nde cada uno de los registros de clima en el volcán Antisana",
              xlab= "%" )

# OJIVAS

#Gráfico No.6 OJIVAS COMBINADAS DE LA FRECUENCIA
plot(LimiteInf,Ni_dsc,main = "Grafica No. 5.6: Ojivas combinadas de humedad \nde cada uno de los registros de clima en el volcán Antisana",
     xlab = "%", ylab="Cantidad", col="black",type = "b")

lines(LimiteSup,Ni_asc,col="blue",type = "b")
legend("left",legend = c("Ojiva descendente", "Ojiva ascendente"),col = c("black", "blue"), pch = 1, lty = 1,cex = 0.7)

#Gráfico No.7 OJIVAS COMBINADAS DE LA FRECUENCIA RELATIVA
plot(LimiteInf,Hi_dsc,main = "Grafica No. 5.6: Ojivas combinadas de humedad \nde cada uno de los registros de clima en el volcán Antisana",
     xlab = "%", ylab="Cantidad (%)", col="black",type = "b")

lines(LimiteSup,Hi_asc,col="blue",type = "b")
legend("left",legend = c("Ojiva descendente", "Ojiva ascendente"),col = c("black", "blue"), pch = 1, lty = 1,cex = 0.7)

INDICADORES

ri<-min(humedad)
rs<-max(humedad)

mediana<-median(humedad)
mediana
## [1] 0.9418875
media_aritmetica<-mean(humedad)
media_aritmetica
## [1] 0.8949008
Mo<-c("[0.95,1]")
desviación_estandar<-sd(humedad)
desviación_estandar
## [1] 0.1085239
coeficiente_variabilidad <- (desviación_estandar/media_aritmetica) * 100 
coeficiente_variabilidad
## [1] 12.12692
As<-skewness(humedad)
As
## [1] -1.179992
curtosis<-kurtosis(humedad)
curtosis
## [1] 0.2023913
#TABLAS
Variable<-c("humedad relativa (%)")
Tabla_indicadores<-data.frame(Variable,round(ri,2),round(rs,2),round(media_aritmetica,2),mediana,Mo,round(desviación_estandar,2), 
                              round(coeficiente_variabilidad,2), round(As,2),round(curtosis,2))
colnames(Tabla_indicadores)<-c("Variable","Mínimo","Máximo","x","Me","Mo","S","Cv (%)","As","K")
library(knitr)
kable(Tabla_indicadores, format = "markdown", caption = "Tabla No. 5.2: Indicadores estadísticos de la variable humedad relativa de cada uno de los registros de clima en el volcán Antisana.")
Tabla No. 5.2: Indicadores estadísticos de la variable humedad relativa de cada uno de los registros de clima en el volcán Antisana.
Variable Mínimo Máximo x Me Mo S Cv (%) As K
humedad relativa (%) 0.56 0.99 0.89 0.9418875 [0.95,1] 0.11 12.13 -1.18 0.2

CONCLUSIONES: CONCLUSIONES: La humedad relativa en el volcán Antisana fluctúa entre 0.56 % y 0.99 % y sus valores están en torno a los 0.9 %, con una desviación estándar de 0.11 siendo un conjunto de valores homogéneos cuyos valores se concentran en la parte derecha baja de la variable y con pronunciado sesgo hacia la izquierda, por lo tanto, el comportamiento de la variable es beneficioso. La humedad alta y constante es fundamental para el equilibrio del Antisana. Ayuda a conservar glaciares, mantener ecosistemas y regular el agua.

.