Benilton,<div><br></div><div>mesmo com --vanilla deu resultados diferentes.</div><div><font face="courier new, monospace"><br></font></div><div><div><font face="courier new, monospace">> sessionInfo()</font></div><div><font face="courier new, monospace">R version 2.15.1 (2012-06-22)</font></div>
<div><font face="courier new, monospace">Platform: x86_64-pc-mingw32/x64 (64-bit)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">locale:</font></div><div><font face="courier new, monospace">[1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C LC_TIME=Portuguese_Brazil.1252 </font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">attached base packages:</font></div><div><font face="courier new, monospace">[1] stats graphics grDevices utils datasets methods base </font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">other attached packages:</font></div><div><font face="courier new, monospace">[1] MASS_7.3-19 copula_0.99-4 pspline_1.0-14</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">loaded via a namespace (and not attached):</font></div><div><font face="courier new, monospace">[1] ADGofTest_0.3 gsl_1.9-9 mvtnorm_0.9-9992 stabledist_0.6-4 stats4_2.15.1 </font></div>
<div><font face="courier new, monospace">> a</font></div><div><font face="courier new, monospace"> [,1] [,2]</font></div><div><font face="courier new, monospace"> [1,] 0.6533333 0.6673587</font></div><div>
<font face="courier new, monospace"> [2,] 0.6956496 0.7000601</font></div><div><font face="courier new, monospace"> [3,] 0.5554635 0.6169322</font></div><div><font face="courier new, monospace"> [4,] 0.4895890 0.5405670</font></div>
<div><font face="courier new, monospace"> [5,] 0.6038044 0.6422919</font></div><div><font face="courier new, monospace"> [6,] 0.7295427 0.7295427</font></div><div><font face="courier new, monospace"> [7,] 0.6548815 0.7080392</font></div>
<div><font face="courier new, monospace"> [8,] 0.5468827 0.6015380</font></div><div><font face="courier new, monospace"> [9,] 0.6630423 0.6955567</font></div><div><font face="courier new, monospace">[10,] 0.7214445 0.7214441</font></div>
<div><font face="courier new, monospace">[11,] 0.6681827 0.6681817</font></div><div><font face="courier new, monospace">[12,] 0.6429204 0.6429205</font></div><div><font face="courier new, monospace">[13,] 0.5678848 0.5918048</font></div>
<div><font face="courier new, monospace">[14,] 0.6829043 0.6985720</font></div><div><font face="courier new, monospace">[15,] 0.7653632 0.7653638</font></div><div><font face="courier new, monospace">[16,] 0.6215542 0.6570017</font></div>
<div><font face="courier new, monospace">[17,] 0.7085038 0.7085053</font></div><div><font face="courier new, monospace">[18,] 0.6620582 0.6958966</font></div><div><font face="courier new, monospace">[19,] 0.6829283 0.6912327</font></div>
<div><font face="courier new, monospace">[20,] 0.6416442 0.6904770</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"><br></font></div><div><div><font face="courier new, monospace">> sessionInfo()</font></div>
<div><font face="courier new, monospace">R version 2.15.1 (2012-06-22)</font></div><div><font face="courier new, monospace">Platform: i386-pc-mingw32/i386 (32-bit)</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">locale:</font></div><div><font face="courier new, monospace">[1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C LC_TIME=Portuguese_Brazil.1252 </font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">attached base packages:</font></div><div><font face="courier new, monospace">[1] stats graphics grDevices utils datasets methods base </font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">other attached packages:</font></div><div><font face="courier new, monospace">[1] MASS_7.3-19 copula_0.99-4 pspline_1.0-14</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">loaded via a namespace (and not attached):</font></div><div><font face="courier new, monospace">[1] ADGofTest_0.3 gsl_1.9-9 mvtnorm_0.9-9992 stabledist_0.6-4 stats4_2.15.1 </font></div>
<div><font face="courier new, monospace">> a</font></div><div><font face="courier new, monospace"> [,1] [,2]</font></div><div><font face="courier new, monospace"> [1,] 0.6673587 0.6673587</font></div><div>
<font face="courier new, monospace"> [2,] 0.7000619 0.7000601</font></div><div><font face="courier new, monospace"> [3,] 0.6169316 0.6169322</font></div><div><font face="courier new, monospace"> [4,] 0.5405662 0.5405670</font></div>
<div><font face="courier new, monospace"> [5,] 0.6423069 0.6422919</font></div><div><font face="courier new, monospace"> [6,] 0.7295427 0.7295427</font></div><div><font face="courier new, monospace"> [7,] 0.7080409 0.7080392</font></div>
<div><font face="courier new, monospace"> [8,] 0.6015378 0.6015380</font></div><div><font face="courier new, monospace"> [9,] 0.6955575 0.6955567</font></div><div><font face="courier new, monospace">[10,] 0.7214445 0.7214441</font></div>
<div><font face="courier new, monospace">[11,] 0.6681827 0.6681817</font></div><div><font face="courier new, monospace">[12,] 0.6429204 0.6429205</font></div><div><font face="courier new, monospace">[13,] 0.5918050 0.5918048</font></div>
<div><font face="courier new, monospace">[14,] 0.6985704 0.6985720</font></div><div><font face="courier new, monospace">[15,] 0.7653632 0.7653638</font></div><div><font face="courier new, monospace">[16,] 0.6569979 0.6570017</font></div>
<div><font face="courier new, monospace">[17,] 0.7085038 0.7085053</font></div><div><font face="courier new, monospace">[18,] 0.6958953 0.6958966</font></div><div><font face="courier new, monospace">[19,] 0.6912321 0.6912327</font></div>
<div><font face="courier new, monospace">[20,] 0.6904751 0.6904770</font></div></div><div><br></div><br><div class="gmail_quote">2012/7/18 Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com" target="_blank">beniltoncarvalho@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Voce consegue executar o R passando argumentos adicionais? (imagino q<br>
vc tenha q fucar na linha de comando dentro do atalho usado pelo<br>
Windows)...<br>
<br>
Se vc estivesse num Linux/Mac, eu recomendaria executar o R como segue:<br>
<br>
R --arch i386 --vanilla<br>
<br>
e<br>
<br>
R --arch x86_64 --vanilla<br>
<br>
E tentar mais uma ultima vez a comparacao entre ambos... (Evite<br>
qualquer firula como RStudio tbm).<br>
<div class="HOEnZb"><div class="h5"><br>
b<br>
<br>
2012/7/18 Rodrigo Coster <<a href="mailto:rcoster@gmail.com">rcoster@gmail.com</a>>:<br>
> Benilton,<br>
><br>
> rodando na mesma maquina o resultado do R32 e R64 é diferente, foi com base<br>
> nisso (R's da mesma maquina) que eu fiz o gráfico que mandei... tentei<br>
> desinstalar o R e reinstalar, mas não resolveu<br>
><br>
><br>
> []'s<br>
><br>
><br>
> 2012/7/18 Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><br>
>><br>
>> Eu comecaria por garantir que as versoes do R e dos demais pacotes<br>
>> sejam as mesmas nos diferentes ambientes.<br>
>><br>
>> Em teoria, a maquina que tem o R 64 bits tambem tem o 32 bits (algum<br>
>> usuario de windows pode te ajudar melhor q eu nesse caso), vc poderia<br>
>> tentar executar o R32 na mesma maquina...<br>
>><br>
>> b<br>
>><br>
>> 2012/7/18 Rodrigo Coster <<a href="mailto:rcoster@gmail.com">rcoster@gmail.com</a>>:<br>
>> > Benilton,<br>
>> ><br>
>> > os teus resultados são iguais aos meus no R 64 bits. Por resultado<br>
>> > antigo tu<br>
>> > diz carregar alguma área de trabalho? Se sim, isso não esta acontecendo.<br>
>> ><br>
>> > Leonard,<br>
>> ><br>
>> > É um i5 3450 com windows 7, R 2.15.1 e os pacotes são as versões mais<br>
>> > recentes. A outra maquina (que é só 32 bits) eu sei que é um Core 2 Quad<br>
>> > com<br>
>> > WinXP, R e pacotes desatualizados<br>
>> ><br>
>> ><br>
>> > []'s<br>
>> ><br>
>> ><br>
>> > 2012/7/17 Leonard de Assis <<a href="mailto:assis.leonard@gmail.com">assis.leonard@gmail.com</a>><br>
>> >><br>
>> >> Eu aqui tenho linux, rodei uma vm windows pra simular o 32 bits e<br>
>> >> comparei<br>
>> >> com o meu linux (64 bits) deu que os resultados não diferem<br>
>> >><br>
>> >> Seria bom apresentar maiores detalhes, como versões do R, do pacote,<br>
>> >> máquina, etc<br>
>> >><br>
>> >> []s<br>
>> >> Leonard de Assis<br>
>> >> assis <dot> leonard <at> gmail <dot> com<br>
>> >><br>
>> >> Em 17/07/2012 12:18, Rodrigo Coster escreveu:<br>
>> >><br>
>> >> Benilton,<br>
>> >><br>
>> >> olha, acho muito pouco provável (até pq os 2 são windows, então acho<br>
>> >> que<br>
>> >> nem existe essa opção). No computador aqui de casa os 2 R divergem no<br>
>> >> resultado (o gráfico que eu mandei foi rodando no R 32 e 64 bits no<br>
>> >> mesmo<br>
>> >> computador, Windows 7).<br>
>> >><br>
>> >> Abri a imagem que tu mandou e rodei o comando, deu que os 2 objetos são<br>
>> >> iguais.<br>
>> >><br>
>> >><br>
>> >> []'s<br>
>> >><br>
>> >> 2012/7/17 Benilton Carvalho <<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>><br>
>> >>><br>
>> >>> Alguma chance de algum dos computadores q vc tenha usado estar<br>
>> >>> utilizando alguma "compilacao propria" do R ou algum esquema de<br>
>> >>> otimizacao?<br>
>> >>><br>
>> >>> Executei o seu codigo no meu computador e tudo parece OK (ie., o<br>
>> >>> objeto 'a' eh identico entre arquiteturas)... O arquivo a seguir<br>
>> >>> possui 2 elements: a32 e a64, resultantes do codigo q vc enviou.<br>
>> >>><br>
>> >>> <a href="https://www.dropbox.com/s/srn0bhtmen4fd6u/R64x32.Rda" target="_blank">https://www.dropbox.com/s/srn0bhtmen4fd6u/R64x32.Rda</a><br>
>> >>><br>
>> >>> Experimente o comando a seguir ao carregar o arquivo:<br>
>> >>><br>
>> >>> all.equal(a32, a64)<br>
>> >>><br>
>> >>> A unica coisa q consigo imaginar e' se alguma "otimizacao" ao compilar<br>
>> >>> o R e qq outro acessorio utilizado nao tiver dado certo.<br>
>> >>><br>
>> >>><br>
>> >>> b<br>
>> >>><br>
>> >>><br>
>> >>> 2012/7/17 Rodrigo Coster <<a href="mailto:rcoster@gmail.com">rcoster@gmail.com</a>>:<br>
>> >>> > Caros,<br>
>> >>> ><br>
>> >>> > programei uma rotina para estimar por máxima verossimilhança os<br>
>> >>> > parâmetros<br>
>> >>> > de uma cópula e para ver se estava certo comparei os resultados com<br>
>> >>> > o<br>
>> >>> > comando do pacote copula. Encontrei diferenças apenas na 5a casa<br>
>> >>> > decimal em<br>
>> >>> > diante, que considerei como sendo por causa do método numérico<br>
>> >>> > utilizado. Só<br>
>> >>> > que, ao fazer a mesma comparação num computador 64 bits os<br>
>> >>> > resultados<br>
>> >>> > são<br>
>> >>> > bastante divergentes (o meu código muda o valor estimado, enquanto o<br>
>> >>> > pacote<br>
>> >>> > copula mantem), mudando na 2a casa decimal (no caso o parâmetro é a<br>
>> >>> > correlação, então a 2a casa decimal é bem importante). Dai me bateu<br>
>> >>> > a<br>
>> >>> > seguinte dúvida: em qual confiar?<br>
>> >>> ><br>
>> >>> > Código:<br>
>> >>> > require(copula)<br>
>> >>> > require(MASS)<br>
>> >>> > set.seed(31415)<br>
>> >>> ><br>
>> >>> > normCop <- function(param,data) {<br>
>> >>> > n <- nrow(data)<br>
>> >>> > if (length(param) != n) { param <- rep(param,nrow(data)) }<br>
>> >>> > cop <- mapply(normalCopula,param=param,MoreArgs=list(dim=2))<br>
>> >>> > datalist <- apply(data,1,list)<br>
>> >>> > for (i in 1:n) { datalist[[i]] <- datalist[[i]][[1]] }<br>
>> >>> > out <- -sum(log(mapply(dcopula,copula=cop,u=datalist)))<br>
>> >>> > if (out == Inf) { out = exp(100) }<br>
>> >>> > return(out)<br>
>> >>> > }<br>
>> >>> > a <- matrix(0,20,2)<br>
>> >>> > n <- 100<br>
>> >>> > Sigma <- matrix(c(10,3,3,2),2,2)<br>
>> >>> > for (j in 1:20) {<br>
>> >>> > data <- mvrnorm(n=n, rep(0, 2), Sigma)<br>
>> >>> > data <- apply(data,2,rank)/(n+1)<br>
>> >>> ><br>
>> >>> > fitNormCop <- function(data) {<br>
>> >>> > optim(cor(data)[2],normCop,data=data, lower = 0, upper =<br>
>> >>> > .9999,method="L-BFGS-B")<br>
>> >>> > }<br>
>> >>> > a[j,1] <- fitNormCop(data)$par # Meu<br>
>> >>> > a[j,2] <- fitCopula(normalCopula(.2,2), data, method="ml")@estimate<br>
>> >>> > #<br>
>> >>> > Pacote copula<br>
>> >>> > }<br>
>> >>> ><br>
>> >>> > E aqui um gráfico de dispersão comparando todos:<br>
>> >>> > <a href="http://img411.imageshack.us/img411/7527/92588280.png" target="_blank">http://img411.imageshack.us/img411/7527/92588280.png</a><br>
>> >>> ><br>
>> >>> ><br>
>> >>> ><br>
>> >>> ><br>
>> >>> ><br>
>> >>> ><br>
>> >>> > []'s<br>
>> >>> ><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<br>
>> >>> > código<br>
>> >>> > mínimo reproduzível.<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<br>
>> >>> código mínimo reproduzível.<br>
>> >><br>
>> >><br>
>> >><br>
>> >><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<br>
>> >> código mínimo reproduzível.<br>
>> >><br>
>> >><br>
>> >><br>
>> >><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<br>
>> >> código mínimo reproduzível.<br>
>> ><br>
>> ><br>
>> ><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<br>
>> > código<br>
>> > mínimo reproduzível.<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<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><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<br>
> mínimo reproduzível.<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>
</div></div></blockquote></div><br></div>