<div><div>Fernando e colegas, boa noite!</div><div><br></div><div>Para um conjunto de dados nominado 'dados', com 4 colunas:</div><div><br></div><div>> head(dados)</div><div>        X1      Y1    MO   ALT</div><div>
1 726023.6 7534573 29.50 387.1</div><div>2 726042.6 7534567 31.60 386.7</div><div>3 726014.5 7534544 33.71 387.1</div><div>4 726033.6 7534538 34.06 386.6</div><div>5 726005.5 7534516 32.65 387.1</div><div>6 726024.6 7534510 34.76 387.3</div>
<div><br></div><div>Eu tenho definido a variável que supostamente promove a tendência como uma covariável dentro do objeto geodata.</div><div>Neste exemplo dG é o objeto geodata, colunas 1 e 2 são as coords, coluna 3 ('MO') como data e coluna 4 ('ALT') como covar:</div>
<div><br></div><div>dG <- as.geodata(dados, coords.col=c(1,2), data.col=3, covar.col=4)<br></div><div><br></div><div>Definida como uma covar, a variável pode ser acessada pelo nome interno que recebe dentro do geodata, neste caso 'ALT'</div>
<div><br></div><div>variog(dG, max.dist=350, trend=~coords)                       ### é a mesma tendência dada por  '1st', usando o nome interno coords</div><div>variog(dG, max.dist=350, trend=~ALT)                           ### tendência promovida por 'ALT' (coluna 4)</div>
<div>variog(dG, max.dist=350, trend=~coords+ALT)               ### combinação de efeitos </div><div>variog(dG, max.dist=350, trend=~coords[,2])                  ### só ordenadas (latitude)</div><div>variog(dG, max.dist=350, trend=~I(coords[,2]^2) +ALT    ### quadrado da latitude + ALT<br>
</div><div>* OBS: não é necessário definir os argumentos coords e data para variog(), pois já são passados diretamente pelo objeto 'dG'.</div><div><br></div><div>Se  você não definir a variável como covar dentro do seu objeto geodata, também dá pra acessar diretamente a partir</div>
<div>do data.frame de dados...</div><div><br></div><div>variog(dG, max.dist=350, trend=~dados$ALT)</div></div><div><br></div><div>Definido o modelo , a interpolação é feita pela função krige.conv() e para validar o procedimento não é jackknifing, mas sim a validação cruzada dada por xvalid().</div>
<div><br></div><div>Sugiro acessar o tutorial:</div><div><a href="http://www.leg.ufpr.br/geoR/geoRdoc/geoRintro.html">http://www.leg.ufpr.br/geoR/geoRdoc/geoRintro.html</a><br></div><div><br></div>Espero que ajude...<br clear="all">
<font face="arial, helvetica, sans-serif"><br><span style="font-size:small">================================================</span><br>Éder Comunello<br><br>Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ)<br>
Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]<br><br>Researcher at Embrapa Western Region Agriculture</font><div><font face="arial, helvetica, sans-serif">Dourados, MS, Brazil [</font>22 16.5'S, 54 49.0'W<span style="font-family:arial,helvetica,sans-serif">]</span></div>
<div><font face="arial, helvetica, sans-serif"><span style="font-size:small">================================================</span><br>UTC-03:00</font></div><br>