1

私は、場所= 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'

コードにもデータにも明らかな問題は見当たりません。前のオプションを削除すると、コードは正常に実行されます。ポインタをいただければ幸いです。

4

1 に答える 1