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:
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
2 matches
Mail list logo