#leer desde excel
library(readxl)
#series de tiempo
library(tseries)
library(forecast)
#Carga de datos de excel
Temper<-read_excel("C:\\Angel\\Programas\\R\\Temperatura.xlsx", col_names=FALSE)
#Transformando los datos a series de tiempO
#Grafica de temperaturas
Temp2<-ts(Temper, start = c(2006,1), frequency = 365)
ts.plot(Temp2, main="TEMPERATURA DEL PERIODO 2006-2017 en C.U. UAS",xlab="Años", ylab="Temperatura Media",col="blue")

#Media temperatura
mean(Temp2)
## [1] 26.29294
#Desviacion estandar temperatura
sd(Temp2)
## [1] 4.026578
#Resumen estadístico
summary(Temp2)
##       X__1      
##  Min.   :10.90  
##  1st Qu.:23.10  
##  Median :26.70  
##  Mean   :26.29  
##  3rd Qu.:29.80  
##  Max.   :35.00
#Varianza
var(Temp2)
##          X__1
## X__1 16.21333
#Cuartiles
quantile(Temp2)
##   0%  25%  50%  75% 100% 
## 10.9 23.1 26.7 29.8 35.0
#Agregando más información (PRECIPITACION).
#Carga de datos de precipitacion
Prec<-read_excel("C:\\Angel\\Programas\\R\\Precipitacion.xlsx", col_names=FALSE)
#Transformando los datos a series de tiempO
Prec2<-ts(Prec, start = c(2006,1), frequency = 365)
## Warning in data.matrix(data): NAs introducidos por coerción
#frequency es para determinar el intervalo entre datos, 365->dias, 12-> mes, 1->años
#Resumen estadístico
summary(Prec2)
##       X__1       
##  Min.   :  0.00  
##  1st Qu.:  1.20  
##  Median :  4.10  
##  Mean   : 11.34  
##  3rd Qu.: 13.43  
##  Max.   :250.30  
##  NA's   :3643
#ANALISIS PREVIOS

#Grafica
plot(Prec,type="l", lwd=1.5, col="blue", xlab="Datos", ylab="Milimetros", main="Precipitación en C.U. UAS 2006-2017")
## Warning in xy.coords(x, y, xlabel, ylabel, log): NAs introducidos por
## coerción

#Histograma
hist(Prec2, prob=T, col="grey", main="Histograma", xlab="Milimetros", ylab="Frecuencia")

#GENERAR MAPAS
library(dismo) #Te carga las bibliotecas necesarias
## Loading required package: raster
## Loading required package: sp
library(ggmap)
## Loading required package: ggplot2
## 
## Attaching package: 'ggmap'
## The following object is masked from 'package:dismo':
## 
##     geocode
library(ggplot2)
library(maptools)
## Checking rgeos availability: TRUE
#Mapa de terreno
Estacion<-geom_point(aes(x=-107.38, y=24.827), col="red", size=5)
Mapa<-get_map(location=c(-107.38, 24.827), source="google", maptype = "terrain", zoom=10)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=24.827,-107.38&zoom=10&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false
ggmap(Mapa)

ggmap(Mapa)+Estacion

#Mapa satelite
Mapa2<-get_map(location=c(-107.38, 24.827), source="google", maptype = "terrain", zoom=18)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=24.827,-107.38&zoom=18&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false
ggmap(Mapa2)

Mapafinal<-ggmap(Mapa2)+Estacion+labs(title="Ubicación de la estación meteorológica", caption="Fuente: ")
Mapafinal