Ubatuba

Marco Arthur
10/03/2010

Chuvas em Ubatuba

  • Autor: Marco Arthur
  • Fonte de dados: Estação metereológica de Ubatuba, bairro Horto
  • Série Temporal: jan/1935 até jul/2001

Dados

rain <- read.table(header=T, '~/data/time_series_rain_ubatuba.txt')
summary(rain)

      Ano            jan             fev             mar       
 Min.   :1935   Min.   : 56.8   Min.   : 16.7   Min.   : 72.5  
 1st Qu.:1952   1st Qu.:252.8   1st Qu.:166.0   1st Qu.:208.8  
 Median :1968   Median :343.4   Median :295.6   Median :280.5  
 Mean   :1968   Mean   :368.5   Mean   :333.5   Mean   :328.1  
 3rd Qu.:1984   3rd Qu.:493.9   3rd Qu.:490.7   3rd Qu.:428.6  
 Max.   :2001   Max.   :787.7   Max.   :898.6   Max.   :958.5  

      abr             mai              jun              jul        
 Min.   : 26.5   Min.   : 10.00   Min.   :  6.60   Min.   :  7.20  
 1st Qu.:136.6   1st Qu.: 73.15   1st Qu.: 41.75   1st Qu.: 43.95  
 Median :180.9   Median :111.80   Median : 70.50   Median : 76.90  
 Mean   :218.5   Mean   :126.23   Mean   : 92.68   Mean   :101.54  
 3rd Qu.:267.4   3rd Qu.:168.85   3rd Qu.:137.30   3rd Qu.:125.70  
 Max.   :628.6   Max.   :367.60   Max.   :349.90   Max.   :433.60  
                                                   NA's   :1       
      ago              set             out             nov       
 Min.   :  7.50   Min.   : 21.9   Min.   : 69.0   Min.   : 59.2  
 1st Qu.: 57.67   1st Qu.:120.3   1st Qu.:152.8   1st Qu.:170.8  
 Median : 83.30   Median :180.4   Median :220.8   Median :244.9  
 Mean   :103.65   Mean   :180.5   Mean   :232.6   Mean   :267.3  
 3rd Qu.:119.42   3rd Qu.:242.9   3rd Qu.:311.9   3rd Qu.:327.9  
 Max.   :383.90   Max.   :376.8   Max.   :500.1   Max.   :717.2  
 NA's   :1        NA's   :1       NA's   :2       NA's   :2      
      dez            Anual     
 Min.   : 95.0   Min.   :1056  
 1st Qu.:236.6   1st Qu.:2223  
 Median :307.8   Median :2735  
 Mean   :327.0   Mean   :2645  
 3rd Qu.:380.7   3rd Qu.:3038  
 Max.   :732.1   Max.   :3854  
 NA's   :2

Médias mensais dados de 1935-2001 em (mm)

plot of chunk unnamed-chunk-2

Série temporal

plot of chunk unnamed-chunk-3

Regressão Linear

t <- seq(from=1935, to=(1935 + length(annualy) -1))
model <- lm(annualy ~ t + 1, na.action = na.exclude)
plot(annualy)
lines(t[1:length(model$fitted.values)], model$fitted.values, col="red")

plot of chunk unnamed-chunk-4