r - Fiting Weibull using R2OpenBUGS -


because estimation of parameters of weibull function using r2openbugs different amounts provided generate data set using rweibull? what's wrong fit?

data<-rweibull(200, 2, 10)  model<-function(){  v ~ dgamma(0.0001,0.0001) lambda ~ dgamma(0.0001,0.0001)  for(i in 1:n){    y[i] ~ dweib(v, lambda)    } }  y<-data n<-length(y) data<-list("y", "n")  inits<-function(){list(v=1, lambda=1)} params<-c("v", "lambda") model.file<-file.path(tempdir(), "model.txt") write.model(model, model.file) weibull<-bugs(data, inits, params, model.file, n.iter = 3000, n.burnin =  2000, n.chains = 3) print(weibull, 4)

the result obtained is:

current: 3 chains, each 3000 iterations (first 2000 discarded) cumulative: n.sims = 3000 iterations saved           mean     sd       2.5%       25%       50%       75%      97.5%    rhat     n.eff v           2.0484 0.1044    1.8450    1.9780    2.0500    2.1180    2.2470 1.0062   780 lambda      0.0097 0.0026    0.0056    0.0078    0.0093    0.0112    0.0159 1.0063   830 deviance 1145.6853 1.8403 1144.0000 1144.0000 1145.0000 1146.0000 1151.0000 1.0047   770  pd = 1.6 , dic = 1147.0

r parameterizes weibull using shape (=2 in case) , scale (=10) default: bugs uses shape , lambda, lambda=(1/scale)^shape. should expect lambda approximately (1/10)^2=0.01, close median of 0.0093.

this question on crossvalidated, , this paper in r journal, compare parameterizations.


Comments

Post a Comment

Popular posts from this blog

asp.net mvc - SSO between MVCForum and Umbraco7 -

i need integrate sso between mvcforum (forum of company) , umbraco7 (manages cms company). however struggling shared authentication work. tried several ways suggested on other threads still not running expected. here did: 1. in mvcforum (forum.mywebsite.com) 1.1. machine key in web.xml: <machinekey validationkey="the same key umbraco7" decryptionkey="the same key umbraco7" validation="sha1" decryption="aes" /> 1.2. authenticate forms in web.xml <forms name="my-sso-auth" protection="all" path="/" timeout="43200" domain="forum.mywebsite.com" loginurl="/members/logon" enablecrossappredirects="true" /> 2. in umbraco7 (news.mywebsite.com) 2.1. machine key in web.xml: <machinekey validationkey="the same key mvcforum" decryptionkey="the same key mvcforum" validation="sha1" decryption="aes" /> 2.2. a
Read more

Python Tkinter keyboard using bind -

i have problem using bind option can click button in keyboard, , call function once key pressed. tried taking other codes have similar purpose , saw i'm doing pretty same, though still have problem must've missed. thankful if me this. here's code: # -*- coding: utf-8 -*- #imports tkinter import * pil import imagetk, image import time import tkmessagebox #===========================================================================================================================================# #tkinter class class app(frame): def __init__(self, master=none): frame.__init__(self, master) self.pack() #===========================================================================================================================================# #game pieces classes class board: def __init__(self, rows = 7, columns = 6, picture_height = none, picture_width = none): self.rows = rows self.columns = columns self.picture
Read more

ubuntu - Selenium Node Not Connecting to Hub, Not Opening Port -

i'm trying set demo selenium grid. have digital ocean ubuntu 14.02 droplet. (i have 2 , tried set using 2 different machines got nowhere.) i'm trying run hub , node both locally. here input both hub, node, , nmap results. https://gist.github.com/mekhami/d5bbecca5ff0de9d2dc9 so, node starting tcp connection reset , port it's supposed operate on closed. i have had no luck in number of different options , configurations work. but... of course, worked randomly once. don't know why, don't know did, know 1 brief moment, node connected , saw on grid console. tried run test , went down again , i've never seen come up. please let me know more information can provide. update: randomly came online again. don't know happened, might related timeout or time_wait or going on. this because vm/server not have enough entropy available. to resolve problem, install haveged (or other similar things). the instructions of how install haveged on ubu
Read more