R 骨架分布

R 骨架分布,r,statistics,R,Statistics,正如你们所知,Skellam包装已从CRAN中移除(不要问为什么)。经过一段时间的互联网研究,我找不到skellam pmf函数,所以我自己编写了。skellamskellam参见VGAM中的dskellam。distributions任务视图有助于查找此类函数。如果你想贡献的话,这个包缺少一个累积质量函数。这很好,但是考虑到这个函数的简单性,有什么意义呢?任何“应该”使用skellam函数的人都可以直接从维基百科中复制它。另请参见VGAM中的dskellam。distributions任务视图

正如你们所知,Skellam包装已从CRAN中移除(不要问为什么)。经过一段时间的互联网研究,我找不到skellam pmf函数,所以我自己编写了。

skellam
skellam参见
VGAM
中的
dskellam
。distributions任务视图有助于查找此类函数。如果你想贡献的话,这个包缺少一个累积质量函数。这很好,但是考虑到这个函数的简单性,有什么意义呢?任何“应该”使用skellam函数的人都可以直接从维基百科中复制它。另请参见
VGAM
中的
dskellam
。distributions任务视图有助于查找此类函数。如果你想贡献的话,这个包缺少一个累积质量函数。这很好,但是考虑到这个函数的简单性,有什么意义呢?任何“应该”使用skellam函数的人都可以直接从维基百科中复制它。当有人说“不要问为什么”时,我很想问“为什么?”。为什么?我不是克朗人,你知道;)mb它没有很好的文档记录或SMTH您总是可以从存档中获取旧版本。然后你可以问维修工杰瑞·W·刘易斯。存档上的版本为我通过了R CMD检查。并安装。当有人说“别问为什么”时,我很想问“为什么?”。为什么?我不是克朗人,你知道;)mb它没有很好的文档记录或SMTH您总是可以从存档中获取旧版本。然后你可以问维修工杰瑞·W·刘易斯。存档上的版本为我通过了R CMD检查。并安装。并运行。。。
skellam <- function(k,mu1,mu2){
  return(exp(-mu1-mu2)*((mu1/mu2)^(k/2))*besselI(2*sqrt(mu1*mu2),k))
}