Re: [R-br] Alteração em função de pacote (novo tópico)
O "script" ao que eu me refiro nada mais é do que uma sequencia de comandos, que pode ser visto a seguir: require(geoComp) require(MASS) require(statmod) require(geoR) data(pivo) dados <- pivo[,c(6,7,8,1,2)] dados <- as.geoComp(dados) bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) estima <- mec(dados) gr <- pred_grid(bor, by=4) md.cov.ck <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados) preditos.gh <- volta.quad(md.cov.ck,n.pontos=7,Variancia=FALSE) preditos.gh <- data.frame(preditos.gh) write.table(preditos.gh,"pred_by4k7ns1000MBM.txt") system.time(preditos.simu <- volta.cokri(md.cov.ck ,num.simu=100,int.conf=0.95)) Informações sobre o computador e sessionInfo: Release: Ubuntu 11.10 (oneiric) Kernel: 3.0.0-19-generic (#33-Ubuntu SMP Thu Apr 19 19:05:14 UTC 2012) OS Type: Linux GCC version: 4.6.1 (x86_64-linux-gnu) Vendor: GenuineIntel CPUs: 8 Model name: Intel(R) Core(TM) i7 CPU Q 720 @ 1.60GHz Frequency: 933.000 MHz L2 cache: 6144 KB R version 2.15.0 (2012-03-30) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] geoR_1.7-2 sp_0.9-99 statmod_1.4.14 MASS_7.3-17 geoComp_1.0 loaded via a namespace (and not attached): [1] grid_2.15.0 lattice_0.20-6 RandomFields_2.0.54 [4] splancs_2.01-31 tools_2.15.0 Att, -- Junior Beleti
O alinhamento dos planetas nao esta' contribuindo para a reproducibilidade dos resultados que postei (e' bem possivel que eu tenha usado um conjunto de dados menor e, por descuido, postedo os resultados dele - desculpas pela falha)... Quando tiver uma folga aqui, volto a olhar os detalhes... de qualquer forma, o codigo q postei e' para ser executado numa fracao do tempo usado na execucao linear... Por exemplo, se vc leva 17h para execucao, minha expectativa e' de q vc precise de cerca de 2h ao usar os recursos de paralelismo e evitar uso desnecessario de recursos gerais. b On 12 June 2012 06:28, Junior Beleti <beleti.junior@gmail.com> wrote:
O "script" ao que eu me refiro nada mais é do que uma sequencia de comandos, que pode ser visto a seguir:
require(geoComp) require(MASS) require(statmod) require(geoR) data(pivo) dados <- pivo[,c(6,7,8,1,2)] dados <- as.geoComp(dados) bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) estima <- mec(dados) gr <- pred_grid(bor, by=4) md.cov.ck <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados) preditos.gh <- volta.quad(md.cov.ck,n.pontos=7,Variancia=FALSE) preditos.gh <- data.frame(preditos.gh) write.table(preditos.gh,"pred_by4k7ns1000MBM.txt") system.time(preditos.simu <- volta.cokri(md.cov.ck,num.simu=100,int.conf=0.95))
Informações sobre o computador e sessionInfo:
Release: Ubuntu 11.10 (oneiric) Kernel: 3.0.0-19-generic (#33-Ubuntu SMP Thu Apr 19 19:05:14 UTC 2012) OS Type: Linux GCC version: 4.6.1 (x86_64-linux-gnu)
Vendor: GenuineIntel CPUs: 8 Model name: Intel(R) Core(TM) i7 CPU Q 720 @ 1.60GHz Frequency: 933.000 MHz L2 cache: 6144 KB
R version 2.15.0 (2012-03-30) Platform: x86_64-unknown-linux-gnu (64-bit)
locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C
attached base packages: [1] stats graphics grDevices utils datasets methods base
other attached packages: [1] geoR_1.7-2 sp_0.9-99 statmod_1.4.14 MASS_7.3-17 geoComp_1.0
loaded via a namespace (and not attached): [1] grid_2.15.0 lattice_0.20-6 RandomFields_2.0.54 [4] splancs_2.01-31 tools_2.15.0
Att,
-- Junior Beleti
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (2)
-
Benilton Carvalho -
Junior Beleti