dongyuexie / gscca Goto Github PK

set.seed(12345)
p = rbeta(100,2,2) # prior mode: 0.5
n = rpois(100,10)
y = rbinom(100,n,p) # simulate Binomial observations
fit=ash(rep(0,length(y)),1,lik=lik_binom(y,n,link='logit'))
round(fit$result$PosteriorMean,2)
[1] -0.50 -0.52 -0.52 -0.52 -0.52 -0.28 -0.52 -0.51 -0.52 -0.50 -0.52 -0.42 -0.51 -0.52 -0.52 -0.52
[17] -0.50 -0.50 -0.37 -0.52 -0.51 -0.39 -0.51 -0.25 0.00 -0.52 -0.52 -0.52 -0.52 -0.52 -0.52 -0.50
[33] -0.52 -0.51 -0.52 -0.46 -0.27 -0.51 -0.50 -0.46 -0.52 -0.52 -0.52 -0.51 -0.02 -0.52 -0.52 -0.43
[49] -0.52 -0.32 -0.52 -0.52 -0.51 -0.52 -0.50 -0.49 -0.02 -0.36 -0.52 -0.28 -0.52 -0.48 -0.52 -0.10
[65] -0.51 -0.52 -0.52 -0.52 -0.51 -0.42 -0.16 -0.50 -0.52 -0.52 -0.36 -0.52 -0.48 -0.51 -0.48 -0.52
[81] -0.06 -0.51 -0.50 -0.52 -0.45 -0.52 -0.51 -0.51 -0.52 -0.52 -0.51 -0.52 -0.51 -0.06 -0.48 -0.43
[97] -0.52 -0.35 -0.50 -0.03