require(geoR)
require(MASS)
##
## Análises geoestatística dos dados geoquimicos
##
teores <- read.table("aguafix.csv", head=T, sep=";")
names(teores)

summary(teores)

ag <- as.geodata(teores, coords.col=2:3, data.col=8)
plot(ag)

al <- as.geodata(teores, coords.col=2:3, data.col=9)
plot(al)
plot(al, lam=0)

b <- as.geodata(teores, coords.col=2:3, data.col=10)
plot(b)

##
## ba   dummy 2 ou 3 niveis?
## indicadora de regioes? em que partição? cortes? onde?
ba <- as.geodata(teores, coords.col=2:3, data.col=11)
plot(ba)
plot(ba, lam=0)

##
## ca 
##
ca <- as.geodata(teores, coords.col=2:3, data.col=12)
plot(ca)
plot(ca, lam=0) ## tem um nitidamente discrepante aqui < -2


cd <- as.geodata(teores, coords.col=2:3, data.col=13)
plot(cd)
co <- as.geodata(teores, coords.col=2:3, data.col=14)
plot(co)
cr <- as.geodata(teores, coords.col=2:3, data.col=15)
plot(cr)
cu <- as.geodata(teores, coords.col=2:3, data.col=16)
plot(cu)

##
## fe : problema com baixos valores (left censoring?)
##
fe <- as.geodata(teores, coords.col=2:3, data.col=17)
plot(fe)
plot(fe, lam=0)

ga <- as.geodata(teores, coords.col=2:3, data.col=18)
plot(ga)

## dummy 2 niveis? (lim minimo e maior que minimo?
ini <- as.geodata(teores, coords.col=2:3, data.col=19)
plot(ini)
plot(ini, lam=0)

## examinar dados atípicos > 20 ou mesmo > 10
k <- as.geodata(teores, coords.col=2:3, data.col=20)
plot(k)
plot(k, lam=0)
ksub1 <- subset(k, (data > 0 & data < 30))
plot(ksub1)
plot(ksub1, lam=0)
summary(ksub1, lam=0)
plot(variog(ksub1, max.dist=400000, lam=0))

li <- as.geodata(teores, coords.col=2:3, data.col=21)
plot(li)

##
## dados atipicos > 50
##
mg <- as.geodata(teores, coords.col=2:3, data.col=22)
plot(mg)
plot(mg, lam=0)

##
## mn valores atipicos > 25 ou > 100
## dummy 2 niveis
mn <- as.geodata(teores, coords.col=2:3, data.col=23)
plot(mn)
plot(mn, lam=0)

mo <- as.geodata(teores, coords.col=2:3, data.col=24)
plot(mo)

## na: 2 valores atípicos > 200
na <- as.geodata(teores, coords.col=2:3, data.col=25)
plot(na)
nasub1 <- subset(na, data < 200)
plot(nasub1)
plot(nasub1, lam=0)

ni <- as.geodata(teores, coords.col=2:3, data.col=26)
plot(ni)

pb <- as.geodata(teores, coords.col=2:3, data.col=27)
plot(pb)

##
## sr: valores atípicos (?) > 0.5
#
sr <- as.geodata(teores, coords.col=2:3, data.col=28)
plot(sr)
plot(sr, lam=0)

tl <- as.geodata(teores, coords.col=2:3, data.col=29)
plot(tl)
v <- as.geodata(teores, coords.col=2:3, data.col=30)
plot(v)
w <- as.geodata(teores, coords.col=2:3, data.col=31)
plot(w)

zn <- as.geodata(teores, coords.col=2:3, data.col=32)
plot(zn)
znsub <- subset(zn, data < 0.5)
plot(znsub)
plot(znsub, lam=0)

##
## f: valore atípioc(?) > 0.8 -- problemas com medições?
##
f <- as.geodata(teores, coords.col=2:3, data.col=33)
plot(f)
plot(f, lam=0)

##
## cl : 2 valores atípicos > 200
##
cl <- as.geodata(teores, coords.col=2:3, data.col=34)
plot(cl)
clsub <- subset(cl, data < 200)
plot(clsub)
plot(clsub, lam=0)

##
## dummy??? 2 níveis 
##
no2 <- as.geodata(teores, coords.col=2:3, data.col=35)
plot(no2)

## br: 2 valores atipicos > 1
br <- as.geodata(teores, coords.col=2:3, data.col=36)
plot(br)
brsub <- subset(br, data < 2)
plot(brsub)
## !!!!! faixas de medição??
plot(brsub, lam=0)

##
## no3: 1 atipico > 20
##
no3 <- as.geodata(teores, coords.col=2:3, data.col=37)
plot(no3)
plot(no3, lam=0)

##
## po4
## permanece assimetrico no log. limtes nas medidas
po4 <- as.geodata(teores, coords.col=2:3, data.col=38)
plot(po4)
plot(po4, lam=0)

## so4: 3 atipicos > 50
so4 <- as.geodata(teores, coords.col=2:3, data.col=39)
plot(so4)
so4sub <- subset(so4, data < 50)
plot(so4sub)
plot(so4sub, lam=0)


## condu: 2 atipicos > 1000
condu <- as.geodata(teores, coords.col=2:3, data.col=40)
plot(condu)
condusub <- subset(condu, data < 2000)
plot(condusub)
plot(condusub, lam=0)

##
## ph: 2 atipicos < 4
##
ph <- as.geodata(teores, coords.col=2:3, data.col=41)
plot(ph)
phsub <- subset(ph, data > 4)
plot(phsub)