牛顿法求解柯西分布参数的最大似然估计（P055）

1.牛顿法求方程根

假设f(x)在x=a处取到0，即f(a)=0。已知c在a的附近，有a=c+Δx。

根据泰勒公式展开，

省略高次项，得

2.简单代码实现

newt_cauchy.R

    

cauchy_mle=function(x,theta){
	sum = 0;
	for(i in 1:length(x)){
		sum = sum+(x[i]-theta)/(1+(x[i]-theta)^2);
	}
	sum;
}
cauchy_mle_der=function(x,theta){
	sum = 0;
	tmp = 0;
	for(i in 1:length(x)){
		tmp