prof.term(gamlss)
prof.term()所属R语言包：gamlss
Plotting the Profile: deviance or information criterion for one of the terms (or hyper-parameters) in a GAMLSS model
在GAMLSS模型绘制的个人资料：偏差或信息标准的条款（或超参数）
译者：生物统计家园网 机器人LoveR
描述----------Description----------
This functions plots the profile deviance for a chosen parameter included in the linear predictor of any of the mu, sigma, nu or tau models so profile confidence intervals can be obtained. In can also be used to plot the profile of a specified information criterion for any hyperparameter.
该函数绘制的档案偏差的线性预测mu，sigma，nu或tau模型，使个人资料的置信区间，可以得到包含在所选的参数。也可以用在，绘制个人资料的任何超参数指定的信息标准。
用法----------Usage----------
prof.term(model = NULL, criterion = "GD", penalty = 2.5, other = NULL,
min = NULL, max = NULL, step = NULL, type = "o", xlabel = NULL,
plot = TRUE, term = TRUE, perc = 95, ... )
参数----------Arguments----------
参数：model
this is a GAMLSS model, e.g. <br> model=gamlss(y~cs(x,df=this),sigma.fo=~cs(x,df=3),data=abdom), where this indicates the (hyper)parameter to be profiled
这是一个GAMLSS的模型，例如参考model=gamlss(y~cs(x,df=this),sigma.fo=~cs(x,df=3),data=abdom)，this表示（超）参数要分析
参数：criterion
whether global deviance ("GD") or information criterion ("IC") is profiled. The default is criterion="GD"
无论是在全球越轨行为（“GD”）或信息准则（“IC”）是异形。默认的是criterion="GD"
参数：penalty
The penalty value if information criterion is used in criterion, default penalty=2.5
处罚信息标准值，如果在criterion，默认penalty=2.5
参数：other
this can be used to evaluate an expression before the actual fitting of the model
这可以用来计算一个表达式前的实际模型的拟合
参数：min
the minimum value for the parameter e.g. min=1
的最小值的参数例如min=1
参数：max
the maximum value for the parameter e.g. max=20
的最大值的参数例如max=20
参数：step
how often to evaluate the global deviance (defines the step length of the grid for the parameter) e.g. step=1
如何经常评估全球越轨（定义网格的步长为参数），例如step=1
参数：type
what type of plot required. This is the same as in type for plot, default value is type="o", that is, both line and points
需要什么类型的图。这是相同的type的plot，默认值是type="o"，就是两条线和点
参数：xlabel
if a label for the axis is required
如果需要用于该轴的标签
参数：plot
whether to plot, plot=TRUE or save the results, plot=FALSE
是否绘制，plot=TRUE或保存结果，plot=FALSE
参数：term
this has the value TRUE and it should be changed to FALSE if a profile global deviance is required for a hyperparameter so the IC are suppressed
此值为TRUE时，它应该被更改为false，如果一个高调的全球性越轨行为的超参数，所以需要抑制的IC
参数：perc
what % confidence interval is required
需要什么％的置信区间
参数：...
for extra arguments
额外的参数
Details
详细信息----------Details----------
This function can be use to provide likelihood based confidence intervals for a parameter involved in terms in the linear predictor(s). These confidence intervals are more accurate than the ones obtained from the parameters' standard errors. The function can also be used to plot a profile information criterion (with a given penalty) against a hyperparameter. This can be used to check the uniqueness in hyperparameter determination using for example find.df.
这个功能可以提供可能性的置信区间为参数的线性预测（S）的使用。这些置信区间参数的标准误差比得到更准确的。针对超参数的功能也可以用来绘制配置文件信息准则（与给定的刑罚）。这可以用来检查的独特性，例如find.df使用的超参数确定。
值----------Value----------
Return a profile plot (if the argument plot=TRUE) or the values of the parameters and the IC or GD values otherwise
如果参数plot=TRUE）的参数值和IC或GD值，否则返回一个剖面图（
警告----------Warning ----------
A dense grid (i.e. small step) evaluation of the global deviance can take a long time, so start with a sparse grid (i.e. large step)
一个密集的网格（即小步骤）评估的全球偏差，可能需要很长的时间，所以开始用稀疏网格（即大步骤）
（作者）----------Author(s)----------
Mikis Stasinopoulos <a href="mailto:d.stasinopoulos@londonmet.ac.uk">d.stasinopoulos@londonmet.ac.uk</a> and Bob Rigby <a href="mailto:r.rigby@londonmet.ac.uk">r.rigby@londonmet.ac.uk</a>
参考文献----------References----------
Appl. Statist., 54, part 3, pp 507-554.
Accompanying documentation in the current GAMLSS help files, (see also http://www.gamlss.org/).
Journal of Statistical Software, Vol. 23, Issue 7, Dec 2007, http://www.jstatsoft.org/v23/i07.
参见----------See Also----------
gamlss, prof.dev
gamlss，prof.dev
实例----------Examples----------
data(aids)
gamlss(y~x+qrt,family=NBI,data=aids)
mod<-quote(gamlss(y ~ offset(this * x) + qrt, data = aids, family = NBI))
prof.term(mod, min=0.06, max=0.11, step=0.001)
mod1<-quote(gamlss(y ~ cs(x,df=this) + qrt, data = aids, family = NBI))
prof.term(mod1, min=1, max=15, step=1, criterion="IC")
mod2 <- quote(gamlss(y ~ x+I((x>this)*(x-this))+qrt,family=NBI,data=aids))
prof.term(mod2, min=1, max=45, step=1, criterion="GD")
rm(mod,mod1,mod2)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-18 07:53:05
