nbolf(VGAM)
nbolf()所属R语言包：VGAM
Negative Binomial-Ordinal Link Function
负二项分布有序链接功能
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Computes the negative binomial-ordinal transformation, including its inverse and the first two derivatives.
计算负二项分布序的转型，包括其逆和前两个衍生物。
用法----------Usage----------
nbolf(theta, earg = stop("argument 'earg' must be given"),
inverse = FALSE, deriv = 0, short = TRUE, tag = FALSE)
参数----------Arguments----------
参数：theta
Numeric or character. See below for further details.
数字或字符。请参阅下面进一步的细节。
参数：earg
Extra argument for passing in additional information. This must be list with components cutpoint and k. Here, k is the k parameter associated with the negative binomial distribution; see negbinomial. The cutpoints should be non-negative integers. If nbolf() is used as the link function in cumulative then one should choose reverse = TRUE, parallel = TRUE, intercept.apply = TRUE.
额外的信息传递额外的参数。这必须是与组件cutpoint和k。在这里，k是k参数与负二项分布，见negbinomial。分割点应该是所有非负整数。如果nbolf()cumulative为纽带功能，那么人们应该选择reverse = TRUE, parallel = TRUE, intercept.apply = TRUE。
参数：inverse
Logical. If TRUE the inverse function is computed.
逻辑。如果TRUE的逆函数计算。
参数：deriv
Order of the derivative. Integer with value 0, 1 or 2.
订购的衍生工具。值0，1或2的整数，带。
参数：short
Used for labelling the blurb slot of a vglmff-class object.
用于标签blurb插槽的vglmff-class对象。
参数：tag
Used for labelling the linear/additive predictor in the initialize slot of a vglmff-class object. Contains a little more information if TRUE.
用于标注线性/对添加剂的预测在initialize插槽的vglmff-class对象。如果TRUE包含了一些信息。
Details
详细信息----------Details----------
The negative binomial-ordinal link function (NBOLF) can be applied to a parameter lying in the unit interval. Its purpose is to link cumulative probabilities associated with an ordinal response coming from an underlying negative binomial distribution.
负二项式序链路的功能（NBOLF）可以应用到趴在单位时间间隔的参数。它的目的是要连结与序响应来自一个潜在的负二项分布的累积概率。
The arguments short and tag are used only if theta is character.
的参数short和tag只有theta是字符。
See Links for general information about VGAM link functions.
见LinksVGAM链接功能的一般信息，。
值----------Value----------
See Yee (2012) for details.
的详细信息，请参阅仪（2012）。
警告----------Warning ----------
Prediction may not work on vglm or vgam etc. objects if this link function is used.
预测vglm或vgam的对象，如果使用这个链接功能可能无法正常工作。
注意----------Note----------
Numerical values of theta too close to 0 or 1 or out of range result in large positive or negative values, or maybe 0 depending on the arguments. Although measures have been taken to handle cases where theta is too close to 1 or 0, numerical instabilities may still arise.
数值theta太接近为0或1，或者在大的正的或负的值的范围的结果，或也许0根据参数的不同。虽然已采取措施，处理情况下，theta是太接近1或0，数值不稳定性仍可能出现。
In terms of the threshold approach with cumulative probabilities for an ordinal response this link function corresponds to the negative binomial distribution (see negbinomial) that has been recorded as an ordinal response using known cutpoints.
在累积概率为序响应这个链接功能对应的负二项分布（negbinomial），已被记录为一个有序的响应，用已知的分割点的阈值方法。
（作者）----------Author(s)----------
Thomas W. Yee
参考文献----------References----------
Ordinal ordination with normalizing link functions for count data, (in preparation).
参见----------See Also----------
Links, negbinomial, polf, golf, nbolf2, cumulative.
Links，negbinomial，polf，golf，nbolf2，cumulative。
实例----------Examples----------
earg = list(cutpoint = 2, k = 1)
nbolf("p", earg = earg, short = FALSE)
nbolf("p", earg = earg, tag = TRUE)
p = seq(0.02, 0.98, by = 0.01)
y = nbolf(p, earg = earg)
y. = nbolf(p, earg = earg, deriv = 1)
max(abs(nbolf(y, earg = earg, inv = TRUE) - p)) # Should be 0[应为0]
## Not run: par(mfrow = c(2, 1), las = 1)[＃不执行：PAR（mfrow = C（2,1），LAS = 1的）]
plot(p, y, type = "l", col = "blue", main = "nbolf()")
abline(h = 0, v = 0.5, col = "red", lty = "dashed")
plot(p, y., type = "l", col = "blue",
main = "(Reciprocal of) first NBOLF derivative")
## End(Not run)[＃（不执行）]
# Another example[另一个例子]
nn = 1000
x2 = sort(runif(nn))
x3 = runif(nn)
mymu = exp( 3 + 1 * x2 - 2 * x3)
k = 4
y1 = rnbinom(nn, mu = mymu, size = k)
cutpoints = c(-Inf, 10, 20, Inf)
cuty = Cut(y1, breaks = cutpoints)
## Not run: plot(x2, x3, col = cuty, pch = as.character(cuty)) [＃未运行图（X2，X3，列= cuty，PCH = as.character（cuty））]
table(cuty) / sum(table(cuty))
fit = vglm(cuty ~ x2 + x3, fam = cumulative(link = "nbolf",
reverse = TRUE, parallel = TRUE, intercept.apply = TRUE,
mv = TRUE, earg = list(cutpoint = cutpoints[2:3], k = k)),
trace = TRUE)
head(fit@y)
head(fitted(fit))
head(predict(fit))
coef(fit)
coef(fit, matrix = TRUE)
constraints(fit)
fit@misc$earg
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-1 15:45:31
