## Instalando/Carregando o Pacote install.packages("RDengue", contriburl = "http://www.leg.ufpr.br/~henrique/RDengue") library(RDengue) rm(list = ls(all = TRUE)) ## 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 = localhost, #"guaja.est.ufpr.br", user = "henrique", pass = "wbhdpj", 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 = "henrique", pass = "wbhdpj", port = "3306", db = "newSaudavel", data=url1) #Plotando gráfico TS armaTs(aed, "EM101", cex.axis = 0.5, n.ahead = 5) plot(rnorm(100), rnorm(100)) legend("bottomleft", legend = "Tipo 1", pch = 1, cex = 1) legend("topleft", legend = "Tipo 1", pch = 1, cex = 2) legend("topright", legend = "Tipo 1", pch = 1, cex = 1.5, bty = 'n') legend("center", legend = c("Tipo 1", "Tipo 2"), pch = 1, cex = 1.8, fill = c("red", "blue")) fac <- factanal(as.formula(paste("~", paste(grep("PRECIP|UMID", names(dados), value = T), collapse = "+"))), 4, data = dados, subset = SITIO == "DI", rotation = "promax") fac.list <- split(apply(fac$loadings, 1, max), apply(fac$loadings, 1, which.max)) fac.list <- fac.list[order(unlist(lapply(fac.list, function(x)as.numeric(gsub(patt="[a-Z]|\\.", rep = "", names(x)[1])))))] names(fac.list) <- 1:4 f <- function(x){ x <- x/sum(x) rowSums(sweep(dados[names(x)], 2, x, FUN="*")) } preumid <- sapply(fac.list, f) head(preumid) ### Temperaturas fac <- factanal(as.formula(paste("~", paste(grep("MAX|MIN", names(dados), value = T), collapse = "+"))), 3, data = dados, subset = SITIO == "DI", rotation = "promax") fac.list <- split(apply(fac$loadings, 1, max), apply(fac$loadings, 1, which.max)) fac.list <- fac.list[order(unlist(lapply(fac.list, function(x)as.numeric(gsub(patt="[a-Z]|\\.", rep = "", names(x)[1])))))] names(fac.list) <- 1:3 f <- function(x){ x <- x/sum(x) rowSums(sweep(dados[names(x)], 2, x, FUN="*")) } fac.temp <- sapply(fac.list, f) head(fac.temp) ## Criando o pacote package.skeleton(name = "RDengue", list = c("armaTs", "bubble.arma", "datachk", "geraPreds", "geraRelatArma", "geraRelatBairro", "is.arma", "kfade", "kfademaps", "ltitle", "makeRelat", "neigh", "plot.arma", "plotDyna", "print.arma", "print.arma.pred", "smooth.arma", "writeonDb"), force = TRUE, namespace = TRUE)