URGENTE - Função de maximização (maxLik)

Ver o tópico anterior Ver o tópico seguinte Ir em baixo

URGENTE - Função de maximização (maxLik)

Mensagem  Orleans Martins em Qui Jan 12, 2012 11:15 am

Prezados,

Nesses últimos dias tenho tentado entender essa função mas ainda me restam algumas dúvidas. Sou leigo em R, estou utilizando porque preciso fazer esta maximização para o modelo que vou utilizar na minha tese e não encontrei outro software que faça esse cálculo.

A questão é a seguinte: preciso calcular uma probabilidade a partir de duas variáveis binárias (B e S). Cada uma contém 50 valores (simulados, apenas para calcular os parâmetros do modelo).

Eis a rotina que estou utilizando:

---INICIAR:

library(maxLik)


---DADOS PARA TESTE (50 observações)

B <- c(1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,1,0,1,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,1,0,1,1,1,1,1,0,0,0,0,0,0,0,0,0)
S <- c(0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0,1,0,1,1,1,0,1,1,1,0,1,1,1,1,1,1,1,0,1,0,0,0,0,0,1,1,1,1,1,1,1,1,1)

---EQUAÇÃO DIGITADA:

loglik <- function(param) {
alpha <- param[1]
mu <- param[2]
Eb <- param[3]
Es <- param[4]
delta <- param[5]
ll <- (1-alpha)*exp(-Eb)*((Eb^B)/factorial(B))*exp(-Es)*((Es^S)/factorial(S)) + alpha*delta*exp(-Eb)*((Eb^B)/factorial(B))*exp(-(mu+Es))*((mu+Es)^S/factorial(S)) + alpha*(1-delta)*exp(-(mu+Eb))*((mu+Eb)^B/factorial(B))*exp(-Es)*((Es^S)/factorial(S))
}
res <- maxLik(loglik, grad=NULL, hess=NULL, start=c(1,1,1,1,1), constraints=NULL)
print(res)
coef(res)


DÚVIDAS:

1) Observação: alpha, mu e delta são probabilidades (deveriam estar no intervalo entre 0 e 1?) Eb e Es são contínuos. Um colega que conhece um pouco mais de R disse que eles seriam "irrestritos" e as probabilidades "restritos". Não estou sabendo limitar isso.
2) start=c devem iniciar de 0 ou de 1? Quando coloco tudo 0 o software passa muito tempo calculando e não externa um resultado. Com 1 os parâmetros são calculados mas todos na casa de centena (por exemplo, alpha = 687.2159). Não estou sabendo analisar.
3) Sugeriram-me utilizar maxBFGS usando um espaço paramétrico. Não encontrei como.

Se puderem me ajudar agradecerei muito! Meu tempo está acabando e não estou encontrando solução.

Aceito qualquer negócio!

Orleans Martins

Mensagens : 3
Data de inscrição : 31/08/2011
Localização : João Pessoa

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Re: URGENTE - Função de maximização (maxLik)

Mensagem  gustavoreis em Qui Jan 12, 2012 10:01 pm

Caro Orleans, estou enviando um e-mail agora, confira sua caixa de entrada do yahoo. (meu e-mail é da hotmail)

att,

Gustavo Mello Reis.

gustavoreis
Admin

Mensagens : 70
Data de inscrição : 01/03/2011
Idade : 32

Ver perfil do usuário http://estatisticanor.forumeiros.com

Voltar ao Topo Ir em baixo

Ver o tópico anterior Ver o tópico seguinte Voltar ao Topo

- Tópicos similares

 
Permissão deste fórum:
Você não pode responder aos tópicos neste fórum