F(x)=(1-F(k) )G(x-k)+F(k)=1-[1+ε((x-k)/σ)]^(-1⁄ε)
ε服从Pareto（a,c）σ服从gamma(a,b)
估计参数ε，σ，在winbugs中应怎么实现啊？不会编程啊，急死人了，大家帮帮忙吧，谢谢了