Prezados,<div>Como geoestatística me interessa tentei rodar o script do Elias, mas está ocorrendo o seguinte:</div><div>require(Rcitrus) ## pediu a versão 2.14 do R,</div><div><br></div><div>Após a instalação da versão 2.14.2 (como a 2.15.2), não carrega o geoR, dá a seguinte msg:</div>
<div><br></div><div><div><b>> require(geoR)</b></div><div><b>Carregando pacotes exigidos: geoR</b></div><div><b>Mensagens de aviso perdidas:</b></div><div><b>In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, :</b></div>
<div><b> there is no package called ‘geoR’</b></div><div><br></div><div><br></div><div>Alguma idéia?</div><div><br></div><div>Hélio</div><br><div class="gmail_quote">Em 1 de dezembro de 2012 17:22, Natalia Martins <span dir="ltr"><<a href="mailto:nsmbarreto@gmail.com" target="_blank">nsmbarreto@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Tentei seguir sua sugestão, a qual por sinal foi muito interessante,<div>no entanto eu gostaria de selecionar apenas as linhas, pois na pratica o problemaa seria:</div>
<div>um trator passa nas ruas(linhas) do talhao (tenho a coordenada de cada linha) e quero verificar se o trator pode pular ruas(linhas) - ou seja- se a linha 1 é = a linha 2, e assim sucessivamente a ssim ver o maximo de linhas q o trator poderia pular.<br>
<br><div class="gmail_quote">Em 1 de dezembro de 2012 15:02, Elias T. Krainski <span dir="ltr"><<a href="mailto:eliaskrainski@yahoo.com.br" target="_blank">eliaskrainski@yahoo.com.br</a>></span> escreveu:<div><div class="h5">
<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div style="font-size:12pt;font-family:times new roman,new york,times,serif"><div><span>Veja <br></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif">
<br><span></span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif"><span>sample.geodata() do pacote 'geoR' para facilitar amostragem de objetos 'geodata'<br>
</span></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif"><span>qcsamples() do pacote 'Rcitrus' (ver <a href="http://leg.ufpr.br/Rcitrus" target="_blank">http://leg.ufpr.br/Rcitrus</a> ) para fazer amostragem em lattice regular de forma aleatoria ou sistematica</span></div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif">
Mais sobre o que qcsamples() pode auxiliar a fazer:<br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif">Suponha um talhao com R linhas e C colunas. </div>
<div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif">comece com n pequeno<br> 1 - Selecione aleatóriamente n quadrats (ou pontos) desse talhão e calcule uma estatística. <br>
2 - Repita o procedimento anterior B vezes <br> 3 - guarde as B estatísticas, faça n = n+1 e volte ao passo 1. pare quando n atingir um valor razoavel
<br>Esse
procedimento termina quando é atingido o N.max ou quando a estatística
"estabiliza". A "estabilização" da estatística é verificada por critério
utilizado em Controle Estatístico de Qualidade. <br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif">
O processo
descrito acima é de seleção com reposição, que pode ser feito utilizando a função
qcsamples(). Esta função implementa duas formas de seleção de quadrats (ou pontos num grid regular) com reposição.</div><br>### exemplo de uso da funcao qcsamples() com dados simulados <br>nr <- 30;<br>nc <- 20;<br>
n <- nc*nr<br>dim(dados <- matrix(rnorm(n), nr, nc))<br><br>require(Rcitrus)<br>### amostras com quadrats 1x1<br>### sem reposicao e quadrats sem interseccao<br>### para selecao de pontos (nao de quadrats) usar dx=1 e dy=1<br>
amr <- qcsamples(dados, dx=1, dy=1, n=20, meth='random')<br>ams <- qcsamples(dados, dx=1, dy=1, nr=5, nc=4, meth='sistematic')<br><br>pr <- sum(amr[,3])/sum(amr[,4])<br>ps <- sum(ams[,3])/sum(ams[,4])<br>
<br>### incidencias estimadas<br>pr<br>ps<br>### verdadeira<br>mean(dados)<br><br>### visualizacao simples<br>image(list(x=1:nr, y=1:nc, z=dados), asp=1)<br>points(amr[,1:2], col=4, pch=19)<br>points(ams[,1:2], col=5, pch=19)<br>
legend("topleft", c("random", "sistematic"), pch=19, col=4:5,
ncol=2, bty='n')<br><br>### outra visualizacao (adequada p situacao de quadrats)<br>par(mfrow=c(1,2), mar=c(3,3,1,1), mgp=c(1.5,.5,0))<br>### vizualizando (aleatorio)<br>image(list(x=1:nr, y=1:nc, z=dados)) <br>
apply(amr, 1, function(x) {<br> mx <- cbind(c(x[1] + c(-.5,0.5,0.5,-.5,-.5)),<br> c(x[2] + c(-.5,-.5,0.5,0.5,-.5)))<br> polygon(mx)<br>})<br><br>### vizualizando (sistematico)<br>image(list(x=1:nr, y=1:nc, z=dados)) <br>
apply(ams, 1, function(x) {<br> mx <- cbind(c(x[1] + c(-.5,.5,.5,-.5,-.5)),<br> c(x[2] + c(-.5,-.5,.5,.5,-.5)))<br> polygon(mx)<br>})<br><br><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:times new roman,new york,times,serif">
Att.<br><span></span></div><div>Elias T.
Krainski<br><blockquote style="border-left:2px solid rgb(16,16,255);margin-left:5px;margin-top:5px;padding-left:5px"> <div style="font-family:times new roman,new york,times,serif;font-size:12pt"> <div style="font-family:times new roman,new york,times,serif;font-size:12pt">
<div dir="ltr"> <font face="Arial"> <hr size="1"> <b><span style="font-weight:bold">De:</span></b> Natalia Martins <<a href="mailto:nsmbarreto@gmail.com" target="_blank">nsmbarreto@gmail.com</a>><br> <b><span style="font-weight:bold">Para:</span></b> Lista R <<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>> <br>
<b><span style="font-weight:bold">Enviadas:</span></b> Sábado, 1 de Dezembro de 2012 9:57<br> <b><span style="font-weight:bold">Assunto:</span></b> [R-br] geoestatistica<br> </font> </div><div><div> <br><div>Prezados bom dia,<div>
estou cursando uma disciplina de amostragem e tenho como tema de seminario a aplicação de amostragem em
geoestatistica.</div><div>Para este estudo tenho as "linhas" que um trator passa em um talhao e o objetivo é verificar a cada quantas linhas o trator pode passar sem que haja prejuízo nas analises de solo.</div>
<div>Acredito ser uma amostragem sistematica na qual devo ir variando os K's e verificando até qual K os resultados sejam iguais quando considero todas as linhas.</div><div>No entanto no momento de fazer a análise geoestatistica na linha não estou conseguindo.</div>
<div>Gostaria de ajuda.</div><div><br></div><div>Estrutura dos dados</div><div><br></div><div><div>LONG<span style="white-space:pre-wrap"> </span><span style="white-space:pre-wrap"> </span>LAT<span style="white-space:pre-wrap"> </span> NDRE<span style="white-space:pre-wrap"> </span> NDVI<span style="white-space:pre-wrap"> </span> linha </div>
<div>-48.0149775<span style="white-space:pre-wrap"> </span>-21.3646713<span style="white-space:pre-wrap"> </span>0.171400000<span style="white-space:pre-wrap"> </span>0.435000000<span style="white-space:pre-wrap"> </span>1</div>
</div><div>.</div><div>.</div><div>.</div><div><br>-48.0149775<span style="white-space:pre-wrap"> </span>-21.3646713<span style="white-space:pre-wrap"> </span>0.171400000<span style="white-space:pre-wrap"> </span>0.435000000<span style="white-space:pre-wrap"> </span>63 <br clear="all">
<div><br></div>Abraços <br><br style="font-family:Tahoma;font-size:13px"><font face="'Times New Roman'" size="4">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: <a href="tel:%2819%29%208306-4743" value="+551983064743" target="_blank">(19) 8306-4743</a></font> <br><div><br></div>
</div>
</div><br></div></div>_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </blockquote></div> </div></div><br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div></div></div><br><br clear="all"><div><br>
</div>-- <br><div class="im">
<br style="font-family:Tahoma;font-size:13px"><font face="'Times New Roman'" size="4">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: (19) 8306-4743</font> <br><div><br></div><br>
</div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>