d****a
发帖数: 193 | 1 I made a graph for a model using winbugs doodle. The graph looked nice when
I
first created it. However, when I opened it the second time, it became very
very small.
I am a winbugs novice. Can someone help tell me how to zoom the graph to
make
it larger. Moreover, can some tell me how to export this graph, to use it in
MS
Word or LaTex? Many thanks. |
|
m***n
发帖数: 127 | 2 用WINBUG能做什么？
在WINBUGS下面写了这么简单的一段，model syntax通过了，可是data就是不能load。
试过list和matrix,都不行。请高手看看为什么。多谢。model {........ |
|
t****x
发帖数: 1650 | 3 有个项目要用winbugs推导missing，从来没用过，谁有coding example没有？谢了。 |
|
x**g
发帖数: 807 | 4 能说具体些吗？
What is the missing mechanism you would consider in the project? Any
specific missing data method you will apply? Why consider WinBugs? There are
tons of package can handle missingness.
I believe your explanation will help the readers provide their suggestion. |
|
s**u
发帖数: 383 | 5 刚开始学WINBUG, 写了个小程序, model 问题, 但是load data 时说 expected
variable name. 请问哪里出错了. 谢谢.
model {
for (i in 1:9) {n[i] <- nwin1[i] + nwin2[i]
winn1[i] ~dbin(p[i],n[i])
logit(p[i]) <- d[dind1[i]]-d[dind2[i]]}
for ( i in 1:2) {d[i] ~ dnorm(0, 0.001)}
for ( i in 4:5) {d[i] ~ dnorm(0, 0.001)}
d[3] <-0
}
list (dind1 = c(1,1,1,1,2,2,2,3,4),
dind2 = c(2,3,4,5,3,4,5,5,5),
nwin1 = c(9,12,6,27,9,12,12,2,2),
nwin2 = c(10,2,3,2,5,3,3,0,4))
list(d =c(0, 0, NA,0,0)) |
|
a*********e
发帖数: 1233 | 6 本人初用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
请问有什么方法能生成这样的数据，或者说能否把生成的几个数据合并成一个变量。
谢谢。 |
|
y*****y
发帖数: 98 | 7 不是很明白你的意思。winbugs是MCMC tool,不是programming language,不能用来
generate random numbers或是做calculation, 你只能把model写进去算MCMC，that's
all. |
|
m**c
发帖数: 88 | 8 model{
for (i in 1:n){
y[i,1:2]~ dmnorm( mu[], T[,] )
}
for (i in 1:P){
mu[i] <- inprod(A[i,],x[,])+ A0[i,1]
}
x[1,1] ~ dnorm(0.0, 1.0E-4)
x[2,1] ~ dnorm(0.0, 1.0E-4)
#x[1,1] ~ dunif(3,7)
#x[2,1] ~ dunif(4,8)
T[1:P,1:P] <- inverse(V[,])
}
INITS
list(x = structure(.Data = c(3,6), .Dim =c(2,1)))
DATA (RECT.)
list( n=10, P=2,
A = structure(.Data = c(1.5,1.2, 2, 0.7),.Dim =c(2,2)),
A0 = structure(.Data... 阅读全帖 |
|
m**c
发帖数: 88 | 9 model{
for (t in 1:n){
y[t,1:5]~ dmnorm( mu[], T[,] )
}
index ~ dcat(Prob[])
for (i in 1:H) {
x_spike[i] ~ dnorm(1,1.0E10)
x_slap[i] ~ dnorm(0.98,100)
lambda[i] <- Index[index,i]
x[i,1] <- (1- lambda[i])* x_spike[i] + lambda[i]* x_slap[i]
}
for (j in 1:P){
mu[j] <- inprod(A[j,],x[,])+ A0[j,1]
}
T[1:P,1:P] <- inverse(V[,])
}
DAtA
.................
.................
.................
END
这个WinBUGS程序我可以在Sampl... 阅读全帖 |
|
G*********o
发帖数: 2045 | 10 安装了winbugs，双击之后看不到GUI界面，请问有人遇到过类似的问题吗？实在是太
weird了，在任务管理器里面能看到winbugs14.exe占用cpu进程，但是就是没有图形显
示。。。 |
|
s******d
发帖数: 501 | 11 本人WinBUGS新手，最近要跑一个simulation，大概要跑50-100组数据，有什么办法能
够自动完成对全部数据的分析，不需要每换一组数据就手动一次吗？
谢谢！ |
|
g*******u
发帖数: 148 | 12 Search "R2WinBUGS"
Manipulate data in R and do analysis iwith WinBUGS |
|
A**P
发帖数: 260 | 13 在WINBUGS下面写了这么简单的一段，model syntax通过了，可是data就是不能load。
试过list和matrix,都不行。请高手看看为什么。多谢。
model {
#model
nI ~ dbin(p, N)
nS ~ dbin(s,nI)
lor <- log(p*s/(1-p*s))
orm ~ dnorm(lor, orp)
#priors
p ~ dbeta(1,1)
s ~ dbeta(1,1)
}
#Data
list (nI=500, nS=200, N=1421, orm= -1.7, orp= 100)
nI [] nS [] N [] orm [] orp []
500 200 1421 -1.7 100
END |
|
p********1
发帖数: 1011 | 14 一共三层：
第一层： individual events,
第二层： urban facilities (e.g. grocery stores, bus stops),
第三层： neighborhood unemployment rate and poverty rate.
要用Hierarchical Bayesian Spatial Model做分析。 求winbugs code. |
|
R*j
发帖数: 248 | 15 刚开始学用Stan, 和Winbugs比起来, 感觉写出的模型不够一目了然。而Bugs就非常清
晰。其它同学觉得呢？ |
|
p********a
发帖数: 5352 | 16 ☆─────────────────────────────────────☆
keke823 (keke823) 于 (Tue Oct 25 20:38:47 2011, 美东) 提到:
一周前投了简历给职位A, 隔了一天就回信问我约个什么时候电面职位B，
晕，看B的要求，最多看懂的只有一半的东西。 于是回信问是不是搞错了， 无回应
隔了一天再问是不是搞错了， 另外通知了available的时间，心想不会就不会吧， 等
了一天依旧无回应
发现A职位的申请已经被拒，B的状态也已经关闭了，于是在发邮件 问到底有没有电面
？已经无回应
本来以为要Move on了， 今天直接打来电话约时间，还是职位B。 彻底晕了
想问问各位有没有什么经验可以参考，本人纯Bayes统计，没有任何CS背景， C和Java
都不会，是不是没有什么希望？
☆─────────────────────────────────────☆
Dreamer1122 (MX) 于 (Tue Oct 25 22:19:46 2011, 美东) 提到:
好牛那～纯bayes～
☆───────────────... 阅读全帖 |
|
o******6
发帖数: 538 | 17 ☆─────────────────────────────────────☆
birspring (bir) 于 (Sat Feb 21 01:19:30 2009) 提到:
最近在看一篇paper.老板打算让我写那篇文章的codes.用的是bayes方法。那篇文章的
作者已经用winbugs做过一些numerical的分析了。但是老板说我们要分析的数据比较大
，而winbugs只是用小数据玩玩。要我用R写codes.所以上来问问。是不是所有的bayes
方法都可以用R写出程序。有大牛给个例子吗？还想问问，用sas能写出bayes方法的
codes吗？多谢！
☆─────────────────────────────────────☆
fanta (Icecode/Ashh/Tallerr) 于 (Sat Feb 21 02:00:19 2009) 提到:
我是菜鸟
我觉得可以吧，用gibbs sampling的话，先把full conditionals算好，再用R照
algorithm写就行了。用metropolis-hasting也是差不多。
我在写一个小的 |
|
o****o
发帖数: 8077 | 18 I think Dr.Burnie wanted to say that because WinBUGS shield users from the
real computational issues and makes Bayesian computing seemly easy, which
has a negative impact on the ppl using WinBUGS to do bayesian statistics |
|
m**c
发帖数: 88 | 19 LZ,你的问题,我不懂啊,帮不了你,你能帮帮我吗?
我在MCMC方面小白啊,只能借助专业软件,比如WinBUGS来做采样,现在碰到一个问题,因
为要调用另外的DLL,这个在WinBUGS里做不到,我必须得用Matlab来做MCMC仿真了,但是
现在让我自写一段MH的Matlab code 很难啊,能不能把你那段"MCMC algorithm using
Metropolis-Hastings algorithm" 让我看看啊? 咱们一起探讨一下啊
LZ是个好人! |
|
a*****i
发帖数: 1045 | 20 第一年，完全没什么概念呢。已经开始选择论文题目了。现在有几个option，不知道有
没有好心人可以帮忙看看哪些方向比较好，或者对将来工作会有帮助。
自己之前定了跟statistical process control 有关的方向。但不怎么样。
1. Title : A comparison of different statistical process monitoring schemes
on real life data （这个是关于statistical process control的）
It is often desirable to provide online monitoring for industrial or
agricultural processes. Faced with the task of performing such monitoring,
practitioner will often choose to apply a control chart. However, there many
options available, and it i... 阅读全帖 |
|
s****e
发帖数: 1180 | 21 面试问起programming ability该怎么回答？有什么样板回答？已经不只一次在面试中
被问到自己的programming ability。已经读到这个阶段了，用过的语言/software也不
少，sas，R，perl，C，minitab，winbugs，但是究竟如何回答这个问题，感觉还是有
点无从下手。大家有什么建议？谁给写个回答样本。谢谢。 |
|
s****e
发帖数: 1180 | 22 【 以下文字转载自 LosAngeles 讨论区 】
发信人: sheide (shei), 信区: LosAngeles
标 题: 面试问起programming ability该怎么回答？有什么样板回答？
发信站: BBS 未名空间站 (Thu Jun 23 16:44:42 2011, 美东)
面试问起programming ability该怎么回答？有什么样板回答？已经不只一次在面试中
被问到自己的programming ability。已经读到这个阶段了，用过的语言/software也不
少，sas，R，perl，C，minitab，winbugs，但是究竟如何回答这个问题，感觉还是有
点无从下手。大家有什么建议？谁给写个回答样本。谢谢。 |
|
|
g*d
发帖数: 6969 | 24 winbugs
和你们的code没发比
那个软件能做的有限。所以要做复杂的很麻烦 |
|
s*****0
发帖数: 357 | 25 【 以下文字转载自 Pharmaceutical 讨论区 】
发信人: HeadhuntingJ (Janelee), 信区: Pharmaceutical
标 题: 谁归？Exploratory Biostatistician&Translational Medicine Leader
关键字: 罗氏 职位 北京&上海
发信站: BBS 未名空间站 (Mon Apr 18 02:10:40 2011, 美东)
Exploratory Biostatistician：PHD+2年博后经验；罗氏 北京；
Translational Medicine Leader：罗氏 上海；
Exploratory Biostatistician
公司：罗氏
地点：北京/上海
招聘要求：
Who you are
You’re someone who wants to influence your own development. You’re looking
for a company where you have the opportunity to pursue your interests
ac... 阅读全帖 |
|
o****o
发帖数: 8077 | 26 是呀，如果bayes那个就是教你如何用WINBUGS，估计一节课也够了 |
|
q**j
发帖数: 10612 | 27 多谢。能够手写出来的还是用手的好。用 Gibbs sampling 还要从matlab里面call
WinBugs。麻烦的很。不容易给别人解释。 |
|
y**g
发帖数: 197 | 28
of
==============================
what kind of your algorithm for Bayesian
analysis? is it different with the one used in WinBug(MCMC+GIBBS SAMPLING)?
implement |
|
i**f
发帖数: 1195 | 29 The theoretical background is the same, combined standard Gibbs and a more
general MH algorithm. However
by writing you own algorithm, it is much more flexible than using the build-
in packages from those
statistical software, ie winbug pkbug,etc. :)
Bayesian
)? |
|
h***t
发帖数: 2540 | 30 there is an item called google |
|
t****x
发帖数: 1650 | 31 用Bayesian statistics， 可能MCMC，不知道一般人用什么
are |
|
x**g
发帖数: 807 | 32 关于Missing Data 的理论和方法有很多。基于不同的Missing Mechanism （比如，
MCAR， MAR，或是MNAR），所采用的处理方法是不一样的。
你可以看看Lee（2007）StructuralEquationModeling：A Bayesian Approach. John
Wiley & Sons Ltd. 里面以一章讲的是Missing Data。 Code应该在网站上找得到。 |
|
c*******7
发帖数: 2506 | 33 我回答你的第二个问题：
你用QQ的截图功能吧，非常好用。截下来paste到work里面。
when
very
in |
|
|
|
f***a
发帖数: 329 | 36 you guys are all frequentists? -_-| |
|
o****o
发帖数: 8077 | 37 I write my own code in R when I conduct Bayesian analysis |
|
r******n
发帖数: 351 | 38 hierarchical linear model analysis，的公式在Gelman 的 "Bayesian Data
Analysis"有较详细的介绍。 很多hierarchical linear models 没有analytic的解，
只能用MCMC求解。 可以用 WinBugs （http://www.mrc-bsu.cam.ac.uk/bugs/） 这个免费软件。 |
|
x**g
发帖数: 807 | 39 这些模型是analytical equivalent的，由于用于不同领域，所以有了不同的名称。一
般用SAS PROC MIXED， PROC NLMIXED， WinBugs， R 都可以做的。 |
|
x**g
发帖数: 807 | 40 不要用SAS。 Proc Calis不能处理Missing Data。 还是用Lisrel， Amos，或是Mplus
好。WinBugs也好。 |
|
b*******g
发帖数: 513 | 41 Hello!
Sometimes, the WinBUGS prompts a warning message saying something like mcmc
can be dangerous.
What does the warning message exactly say?
Can someone tell me?
Thank you very much! |
|
|
D*********2
发帖数: 535 | 43 The non-clinical statistics group of MedImmune (Gaithersburg, MD) is seeking
an summer intern. The position is available immediately. The intern will
work under the supervision of a mentor to carry out statistical methodology
research. The responsibilities include implementing of novel methods,
performing simulation studies, and assisting in manuscript preparation.
Excelling computing skills using R and/or WinBUGS required. Knowledge of
longitudinal data analysis, survival analysis and Bayesia |
|
W*****r
发帖数: 193 | 44 model {
for (j in 1:J){
y[j] ~ dnorm(theta[j], tau.y[j])
theta[j] ~dnorm(mu.theta, tau.theta)
tau.y[j] <- pow(sigma.y[j], -2)
}
mu.theta ~ dnorm(0, 1.0E-6)
tau.theta <- pow(sigma.theta, -2)
sigma.theta ~ dunif (0, 1000)
either<-max(theta1,theta2，theta3,theta4, theta5, theta6, the7, theta8)
问题在这一行，我共有8个theta值。这样输入说“没有定义theta1...“
我想请教一下：应该怎么输入？
多谢多谢。会发包子的。
}
DATA
list(J=8, y = c(28, 8, -3, 7, -1,1,18,12), sigma.y=c(15,10,16,11,9,11,
10,18))
INITIAL VALUES
list(theta= c(0,0,0,0,0,0,0,0), mu.th... 阅读全帖 |
|
W*****r
发帖数: 193 | 45 我知道max(a,b)
但这个（a,b)怎么表示8个theta呢？
难道都两个两个过一遍？
我试过 max(theta[j],theta[j+1])
还是不行撒 |
|
l*********s
发帖数: 5409 | 46 1. theta[1] not theta1
2. max() takes only 2 arguments. |
|
l*********s
发帖数: 5409 | 47 use sort function, which takes vector input |
|
|
W*****r
发帖数: 193 | 49 model {
for (j in 1:J){
y[j] ~ dnorm(theta[j], tau.y[j])
theta[j] ~dnorm(mu.theta, tau.theta)
tau.y[j] <- pow(sigma.y[j], -2)
}
mu.theta ~ dnorm(0, 1.0E-6)
tau.theta <- pow(sigma.theta, -2)
sigma.theta ~ dunif (0, 1000)
either<-max(sort[theta[1], theta[2],theta[3],theta[4]],sort[theta[5],
theta[6],theta[7],theta[8]]))
}
我又试了一下，说有variable未定义。BTW，我只用两个theta值的时候，either可以
set成node，density也可以显示。
多谢多谢。
DATA
list(J=8, y = c(28, 8, -3, 7, -1,1,18,12), sigma.y=c(15,10,16,11,9,11,
10,18))
INITIAL VALUES... 阅读全帖 |
|
l*********s
发帖数: 5409 | 50 do this:
thetaorder[1:8] <- sort(theta)
either <- tetaorder[8] |
|
