Temperatura mínima

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
TempMin<-datos$Min.Temperature

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

# DETERMINACIÓN DE INTERVALOS CON R

histo_TempMin<-hist(TempMin,main= "Grafica No. 2.1:Distribución de la cantidad de temperatura mínima \nde cada uno de los registros de clima en el volcán Antisana",xlab= "°C", ylab= "Cantidad", col="blue", las=2)

TABLA DE FRECUENCIAS

# CREAR LOS LIMITES
Limites<-histo_TempMin$breaks
Limites
##  [1]  2  3  4  5  6  7  8  9 10 11
LimiteSup<-Limites[2:10]
LimiteSup
## [1]  3  4  5  6  7  8  9 10 11
LimiteInf<-Limites[1:9]
LimiteInf
## [1]  2  3  4  5  6  7  8  9 10
# MARCA DE CLASE
MC<-histo_TempMin$mids
MC
## [1]  2.5  3.5  4.5  5.5  6.5  7.5  8.5  9.5 10.5
#Tabla
ni<-histo_TempMin$counts
ni
## [1]   1   2   7  16  47 110  90  66  27
sum(ni)
## [1] 366
hi<-ni/sum(ni)*100
hi
## [1]  0.2732240  0.5464481  1.9125683  4.3715847 12.8415301 30.0546448 24.5901639
## [8] 18.0327869  7.3770492
sum(hi)
## [1] 100
Ni_asc<-cumsum(ni)
Ni_asc
## [1]   1   3  10  26  73 183 273 339 366
Hi_asc<- cumsum(hi)
Hi_asc
## [1]   0.2732240   0.8196721   2.7322404   7.1038251  19.9453552  50.0000000
## [7]  74.5901639  92.6229508 100.0000000
Ni_dsc<- rev(cumsum(rev(ni)))  
Ni_dsc  
## [1] 366 365 363 356 340 293 183  93  27
Hi_dsc<- rev(cumsum(rev(hi)))  
Hi_dsc 
## [1] 100.000000  99.726776  99.180328  97.267760  92.896175  80.054645  50.000000
## [8]  25.409836   7.377049
Tabla_TemMin <- data.frame(LimiteInf,LimiteSup,MC,ni,round(hi,2),Ni_asc,round(Hi_asc,2),
                           Ni_dsc,round(Hi_dsc,2))  
colnames(Tabla_TemMin)<- c("LimiteInf","LimiteSup","MC","ni","hi (%)","Ni_asc","Hi_asc(%)","Ni_dsc","Hi_dsc(%)")
Tabla_TemMin
##   LimiteInf LimiteSup   MC  ni hi (%) Ni_asc Hi_asc(%) Ni_dsc Hi_dsc(%)
## 1         2         3  2.5   1   0.27      1      0.27    366    100.00
## 2         3         4  3.5   2   0.55      3      0.82    365     99.73
## 3         4         5  4.5   7   1.91     10      2.73    363     99.18
## 4         5         6  5.5  16   4.37     26      7.10    356     97.27
## 5         6         7  6.5  47  12.84     73     19.95    340     92.90
## 6         7         8  7.5 110  30.05    183     50.00    293     80.05
## 7         8         9  8.5  90  24.59    273     74.59    183     50.00
## 8         9        10  9.5  66  18.03    339     92.62     93     25.41
## 9        10        11 10.5  27   7.38    366    100.00     27      7.38
library(knitr)
kable(Tabla_TemMin, format = "markdown", caption = "Tabla 2.1:Tabla de frecuencias agrupadas de la temperatura mínima (°C) de los registros de clima")
Tabla 2.1:Tabla de frecuencias agrupadas de la temperatura mínima (°C) de los registros de clima
LimiteInf LimiteSup MC ni hi (%) Ni_asc Hi_asc(%) Ni_dsc Hi_dsc(%)
2 3 2.5 1 0.27 1 0.27 366 100.00
3 4 3.5 2 0.55 3 0.82 365 99.73
4 5 4.5 7 1.91 10 2.73 363 99.18
5 6 5.5 16 4.37 26 7.10 356 97.27
6 7 6.5 47 12.84 73 19.95 340 92.90
7 8 7.5 110 30.05 183 50.00 293 80.05
8 9 8.5 90 24.59 273 74.59 183 50.00
9 10 9.5 66 18.03 339 92.62 93 25.41
10 11 10.5 27 7.38 366 100.00 27 7.38

GRÁFICAS

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

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

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

#Gráfico No. 5
# DIAGRMA DE CAJA Y BIGOTE
Caja<-boxplot(TempMin,horizontal = T, col = "blue",main="Gráfica No. 2.5:Distribución de la cantidad de la temperatura mínima \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. 2.6: Ojivas combinadas de la temperatura mínima \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("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. 2.6: Ojivas combinadas de la temperatura mínima \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("left",legend = c("Ojiva descendente", "Ojiva ascendente"),col = c("black", "blue"), pch = 1, lty = 1,cex = 0.7)

INDICADORES

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

mediana<-median(TempMin)
mediana
## [1] 8.001
media_aritmetica<-mean(TempMin)
media_aritmetica
## [1] 8.045505
Mo<-c("[7,8]")
desviación_estandar<-sd(TempMin)
desviación_estandar
## [1] 1.374138
coeficiente_variabilidad <- (desviación_estandar/media_aritmetica) * 100 
coeficiente_variabilidad
## [1] 17.07957
As<-skewness(TempMin)
As
## [1] -0.4343909
curtosis<-kurtosis(TempMin)
curtosis
## [1] 0.56792
#TABLAS
Variable<-c("Temperatura mínima (°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. 2.2: Indicadores estadísticos de la variable Temperatura mínima de cada uno de los registros de clima en el volcán Antisana.")
Tabla No. 2.2: Indicadores estadísticos de la variable Temperatura mínima 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 mínima (°C) 2.65 10.853 8.05 8.001 [7,8] 1.37 17.08 -0.43 0.57

CONCLUSIONES: CONCLUSIONES: La temperatura mínima en el volcán Antisana fluctúa entre 2.65 °C y 10.85 °C y sus valores están en torno a los 8.05 °C, con una desviación estándar de 1.37 siendo un conjunto de valores ligeramente homogéneos cuyos valores se concentran en la parte media alta de la variable y con un sesgo hacia la izquierda, por lo tanto, el comportamiento de la variable es beneficiosa. La mayoría de tiempo la temperatura mínima se encuentra a más de 8°C, estas temperaturas no perjudican a la prosperidad del ecosistema del lugar.

.