====== Mini-curso: métodos computacionais para inferência estatística ======
===== Escola de Verão, UFRN, Natal, 28-30/01/2014 =====
O //download// de arquivos abaixo **PODE NÃO FUNCIONAR** com o //internet explorer//. \\
Utilize qualquer outro navegador (firefox, opera, Chromium, etc)
- [[http://www.leg.ufpr.br/~paulojus/mcie/masterSINAPE-2012-08-13.pdf|Texto base do curso]] (cap 1-3)
- Dia 1:
- {{:cursos:natal2014:binom.r|Espaço do modelo}}
- Conceitos básicos: {{:cursos:natal2014:binomvero.r|exemplo com a distribuição binomial}}
- Dia 2
- {{:cursos:mcie:exponencial.r|Inferência na distribuição exponencial}}
- {{:cursos:mcie:infpois.r|Inferência na Poisson usando diferentes algorítimos/funções}}
- Dia 3
*- {{:cursos:mcie:ar1.r|Exemplo AR1}}
- {{:disciplinas:ce089-2013-02:misc.r|Miscelânea de tópicos}}
====== Alguns Códigos ======
* {{:cursos:mcie:ar1.r|Exemplo AR1}}
* {{:disciplinas:ce089-2013-02:misc.r|Miscelânea de tópicos}}
* {{:disciplinas:ce089-2013-02:integranum.r|Integração numérica}}
* Atividades:
- Considere a seguinte amostra de uma distribuição G(\alpha, \beta = 2) com densidade f(x) = \frac{\beta^\alpha}{\Gamma(\alpha)} x^{\alpha-1} \exp\{-\beta x\}
7,0 2,6 8,8 6,0 5,7 4,6 3,0 3,5 5,3 5,8 4,8 4,3
* Obtenha a expressão das funções de verossimilhança, log-verossimilhança, escore e hessiana
* Faça gráficos para cada uma dessas funções
* Escreva códigos para obtenção da estimativa de verossimilhança
* por maximização numérica
* por solução de sistema
* Obtenha um intervalos de confiança
* baseado na aproximação quadrática
* baseado na função de verossimilhança
- Seja a seguinte amostra de uma distribuição X \sim {\rm N}(\mu, \sigma^2 = 9)
27 <25 22 20 [20-24] 22 26 23 >23 18 26 <32 22 23 >20
* resolva os mesmos itens da questão anterior
* {{:disciplinas:ce089-2013-02:gamma.pdf|Inferência na distribuição Gama}}
* Refazer os exemplos do arquivo, complementando para cada parametrização
* Montar/escrever o algoritmo Newton-Rapson para este problema
* Montar um material similar para a distribuição Beta
* Verificar possível(eis) (ao menos 1) reparametrização(ções)
* Verificar se é possível reduzir a dimensão do problema numérico (isolar 1 parâmetro)
* Montar soluções similares às do material da distribuição Gamma
====== Dados ======
- Capítulo 1:
- Capítulo 2:
- Dados simulados Exemplo 5 - Distribuição Gaussiana
set.seed(123)
y10 <- rnorm(10,10,1.5)
y50 <- rep(y10,5)
y100 <- rep(y10,10)
y1000 <- rep(y10,100)
- Capítulo 3:
* Exemplo: Regressão Beta - {{:cursos:mcie:dadosbeta.txt|}}