cennormal1(VGAM)
cennormal1()所属R语言包：VGAM
Censored Normal Distribution
定数截尾正态分布
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Maximum likelihood estimation for the normal distribution with left and right censoring.
左，右设限的正态分布的最大似然估计。
用法----------Usage----------
cennormal1(lmu = "identity", lsd = "loge",
emu = list(), esd = list(), imethod = 1, zero = 2)
参数----------Arguments----------
参数：lmu, lsd, emu, esd
Parameter link functions and its extra arguments, applied to the mean and standard deviation parameters. See Links for more choices. The standard deviation is a positive quantity, therefore a log link is the default.
参数链接的功能和额外的参数，应用的均值和标准差的参数。见Links更多的选择。标准差是一个积极的数量，因此，log链接是默认的。
参数：imethod
Initialization method. Either 1 or 2, this specifies two methods for obtaining initial values for the parameters.
初始化方法。为1或2，这指定了两种方法获得的参数的初始值。
参数：zero
An integer vector, containing the value 1 or 2. If so, the mean or standard deviation respectively are modelled as an intercept only. Setting zero = NULL means both linear/additive predictors are modelled as functions of the explanatory variables.
，含有1或2的值的整数向量。如果是这样，平均值或标准偏差分别参照仅作为截距。设置zero = NULL是指线性/添加剂的预测建模为解释变量的函数。
Details
详细信息----------Details----------
This function is like normal1 but handles observations that are left-censored (so that the true value would be less than the observed value) else right-censored (so that the true value would be greater than the observed value). To indicate which type of censoring, input extra = list(leftcensored = vec1, rightcensored = vec2) where vec1 and vec2 are logical vectors the same length as the response. If the two components of this list are missing then the logical values are taken to be FALSE. The fitted object has these two components stored in the extra slot.
此功能是一样normal1但处理左审查的意见（这样的真正价值将低于观测值），其他合适的审查（这样的真正价值将大于观测值）。为了表明哪种类型的审查，输入extra = list(leftcensored = vec1, rightcensored = vec2)其中vec1和vec2是逻辑的响应向量的长度相同。如果此列表中的两个组成部分丢失，则对应的逻辑值分别取为FALSE。设备齐全的对象有这两个组件存储在extra插槽。
值----------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 function is an alternative to tobit but cannot handle a matrix response and uses different working weights. If there are no censored observations then normal1 is recommended instead.
此功能的替代tobit但不能处理矩阵响应，并使用不同的权重。如果有没有审查意见，然后normal1建议。
（作者）----------Author(s)----------
T. W. Yee
参见----------See Also----------
tobit, normal1, dcennormal1.
tobit，normal1，dcennormal1。
实例----------Examples----------
cdata = data.frame(x2 = runif(nn <- 1000)) # ystar are true values[ystar是真正的价值]
cdata = transform(cdata, ystar = rnorm(nn, m = 100 + 15 * x2, sd = exp(3)))
## Not run: with(cdata, hist(ystar))[＃不运行：与（CDATA ystar），hist()]
cdata = transform(cdata, L = runif(nn, 80, 90), # Lower censoring points[下审查]
U = runif(nn, 130, 140)) # Upper censoring points[上审查点]
cdata = transform(cdata, y = pmax(L, ystar)) # Left censored[左审查]
cdata = transform(cdata, y = pmin(U, y)) # Right censored[右删失]
## Not run: with(cdata, hist(y))[＃不运行：（CDATA，历史（Y））]
Extra = list(leftcensored = with(cdata, ystar < L),
rightcensored = with(cdata, ystar > U))
fit1 = vglm(y ~ x2, cennormal1, cdata, crit = "c", extra = Extra, trace = TRUE)
fit2 = vglm(y ~ x2, tobit(Lower = with(cdata, L), Upper = with(cdata, U)),
cdata, crit = "c", trace = TRUE)
coef(fit1, matrix = TRUE)
max(abs(coef(fit1, matrix = TRUE) - coef(fit2, matrix = TRUE))) # Should be 0[应为0]
names(fit1@extra)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-10-1 15:29:15
