{

  p ~ dbeta( 1.0, 1.0 )
  r ~ dgamma( 0.001, 0.001 )

  for ( i in 1:n ) {

    y[ i ] ~ dnegbin( p, r )

  }

  y.new ~ dnegbin( p, r )

}