原因可能在内不在外：loss中有可能有Inf值

但是我的数据是自己输入的我觉得不可能存在缺失值，而且现在又出现了Error in if (min(y) <= 0) l <- 10^6 else { :
  需要TRUE/FALSE值的地方不可以用缺少值。能帮我看下这是什么原因，万分感谢！

最后那个是因为你的取值中出现了负值，我之前也是，因为取了对数，所以在0-1直接的数就变成了负数，跟你的错误一样，所以要剔除掉负数就可以了

可以把用到的两个函数拿出来简单改写， 请参考
#问题一，mrl.plot, 这个函数写的比较稳健，直接把源代码拿出来绘图


## 问题二：在gdp.fit中调用了optim，在某些情况下无最优质，思路屏蔽掉无最优值的项，代码如下