fn<-sum(log(clayton(lamda1,U))*s1+log(gumbel(lamda2,U))*s2+log(frank(lamda3,U))*s3)
}
a1[k]<-mean(a1[k-1]*clayton(theta1[k-1],U)/(a1[k-1]*clayton(theta1[k-1],U)+a2[k-1]*gumbel(theta2[k-1],U)+a3[k-1]*frank(theta3[k-1],U)))
a2[k]<-mean(a2[k-1]*gumbel(theta2[k-1],U)/(a1[k-1]*clayton(theta1[k-1],U)+a2[k-1]*gumbel(theta2[k-1],U)+a3[k-1]*frank(theta3[k-1],U)))
a3[k]<-mean(a3[k-1]*frank(theta3[k-1],U)/(a1[k-1]*clayton(theta1[k-1],U)+a2[k-1]*gumbel(theta2[k-1],U)+a3[k-1]*frank(theta3[k-1],U)))

lamda<-optim(c(theta1[k-1],theta2[k-1],theta3[k-1]),f,method="L-BFGS-B",lower=c(0,1,-Inf),upper=c(Inf,Inf,Inf),control=list(fnscale=-1))$par


运行报错 “L-BFGS-B不能有'fn'的无限值”是怎么回事呢？？求大神帮忙