<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">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 class="moz-signature" 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>
<blockquote
cite="mid:CAKU4wovZN=dbbcHYeY44+X6fWawbZfQ5FEtQjiTfVr6v+1Zcnw@mail.gmail.com"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:rcoster@gmail.com">rcoster@gmail.com</a>>:<br>
<div>
<div class="h5">> 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 moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a moz-do-not-send="true"
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 moz-do-not-send="true"
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 moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
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 moz-do-not-send="true"
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 class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<br>
</body>
</html>