作者celestialgod (攸蓝)
看板Statistics
标题Re: [程式] R用数值分析方法找方程式的根
时间Sat Jun 8 16:39:14 2013
前面恕删
X_i服从ZIP(π,θ),π和θ的MLE分别为
_ θhat n-n_0
x = ------------- ; πhat = ------------------
1-exp(-θhat) n*(1-exp(-θhat))
其中n_0: x_i for i=1,...,n 中0的个数
_
x : (不为0的x_i的总合)/(n-n_0)
[程式修正]:
n = 2000
p = c(0.3,2)
y = sapply(1:n,function(v) ifelse(rbinom(1,1,p[1])==0, 0, rpois(1,p[2]))
mean.num0 = mean(y[y!=0])
eq = function(x) x/(1-exp(-x))-mean.num0
eq1 = function(x) (1-exp(-x)-x*exp(-x))/((1-exp(-x))^2)
theta = 3
tol = 1e-6
repeat{
theta.hat = theta - eq(theta) / eq1(theta)
eps = abs(theta.hat - theta)
if(eps <= tol) break
theta = theta.hat
}
pi.hat = length(y[y!=0]) / (1-exp(-theta.hat)) / length(y)
嵇老师在你後面
※ 编辑: celestialgod 来自: 218.166.8.214 (06/08 20:08)
1F：推 tokyo291:老师在4楼 06/10 17:17
2F：→ celestialgod:楼上幽默喔XDD 06/10 19:09