i new user r , perform exact logistic regression. have chosen method because dependent variable binary (i.e. 0 failure , 1 success), , have continuous independent variables. have small sample size of n=10. following sample of code have written. library(elrm) baselinesteps = c(4202, 7244, 4374.6, 2965.6, 4263.5, 1814.1, 3243.1, 3102.9, 4652.5, 6324.9) stepsclass = c(1, 0, 1, 1, 1, 1, 1, 1, 1, 1) const = rep(1,length(stepsclass)) dat = data.frame(pred = baselinesteps, trials = const, success = stepsclass) m.stepsfromsteps = elrm(formula = success/trials ~ pred, interest = ~pred, iter = 22000, dataset = dat, burnin = 2000) i trying use baseline average steps/day taken subject determine class of recovery pattern fall (0 or 1). however, when run last line elrm implemented (changing nothing), obtain drastically different p-values every time run it. see following: > m.stepsfromsteps = elrm(formula = success/trials ~pred, interest = ~pred, iter=30000, dataset = dat, burnin=5000) ...