用R码拟合新的概率模型

我使用混合的广义技术开发了新的概率模型。 现在我希望它适合离散数据集。 但我得到了一个错误

序列默认值（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

#rm（list=ls（all=TRUE））
obs=代表（序列（0,6），c（260,87,32,4,1,0,0））
NBWED
seq（0，n）
创建一个从
0
到
n
的向量。很可能，您的
x
是一个向量或类似的东西，因此它会抛出一个错误

只需尝试：
seq（0,5）
并查看结果。这会有帮助的