invbinomial(VGAM)
invbinomial()所属R语言包:VGAM
Inverse Binomial Distribution Family Function
逆二项分布家庭功能
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Estimates the two parameters of an inverse binomial distribution by maximum likelihood estimation.
估计这两个参数的反二项式分布的最大似然估计。
用法----------Usage----------
invbinomial(lrho="elogit", llambda="loge",
erho=if(lrho=="elogit") list(min = 0.5, max = 1) else list(),
elambda=list(), irho=NULL, ilambda=NULL, zero=NULL)
参数----------Arguments----------
参数:lrho, llambda
Link function for the rho and lambda parameters. See Links for more choices.
rho和lambda参数的链接功能。见Links更多的选择。
参数:erho, elambda
List. Extra argument for each of the links. See earg in Links for general information.
列表。每个环节的额外参数。见earg中Links的一般信息。
参数:irho, ilambda
Numeric. Optional initial values for rho and lambda.
数字。可选的初始值rho和lambda。
参数:zero
See CommonVGAMffArguments.
见CommonVGAMffArguments。
Details
详细信息----------Details----------
The inverse binomial distribution of Yanagimoto (1989) has density function
柳本逆二项分布(1989)具有密度函数
where y=0,1,2,... and 0.5 < rho < 1, and lambda > 0. The first two moments exist for rho>0.5; then the mean is lambda*(1-rho)/(2*rho-1) (returned as the fitted values) and the variance is lambda*rho*(1-rho)/(2*rho-1)^3. The inverse binomial distribution is a special case of the generalized negative binomial distribution of Jain and Consul (1971). It holds that Var(Y) > E(Y) so that the inverse binomial distribution is overdispersed compared with the Poisson distribution.
在那里y=0,1,2,...和0.5 < rho < 1和lambda > 0“。前两个时刻存在rho>0.5,意思是lambda*(1-rho)/(2*rho-1)(拟合值)返回,方差为lambda*rho*(1-rho)/(2*rho-1)^3。逆二项分布是一种特殊的情况,广义负二项分布耆那教和领事(1971年)的。这认为Var(Y) > E(Y),使得逆二项式分布相比,泊松分布overdispersed。
值----------Value----------
An object of class "vglmff" (see vglmff-class). The object is used by modelling functions such as vglm and vgam.
类的一个对象"vglmff"(见vglmff-class)。该对象被用于建模功能如vglm和vgam。
注意----------Note----------
This VGAM family function only works reasonably well with intercept-only models. Good initial values are needed; if convergence failure occurs use irho and/or ilambda.
这VGAM家庭功能工作得相当好,仅截距模型。需要良好的初始值;,如果融合发生故障时使用“irho和/或ilambda。
Some elements of the working weight matrices use the expected information matrix while other elements use the observed information matrix. Yet to do: using the mean and the reciprocal of lambda results in a EIM that is diagonal.
工作的权重矩阵的某些元素的使用预期的信息矩阵,而其他元素使用的观测信息矩阵。然而,做的均值和lambda结果在EIM是对角线的倒数。
(作者)----------Author(s)----------
T. W. Yee
参考文献----------References----------
The inverse binomial distribution as a statistical model. Communications in Statistics: Theory and Methods, 18, 3625–3633.
A generalized negative binomial distribution. SIAM Journal on Applied Mathematics, 21, 501–513.
The Theory of Dispersion Models. London: Chapman & Hall
参见----------See Also----------
negbinomial, poissonff.
negbinomial,poissonff。
实例----------Examples----------
idata = data.frame(y = rnbinom(n <- 1000, mu=exp(3), size=exp(1)))
fit <- vglm(y ~ 1, invbinomial, idata, trace=TRUE)
with(idata, c(mean(y), head(fitted(fit), 1)))
summary(fit)
coef(fit, matrix=TRUE)
Coef(fit)
sum(weights(fit)) # sum of the prior weights[现有的权重的总和]
sum(weights(fit, type="w")) # sum of the working weights[工作权重的总和]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|