私は、場所= 0およびスケール= 2.5の情報量の少ないCauchy事前分布を使用して、ベイジアンロジスティックモデルを推定しようとしています( Gelman et al、2008 )。私はパッケージで使用MCMClogit
しています:MCMCpack
df = read.csv("http://dl.dropbox.com/u/1791181/MCMC.csv")#Read in data
library(MCMCpack)
logpriorfun <- function(beta, location, scale){
sum(dcauchy(beta, location, scale, log=T))
} #From the MCMCpack documentation
mcmc.logit = MCMClogit(SECONDARY.LEVEL ~ AGE + SEX + as.factor(DISTRICT),
tune=.25,
user.prior.density=logpriorfun, location=0, scale=2.5, log=T,
mcmc=10000,
data=df)
しかし、私は以下のエラーメッセージを受け取り続けます:
Error in MCMClogit(SECONDARY.LEVEL ~ AGE + SEX + as.factor(DISTRICT), :
REAL() can only be applied to a 'numeric', not a 'NULL'
コードにもデータにも明らかな問題は見当たりません。前のオプションを削除すると、コードは正常に実行されます。ポインタをいただければ幸いです。