你的上一个帖子里我给你写了一个求Hill统计量的函数，这里借助函数Hill()来求：
<br />
Hill=function(x,k){<br />
ox=sort(x);<br />
n=length(x);<br />
Hill=sum(log(ox[(n-k+1):n]))/k-log(ox[n-k]);<br />
Hill<br />
}</p>
<p>y=numeric((length(x)-1))<br />
for(i in 1:(length(x)-1)){<br />
y[i]=Hill(x,i)<br />
}<br />
which.min(y)<br />
</p>
