Fatal编程技术网
  • r/
  • optim中的错误（par=Tm1，fn=.logLgam，x=x）：对象'；Tm1和x27；找不到

optim中的错误（par=Tm1，fn=.logLgam，x=x）：对象'；Tm1和x27；找不到

r statistics

optim中的错误（par=Tm1，fn=.logLgam，x=x）：对象'；Tm1和x27；找不到,r,statistics,R,Statistics,我正在尝试在R中运行nsRFA库的MSClaio2008功能。根据文档，我已成功运行网站上给出的示例。但是，当我创建自己的数据集以使用MSClaio2008函数运行时，如下所示： library(nsRFA) data <- c(91,84,42,30,66,95,65,12,27,61,31,101,48,52,53,55,80,23,87,46,50,33,75,88,54,17,57,39,10,89,59,24,11,43,13,93,105,28,104,18,103,62,22

我正在尝试在R中运行
nsRFA
库的
MSClaio2008
功能。根据文档，我已成功运行网站上给出的示例。但是，当我创建自己的数据集以使用
MSClaio2008
函数运行时，如下所示：

library(nsRFA)
data <- c(91,84,42,30,66,95,65,12,27,61,31,101,48,52,53,55,80,23,87,46,50,33,75,88,54,17,57,39,10,89,59,24,11,43,13,93,105,28,104,18,103,62,22,58,15,34,74,51,97,44,99,76,14,16,109,92,110,40,78,83,60,49,96,36,32,81,68,20,56,25,63,47,37,29,100,71,106,41,90,70,85,38,19,108,73,102,26,82,98,45,77,35,94,79,86,72,107,21,67,69,64)
MSC <- MSClaio2008(data)

库（nsRFA）

data这看起来像是包中的一个bug，但它与您的数据很奇怪这一事实有关。这个函数被认为是将分布拟合到水文极值的观测值：你给它一个从10到110的整数序列

事实证明，当计算某些可能分布的拟合时（
P3
和
GEV
，我认为），函数会在内部计算数据的偏度，并测试偏度为正或负时该做什么。它不考虑偏度将正好是0的可能性！
内部
nsRFA:：：ML\u估算
：


skx0）{

Tm1向我们展示所有的代码。请，这就是我的代码工作的全部内容。这是他们在网站上给出的示例。非常感谢编辑和回复。你的答案对我有效，我从解释中得到了很多。但我看到了一些奇怪的东西。我的数据长度是101，当我将其减少到100时，它工作得非常好。我不知道为什么，但我知道t有效。因为您创建的数据集不是完全对称/不倾斜的。如果这解决了您的问题，建议您单击复选标记以接受它