[R-br] Resíduos e MCMCglmm

[Elias Teixeira Krainski]

Embora hoje em dia seja possível  encontrar o termo "teste de hipóteses 
Bayesiano", isso não é comum. O mais comum é estimar os modelos com e 
sem a componente e compará-los. A opção mais simples é usar fator de 
Bayes. A PC-prior não é contruida para fazer teste de hipóteses, mas 
funciona nesse caso, veja teorema 2 do paper.

Como é uma lista de R, segue um código:

n= 100; a = 2; tau1 = 1
x = rnorm(n, 0, sqrt(1/tau1))
y = rpois(n, exp(a + x))

require(INLA)
res0 = inla(resp ~ 1, family='poisson', data=list(id=1:n, resp=y))
### priori para a precisao do efeito aleatorio IID, x_1, ..., x_n
###    P( std(x) >1 ) = 0.01
pp = list(theta=list(prior='pc.prec', param=c(1,0.01)))
res1 = inla(resp ~ f(id, model='iid', hyper=pp),
    family='poisson', data=list(id=1:n, resp=y))
res1$mlik/res0$mlik

Elias