Temperatura máxima

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) Temperatura máxima
TempMax<-datos$Max.Temperature

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

# DETERMINACIÓN DE INTERVALOS CON R

histo_TempMax<-hist(TempMax,main= "Grafica No. 1.1:Distribución de la cantidad de temperatura máxima \nde cada uno de los registros de clima en el volcán Antisana",xlab= "°C", ylab= "Cantidad", col="red2", las=2)

TABLA DE FRECUENCIAS

# CREAR LOS LIMITES
Limites<-histo_TempMax$breaks
Limites
##  [1] 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
LimiteSup<-Limites[2:15]
LimiteSup
##  [1] 11 12 13 14 15 16 17 18 19 20 21 22 23 24
LimiteInf<-Limites[1:14]
LimiteInf
##  [1] 10 11 12 13 14 15 16 17 18 19 20 21 22 23
# MARCA DE CLASE
MC<-histo_TempMax$mids
MC
##  [1] 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5 18.5 19.5 20.5 21.5 22.5 23.5
#Tabla
ni<-histo_TempMax$counts
ni
##  [1]  7 23 41 53 37 40 41 42 28 22 16  9  4  3
sum(ni)
## [1] 366
hi<-ni/sum(ni)*100
hi
##  [1]  1.9125683  6.2841530 11.2021858 14.4808743 10.1092896 10.9289617
##  [7] 11.2021858 11.4754098  7.6502732  6.0109290  4.3715847  2.4590164
## [13]  1.0928962  0.8196721
sum(hi)
## [1] 100
Ni_asc<-cumsum(ni)
Ni_asc
##  [1]   7  30  71 124 161 201 242 284 312 334 350 359 363 366
Hi_asc<- cumsum(hi)
Hi_asc
##  [1]   1.912568   8.196721  19.398907  33.879781  43.989071  54.918033
##  [7]  66.120219  77.595628  85.245902  91.256831  95.628415  98.087432
## [13]  99.180328 100.000000
Ni_dsc<- rev(cumsum(rev(ni)))  
Ni_dsc  
##  [1] 366 359 336 295 242 205 165 124  82  54  32  16   7   3
Hi_dsc<- rev(cumsum(rev(hi)))  
Hi_dsc 
##  [1] 100.0000000  98.0874317  91.8032787  80.6010929  66.1202186  56.0109290
##  [7]  45.0819672  33.8797814  22.4043716  14.7540984   8.7431694   4.3715847
## [13]   1.9125683   0.8196721
Tabla_TemMax <- data.frame(LimiteInf,LimiteSup,MC,ni,round(hi,2),Ni_asc,round(Hi_asc,2),
                               Ni_dsc,round(Hi_dsc,2))  
colnames(Tabla_TemMax)<- c("LimiteInf","LimiteSup","MC","ni","hi (%)","Ni_asc","Hi_asc(%)","Ni_dsc","Hi_dsc(%)")
library(knitr)
kable(Tabla_TemMax, format = "markdown", caption = "Tabla 1.1:Tabla de frecuencias agrupadas de la temperatura máxima (°C) de los registros de clima")
Tabla 1.1:Tabla de frecuencias agrupadas de la temperatura máxima (°C) de los registros de clima
LimiteInf LimiteSup MC ni hi (%) Ni_asc Hi_asc(%) Ni_dsc Hi_dsc(%)
10 11 10.5 7 1.91 7 1.91 366 100.00
11 12 11.5 23 6.28 30 8.20 359 98.09
12 13 12.5 41 11.20 71 19.40 336 91.80
13 14 13.5 53 14.48 124 33.88 295 80.60
14 15 14.5 37 10.11 161 43.99 242 66.12
15 16 15.5 40 10.93 201 54.92 205 56.01
16 17 16.5 41 11.20 242 66.12 165 45.08
17 18 17.5 42 11.48 284 77.60 124 33.88
18 19 18.5 28 7.65 312 85.25 82 22.40
19 20 19.5 22 6.01 334 91.26 54 14.75
20 21 20.5 16 4.37 350 95.63 32 8.74
21 22 21.5 9 2.46 359 98.09 16 4.37
22 23 22.5 4 1.09 363 99.18 7 1.91
23 24 23.5 3 0.82 366 100.00 3 0.82

GRÁFICAS

#Gráfico No. 2
# DIAGRAMA DE BARRAS LOCAL
histograma_TemMAx<-hist(TempMax,main= "Grafica No.1.2: Distribución de la cantidad de temperatura máxima \nde cada uno de los registros de clima en el volcán Antisana",
                        xlab= "°C", ylab= "Cantidad", col="red", las=2)

#Gráfico No. 3
# DIAGRAMA DE BARRAS LOCAL
histograma_TemMAx<-hist(TempMax,main= "Grafica No.1.3: Distribución de la cantidad de temperatura máxima \nde cada uno de los registros de clima en el volcán Antisana",
                        xlab= "°C", ylab= "Cantidad", col="red", las=2,ylim = c(0,length(TempMax)))

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

#Gráfico No. 5
# DIAGRMA DE CAJA Y BIGOTE
Caja<-boxplot(TempMax,horizontal = T, col = "Salmon",main="Gráfica No. 1.5:Distribución de la cantidad de la temperatura máxima \nde cada uno de los registros de clima en el volcán Antisana",
              xlab= "°C" )

# OJIVAS

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

lines(LimiteSup,Ni_asc,col="blue",type = "b")
legend("right",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. 1.6: Ojivas combinadas de la temperatura máxima \nde cada uno de los registros de clima en el volcán Antisana",
     xlab = "°C", ylab="Cantidad (%)", col="black",type = "b")

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

INDICADORES

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

mediana<-median(TempMax)
mediana
## [1] 15.5105
media_aritmetica<-mean(TempMax)
media_aritmetica
## [1] 15.74001
Mo<-c("[13,14]")
desviación_estandar<-sd(TempMax)
desviación_estandar
## [1] 2.867474
coeficiente_variabilidad <- (desviación_estandar/media_aritmetica) * 100 
coeficiente_variabilidad
## [1] 18.21775
As<-skewness(TempMax)
As
## [1] 0.387115
curtosis<-kurtosis(TempMax)
curtosis
## [1] -0.5588095
#TABLAS
Variable<-c("Temperatura maxima (°C)")
Tabla_indicadores<-data.frame(Variable,ri,rs,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. 1.2: Indicadores estadísticos de la variable Temperatura máxima de cada uno de los registros de clima en el volcán Antisana.")
Tabla No. 1.2: Indicadores estadísticos de la variable Temperatura máxima 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
Temperatura maxima (°C) 10.322 23.788 15.74 15.5105 [13,14] 2.87 18.22 0.39 -0.56

CONCLUSIONES: CONCLUSIONES: La temperatura máxima en el volcán Antisana fluctúa entre 10.32°C y 23.71 °C y sus valores están en torno a los 15.74 °C, con una desviación estándar de 2.87 siendo un conjunto de valores poco heterogéneo cuyos valores se concentran en la parte media de la variable y con un ligero sesgo hacia la derecha, por lo tanto el comportamiento de la variable es beneficiosa ya que la mayoría de tiempo la temperatura máxima se encuentra a menos de 18 °C, lo que favorece a la flora y fauna del lugar, al igual que la conservación de los glaciares por ende su hidrología.

.