Caracterização dos Isolados

Session settings »

pkg <- c("ggplot2", "knitr","pander", "ggmap")
sapply(pkg, library, character.only=TRUE, logical.return=TRUE)

## ggplot2   knitr  pander   ggmap 
##    TRUE    TRUE    TRUE    TRUE

opts_chunk$set(
    tidy=FALSE,
    fig.width=6,
    fig.height=5,
    fig.align="center",
    dpi=150,
    dev="png",
    dev.args=list(png=list(family="Ubuntu Light", type="cairo")))
options(width=80)

Dataset »

#str(dat)
#head(dat)
dat$inc =round(dat[,"posit"]/dat[,"tot"],1)
dat$sev = round(rowMeans(dat[,12:20], na.rm=TRUE),2)
dat$int= dat$inc * dat$sev
attach(dat)
df= data.frame(uf, iso, inc, sev, int)
df=df[with(df, order(-int)), ]
rownames(df) <- NULL
names(df) = c("UF","Isolados", "Incidência",  "Severidade", "Intensidade"  )
#pandoc.table(df, split.tables=Inf, style='rmarkdown')

Tabela resumo » 7 dias após inoculação

UF	Isolados	Incidência	Severidade	Intensidade
MT	Querência	1	7.78	7.78
MT	Campo novo Parecis	1	7.44	7.44
MT	Sinop	1	5.56	5.56
GO	Río Verde	1	5	5
GO	Goiatuba	1	4.89	4.89
MT	Nova Xavantina	1	3.89	3.89
MA	Tasso Fragoso	1	3.22	3.22
TO	Porto Nacional	0.7	4.5	3.15
MT	Campo Verde	1	2.28	2.28
PR	Mauã da Serra	1	1.89	1.89
TO	Afonso	0.8	2.28	1.824
PR	Londrina	0.9	1.33	1.197
MT	Itiquira	0.6	1.5	0.9
PR	Sarandí	0.7	0.83	0.581
MT	Jaciara	0.6	0.69	0.414
MT	Nova Mutum	0.6	0.64	0.384
MA	Riachao	0.5	0.58	0.29
MA	Mangabeiras	0.4	0.64	0.256
GO	Morrinhos	0.3	0.62	0.186
TO	Brejinho de Nazaré	0.1	1	0.1

dat = subset(dat, !iso=="Sorriso")
location = c(-50, -8,-49,-28) #location(esqu izq/inf, esq der/sup)
ma_map <- qmap(location=location, zoom = 5)
#ma_map + labs(title="c(-50, -8,-49,-28)zoom = 5")

ma_map=ma_map + 
            geom_point(data = dat, aes(x = long, y = lat, color=int), # size=inci, 
                       alpha=I(7/10), size = 5) + 
  scale_colour_gradient(low="yellow", high="red") + 
  #labs(title="Intensidade mancha alvo")   + 
  scale_size_continuous(range= c(1, 25))  
  
ma_map

Caracterização dos Isolados

Juan Pablo Edwards

21/12/2014