Skip to main content
título editado
Enlace

Estadística bayesiana, percentiles 5 % y 95 % en la distribución beta con r y R2OPENBGSR2OPENBUGS

se añadieron 3 caracteres en el cuerpo
Origen Enlace
Alfabravo
  • 9.9k
  • 5
  • 24
  • 32

Necesito obtener los percentiles 5 % y 95% en la salida de la instrucción bugs en R2openBUGS, para poder hacer un contraste de hipótesis unilateral, con un intervalo de credibilidad bayesiano del 95%. Les adjunto el código en r y openBUGS. El problema es que no encuentro la instrucción dentro de Bugs en R2openBUGS para poder indicar o setear más percentiles o cuantiles. ¿Pueden ustedes por favor ayudarme en esta tarea? Muchas gracias por la ayuda.

```        library("R2OpenBUGS") 
        library("coda") 
        library("lattice")

        n4 = 10000  
        n5 = 1000
        n6 = 1
        n7 = 1
        N= 10
        r = 8
        a = 1
        b = 1

        mod<-"mbp2022.txt"

        data <- list( N = N,
              r = r,
              a = a,
              b = b          
        )
        sink(mod)
        cat("
       
        model{

       pi1 ~ dbeta(a,b)    
       r ~ dbin(pi1, N)

      }

      ",fill=TRUE) 
      sink()

      modelfile <- mod


      vars <-list("pi1")

      inits <-function(){list(pi1=0.5)} 

 

      out <-bugs(
  
      model.file=file.path(modelfile),
  
      data=data,
  
      inits=inits,
  
      parameters.to.save=vars, 
   
      n.iter=n4, 
      n.burnin=n5, 
      n.chains=n6,
      n.thin=n7,
      summary.only = F,
      debug=F, 
      DIC=F, 
      working.directory = getwd(),
  
      )
      cat("\n")  
      print (out, digits = 03)
  

Necesito obtener los percentiles 5 % y 95% en la salida de la instrucción bugs en R2openBUGS, para poder hacer un contraste de hipótesis unilateral, con un intervalo de credibilidad bayesiano del 95%. Les adjunto el código en r y openBUGS. El problema es que no encuentro la instrucción dentro de Bugs en R2openBUGS para poder indicar o setear más percentiles o cuantiles. ¿Pueden ustedes por favor ayudarme en esta tarea? Muchas gracias por la ayuda.

``` library("R2OpenBUGS") 
    library("coda") 
    library("lattice")

    n4 = 10000  
    n5 = 1000
    n6 = 1
    n7 = 1
    N= 10
    r = 8
    a = 1
    b = 1

    mod<-"mbp2022.txt"

    data <- list( N = N,
          r = r,
          a = a,
          b = b          
    )
    sink(mod)
    cat("
   
    model{

   pi1 ~ dbeta(a,b)    
   r ~ dbin(pi1, N)

  }

  ",fill=TRUE) 
  sink()

  modelfile <- mod


  vars <-list("pi1")

  inits <-function(){list(pi1=0.5)} 



  out <-bugs(

  model.file=file.path(modelfile),

  data=data,

  inits=inits,

  parameters.to.save=vars, 

  n.iter=n4, 
  n.burnin=n5, 
  n.chains=n6,
  n.thin=n7,
  summary.only = F,
  debug=F, 
  DIC=F, 
  working.directory = getwd(),

  )
  cat("\n")  
  print (out, digits = 03)
 

Necesito obtener los percentiles 5 % y 95% en la salida de la instrucción bugs en R2openBUGS, para poder hacer un contraste de hipótesis unilateral, con un intervalo de credibilidad bayesiano del 95%. Les adjunto el código en r y openBUGS. El problema es que no encuentro la instrucción dentro de Bugs en R2openBUGS para poder indicar o setear más percentiles o cuantiles. ¿Pueden ustedes por favor ayudarme en esta tarea? Muchas gracias por la ayuda.

        library("R2OpenBUGS") 
        library("coda") 
        library("lattice")

        n4 = 10000  
        n5 = 1000
        n6 = 1
        n7 = 1
        N= 10
        r = 8
        a = 1
        b = 1

        mod<-"mbp2022.txt"

        data <- list( N = N,
              r = r,
              a = a,
              b = b          
        )
        sink(mod)
        cat("
       
        model{

       pi1 ~ dbeta(a,b)    
       r ~ dbin(pi1, N)

      }

      ",fill=TRUE) 
      sink()

      modelfile <- mod


      vars <-list("pi1")

      inits <-function(){list(pi1=0.5)} 

 

      out <-bugs(
  
      model.file=file.path(modelfile),
  
      data=data,
  
      inits=inits,
  
      parameters.to.save=vars, 
   
      n.iter=n4, 
      n.burnin=n5, 
      n.chains=n6,
      n.thin=n7,
      summary.only = F,
      debug=F, 
      DIC=F, 
      working.directory = getwd(),
  
      )
      cat("\n")  
      print (out, digits = 03)
 
Origen Enlace

Estadística bayesiana, percentiles 5 % y 95 % en la distribución beta con r y R2OPENBGS

Necesito obtener los percentiles 5 % y 95% en la salida de la instrucción bugs en R2openBUGS, para poder hacer un contraste de hipótesis unilateral, con un intervalo de credibilidad bayesiano del 95%. Les adjunto el código en r y openBUGS. El problema es que no encuentro la instrucción dentro de Bugs en R2openBUGS para poder indicar o setear más percentiles o cuantiles. ¿Pueden ustedes por favor ayudarme en esta tarea? Muchas gracias por la ayuda.

``` library("R2OpenBUGS") 
    library("coda") 
    library("lattice")

    n4 = 10000  
    n5 = 1000
    n6 = 1
    n7 = 1
    N= 10
    r = 8
    a = 1
    b = 1

    mod<-"mbp2022.txt"

    data <- list( N = N,
          r = r,
          a = a,
          b = b          
    )
    sink(mod)
    cat("
   
    model{

   pi1 ~ dbeta(a,b)    
   r ~ dbin(pi1, N)

  }

  ",fill=TRUE) 
  sink()

  modelfile <- mod


  vars <-list("pi1")

  inits <-function(){list(pi1=0.5)} 



  out <-bugs(

  model.file=file.path(modelfile),

  data=data,

  inits=inits,

  parameters.to.save=vars, 

  n.iter=n4, 
  n.burnin=n5, 
  n.chains=n6,
  n.thin=n7,
  summary.only = F,
  debug=F, 
  DIC=F, 
  working.directory = getwd(),

  )
  cat("\n")  
  print (out, digits = 03)