R是否存在数值稳定性问题？

Question

我正在努力在R中编写一个排名模型的模拟。我似乎遇到了一些可能与数字稳定性有关的事情＆＃34;。我想知道是否有人可以建议一些解决方法或方法来处理它。

ranking_algo_fn<- function(old_rank, website_choice,kappa)
 {  
   v<- (kappa/ (1+ kappa))
   new_rank<- (v*old_rank)+((1-v)*website_choice)
   return(new_rank)
 }

我已将此问题追溯到此功能。基本上这是一个更新函数，它迭代执行某些操作的代理 - website choice

等级是概率的向量（总和为1 - 应该如此）。从均匀分布开始 - 如果有20个网站，则每个网站为0.05。即[0.05,0.05,0.05,0.05,0.05，....]

问题在于，直观地说，每个后续的new_rank应该总和为1，但它并不是。绘制它基本上需要它自己的生命，（不同的行为基于不同的参数，但从来没有它应该）。有关此类实例，请参见附图。

website_choice和以前所有其他成分之和的总和，总是一条直线。任何解决方法？

我被告知使用MATLAB或Mathematica等其他软件来做这类事情会更好。这是真的吗？