过程:
设置初始点x1,...,xn+1
对点进行排序,使得f(x1)≤f(x2)≤⋯≤f(xn+1)。
计算xo作为x1,...,xn的中心点。
反射
计算反射点xr=xo+α(xo-xn+1)。...样本的随机生成
#(1) beta分布n x 4)lnl(c(3, 4), x) #检验
hist(x, prob=TRUE)
拟合Beta分布
定义控制参数...plot.np=50, data=x, distr="beta", back.col = FALSE)points(unconstropt[1,"BFGS"], unconstropt[2,"BFGS"...], pch="+", col="red")points(3, 3/4, pch="x", col="green")
我们可以用bootdist函数来模拟bootstrap 复制的情况。...boot(fit(x, "beta", method="mle", optim.method="BFGS"))
plot(b1)abline(v=3, h=3/4, col="red", lwd=1.5