我试图从0到7(随替换)随机抽样7个数字,但受限于所选数字加起来为7.因此,例如,输出0 1 1 2 3 0 0是可以的,但输出1 2 3 4 5 6 7不是.有没有办法使用带有附加约束的sample命令？
我打算使用带有sample命令的replicate()函数作为参数,从sample命令返回N个不同向量的列表.我目前使用示例命令的方式(没有任何约束),我需要N非常大,以便获得尽可能多的可能向量,总和精确到7.我认为必须有一个更简单的方法来做到这一点!
这是我的部分代码:
x <- replicate(100000, sample(0:7, 7, replace=T))
理想情况下,我希望x中的10,000或100,000个向量总和为7,但是需要一个巨大的N值才能做到这一点.谢谢你的帮助.
我试图在JAGS中使用逻辑回归模型,但我有(#success y,#attempts n)形式的数据,而不是二进制变量.在R中,可以通过使用glm(y/n~)和"weights"参数将模型拟合到这些数据,但我不确定如何在JAGS中使用它.
这是一个简单的例子,我希望解决我想要问的问题.请注意,我使用的是rjags包.谢谢你的帮助!
y <- rbinom(10, 500, 0.2)
n <- sample(500:600, 10)
p <- y/n
x <- sample(0:100, 10) # some covariate
data <- data.frame(y, n, p, x)
model <- "model{
# Specify likelihood
for(i in 1:10){
y[i] ~ dbin(p[i], n[i])
logit(p[i]) <- b0 + b1*x
}
# Specify priors
b0 ~ dnorm(0, 0.0001)
b1 ~ dnorm(0, 0.0001)
}"