### GWR dengue - modelar os casos de dengue de 2005_6 em relacao as medidas de controle e IB deste mesmo per?odo

# carregar os pacotes

library(maptools)
library (spdep)
library(sf)
library(lattice)
library(lctools) # Poisson GWR

source("HighstatLibV9.R")

#getwd()

## lendo shape

dengue <- st_read("./bc_aula10/deng_med_control_IB_2000_6.shp") # lendo mapa

plot(st_geometry(dengue))
names(dengue)
str(dengue)

## Análise exploratória 

# transformar banco com geometria em dataframe

dengue.df <- as.data.frame(dengue)

# outliers

str(dengue.df)

dotchart(dengue$INC_05_06) #OK!

MyVar <- c("A4COB05_06","A5COB05_06","A6COB05_06","CC_5_6","IB_5_6")
Mydotplot(dengue.df[,MyVar])

# outlier em 'A6COB05_06'

# raiz quadrada

rz_A6COB05_06 <- sqrt(dengue.df$A5COB05_06)
dotchart(rz_A6COB05_06) # ok!

# incluindo covar nos banco com geometria e sem geometria

dengue$rz_A6COB05_06 <-  rz_A6COB05_06
#str(dengue) # ok!

dengue.df$rz_A6COB05_06 <-  rz_A6COB05_06
#str(dengue.df) # ok!

# colinearidade

MyVar <- c("A4COB05_06","A5COB05_06","rz_A6COB05_06","CC_5_6","IB_5_6")
Mypairs(dengue.df[,MyVar])

#VIF

corvif(dengue.df[,MyVar])

# remover A5COB05_06 ou rz_A6cof05_06 ?

MyVar <- c("A4COB05_06","rz_A6COB05_06","CC_5_6","IB_5_6")
corvif(dengue.df[,MyVar]) # ok!

MyVar <- c("A4COB05_06","A5COB05_06","CC_5_6","IB_5_6")
corvif(dengue.df[,MyVar]) # ok!

# deixar para decidir olhando relacao ente X e covars

# relacao x versus y - usar tx de incid

MyX  <- c("A4COB05_06","A5COB05_06","rz_A6COB05_06","CC_5_6","IB_5_6")
Myxyplot(dengue.df, MyX, "INC_05_06", MyYlab = "Inc 05 06")

# retiramos A5COB05_06!

### Rodando modelo GWR

## coordenadas dos centroides

deng.centr <- st_geometry(st_centroid(dengue))
class(deng.centr)

# obtendo a matriz com as coord dos centroides

deng.coord <- st_coordinates(deng.centr)
class(deng.coord)

#str(dengue)

?gw.glm.bw # so para Poisson

bw <- gw.glm.bw(TOTAL0506 ~ A4COB05_06 + A5COB05_06 + rz_A6COB05_06 + CC_5_6 + IB_5_6, "poisson", dengue, deng.coord, kernel = 'adaptive', algorithm="exhaustive",optim.method="Nelder-Mead")

## usar 67 viznhos - ver gráfico!

?gw.glm

gwpr <- gw.glm(TOTAL05 ~ A4COB05 + A5COB05 + r_A6COB + CC_5_6 + IB_5_6, "poisson", dengue, 67, kernel = 'adaptive', deng.coord)

summary(gwpr$GGLM_LEst)
summary(gwpr$GGLM_LPvalues)
summary(gwpr$GGLM_GofFit)


#Exportando os resultados
write.table(gwpr$GGLM_LEst, file = "gwpr.Lest.txt", sep = ",", col.names = NA, qmethod = "double")

write.table(gwpr$GGLM_LPvalues, file = "gwpr.LPvalues.txt", sep = ",", col.names = NA, qmethod = "double")

write.table(gwpr$GGLM_GofFit, file = "gwpr.GofFit.txt", sep = ",", col.names = NA, qmethod = "double")