R

时间:2016-10-21 15:51:59

标签: r function package

我试图找到方程的根并使用这个根,我试图通过运行10000次迭代(与统计相关)来找到一个名为“偏差”的值。

我在R中使用uniroot函数来查找根,但问题是在10000次迭代之后多次,这些值没有收敛。

我告诉我的教授这个,他告诉我在R中找到unirootuniroot.all函数的替代函数。

任何人都可以帮我在R中找到这样的替代功能吗?

编辑:我的功能如下:

  f=function(p)
  {
    g=NULL
    h=0
    for(i in 1:length(a1))
    {
      h=mat1[i,2]/dpois(mat1[i,1],p)
      g[i]=((mat1[i,1]/p)-1)*(dpois(mat1[i,1],p)^(1+alpha))*(((h^A)-1)/A)
      h=0
    }
    sum(g)
  }
  msde[k]=uniroot(f,c(1,10))$root

我需要找到uniroot的替代函数,因为msde [k]的值在不同的运行期间会显着不同。不同运行期间msde[k]的值不会收敛。

0 个答案:

没有答案