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/04/15 18:48] paulojus |
disciplinas:ce097-2015-01:historico2015 [2015/04/15 19:25] paulojus |
||
---|---|---|---|
Linha 74: | Linha 74: | ||
- Obtenha um mapa do 1o quartil, um do 3o quantil e dos percentis 10 e 90. | - Obtenha um mapa do 1o quartil, um do 3o quantil e dos percentis 10 e 90. | ||
- Como voce faria para obter a predição da proporção da área acima de um determinado limiar (por exemplo 1,2). | - Como voce faria para obter a predição da proporção da área acima de um determinado limiar (por exemplo 1,2). | ||
+ | * Explore um segundo exemplo de análises | ||
+ | - {{:disciplinas:ce097-2015-01:wolf.r|Comandos do R}} | ||
+ | - {{:disciplinas:ce097-2015-01:wolfcamp.txt|Arquivo de dados}} | ||
+ | * Mais um exemplo:<code R> | ||
+ | class(Ksat) | ||
+ | names(Ksat) | ||
+ | summary(Ksat) | ||
+ | plot(Ksat) | ||
+ | plot(Ksat, lam=0) | ||
+ | |||
+ | par(mfrow=c(1,1)) | ||
+ | K.v <- variog(Ksat, max.dist=9, lam=0) | ||
+ | plot(K.v) | ||
+ | |||
+ | var(log(Ksat$dat)) | ||
+ | |||
+ | K.eye <- eyefit(K.v) | ||
+ | K.eye | ||
+ | |||
+ | gr <- expand.grid(seq(0, 22.5, by=0.2), seq(0,13, by=0.2)) | ||
+ | kr <- krige.conv(Ksat, loc=gr, krige=krige.control(cov.pars=c(2.3,2.35), lam=0)) | ||
+ | image(kr, col=gray(seq(1,0,l=21))) | ||
+ | image(kr, col=terrain.colors(16)) | ||
+ | </code> | ||
+ | Como ficariam perguntas semelhantes às feitas para o 1o exemplo neste último caso? | ||
+ | Por exemplo: mapa da P[Y > 1] , mapa de quantis etc | ||