statistiques bayésiennes
Répondre à la discussion
Affichage des résultats 1 à 1 sur 1

statistiques bayésiennes



  1. #1
    invitece4ad8a9

    statistiques bayésiennes


    ------

    Bonjour,
    J'ai un souci avec un modèle que j'ai écrit dans le cadre d'un travail en statistiques bayésiennes. Je n'arrive pas à le compiler sur openbugs pourtant le fichier
    log montre que le modèle est correct.
    voici le modèle que je cherche à compiler sur openbugs
    Code:
    model { 
            for(i in 1:N){ 
               y[i] ~ dbin(psi[i],1)   #binomial distribution   
               logit(psi[i]) <- b0 + b1*foot[i] + b2*sex[i]             
                
             } 
            #prior
            p ~ dbeta(1,1)                
            b0 ~ dnorm(0,.001) 
            b1 ~ dnorm(0,.001)
            b2 ~ dnorm(0,.001)  
            logit(psi0)<-b0
          }
    et voici celui que j'ai fait sur R en premier lieu
    library(R2OpenBUGS)
    library(coda)
    
    #Extraction des données
    data=read.table("données.txt",h=T)
    attach(data)
    names(data)
    dim(data)
    
    N=nrow(data)-1 #nombre d'individus
    
    foot=data[,1] 
    y=data[,2]
    sex=data[,3]
    
    #on  fabrique les données
    data=list("foot","y","N","sex")
    
    #
    ni=5000 #nbre d'iteration
    nb=2000 #nbre d'iteration qu'on jette au début
    nt=2
    nc=2 #nbre de chaines
    
    inits <- function()
      list (b0=rnorm(1),p=runif(1),b1=rnorm(1),b2=rnorm(1))
    
    #Paramètres à monitorer
    parameters <- c("p","b0","b1","b2","psi0") 
    
    out<-bugs(data, inits, parameters, "modèle.txt", n.thin=nt,n.chains=nc, 
              n.burnin=nb,n.iter=ni,debug=TRUE, coda=T,working.directory="D:/COURS/Statistiques bayésiennes/projet",
              OpenBUGS.pgm = NULL)
    modèle.txt est celui que j'ai écrit plus haut

    J'ai vraiment besoin d'un coup de main. Merci d'avance!

    -----
    Dernière modification par JPL ; 02/02/2015 à 17h33. Motif: Ajout de la balise Code (#) pour garder l'indentation

Discussions similaires

  1. [Biotechnologie] Statistiques
    Par invite701a4543 dans le forum Biologie
    Réponses: 2
    Dernier message: 27/01/2011, 12h59
  2. Statistiques un jour, statistiques toujours...
    Par invite15f14b48 dans le forum Mathématiques du supérieur
    Réponses: 7
    Dernier message: 04/01/2009, 12h51
  3. Statistiques
    Par invite90128b89 dans le forum Mathématiques du supérieur
    Réponses: 1
    Dernier message: 10/08/2007, 20h32
  4. statistiques
    Par invite7a5c5550 dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 26/07/2007, 17h28
  5. Probas bayesiennes en phylogénie
    Par invite533a42a8 dans le forum Biologie
    Réponses: 0
    Dernier message: 16/12/2006, 19h48