weight=c(11260,15960,170304,17220,19470,20900,21620,21950,41410, 50760,53060,58600,65150,75130,85470,86240, 98580,130690,143630,153310, 158330)

LL<-function(params,weight)
{
f=exp(-exp(-(data-params[1])/params[2])-(data-params[1])/params[2])/params[2]
ll<-sum(log(f))
return(-ll)
}


params<-c(30000,40000)
geyser.res<-optim(c(30000,40000),LL, hessian = TRUE,method= "Nelder-Mead",control=list(fnscale=-1))
运行后：Error in data - params[1] : non-numeric argument to binary operator
请问这是为什么，是哪出错了，我就是想用这种估计的方式求出其系数，因为后还想做其它的 加权拟合 ，极值与对数正交分布的加权，请大神赐教。。