a*********e 发帖数: 1233 | 1 本人初用WinBUGS。
假设有模型 mu[i] <- beta0+ beta1*p[i]
其中i=100,
mu[i] 是100个已知数据,
p[i]是生成的数据。
问题是,p[i]分成两部分,例如前50个是均匀分布,后50个是正态分布。
下面的代码显然不能运行
for (i in 1:m){
y[i] = dunif(0,10)
}
for(i in m+1:n){
y[i] = dnorm(0,1)
}
m=50, n=100
请问有什么方法能生成这样的数据,或者说能否把生成的几个数据合并成一个变量。
谢谢。 | a********s 发帖数: 188 | 2 You could use "R2WinBUGS" package in R, if you knew how to use R. Data is
generated quite easily in R. | y*****y 发帖数: 98 | 3 不是很明白你的意思。winbugs是MCMC tool,不是programming language,不能用来
generate random numbers或是做calculation, 你只能把model写进去算MCMC,that's
all. | t******l 发帖数: 32 | |
|