Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior | Próxima revisão Ambos lados da revisão seguinte | ||
disciplinas:ce097-2015-01:historico2015 [2015/03/25 10:59] paulojus |
disciplinas:ce097-2015-01:historico2015 [2015/04/13 22:19] paulojus |
||
---|---|---|---|
Linha 15: | Linha 15: | ||
| 01/04 |(Venessa e Bruno) | | | | 01/04 |(Venessa e Bruno) | | | ||
+ | |||
+ | === 13/04 === | ||
+ | |||
+ | <code R> | ||
+ | require(geoR) | ||
+ | data(s100) | ||
+ | plot(s100) | ||
+ | plot(s100, low=T) | ||
+ | plot(Ksat) | ||
+ | plot(ca20) | ||
+ | plot(parana) | ||
+ | plot(s100, low=T) | ||
+ | ml <- likfit(s100, ini=c(var(s100$data), 0.1)) | ||
+ | ml | ||
+ | g <- expand.grid(seq(0,1, l=11), seq(0,1, l=11)) | ||
+ | points(s100) | ||
+ | points(g) | ||
+ | points(g, col=2, pch=19) | ||
+ | g <- expand.grid(seq(0,1, l=51), seq(0,1, l=51)) | ||
+ | points(g, col=2, pch=19) | ||
+ | g1 <- expand.grid(seq(0,1, l=11), seq(0,1, l=11)) | ||
+ | g2 <- expand.grid(seq(0,1, l=51), seq(0,1, l=51)) | ||
+ | ml | ||
+ | kr1 <- krige.conv(s100, loc=gr1, krige=krige.control(obj.m=ml)) | ||
+ | kr1 <- krige.conv(s100, loc=g1, krige=krige.control(obj.m=ml)) | ||
+ | image(g1) | ||
+ | image(kr1) | ||
+ | points(s100, add=T) | ||
+ | kr2 <- krige.conv(s100, loc=g2, krige=krige.control(obj.m=ml)) | ||
+ | image(kr2) | ||
+ | image(kr2, col=gray(seq(1, 0, l=50))) | ||
+ | image(kr2, col=gray(seq(1, 0.2, l=50))) | ||
+ | points(s100, add=T) | ||
+ | image(kr2, col=gray(seq(0, 1, l=50))) | ||
+ | points(s100, add=T) | ||
+ | image(kr2, col=rev(gray(seq(0, 1, l=50)))) | ||
+ | image(kr2, col=terrain.colors(21)) | ||
+ | image(kr2, col=terrain.colors(21), val = sqrt(kr2$var)) | ||
+ | image(kr2, col=terrain.colors(21), val = sqrt(kr2$krige.var)) | ||
+ | points(s100, add=T) | ||
+ | </code> | ||