Análisis de Temperatura Superficial del Mar de cada Estación Oceanográfica con respecto a datos de actividad sismicas en Perú

Autor: Morante Suarez, Diego Humberto

Paquete utilizado para la la limitación de datos:

library(tidyverse)

## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --

## v ggplot2 3.3.5     v purrr   0.3.4
## v tibble  3.1.6     v dplyr   1.0.7
## v tidyr   1.1.4     v stringr 1.4.0
## v readr   2.1.1     v forcats 0.5.1

## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

Datos utilizados: https://drive.google.com/drive/folders/1qIFXwENnjtKkbfrO20RgBK8Z7liXoCyJ?usp=sharing

Lectura de data:

sismos <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Sismos IGP/IGP_datos_sismicos 2017-2021.csv", sep = ";")
view(sismos)

Delimitación de datos sismicos con respecto a cada estación oceanográfica de IMARPE:

Coordenadas de la estación de Pisco de latitud es Y=-13.70927778

sismos %>% filter(Y <= -12.88766389, Y >= -14.53089167, X < -76.0, magnitud..M. >= 5)

##   ï..fecha.UTC hora.UTC        Y        X profundidad..km. magnitud..M. Dia Mes
## 1   27/01/2017  27:56.0 -14.0856 -76.3863               58          5.5  27   1
## 2   28/01/2017  39:38.0 -14.5102 -76.5473               21          5.2  28   1
## 3   12/08/2017  30:43.0 -13.7946 -76.4816               46          5.1  12   8
## 4    8/11/2017  13:12.0 -13.5205 -77.0244               48          5.1   8  11
## 5    3/09/2018  24:19.0 -14.0271 -76.2349               35          5.0   3   9
## 6   13/10/2018  13:19.0 -13.9184 -76.6090               26          5.5  13  10
## 7   11/04/2019  50:31.5 -13.6574 -76.3434               64          5.1  11   4
## 8    5/12/2019  23:18.4 -13.8251 -76.2927               71          5.3   5  12
##   AÃ.o
## 1 2017
## 2 2017
## 3 2017
## 4 2017
## 5 2018
## 6 2018
## 7 2019
## 8 2019

Coordenadas de la estación de Callao de latitud es Y=-12.06605

sismos %>% filter(Y <= -11.24443611, Y >= -12.88766389, X < -76.0, magnitud..M. >= 5)

##    ï..fecha.UTC hora.UTC        Y        X profundidad..km. magnitud..M. Dia
## 1     6/01/2017  55:27.0 -11.4853 -78.5790               28          5.4   6
## 2     9/07/2017  07:24.0 -11.8245 -77.5660               54          5.2   9
## 3    14/09/2017  19:23.0 -11.9032 -76.3349               11          5.2  14
## 4    23/10/2017  09:11.0 -12.4138 -77.4528               34          5.1  23
## 5    24/11/2017  15:52.0 -12.1577 -76.2506              110          5.1  24
## 6     5/02/2018  27:06.0 -12.7538 -78.1940               31          5.2   5
## 7    17/05/2018  07:12.0 -12.7852 -76.8737               49          5.8  17
## 8    22/05/2018  08:36.0 -11.2590 -78.1469               48          5.3  22
## 9    27/09/2018  00:07.0 -11.5097 -78.1258               52          5.3  27
## 10   20/04/2021 04:53:29 -12.6480 -76.8788               44          5.0  20
## 11   23/06/2021 02:54:17 -12.7668 -77.0052               45          6.0  23
## 12   25/08/2021 07:15:09 -12.4127 -77.5699               31          5.0  25
##    Mes AÃ.o
## 1    1 2017
## 2    7 2017
## 3    9 2017
## 4   10 2017
## 5   11 2017
## 6    2 2018
## 7    5 2018
## 8    5 2018
## 9    9 2018
## 10   4 2021
## 11   6 2021
## 12   8 2021

Coordenadas de la estación de Huacho de latitud es Y=-11.12189722

sismos %>% filter(Y <= -10.99935833, Y >= -11.24443611, X < -77.4, magnitud..M. >= 4.5)

