<div dir="ltr">Romero,<div><br></div><div>O princípio não seria "Com base numa amostra causal retirada de uma população, as estimativas dos parâmetros da função de distribuição da população são aqueles valores dos parâmetros que tornam máxima a probabilidade da <u>amostra</u> observada."? (grifo meu).</div><div><br></div><div>Como a estimativa é baseada em uma amostragem, por óbvio ela é afetada de um erro que geralmente controla-se via um IC.</div><div><br></div><div>Por outro lado como você começa de uma "amostra controlada" ou seja você acredita que já sabe o(s) parâmetro(s) da distribuição de onde amostra, então você tem dizer qual é  sua métrica para dizer se o parâmetro estimado por <span style="font-size:12.8px"> </span><span style="font-size:12.8px">log-verossimilhança está "próximo" ou "afastado"?</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-03-11 13:31 GMT-03:00 Romero Luiz M. Sales Filho <span dir="ltr"><<a href="mailto:romero.sfilho@gmail.com" target="_blank">romero.sfilho@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div><span style="font-size:12.8px">Oi Cesar,</span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">esse é o princípio da máxima verossimilhança... maximizando a função de log-verossimilhança (que é mais simples que a verossimilhança) o valor máximo é a estimativa do parâmetro estudado. Como estou simulando dados, ou seja, tenho uma amostra controlada, gerada computacionalmente e que sei a priori qual o valor do parâmetro que foi utilizado para gerar essa amostra, então espero que essa estimativa de máxima verossimilhança seja bastante próximas do valor do parâmetro. Na medida que o tamanho da amostra aumenta, é esperado que essa estimativa torne-se cada vez mais próxima desse valor do parâmetro devido a propriedade de consistência do estimador de máxima verossimilhança.</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Eu percebi que o código que eu estava utilizando estava incorreto... fiz algumas correções e os resultados melhoraram significativamente. A seguir o código correto:</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><div><span class=""><span><div>### Gerando a amostra baseando-se nos parâmetros iniciais</div><div>### Considerando um parâmetro q positivo.</div><div>q=1.5</div><div>eta=5</div><div>n=100</div></span></span><span class=""><span><div>u<-runif(n)<span style="white-space:pre-wrap">        </span></div><div>y = eta*((1-(u)^((1-q)/(2-q)))/(1-q))</div><div>### Obtendo as estimativas dos parâmetros a partir da amostra gerada anteriormente</div><div>vero <- function(par,x){</div><div>q = par[1]</div><div>eta = par[2]</div><div>   saida<-((sum(log(1 - ((1 - q)*x*(1/eta)))))/(1 - q)) + ((n)*log(2 - q)) + (n*log(1/eta))</div><div>return(-saida)</div><div>}</div><div>saiday<-optim(par=c(1.8,4),fn=vero,  method= "Nelder-Mead",x=y  )</div><div>q_chap<-saiday[1]$par[1]</div><div>eta_chap<-saiday[1]$par[2]</div></span></span><span class=""><span>### Construindo os gráficos de log-verossimilhança</span></span><div>par.vals <- expand.grid(q=seq(1.1, 1.9, l = 25), eta=seq(2, 10, l = 25)) </div><div>dim(par.vals)</div><div>f <- function(pars,dados) { ((sum(log(1 - ((1 - pars[1])*dados*(1/pars[2])))))/(1 - pars[1])) + (n*log(2 - pars[1])) + (n*log(1/pars[2]))}</div><div>par.vals$L <- apply(par.vals, 1, f, dados = y) </div><div><br></div><div>with(par.vals, persp(unique(q), unique(eta), matrix(L, ncol = length(unique(eta))), xlab = expression(q), ylab = expression(eta), zlab = expression(l(q, eta)), theta= 30, phi = 30)) </div><div>z<-as.matrix(par.vals)</div><div>z1 <- z[order(z[,3],decreasing=T),]</div><div>result<-c(z1[1,1],z1[1,2],z1[1,3])</div><div>###Valores máximos para o gráfico da função de log-verossimilhança:</div><div>result</div><div>###Estimativas de Máxima verossimilhança:</div><div>q_chap</div><div>eta_chap</div></div><div><br></div><div>Apesar disso, o caso em que o q é negativo ainda não me da bons resultados ... alguém sabe o que pode está acontecendo, segue o código de um exemplo com o q negativo:<br></div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><span class=""><span><div>### Gerando a amostra baseando-se nos parâmetros iniciais</div><div>### Considerando um parâmetro q positivo.</div></span></span><div>q=-1.5</div><div>eta=5</div><div>n=100</div><span class=""><span><div>u<-runif(n)<span style="white-space:pre-wrap">     </span></div><div>y = eta*((1-(u)^((1-q)/(2-q)))/(1-q))</div><div>### Obtendo as estimativas dos parâmetros a partir da amostra gerada anteriormente</div><div>vero <- function(par,x){</div><div>q = par[1]</div><div>eta = par[2]</div><div>   saida<-((sum(log(1 - ((1 - q)*x*(1/eta)))))/(1 - q)) + ((n)*log(2 - q)) + (n*log(1/eta))</div><div>return(-saida)</div><div>}</div><div>saiday<-optim(par=c(1.8,4),fn=vero,  method= "Nelder-Mead",x=y  )</div><div>q_chap<-saiday[1]$par[1]</div><div>eta_chap<-saiday[1]$par[2]</div></span></span><span class=""><span>### Construindo os gráficos de log-verossimilhança</span></span><div>par.vals <- expand.grid(q=seq(-2, 0, l = 25), eta=seq(2, 10, l = 25)) </div><div>dim(par.vals)</div><div>f <- function(pars,dados) { ((sum(log(1 - ((1 - pars[1])*dados*(1/pars[2])))))/(1 - pars[1])) + (n*log(2 - pars[1])) + (n*log(1/pars[2]))}</div><div>par.vals$L <- apply(par.vals, 1, f, dados = y) </div><div><br></div><div>with(par.vals, persp(unique(q), unique(eta), matrix(L, ncol = length(unique(eta))), xlab = expression(q), ylab = expression(eta), zlab = expression(l(q, eta)), theta= 30, phi = 30)) </div><div>z<-as.matrix(par.vals)</div><div>z1 <- z[order(z[,3],decreasing=T),]</div><div>result<-c(z1[1,1],z1[1,2],z1[1,3])</div><div>###Valores máximos para o gráfico da função de log-verossimilhança:</div><div>result</div><div>###Estimativas de Máxima verossimilhança:</div><div>q_chap</div><div>eta_chap</div></div></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" rel="noreferrer" 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" rel="noreferrer" 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>