seq2binomial(VGAM)
seq2binomial()所属R语言包：VGAM
The Two-stage Sequential Binomial Distribution Family Function
这两个阶段的连续二项分布家庭功能
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Estimation of the probabilities of a two-stage binomial distribution.
一个两阶段的二项式分布的概率估计。
用法----------Usage----------
seq2binomial(lprob1 = "logit", lprob2 = "logit", eprob1 = list(),
eprob2 = list(), iprob1 = NULL, iprob2 = NULL, zero = NULL)
参数----------Arguments----------
参数：lprob1, lprob2
Parameter link functions applied to the two probabilities, called p and q below. See Links for more choices.
参数链接功能，被称为p和q下面的两个概率。见Links更多的选择。
参数：eprob1, eprob2
Lists. Extra arguments for the links. See earg in Links for general information.
列表。额外的参数的链接。见earg中Links的一般信息。
参数：iprob1, iprob2
Optional initial value for the first and second probabilities respectively. A NULL means a value is obtained in the initialize slot.
可选的初始值分别为第一和第二的概率。 ANULL是指initialize插槽中获得的值。
参数：zero
An integer-valued vector specifying which linear/additive predictors are modelled as intercepts only. If used, the value must be from the set {1,2} which correspond to the first and second probabilities respectively. A NULL value means none.
指定一个整数值向量线性/添加剂的预测模型仅作为拦截。如果使用的话，该值必须是从集合{1,2}，其分别对应于所述第一和第二的概率。 ANULL这个值意味着没有。
Details
详细信息----------Details----------
This VGAM family function fits the model described by Crowder and Sweeting (1989) which is described as follows. Each of m spores has a probability p of germinating. Of the y1 spores that germinate, each has a probability q of bending in a particular direction. Let y2 be the number that bend in the specified direction. The probability model for this data is P(y1,y2) =
该VGAM家庭功能适合Crowder和斯威汀（1989）被描述为如下所描述的模型。 m孢子的概率p发芽的。的y1孢子萌发，每个有一个概率q在一个特定的方向弯曲。让我们y2是，在指定的方向弯曲。此数据是基于概率模型P(y1,y2) =
for 0 < p < 1, 0 < q < 1, y1=1,…,m and y2=1,…,y1. Here, p is prob1, q is prob2.
0 < p < 1，0 < q < 1，y1=1,…,m和y2=1,…,y1。在这里，p是prob1，q是prob2。
Although the Authors refer to this as the bivariate binomial model, I have named it the (two-stage) sequential binomial model. Fisher scoring is used.
虽然作者称此为二元二项式模型，我把它命名为（两阶段）顺序二项式模型。使用Fisher评分。
值----------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 response must be a two-column matrix of sample proportions corresponding to y1 and y2. The m values should be inputted with the weights argument of vglm and vgam. The fitted value is a two-column matrix of estimated probabilities p and q. A common form of error is when there are no trials for y1, e.g., if mvector below has some values which are zero.
响应必须是一个两列的矩阵对应的样本比例y1和y2。 m值应输入weightsvglm和vgam参数。的拟合值，估计概率p和q是一个两列的矩阵。一种常见的错误是，当有没有y1，例如试验，如果mvector有一定的值是零。
（作者）----------Author(s)----------
Thomas W. Yee
参考文献----------References----------
Bayesian inference for a bivariate binomial distribution. Biometrika, 76, 599–603.
参见----------See Also----------
binomialff.
binomialff。
实例----------Examples----------
sdata = data.frame(mvector = round(rnorm(nn <- 100, m = 10, sd = 2)),
x = runif(nn))
sdata = transform(sdata, prob1 = logit(+2 - x, inverse = TRUE),
prob2 = logit(-2 + x, inverse = TRUE))
sdata = transform(sdata, successes1 = rbinom(nn, size=mvector, prob=prob1))
sdata = transform(sdata, successes2 = rbinom(nn, size=successes1, prob=prob2))
sdata = transform(sdata, y1 = successes1 / mvector)
sdata = transform(sdata, y2 = successes2 / successes1)
fit = vglm(cbind(y1,y2) ~ x, seq2binomial, weight=mvector,
data = sdata, trace=TRUE)
coef(fit)
coef(fit, matrix = TRUE)
head(fitted(fit))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-1 15:52:53
