不知道有没更简单的，自己编了个，先安装一个quantmod包
library(quantmod)
x=c(1,2,3,2,1,-1,2,5,3,7)
(X=mean(x))#计算均值X及峰值K
(K=findPeaks(x))# 计算峰值
while (length(K[K>=3])>0)
{ (x=x[-which(abs((x-X))==max(abs(x-X)))])
K=findPeaks(x)
X=mean(x)}
#对K进行判断，若K>=3,则选择|Xi-X|最大的Xi，将其从样本中除去,重复以上操作,直到峰值小于3
max(x)
#从剩余的样本中输出最大的Xi

运行过了，出现下边的错误，还有，我刚开始写错了，K应该是峰度，不好意思
> while (length(K[K>=3])>0)
+ {(zt=zt[-which(abs((zt-X))==max(abs(zt-X)))]
+ K=kurtosis(zt),X=mean(zt)}
错误: unexpected symbol in:
"{(zt=zt[-which(abs((zt-X))==max(abs(zt-X)))]
K"

是我的问题，少写个括号，现在出来结果了