##   ï..fecha.UTC hora.UTC        Y        X profundidad..km. magnitud..M. Dia Mes
## 1    2/04/2017  48:54.0 -11.1003 -78.3529               39          4.6   2   4
## 2   28/06/2017  50:30.0 -11.0131 -78.0862               66          4.8  28   6
## 3    1/11/2017  19:50.0 -11.1984 -78.1483               42          4.6   1  11
## 4    5/02/2021 23:55:15 -11.0763 -78.2306               36          4.8   5   2
##   AÃ.o
## 1 2017
## 2 2017
## 3 2017
## 4 2021

Coordenadas de la estación de Huanchaco de latitud es Y=-8.078111111

sismos %>% filter(Y <= -7.424472222, Y >= -8.73175, X < -79.4, magnitud..M. >= 4.5, AÃ.o == 2018)

##    ï..fecha.UTC hora.UTC       Y        X profundidad..km. magnitud..M. Dia Mes
## 1     2/01/2018  13:23.0 -8.5175 -80.0582               43          4.6   2   1
## 2    17/01/2018  29:22.0 -8.4244 -79.4254               52          4.7  17   1
## 3    22/01/2018  11:35.0 -8.5037 -80.7142                4          4.6  22   1
## 4    27/02/2018  45:14.0 -8.6614 -80.9012               26          4.5  27   2
## 5     7/03/2018  46:53.0 -8.1456 -79.8195               31          4.6   7   3
## 6    31/03/2018  52:19.0 -8.4371 -80.2103                6          4.6  31   3
## 7     8/05/2018  32:32.0 -8.0692 -79.4353               56          4.9   8   5
## 8    29/06/2018  44:19.0 -8.5748 -79.7173               34          4.6  29   6
## 9    30/07/2018  46:39.0 -8.6246 -79.7349               36          4.6  30   7
## 10    2/10/2018  40:51.0 -7.6881 -79.5124               55          4.7   2  10
## 11   18/10/2018  25:53.0 -7.5070 -79.9704               37          4.8  18  10
##    AÃ.o
## 1  2018
## 2  2018
## 3  2018
## 4  2018
## 5  2018
## 6  2018
## 7  2018
## 8  2018
## 9  2018
## 10 2018
## 11 2018

Coordenadas de la estación de San Jose de latitud es Y=-6.770833333

sismos %>% filter(Y <= -6.117194444, Y >= -7.424472222, X < -79.4, magnitud..M. >= 4.5, AÃ.o == 2018)

##    ï..fecha.UTC hora.UTC       Y        X profundidad..km. magnitud..M. Dia Mes
## 1    21/01/2018  13:00.0 -6.8134 -81.0015               37          4.5  21   1
## 2    24/01/2018  52:50.0 -6.9624 -81.1062               29          4.5  24   1
## 3    10/02/2018  54:40.0 -7.2179 -80.4496               33          4.6  10   2
## 4     6/05/2018  01:11.0 -6.4461 -80.8868               24          4.9   6   5
## 5    25/05/2018  21:22.0 -6.5523 -80.9063               11          4.9  25   5
## 6    25/05/2018  25:19.0 -6.5493 -81.0147               27          4.9  25   5
## 7    25/05/2018  20:57.0 -6.5554 -80.9374               26          4.7  25   5
## 8     3/06/2018  01:29.0 -6.3607 -81.1835               43          4.6   3   6
## 9    26/07/2018  30:58.0 -7.1972 -80.6060               10          4.7  26   7
## 10    2/10/2018  40:11.0 -6.5238 -80.6658              191          5.1   2  10
## 11   10/10/2018  48:45.0 -6.5226 -81.2228               22          4.9  10  10
## 12   17/10/2018  27:25.0 -6.9735 -80.0897               50          4.5  17  10
##    AÃ.o
## 1  2018
## 2  2018
## 3  2018
## 4  2018
## 5  2018
## 6  2018
## 7  2018
## 8  2018
## 9  2018
## 10 2018
## 11 2018
## 12 2018

Coordenadas de la estación de Tumbes de latitud es Y=-12.06605

sismos %>% filter( Y >= -4.5, X < -80.0, magnitud..M. >= 5)

