====== Reunião 23 ======

=== PI, 14 de outubro de 2016 ===

**Pautas: dúvidas sobre o trabalho de Modelagem Estatística**

Teóricas
  * Dúvida fórum 
      * "A questão que me incomoda em (necessariamente) modelar os parâmetros é: suponha que eu vá fazer vários modelos, cada um com uma distribuição diferente, e que eu esteja na fase de explorar a relação entre minha variável-resposta e a covariável. Suponha que eu avalie que a função logística possa explicar bem a relação entre elas, e que eu queira então incluir a covariável em todos os meus modelos. Me parece estranho eu assumir que será JUSTAMENTE a função LOGÍSTICA que irá explicar a relação entre minha covariável e algum PARÂMETRO de QUALQUER distribuição que eu usar! A relação eu eu vi anteriormente foi entre a covariável e minha variável de interesse! Entendeu? Neste caso, existem outras duas opções que me parecem fazer mais sentido: (1) em todos os modelos, usar a logística pra modelar a ESPERANÇA da distribuição em função da covariável e a partir da esperança deduzir o valor do parâmetro; (2) escolher uma função específica que explique a relação entre a covariável e algum PARÂMETRO de CADA distribuição que eu vou usar (ou seja, esquece a relação entre a covariável e a variável-resposta, esquece a logística, e olhe pra relação entre a covariável e um determinado parâmetro de alguma distribuição DE CADA VEZ)"
  * E-mail Pi do dia 07/10/16

Entre Teórica e Prática
  * Gostaria de pensar em uma distribuição que fizesse sentido "mecanístico" segundo a natureza dos meus dados. Pensei em: cada estratégia de vida pode estar associada a uma probabilidade de "estabelecimento" na comunidade, e é este parâmetro que a intensidade e a frequência do distúrbio influenciam. Isso me parece uma distribuição categórica (que é análoga à Bernoulli mas para multinomial), mas para distribuições contínuas. Tinha encontrado a Dirichlet, que parecia adequada à primeira vista, mas acho que ela gera as probabilidades, e não os valores em si (é como se ela estivesse por trás dos meus dados). O que vc acha, PI?

Práticas
  * Por que (usando a função mle2) o erro "Error in optim(par = c(..) : initial value in 'vmmin' is not finite" é recorrente e como me livrar dele?
  * Devo unir intensidade e frequência do distúrbio em um índice?

----

<- [[equipe:luisa:restrito:start|Voltar para Área Restrita]]