Coleção de Materiais sobre Métodos Computacionalmente Intensivos

• Notas sobre MCI
• arquivo compatível com R com dados das notas sobre MCI acima. Carregue no R com o comando:

source("cimdata.r")
ls(pat=".dat")

• A Compendium of Conjugate Priors

• Máxima Verossimilhança
• Inferência Bayesiana (em desenvolvimento)
  • Entendendo Inferência Bayesiana (Rnw/Tex)(28/08/2007) Para compilar este arquivo há a necessidade de estar instalado o pacote 'mtable' no R.
  • Entendendo Inferência Bayesiana (PDF)(28/08/2007)

• Um exemplo de Regressão Linear Simples
  • Via Método de Reamostragem Ponderada (30/10/2007)
    • Dados da Tese de Edson A. A. Silva
    • Regressão Linear Simples via Reamostragem Ponderada (rnw/tex)
    • Regressão Linear Simples via Reamostragem Ponderada (PDF)
    • Via JAGS (12/11/07)
    • Dados da Tese de Edson A. A. Silva (txt)
    • reglijags.r - Regressão Linear Simples, Arquivo mãe (.r)
    • jags.r - Arquivo que deve estar no mesmo diretório que reglijags.r
    • reglijags.bug - Arquivo contendo o modelo (.bug)
    • reglijags.cmd - Arquivo criado por reglijags.r (.cmd)
    • reglijags-data.r - Arquivo criado por reglijags.r (.r)
    • reglijags-inits.r - Arquivo criado por reglijags.r (.r)
    • codachain1.txt - Arquivo de saída, criado por reglijags.r (.txt)
    • codaindex.txt - Arquivo de saída, criado por reglijags.r (.txt)
  • Via UMacs (22/11/07)
    • Dados da Tese de Edson A. A. Silva (txt)
    • Regressão Linear Simples-UMacs (.R)
• Exemplos de Inferência Bayesiana
  • Modelos Poisson e Espaciais (PDF)-Luiz (08/11/07)
• Estimação da precisão de um estimador via Bootstrap - Ana e Guerino (14/01/08)
  • boot.r - Precisão do estimador (.r)
    • theta1.r - Arquivo utilizado em boot.r se o estimador é a média (.r)
    • theta2.r - Arquivo utilizado em boot.r se o estimador é a mediana (.r)
    • theta3.r - Arquivo utilizado em boot.r se o estimador é a variância (.r)
    • theta4.r - Arquivo utilizado em boot.r se o estimador é a correlação (.r)
    • theta5.r - Arquivo utilizado em boot.r se o estimador é a regressão (.r)
  • boot1.r - Precisão do estimador - uma versão melhorada da anterior (.r)
  • bootcor.r - Precisão do estimador correlação-Ana (.r)
  • bootcorg.r - Precisão do estimador correlação-Guerino (.r)
  • regliboot.r - Precisão do estimador regressão-Ana (.r)
  • bootregg.r - Precisão do estimador regressão-Guerino (.r)
  • Seminário sobre Boostrap (.tex)
• Algoritmo EM
  • distmultiem.r - Algoritmo EM, Distribuição Multinomial (.r)
  • emseminario.tex - Seminário sobre Algoritmo EM
  • Um exemplo de regressão linear via Algoritmo EM -Guerino 22/01/08
    • dadosem.txt
    • emh.r
    • ordem.r
    • regressao_em.r

Dei uma melhorada no algoritmo de regressão para dados censurados. Na apostila do Stuart, o mesmo está limitado a um conjunto

    específico de dados (40 observações, sendo 23 censurados e 17 não censurados).
         Além de algumas alterações no algoritmo inicial, criei uma função "ordena" (dentro do script ordem.r), a qual toma qualquer conjunto 
    de dados (contendo valores censurados ou não) e ordena o mesmo de forma que as informações censuradas venham primeiro e sejam contadas. 
       Assim o algoritmo fica generalizado, fazendo a regressão p/ quaisquer dados (censurados e não censurados).
       A unica recomendação é que a matriz contendo os dados tenha 3 colunas:
          - 1a coluna: contem zero's e um's, indicando se o dado é não censurado ou censurado
          - 2a coluna: contem a variável independente
          - 3a coluna: contem a variável dependente
 

1. Exemplos utilizando funcao MCMCmetrop1R() do pacote 'MCMCpack'

• Modelo Normal com efeito aleatorio
  • Modelo normal misto (.Rnw)
  • Modelo normal misto (.pdf)
• Poisson
  • Poisson (.R)
• CAR
  • CAR (.R)
• Geoestatistico
  • Modelo geoestatistico (.R)

1. Simple linear regression using UMacs package
   • Simple linear regression using UMacs package (.R) (08/11/07)

1. Modelos de Inferência Bayesiana via WinBugs
   • Modelos Poisson e Espaciais (.txt)-Luiz (08/11/07)