[R-br] Script de Bernouilli

12 Out 2013

      Olá!!

Alguém, por gentileza, poderia me dizer por que o modelo está dando erro?

Segue:

# Gerando dados
tabela<-read.csv("morfo.csv",sep=";",head=T)
tabela$SX<-as.factor(tabela$SX)

#criar a lista de dados para o jags
lddados <- list(y=tabela$SX,
    N=length(tabela$SX),
    CT=tabela$CT,
    CF=tabela$CF,
    CP=tabela$CP,
    CC=tabela$CC,
    FO=tabela$FO,
    AC=tabela$AC,
    PD=tabela$PD,
    CD=tabela$CD,
    BD=tabela$BD,
    PP=tabela$PP,
    CPT=tabela$CPT,
    BP=tabela$BP,
    PV=tabela$PV,
    BV=tabela$BV,
    LC=tabela$LC,
    CV=tabela$CV,
    PA=tabela$PA,
    BA=tabela$BA,
    CA=tabela$CA,
    AP=tabela$AP,
    LP=tabela$LP,
    CM=tabela$CM,
    DOH=tabela$DOH,
    DOV=tabela$DOV)

#criar modelo
sink('morfo.txt')
cat('
model{
for(i in 1:N){                   # loop sobre as observações
y[i] ~ dbern(p[i])         # saída binária
logit(p[i]) <- beta0       # regressão multivariada
 + beta1*CT[i]
 + beta2*CF[i]
 + beta3*CP[i]
 + beta4*CC[i]
 + beta5*FO[i]
 + beta6*AC[i]
 + beta7*PD[i]
 + beta8*CD[i]
 + beta9*BD[i]
 + beta10*PP[i]
 + beta11*CPT[i]
 + beta12*BP[i]
 + beta13*PV[i]
 + beta14*BV[i]
 + beta15*LC[i]
 + beta16*CV[i]
 + beta17*PA[i]
 + beta18*BA[i]
 + beta19*CA[i]
 + beta20*AP[i]
 + beta21*LP[i]
 + beta22*CM[i]
 + beta23*DOH[i]
 + beta24*DOV[i]
}
# Prioris
beta0 ~ dnorm(0,1.0E-6)
beta1 ~ dnorm(0,1.0E-6)
beta2 ~ dnorm(0,1.0E-6)
beta3 ~ dnorm(0,1.0E-6)
beta4 ~ dnorm(0,1.0E-6)
beta5 ~ dnorm(0,1.0E-6)
beta6 ~ dnorm(0,1.0E-6)
beta7 ~ dnorm(0,1.0E-6)
beta8 ~ dnorm(0,1.0E-6)
beta9 ~ dnorm(0,1.0E-6)
beta10 ~ dnorm(0,1.0E-6)
beta11 ~ dnorm(0,1.0E-6)
beta12 ~ dnorm(0,1.0E-6)
beta13 ~ dnorm(0,1.0E-6)
beta14 ~ dnorm(0,1.0E-6)
beta15 ~ dnorm(0,1.0E-6)
beta16 ~ dnorm(0,1.0E-6)
beta17 ~ dnorm(0,1.0E-6)
beta18 ~ dnorm(0,1.0E-6)
beta19 ~ dnorm(0,1.0E-6)
beta20 ~ dnorm(0,1.0E-6)
beta21 ~ dnorm(0,1.0E-6)
beta22 ~ dnorm(0,1.0E-6)
beta23 ~ dnorm(0,1.0E-6)
beta24 ~ dnorm(0,1.0E-6)
}
')

sink()

#carregar a biblioteca jags
library(rjags)

#definir os parâmetros para as posteriores

params <- c("beta0",
            "beta1",
  "beta2",
            "beta3",
"beta4",
            "beta5",
            "beta6",
  "beta7",
            "beta8",
"beta9",
            "beta10",
            "beta11",
  "beta12",
            "beta13",
"beta14",
            "beta15",
            "beta16",
  "beta17",
            "beta18",
"beta19",
            "beta20",
            "beta21",
  "beta22",
            "beta23",
"beta24")

#inicializar cadeia
inicio <- list(beta0=0,
   beta1=0,
   beta2=0,
   beta3=0,
   beta4=0,
   beta5=0,
   beta6=0,
   beta7=0,
   beta8=0,
   beta9=0,
   beta10=0,
   beta11=0,
   beta12=0,
   beta13=0,
   beta14=0,
   beta15=0,
   beta16=0,
   beta17=0,
   beta18=0,
   beta19=0,
   beta20=0,
   beta21=0,
   beta22=0,
   beta23=0,
   beta24=0)

#chamar o modelo no jags
morfo.m <-
jags.model("morfo.txt",data=lddados,inicio,n.chains=1,n.adapt=1000)

Abraço,
Gisele

[R-br] Script de Bernouilli

Gisele Correia