[R] MCMClogit: normal intercept priors don't work

2009-10-07 Thread AlexanderDueDo
Dear R users, I think that the question I posted yesterday was not specific enough. Apologies for this. The question should have been: why does MCMClogit stop working if I a) set a normal prior for the intercept and b) remove all variables except for the intercept as explanatory? Example code:

[R] MCMClogit confusion

2009-10-06 Thread AlexanderDueDo
Hello, running the code simulatedCase <- rbinom(100,1,0.5) simDf <- data.frame(CASE = simulatedCase) posterior_m0 <<- MCMClogit(CASE ~ 1, data = simDf, b0 = 0, B0 = 1) always results in an acceptance ratio of 0. If I do not specify b0 and B0, this does not occur. I do not understand the reasons