zetaff(VGAM)
zetaff()所属R语言包:VGAM
Zeta Distribution Family Function
泽塔分配家庭功能
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Estimates the parameter of the zeta distribution.
估计参数的ζ分布。
用法----------Usage----------
zetaff(link = "loge", earg = list(), init.p = NULL)
参数----------Arguments----------
参数:link, earg, init.p
See CommonVGAMffArguments for more information. These arguments apply to the (positive) parameter p. See Links for more choices. Choosing loglog constrains p>1, but may fail if the maximum likelihood estimate is less than one.
见CommonVGAMffArguments更多信息。这些论点也适用于(正)参数p。见Links更多的选择。选择loglog限制p>1,但可能会失败,如果最大似然估计是小于1的。
Details
详细信息----------Details----------
In this long tailed distribution the response must be a positive integer. The probability function for a response Y is
在这漫长的尾分布的反应必须是一个正整数。的概率函数的响应Y是
where zeta is Riemann's zeta function. The parameter p is positive, therefore a log link is the default. The mean of Y is mu = zeta(p)/zeta(p+1) (provided p>1) and these are the fitted values. The variance of Y is zeta(p-1) / zeta(p+1) - mu^2 provided p>2.
其中zeta是黎曼zeta函数。参数p是积极的,因此,log链接是默认的。 Y的平均mu = zeta(p)/zeta(p+1)(提供p>1),这些都是拟合值。 是Y的方差zeta(p-1) / zeta(p+1) - mu^2提供p>2。
It appears that good initial values are needed for successful convergence. If convergence is not obtained, try several values ranging from values near 0 to values about 10 or more.
看来,良好的初始值所需要的成功融合。如果还没有获得收敛,尝试从接近0的值至约10个或更多的值的几个值。
值----------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----------
The zeta function may be used to compute values of the zeta function.
zeta函数可被用于计算zeta函数的值。
(作者)----------Author(s)----------
T. W. Yee
参考文献----------References----------
Johnson N. L., Kemp, A. W. and Kotz S. (2005) Univariate Discrete Distributions, 3rd edition, Hoboken, New Jersey: Wiley.
Mathematical Statistics. Boca Raton: Chapman & Hall/CRC Press.
参见----------See Also----------
zeta, dzeta, hzeta, zipf.
zeta,dzeta,hzeta,zipf。
实例----------Examples----------
zdata = data.frame(y = 1:5, w = c(63, 14, 5, 1, 2)) # Knight, p.304[骑士,p.304]
fit = vglm(y ~ 1, zetaff, zdata, trace = TRUE, weight = w, crit = "c")
(phat = Coef(fit)) # 1.682557[1.682557]
with(zdata, cbind(round(dzeta(y, phat) * sum(w), 1), w))
with(zdata, weighted.mean(y, w))
fitted(fit, matrix = FALSE)
predict(fit)
# The following should be zero at the MLE:[下面列出的是零的极大似然估计:]
with(zdata, mean(log(rep(y, w))) + zeta(1+phat, deriv = 1) / zeta(1+phat))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|