dirmul.old(VGAM)
dirmul.old()所属R语言包：VGAM
Fitting a Dirichlet-Multinomial Distribution
一类Dirichlet多项分布拟合
译者：生物统计家园网 机器人LoveR
----------Description----------
Fits a Dirichlet-multinomial distribution to a matrix of non-negative integers.
适合的Dirichlet多项式分布矩阵的非负整数。
用法----------Usage----------
dirmul.old(link = "loge", earg = list(),
init.alpha = 0.01, parallel = FALSE, zero = NULL)
参数----------Arguments----------
参数：link
Link function applied to each of the M (positive) shape parameters alpha_j for j=1,…,M. See Links for more choices. Here, M is the number of columns of the response matrix.
Link功能M的alpha_jj=1,…,M（正）形状参数。见Links更多的选择。在这里，M是响应矩阵的列的数量。
参数：earg
List. Extra argument for link. See earg in Links for general information.
列表。额外参数link。见earg中Links的一般信息。
参数：init.alpha
Numeric vector. Initial values for the alpha vector. Must be positive. Recycled to length M.
数字矢量。初始值为alpha向量。必须为正数。回收长度M。
参数：parallel
A logical, or formula specifying which terms have equal/unequal coefficients.
一个逻辑，或指定的条款有平等/不平等系数，公式。
参数：zero
An integer-valued vector specifying which linear/additive predictors are modelled as intercepts only. The values must be from the set {1,2,...,M}.
指定一个整数值向量线性/添加剂的预测模型仅作为拦截。这些值必须是集合{1,2，...，M}。
Details
详细信息----------Details----------
The Dirichlet-multinomial distribution, which is somewhat similar to a Dirichlet distribution, has probability function
的Dirichlet多项式分布，这是有点类似的Dirichlet分布的概率函数
for alpha_j > 0, alpha_+ = alpha_1 + \cdots + alpha_M, and 2y_* = y_1 + \cdots + y_M. Here, C_b^a means “a choose b” and refers to combinations (see choose). The (posterior) mean is
alpha_j > 0，alpha_+ = alpha_1 + \cdots + alpha_M和2y_* = y_1 + \cdots + y_M。在这里，C_b^a的意思是“a选择b”是指组合（见choose）。 （后）平均
for j=1,…,M, and these are returned as the fitted values as a M-column matrix.
j=1,…,M，这些传回的拟合值作为M列的矩阵。
值----------Value----------
An object of class "vglmff" (see vglmff-class). The object is used by modelling functions such as vglm, rrvglm and vgam.
类的一个对象"vglmff"（见vglmff-class）。该对象被用于建模功能，如vglm，rrvglm和vgam。
注意----------Note----------
The response should be a matrix of non-negative values. Convergence seems to slow down if there are zero values. Currently, initial values can be improved upon.
的反应应该是一个非负的值的矩阵。收敛似乎放慢是否存在零值。目前，初始值可以改进。
This function is almost defunct and may be withdrawn soon. Use dirmultinomial instead.
此功能几乎已不存在，并可能被撤销。使用dirmultinomial“。
----------Author(s)----------
Thomas W. Yee
参考文献----------References----------
Mathematical and Statistical Methods for Genetic Analysis, 2nd ed. New York: Springer-Verlag.
Statistical Distributions, New York: Wiley-Interscience, Third edition.
Fisher information matrix of the Dirichlet-multinomial distribution. Biometrical Journal, 47, 230–236.
Overdispersion in allelic counts and <code>θ</code>-correction in forensic genetics. Theoretical Population Biology, 78, 200–210.
http://www.stat.auckland.ac.nz/~yee contains further information and examples.
参见----------See Also----------
dirmultinomial, dirichlet, betabinomial.ab, multinomial.
dirmultinomial，dirichlet，betabinomial.ab，multinomial。
实例----------Examples----------
# Data from p.50 of Lange (2002)[从第50页的数据兰格（2002）]
alleleCounts <- c(2, 84, 59, 41, 53, 131, 2, 0,
0, 50, 137, 78, 54, 51, 0, 0,
0, 80, 128, 26, 55, 95, 0, 0,
0, 16, 40, 8, 68, 14, 7, 1)
dim(alleleCounts) <- c(8, 4)
alleleCounts <- data.frame(t(alleleCounts))
dimnames(alleleCounts) <- list(c("White","Black","Chicano","Asian"),
paste("Allele", 5:12, sep = ""))
set.seed(123) # @initialize uses random numbers[@初始化使用随机数]
fit <- vglm(cbind(Allele5,Allele6,Allele7,Allele8,Allele9,
Allele10,Allele11,Allele12) ~ 1, dirmul.old,
trace = TRUE, crit = "c", data = alleleCounts)
(sfit <- summary(fit))
vcov(sfit)
round(eta2theta(coef(fit), fit@misc$link), dig = 2) # not preferred[不喜欢]
round(Coef(fit), dig = 2) # preferred # preferred[首选＃首选]
round(t(fitted(fit)), dig = 4) # 2nd row of Table 3.5 of Lange (2002)[第二排的朗格表3.5（2002年）]
coef(fit, matrix = TRUE)
pfit <- vglm(cbind(Allele5,Allele6,Allele7,Allele8,Allele9,
Allele10,Allele11,Allele12) ~ 1,
dirmul.old(parallel = TRUE), trace = TRUE,
data = alleleCounts)
round(eta2theta(coef(pfit), pfit@misc$link), dig = 2) # not preferred[不喜欢]
round(Coef(pfit), dig = 2) # preferred[首选]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-1 15:32:21
