我试图找到方程的根并使用这个根,我试图通过运行10000次迭代(与统计相关)来找到一个名为“偏差”的值。
我在R中使用uniroot
函数来查找根,但问题是在10000次迭代之后多次,这些值没有收敛。
我告诉我的教授这个,他告诉我在R中找到uniroot
和uniroot.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]
的值不会收敛。