<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Veja <br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>sample.geodata() do pacote 'geoR' para facilitar amostragem de objetos 'geodata'<br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>qcsamples() do pacote 'Rcitrus' (ver http://leg.ufpr.br/Rcitrus ) para fazer amostragem em lattice regular de forma aleatoria ou sistematica</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family:
times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">Mais sobre o que qcsamples() pode auxiliar a fazer:<br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">Suponha um talhao com R linhas e C colunas. </div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">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="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">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="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">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" size="2"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Natalia Martins <nsmbarreto@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> Lista R <R-br@listas.c3sl.ufpr.br> <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> <br><meta http-equiv="x-dns-prefetch-control" content="off"><div id="yiv1739376102">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 class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span><span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>LAT<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span> NDRE<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span> NDVI<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span> linha </div>
<div>-48.0149775<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>-21.3646713<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>0.171400000<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>0.435000000<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>1</div>
</div><div>.</div><div>.</div><div>.</div><div><br>-48.0149775<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>-21.3646713<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>0.171400000<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </span>0.435000000<span class="yiv1739376102Apple-tab-span" style="white-space:pre;"> </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: (19) 8306-4743</font> <br><div><br></div>
</div>
</div><meta http-equiv="x-dns-prefetch-control" content="on"><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" 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><br> </div> </div> </blockquote></div> </div></body></html>