##    ï..fecha.UTC hora.UTC       Y        X profundidad..km. magnitud..M. Dia Mes
## 1    21/03/2017  29:08.0 -3.5874 -81.1816               47          5.2  21   3
## 2     5/06/2017  34:08.0 -4.1623 -80.6432               52          6.2   5   6
## 3     6/10/2017  32:08.0 -4.2772 -80.8156               47          5.1   6  10
## 4     7/12/2017  37:44.0 -4.3129 -80.2310               28          5.1   7  12
## 5    16/08/2018  09:14.0 -3.1327 -81.7365               48          5.0  16   8
## 6    29/10/2018  20:50.0 -3.8091 -80.2878               53          5.2  29  10
## 7     4/07/2019  20:57.8 -4.3062 -80.5699               50          5.0   4   7
## 8     3/01/2020 14:16:39 -3.3280 -80.2922               57          5.0   3   1
## 9    25/08/2020 05:18:09 -4.1455 -80.9244               31          5.5  25   8
## 10   18/10/2020 06:58:26 -3.7723 -82.0227               27          5.4  18  10
## 11   13/02/2021 05:18:31 -4.1493 -80.4579               75          5.4  13   2
##    AÃ.o
## 1  2017
## 2  2017
## 3  2017
## 4  2017
## 5  2018
## 6  2018
## 7  2019
## 8  2020
## 9  2020
## 10 2020
## 11 2021

Gráficas de Temperatura Superficial del Mar de cada Estación

Leyenda de gráficas:

linea verdes: Variacion de Temperatura superficial del Mar
Linea rojas: Dia del sismo
Linea negras: Desviación estándar

Estación de Tumbes:

Analisis con respecto al sismo 16 de Agosto 2018 - Tumbes

anl_tumb1 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/TUMBES/tumbes1.csv", sep = ";")

Media:

med_tumb1 <- mean(anl_tumb1$TSM)

## Warning in mean.default(anl_tumb1$TSM): argument is not numeric or logical:
## returning NA

Desviación estándar:

ds_tumb1 <- sd(anl_tumb1$TSM)

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introducidos por coerción

Máximo:

max_tumb1 <- med_tumb1 + ds_tumb1

Mínimo:

min_tumb1 <- med_tumb1 - ds_tumb1

Gráfica:

