Licenciamento

d <- read.csv("licenciamentoS.csv")

Modelo autoregressivo AR3:

d2 <- cbind(d[5:150,]$LICENCIAMENTO, 
            d[4:149,]$LICENCIAMENTO, 
            d[3:148,]$LICENCIAMENTO, 
            d[2:147,]$LICENCIAMENTO)

colnames(d2) <- c("LICENCIAMENTO", "LICENCIAMENTO-1", "LICENCIAMENTO-2", "LICENCIAMENTO-3")
d2 <- as.data.frame(d2)
m1 <- lm(LICENCIAMENTO ~ ., data=d2)
summary(m1)

## 
## Call:
## lm(formula = LICENCIAMENTO ~ ., data = d2)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -84747 -20070  -2435  13662 130234 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)    
## (Intercept)       8.001e+03  9.394e+03   0.852    0.396    
## `LICENCIAMENTO-1` 5.230e-01  7.810e-02   6.697  4.6e-10 ***
## `LICENCIAMENTO-2` 6.848e-02  8.945e-02   0.766    0.445    
## `LICENCIAMENTO-3` 3.679e-01  7.883e-02   4.667  7.0e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 34380 on 142 degrees of freedom
## Multiple R-squared:  0.8158, Adjusted R-squared:  0.812 
## F-statistic: 209.7 on 3 and 142 DF,  p-value: < 2.2e-16

Yjun <- d[1, c("LICENCIAMENTO")]
Ymai <- d[2, c("LICENCIAMENTO")]
Yabr <- d[3, c("LICENCIAMENTO")]

alf1 <- m1$coefficients[2]
alf2 <- m1$coefficients[3]
alf3 <- m1$coefficients[4]
beta <- m1$coefficients[1]

Yjul <- alf1*Yjun + alf2*Ymai + alf3*Yabr + beta

Yago <- alf1*Yjul + alf2*Yjun + alf3*Ymai + beta

Yset <- alf1*Yago + alf2*Yjul + alf3*Yjun + beta

Julho:

## `LICENCIAMENTO-1` 
##          214415.7

Agosto:

## `LICENCIAMENTO-1` 
##          212944.9

Setembro:

## `LICENCIAMENTO-1` 
##          212241.9