有winbugs的高手，帮忙看看这个文献，感激不尽

epoh 发表于 2011-12-27 10:49
你说的问题
就compressor而言只要更改initia values
inits = function() {list(alpha1 =14, alpha2 = ...

老师您好，initial values 的设置对结果影响会这么大啊，那怎么来选择initial values呢

博公 发表于 2011-12-27 14:55
老师您好，initial values 的设置对结果影响会这么大啊，那怎么来选择initial values呢

1.对模型的了解
2.try and error

epoh 发表于 2011-12-27 15:23
1.对模型的了解
2.try and error

那结果的重现性和 initial values 有很大的关系了？可不可以这样认为：只要选到合适的 initial values，文献的结果完全可以重现呢？

博公 发表于 2011-12-27 16:56
那结果的重现性和 initial values 有很大的关系了？可不可以这样认为：只要选到合适的 initial values，文 ...

对,可以这样说
different components use the same model
different datas use the same model
hydro_revised(need pw)
  include compressor,flanges,pipes.
   hydro_revised.rar (2.56 KB)

2011-12-28 09:59:50 上传

老师，您好，非常感谢您的解答，请问文献中Figure 4-1，以及后面的Figure 4-2 Figure 4-3，······都是用什么软件画出来的？

博公 发表于 2011-12-28 17:06
老师，您好，非常感谢您的解答，请问文献中Figure 4-1，以及后面的Figure 4-2 Figure 4-3，······都是 ...

Figure 2-1.Figure 2-2.Figure 2-3.
用Microsoft Office 就可以画出来了
而且比一般计量软件漂亮

epoh 发表于 2011-12-28 19:00
Figure 2-1.Figure 2-2.Figure 2-3.
用Microsoft Office 就可以画出来了
而且比一般计量软件漂亮

POT模型在巨灾损失预测中的应用——基于MCMC方法的估计.pdf (312.18 KB)

2011-12-30 13:37:33 上传

epoh老师，您好！看了这帖子，感觉您对winbugs软件很了解，我最近看到这篇文献中用winbugs软件进行参数估计，所有想要利用它实现mcmc对pot模型参数估计，但是因为自己从来没有用过winbugs软件，所以不能很顺利地完成，如果您有时间能不能帮我看看呢，实在很感谢

model
        {
                for(i in 1:N)  
                {
                        p <-N/6000*tao*pow(1+cigam*tao*(x-k),-1/cigam-1)
                 }
                tao~ dgamma(0.01,0.01)
                cigam~ dpar(0.01,0.1)
         }
        data
        list(N=53,k=308,x=c(309,310,311,313,313,316,319,320,328,330,332,332,338,341,341,349,351,354,356,365,368,370,371,375,384,384,385,389,389,391,396,408,417,418,418,435,447,448,460,463,465,473,500,507,521,522,554,618,702,754,842,858,1152))
        initial values
         list(tao=1,cigam=1)

sky-xiaotian 发表于 2011-12-30 13:41
epoh老师，您好！看了这帖子，感觉您对winbugs软件很了解，我最近看到这篇文献中用winbugs软件进行参数 ...

Openbugs 有
Generalized Pareto Distribution : z~dgpar(mu,sigma,xi)
http://mathstat.helsinki.fi/open ... lSpecification.html
你给数据及mu,就可算出sigma,xi
你先安装试试
有问题再说

sky-xiaotian 发表于 2011-12-30 13:41
epoh老师，您好！看了这帖子，感觉您对winbugs软件很了解，我最近看到这篇文献中用winbugs软件进行参数 ...

R package "POT"
提供了多种估计 GPD 参数的方法
method of moments moments,
maximum likelihood mle,
biased and unbiased probability weighted moments pwmb, pwmu,
mean power density divergence mdpd,
median med,
pickands' pickands,
maximum penalized likelihood mple and
maximum goodness-of-fit mgf estimators.
我用function rgpd()
产生了一组数据
obs  : 200
mu   : 1
sigma: 2
xi     : 0.25
分别在R及openbugs做比较,结果如下,请参考.
#####in R
library(POT)
x <- rgpd(200, 1, 2, 0.25)
mom <- fitgpd(x, 1, "moments")$param
mle <- fitgpd(x, 1, "mle")$param
pwmu <- fitgpd(x, 1, "pwmu")$param
pwmb <- fitgpd(x, 1, "pwmb")$param
pickands <- fitgpd(x, 1, "pickands")$param
med <- fitgpd(x, 1, "med", start = list(scale = 2, shape = 0.25))$param
mdpd <- fitgpd(x, 1, "mdpd")$param
mple <- fitgpd(x, 1, "mple")$param
ad2r <- fitgpd(x, 1, "mgf", stat = "AD2R")$param
print(rbind(mom, mle, pwmu, pwmb, pickands, med, mdpd, mple, ad2r))

                scale     shape
mom       2.386221 0.2162865
mle         2.225754 0.2759815
pwmu     2.199347 0.2776621
pwmb     2.212050 0.2734900
pickands  2.088740 0.4066740
med       2.142015 0.4706414
mdpd     2.199591 0.2924255
mple      2.254318 0.2603019
ad2r       2.263311 0.2718832

######in openbugs
                 mean        sd              MC_error val2.5pc median val97.5pc   start        sample
sigma        2.246        0.267        0.004287 1.755   2.236          2.8       10001        10000
xi              0.294        0.0992       0.001591 0.119  0.2859  0.5055       10001        10000

######in openbugs
                 mean        sd              MC_error val2.5pc median val97.5pc   start        sample
sigma        2.246        0.267        0.004287 1.755   2.236          2.8       10001        10000
xi              0.294        0.0992       0.001591 0.119  0.2859  0.5055       10001        10000

epoh老师，您好！
以上部分的openbugs程序代码是什么？
非常感谢！