
## Instalando/Carregando o Pacote
install.packages("RDengue",
                 contriburl = "http://www.leg.ufpr.br/~henrique/RDengue")
library(RDengue)



## Criando os objetos necessários
codbairro <- list(BT = 94,        # Brasilia Teimosa
                  CFP = c(54,55), # Casa Forte/Parnamirim
                  DI = c(39,40),  # Dois Irmaos
                  EM = 73,        # Engenho do Meio 
                  MCP = c(51,52)) # Morro da Conceição

## Lendo os dados do Banco
aed <- datachk(host = "guaja.est.ufpr.br",
               user = "usuario",
               pass = "senha",
               codBairros = codbairro)

## Criando a estrutura de Vizinhança
neigh(aed)


## Plotando armadilhas
plot(aed,
     "EM101",
     log = FALSE,
     leg = T)


## Plotando uma superfície Suavizada
smooth.arma(aed, "DI",
            data = 1, method = "lm",
            log = F, points = T, plot = T,
            statistics = F,
            contour = T, showTitle = F)

X11()
smooth.arma(aed, "EM",
            data = "2004-03-02",
            method = "loess",
            points = T, log = F,contour = T)


##Gera Objetos necessários para animação com a kfade:
preds <- geraPreds(aed, "MCP", method="gam", idx=1:108)
#   user  system elapsed 
# 95.958   0.356  98.992



## Gerando o filme(Necessita do mplayer e mencoder - $ sudo apt-get install mplayer mencoder)
## O filme e as imagens são salvos no diretorio especificado em 'loc'
system.time(plotDyna(preds, aed, loc = "~/imgMCP"))
#   user  system elapsed 
# 90.366  19.121 132.410 

##Demonstracao da Vizinhança por Delaunay
require(deldir)
dd <- deldir(coordinates(aed$pts)[1:20,][,1],
             (coordinates(aed$pts)[1:20,][,2]))
tt <- tile.list(dd)
plot(tt, polycol=colors()[1:20], close=T)
title("Bairro BT")


##Gerando os Relatórios em HTML
makeRelat(aed)
url1 <- data.frame(object_id  = sort(unique(aed$dados$COD_ARMADILHA)),
                   media_name = paste("http://www.leg.ufpr.br/~henrique/dengue/Relat/",
                                         sort(unique(aed$dados$COD_ARMADILHA)), sep = ""),
                   stringsAsFactors = FALSE)

url2 <- data.frame(object_id  = as.numeric(paste(2611606050, sprintf("%02d", 1:94), sep="")),
                   media_name = paste("http://www.leg.ufpr.br/~henrique/dengue/Relat/Bairros/",
                                         1:94, ".avi", sep = ""),
                   stringsAsFactors = FALSE)

writeonDb(host = 'guaja.est.ufpr.br',
          user = "usuario",
          pass = "senha",
          port = "3306",
          db = "newSaudavel", data=url1)


#Plotando gráfico TS
armaTs(aed, "EM101", cex.axis = 0.5, n.ahead = 5)