O "script" ao que eu me refiro nada mais é do que uma sequencia de comandos, que pode ser visto a seguir:<div><br></div><div><div>require(geoComp)</div><div>require(MASS)</div><div>require(statmod)</div><div>require(geoR)</div>
<div>data(pivo)</div><div>dados <- pivo[,c(6,7,8,1,2)]</div><div>dados <- as.geoComp(dados)</div><div>bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) </div><div>estima <- mec(dados)</div>
<div>gr <- pred_grid(bor, by=4)</div><div><a href="http://md.cov.ck">md.cov.ck</a> <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados)</div><div><a href="http://preditos.gh">preditos.gh</a> <- volta.quad(<a href="http://md.cov.ck">md.cov.ck</a>,n.pontos=7,Variancia=FALSE)</div>
<div><a href="http://preditos.gh">preditos.gh</a> <- data.frame(<a href="http://preditos.gh">preditos.gh</a>)</div><div>write.table(<a href="http://preditos.gh">preditos.gh</a>,"pred_by4k7ns1000MBM.txt")</div>
<div>system.time(preditos.simu <- volta.cokri(<a href="http://md.cov.ck">md.cov.ck</a>,num.simu=100,int.conf=0.95))</div><div><br></div><div><br></div><div>Informações sobre o computador e sessionInfo:</div><div><br></div>
<div>Release: Ubuntu 11.10 (oneiric)</div><div>Kernel: 3.0.0-19-generic (#33-Ubuntu SMP Thu Apr 19 19:05:14 UTC 2012)</div><div>OS Type: Linux</div><div>GCC version: 4.6.1 (x86_64-linux-gnu)</div><div><br></div><div>Vendor: GenuineIntel</div>
<div>CPUs: 8</div><div>Model name: Intel(R) Core(TM) i7 CPU       Q 720  @ 1.60GHz</div><div>Frequency: 933.000 MHz</div><div>L2 cache: 6144 KB</div><div><br></div><div><br></div><div><div>R version 2.15.0 (2012-03-30)</div>
<div>Platform: x86_64-unknown-linux-gnu (64-bit)</div><div><br></div><div>locale:</div><div> [1] LC_CTYPE=pt_BR.UTF-8       LC_NUMERIC=C              </div><div> [3] LC_TIME=pt_BR.UTF-8        LC_COLLATE=pt_BR.UTF-8    </div>
<div> [5] LC_MONETARY=pt_BR.UTF-8    LC_MESSAGES=pt_BR.UTF-8   </div><div> [7] LC_PAPER=C                 LC_NAME=C                 </div><div> [9] LC_ADDRESS=C               LC_TELEPHONE=C            </div><div>[11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C       </div>
<div><br></div><div>attached base packages:</div><div>[1] stats     graphics  grDevices utils     datasets  methods   base     </div><div><br></div><div>other attached packages:</div><div>[1] geoR_1.7-2     sp_0.9-99      statmod_1.4.14 MASS_7.3-17    geoComp_1.0   </div>
<div><br></div><div>loaded via a namespace (and not attached):</div><div>[1] grid_2.15.0         lattice_0.20-6      RandomFields_2.0.54</div><div>[4] splancs_2.01-31     tools_2.15.0</div></div><div><br></div><div><br></div>
<div>Att,</div><div><br></div>-- <br>Junior Beleti
</div>