[英]New Probability Model fitting using R Code
我使用广义混合技术开发了新的概率模型。 现在我希望它适合离散数据集。 但我得到的错误是
seq.default(0, x) 中的错误:'to' 的长度必须为 1 我不明白,如何处理这个问题。 R代码如下:
# rm(list=ls(all=TRUE))
obs=rep(seq(0,6),c(260,87,32,4,1,0,0))
NBWED<-function(x,r,alpha,beta){
j=seq(0,x)
C=function(n,x){
factorial(n)/(factorial(n-x)*factorial(x))
}
C(x+r-1,x)*sum(C(x,j)*(-1)^j*(alpha^2/(alpha+beta))*((r+j+alpha)+beta)/(r+j+alpha)^2)
}
library(MASS)
fit09=fitdistr(x = obs,densfun = NBWED,start = list(r=1,alpha=0.5,beta=9.4),lower = list(a = 0.1,0.001,0.001),upper=c(Inf,Inf,Inf))
fit09
seq(0,n)
创建一个从0
到n
的向量。 可能你的x
是一个向量或类似的东西,因此它会引发错误。
只需尝试: seq(0,5)
并查看结果。 它会有所帮助。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.