Benilton,<div><br></div><div>os teus resultados são iguais aos meus no R 64 bits. Por resultado antigo tu diz carregar alguma área de trabalho? Se sim, isso não esta acontecendo.</div><div><br></div><div>Leonard,</div><div>
<br></div><div>É um i5 3450 com windows 7, R 2.15.1 e os pacotes são as versões mais recentes. A outra maquina (que é só 32 bits) eu sei que é um Core 2 Quad com WinXP, R e pacotes desatualizados</div><div><br></div><div>
<br></div><div>[]'s<br><div><br><div class="gmail_quote">2012/7/17 Leonard de Assis <span dir="ltr"><<a href="mailto:assis.leonard@gmail.com" target="_blank">assis.leonard@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

  
    
  
  <div text="#000000" bgcolor="#FFFFFF">
    <div>Eu aqui tenho linux, rodei uma vm
      windows pra simular o 32 bits e comparei 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, máquina, etc<br>
      <pre cols="72">[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
</pre>
      Em 17/07/2012 12:18, Rodrigo Coster escreveu:<br>
    </div><div><div class="h5">
    <blockquote type="cite">Benilton,
      <div><br>
      </div>
      <div>olha, acho muito pouco provável (até pq os 2 são windows,
        então acho que nem existe essa opção). No computador aqui de
        casa os 2 R divergem no resultado (o gráfico que eu mandei foi
        rodando no R 32 e 64 bits no mesmo computador, Windows 7).</div>
      <div><br>
      </div>
      <div>Abri a imagem que tu mandou e rodei o comando, deu que os 2
        objetos são iguais.</div>
      <div><br>
      </div>
      <div><br>
      </div>
      <div>[]'s<br>
        <br>
        <div class="gmail_quote">2012/7/17 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">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" target="_blank">rcoster@gmail.com</a>>:<br>
            <div>
              <div>> Caros,<br>
                ><br>
                > programei uma rotina para estimar por máxima
                verossimilhança os parâmetros<br>
                > de uma cópula e para ver se estava certo comparei
                os resultados com o<br>
                > comando do pacote copula. Encontrei diferenças
                apenas na 5a casa decimal em<br>
                > diante, que considerei como sendo por causa do
                método numérico utilizado. Só<br>
                > que, ao fazer a mesma comparação num computador 64
                bits os resultados são<br>
                > bastante divergentes (o meu código muda o valor
                estimado, enquanto o 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 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>
                > 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>
              </div>
            </div>
            > _______________________________________________<br>
            > R-br mailing list<br>
            > <a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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" target="_blank">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>
          </blockquote>
        </div>
        <br>
      </div>
      <br>
      <fieldset></fieldset>
      <br>
      <pre>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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.</pre>
    </blockquote>
    <br>
    <br>
  </div></div></div>

<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></blockquote></div><br></div></div>