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