model { X ~ dbin(p, n); p <- 0.2; n <- round(u); u ~ dflat(); } DATI list( X=16 ) INITS list( u=100 )