plot(anl_tumb1$Dias, anl_tumb1$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Analisis con respecto al sismo 16 de Agosto 2018 - Tumbes")

## Warning in xy.coords(x, y, xlabel, ylabel, log): NAs introducidos por coerción

## Warning in xy.coords(x, y, xlabel, ylabel, log): NAs introducidos por coerción

abline(h =c(med_tumb1,max_tumb1,min_tumb1), 
       lty = 10, 
       col = "black")

abline(v = c(16), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 29 de Octubre 2018 - Tumbes

anl_tumb2 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/TUMBES/tumbes2.csv", sep = ";")

Media:

med_tumb2 <- mean(anl_tumb2$TSM)

Desviación estándar:

ds_tumb2 <- sd(anl_tumb2$TSM)

Máximo:

max_tumb2 <- med_tumb2 + ds_tumb2

Mínimo:

min_tumb2 <- med_tumb2 - ds_tumb2

Gráfica:

plot(anl_tumb2$Dias, anl_tumb2$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 29 de Octubre 2018 - Tumbes")

abline(h =c(med_tumb2,max_tumb2,min_tumb2), 
       lty = 10, 
       col = "black")

abline(v = c(29), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 4 de Julio 2019 - Tumbes

anl_tumb3 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/TUMBES/tumbes3.csv", sep = ";")

Media:

med_tumb3 <- mean(anl_tumb3$TSM)

Desviación estándar:

ds_tumb3 <- sd(anl_tumb3$TSM)

Máximo:

max_tumb3 <- med_tumb3 + ds_tumb3

Mínimo:

min_tumb3 <- med_tumb3 - ds_tumb3

Gráfica:

plot(anl_tumb3$Dias, anl_tumb3$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 4 de Julio 2019 - Tumbes")

abline(h =c(med_tumb3,max_tumb3,min_tumb3), 
       lty = 10, 
       col = "black")

abline(v = c(4), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 3 de Enero 2020 - Tumbes

anl_tumb4 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/TUMBES/tumbes4.csv", sep = ";")

Media:

med_tumb4 <- mean(anl_tumb4$TSM)

Desviación estándar:

ds_tumb4 <- sd(anl_tumb4$TSM)

Máximo:

max_tumb4 <- med_tumb4 + ds_tumb4

Mínimo:

min_tumb4 <- med_tumb4 - ds_tumb4

Gráfica:

plot(anl_tumb4$Dias, anl_tumb4$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 3 de Enero 2020 - Tumbes")

abline(h =c(med_tumb4,max_tumb4,min_tumb4), 
       lty = 10, 
       col = "black")

abline(v = c(3), 
       lty = 10, 
       col = "red")

legend("bottomright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Estación de Pisco:

Análisis con respecto al sismo 3 de Setiembre 2018 - Pisco

anl_pisc1 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/PISCO/pisco1.csv", sep = ";")

Media:

med_pisc1 <- mean(anl_pisc1$TSM)

Desviación estándar:

ds_pisc1 <- sd(anl_pisc1$TSM)

Máximo:

max_pisc1 <- med_pisc1 + ds_pisc1

Mínimo:

min_pisc1 <- med_pisc1 - ds_pisc1

Gráfica:

plot(anl_pisc1$Dias, anl_pisc1$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 3 de Setiembre 2018 - Pisco")

abline(h =c(med_pisc1,max_pisc1,min_pisc1), 
       lty = 10, 
       col = "black")

abline(v = c(3), 
       lty = 10, 
       col = "red")

legend("bottomright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 13 de Octubre 2018 - Pisco

anl_pisc2 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/PISCO/pisco2.csv", sep = ";")

Media:

med_pisc2 <- mean(anl_pisc2$TSM)

Desviación estándar:

ds_pisc2 <- sd(anl_pisc2$TSM)

Máximo:

max_pisc2 <- med_pisc2 + ds_pisc2

Mínimo:

min_pisc2 <- med_pisc2 - ds_pisc2

Gráfica:

plot(anl_pisc2$Dias, anl_pisc2$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 13 de Octubre 2018 - Pisco")

abline(h =c(med_pisc2,max_pisc2,min_pisc2), 
       lty = 10, 
       col = "black")

abline(v = c(13), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Estación de Callao

Análisis con respecto al sismo 5 de Febrero 2018 - Pisco

anl_calla1 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/CALLAO/callao1.csv", sep = ";")

Media:

med_calla1 <- mean(anl_calla1$TSM)

Desviación estándar:

ds_calla1 <- sd(anl_calla1$TSM)

Máximo:

max_calla1 <- med_calla1 + ds_calla1

Mínimo:

min_calla1 <- med_calla1 - ds_calla1

Gráfica:

plot(anl_calla1$Dias, anl_calla1$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 5 de Febrero 2018 - Pisco")

abline(h =c(med_calla1,max_calla1,min_calla1), 
       lty = 10, 
       col = "black")

abline(v = c(5), 
       lty = 10, 
       col = "red")

legend("bottomright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 17 y 22 de Mayo 2018 - Pisco

anl_calla2 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/CALLAO/callao2.csv", sep = ";")

Media:

med_calla2 <- mean(anl_calla2$TSM)

Desviación estándar:

ds_calla2 <- sd(anl_calla2$TSM)

Máximo:

max_calla2 <- med_calla2 + ds_calla2

Mínimo:

min_calla2 <- med_calla2 - ds_calla2

Gráfica:

plot(anl_calla2$Dias, anl_calla2$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 17 y 22 de Mayo 2018 - Pisco")

abline(h =c(med_calla2,max_calla2,min_calla2), 
       lty = 10, 
       col = "black")

abline(v = c(17, 22), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 27 de Setiembre 2018 - Pisco

anl_calla3 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/CALLAO/callao3.csv", sep = ";")

Media:

med_calla3 <- mean(anl_calla3$TSM)

Desviación estándar:

ds_calla3 <- sd(anl_calla3$TSM)

Máximo:

max_calla3 <- med_calla3 + ds_calla3

Mínimo:

min_calla3 <- med_calla3 - ds_calla3

Gráfica:

plot(anl_calla3$Dias, anl_calla3$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 27 de Setiembre 2018 - Pisco")

abline(h =c(med_calla3,max_calla3,min_calla3), 
       lty = 10, 
       col = "black")

abline(v = c(27), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Estación de Huanchaco

Análisis con respecto al sismo 2, 17 y 22 de Enero 2018 - Huanchaco

anl_hnchac1 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/HUANCHACO/huanchaco1.csv", sep = ";")

Media:

med_hnchac1 <- mean(anl_hnchac1$TSM)

Desviación estándar:

ds_hnchac1 <- sd(anl_hnchac1$TSM)

Máximo:

max_hnchac1 <- med_hnchac1 + ds_hnchac1

Mínimo:

min_hnchac1 <- med_hnchac1 - ds_hnchac1

Gráfica:

plot(anl_hnchac1$Dias, anl_hnchac1$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 2, 17 y 22 de Enero 2018 - Huanchaco")

abline(h =c(med_hnchac1,max_hnchac1,min_hnchac1), 
       lty = 10, 
       col = "black")

abline(v = c(2, 17, 22), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

**Análisis con respecto al sismo 27(-1) de Febrero, 7 y 31 de Marzo 2018

Huanchaco**

anl_hnchac2 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/HUANCHACO/huanchaco2.csv", sep = ";")

Media:

med_hnchac2 <- mean(anl_hnchac2$TSM)

Desviación estándar:

ds_hnchac2 <- sd(anl_hnchac2$TSM)

Máximo:

max_hnchac2 <- med_hnchac2 + ds_hnchac2

Mínimo:

min_hnchac2 <- med_hnchac2 - ds_hnchac2

Gráfica:

plot(anl_hnchac2$Dias, anl_hnchac2$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 27(-1) de Febrero, 7 y 31 de Marzo 2018 - Huanchaco")

abline(h =c(med_hnchac2,max_hnchac2,min_hnchac2), 
       lty = 10, 
       col = "black")

abline(v = c(-1, 7, 31), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 8 de Mayo 2018 - Huanchaco

anl_hnchac3 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/HUANCHACO/huanchaco3.csv", sep = ";")

Media:

med_hnchac3 <- mean(anl_hnchac3$TSM)

Desviación estándar:

ds_hnchac3 <- sd(anl_hnchac3$TSM)

Máximo:

max_hnchac3 <- med_hnchac3 + ds_hnchac3

Mínimo:

min_hnchac3 <- med_hnchac3 - ds_hnchac3

Gráfica:

plot(anl_hnchac3$Dias, anl_hnchac3$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 8 de Mayo 2018 - Huanchaco")

abline(h =c(med_hnchac3,max_hnchac3,min_hnchac3), 
       lty = 10, 
       col = "black")

abline(v = c(8), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 29 de Junio 2018 - Huanchaco

anl_hnchac4 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/HUANCHACO/huanchaco4.csv", sep = ";")

Media:

med_hnchac4 <- mean(anl_hnchac4$TSM)

Desviación estándar:

ds_hnchac4 <- sd(anl_hnchac4$TSM)

Máximo:

max_hnchac4 <- med_hnchac4 + ds_hnchac4

Mínimo:

min_hnchac4 <- med_hnchac4 - ds_hnchac4

Gráfica:

plot(anl_hnchac4$Dias, anl_hnchac4$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 29 de Junio 2018 - Huanchaco")

abline(h =c(med_hnchac4,max_hnchac4,min_hnchac4), 
       lty = 10, 
       col = "black")

abline(v = c(29), 
       lty = 10, 
       col = "red")

legend("bottomright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 30 de Julio 2018 - Huanchaco

anl_hnchac5 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/HUANCHACO/huanchaco5.csv", sep = ";")

Media:

med_hnchac5 <- mean(anl_hnchac5$TSM)

Desviación estándar:

ds_hnchac5 <- sd(anl_hnchac5$TSM)

Máximo:

max_hnchac5 <- med_hnchac5 + ds_hnchac5

Mínimo:

min_hnchac5 <- med_hnchac5 - ds_hnchac5

Gráfica:

plot(anl_hnchac5$Dias, anl_hnchac5$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 30 de Julio 2018 - Huanchaco")

abline(h =c(med_hnchac5,max_hnchac5,min_hnchac5), 
       lty = 10, 
       col = "black")

abline(v = c(30), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 2, 18 de Octubre 2018 - Huanchaco

anl_hnchac6 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/HUANCHACO/huanchaco6.csv", sep = ";")

Media:

med_hnchac6 <- mean(anl_hnchac6$TSM)

Desviación estándar:

ds_hnchac6 <- sd(anl_hnchac6$TSM)

Máximo:

max_hnchac6 <- med_hnchac6 + ds_hnchac6

Mínimo:

min_hnchac6 <- med_hnchac6 - ds_hnchac6

Gráfica:

plot(anl_hnchac6$Dias, anl_hnchac6$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 2, 18 de Octubre 2018 - Huanchaco")

abline(h =c(med_hnchac6,max_hnchac6,min_hnchac6), 
       lty = 10, 
       col = "black")

abline(v = c(2, 18), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Estación de San José

Análisis con respecto al sismo 21 y 24 de Enero, 10 de Febrero 2018 - San josé

anl_sanjose1 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/SAN JOSE/sanjose1.csv", sep = ";")

Media:

med_sanjose1 <- mean(anl_sanjose1$TSM)

Desviación estándar:

ds_sanjose1 <- sd(anl_sanjose1$TSM)

Máximo:

max_sanjose1 <- med_sanjose1 + ds_sanjose1

Mínimo:

min_sanjose1 <- med_sanjose1 - ds_sanjose1

Gráfica:

plot(anl_sanjose1$Dias, anl_sanjose1$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 21 y 24 de Enero, 10 de Febrero 2018 - San josé")

abline(h =c(med_sanjose1,max_sanjose1,min_sanjose1), 
       lty = 10, 
       col = "black")

abline(v = c(21,24,41), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 6 y 25 de Mayo, 3(34) de Junio 2018 - San josé

anl_sanjose2 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/SAN JOSE/sanjose2.csv", sep = ";")

Media:

med_sanjose2 <- mean(anl_sanjose2$TSM)

Desviación estándar:

ds_sanjose2 <- sd(anl_sanjose2$TSM)

Máximo:

max_sanjose2 <- med_sanjose2 + ds_sanjose2

Mínimo:

min_sanjose2 <- med_sanjose2 - ds_sanjose2

Gráfica:

plot(anl_sanjose2$Dias, anl_sanjose2$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 6 y 25 de Mayo, 3(34) de Junio 2018 - San josé")

abline(h =c(med_sanjose2,max_sanjose2,min_sanjose2), 
       lty = 10, 
       col = "black")

abline(v = c(6, 25,34), 
       lty = 10, 
       col = "red")

legend("top",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 26 de Julio 2018 - San josé

anl_sanjose3 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/SAN JOSE/sanjose3.csv", sep = ";")

Media:

med_sanjose3 <- mean(anl_sanjose3$TSM)

Desviación estándar:

ds_sanjose3 <- sd(anl_sanjose3$TSM)

Máximo:

max_sanjose3 <- med_sanjose3 + ds_sanjose3

Mínimo:

min_sanjose3 <- med_sanjose3 - ds_sanjose3

Gráfica:

plot(anl_sanjose3$Dias, anl_sanjose3$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 26 de Julio 2018 - San josé")

abline(h =c(med_sanjose3,max_sanjose3,min_sanjose3), 
       lty = 10, 
       col = "black")

abline(v = c(26), 
       lty = 10, 
       col = "red")

legend("topleft",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)

Análisis con respecto al sismo 2, 10 y 17 de Octubre 2018 - San José

anl_sanjose4 <- read.csv("C:/Users/diego/Documents/climato/Material para investigacion/Temperatura Superficial del Mar IMARPE/analisis/SAN JOSE/sanjose4.csv", sep = ";")

Media:

med_sanjose4 <- mean(anl_sanjose4$TSM)

Desviación estándar:

ds_sanjose4 <- sd(anl_sanjose4$TSM)

Máximo:

max_sanjose4 <- med_sanjose4 + ds_sanjose4

Mínimo:

min_sanjose4 <- med_sanjose4 - ds_sanjose4

Gráfica:

plot(anl_sanjose4$Dias, anl_sanjose4$TSM, 
     type="b", pch=17, col="forestgreen", 
     xlab="Día", ylab="TSM", 
     main = "Análisis con respecto al sismo 2, 10 y 17 de Octubre 2018 - San josé")

abline(h =c(med_sanjose4,max_sanjose4,min_sanjose4), 
       lty = 10, 
       col = "black")

abline(v = c(2, 10, 17), 
       lty = 10, 
       col = "red")

legend("topright",
       legend = c("Sismo", "TSM", "Desviacion estandar"),
       title = "LEYENDA",
       inset = 0.05,
       lty = c(1, 2),
       col = c(2, 3, 1),
       cex = 0.75, 
       lwd = 2)