cred.set(dks)
cred.set()所属R语言包：dks
Calculate a credible set for the posterior distribution on the Beta hyperparameters.
计算后验分布为测试hyperparameters的可信。
译者：生物统计家园网 机器人LoveR
描述----------Description----------
This function accepts a distribution calculated with pprob.dist and calculates a credible set of the specified level for the hyperparameters. If the credible set includes the value (1,1) the sample is likely to be uniform.
这个函数接受分布与pprob.dist计算，并计算出一套可信的为hyperparameters的指定级别。如果可信集包含的值（1,1），样本可能要均匀。
用法----------Usage----------
cred.set(dist,delta=NULL,level=0.95)
参数----------Arguments----------
参数：dist
The posterior distribution for the hyperparameters computed with pprob.dist.
与pprob.dist计算的hyperparameters，后验分布。
参数：delta
The grid size, must match the grid size from pprob.dist.
网格尺寸，必须符合从pprob.dist电网大小。
参数：level
The level of the credible set.
可信集的水平。
Details
详情----------Details----------
The cred.set function calculates a credible set of the specified level based on the distribution calculated with pprob.dist. The grid size, delta, should match the grid size from the call to pprob.dist. The result is a matrix of the same size as dist which indicates whether each point is in the credible set.
cred.set功能计算出一套可信的指定水平上分布计算与pprob.dist的基础。网格的大小，Delta，网格的大小相匹配从通话pprob.dist。其结果是一个区相同大小的矩阵，这表明每个点是否是可信的一套。
值----------Value----------
参数：cred
The credible set for the hyper-parameters of the beta distribution.
β分布的超参数的可信集。
参数：level
The user specified level of the set.
用户指定集的水平。
参数：elevel
The empirical level of the set, the smaller delta is, the closer elevel will be to level.
集的经验水平，较小的增量将接近elevel水平。
作者（S）----------Author(s)----------
Jeffrey T. Leek <a href="mailto:jleek@jhsph.edu">jleek@jhsph.edu</a>
参考文献----------References----------
参见----------See Also----------
dks, dks.pvalue, pprob.dist,cred.set
dks，dks.pvalue，pprob.dist，cred.set
举例----------Examples----------
## Load data[＃将数据]
data(dksdata)
## Calculate the posterior distribution[＃计算后验分布]
dist1 <- pprob.dist(P[,1])
delta = 0.1
## Calculate a 95% credible set[＃计算95％的可信集。]
cred1 <- cred.set(dist1,delta=0.1)
## Plot the posterior and the credible set[＃绘制后和可信的一套。]
alpha <- seq(0.1,10,by=delta)
beta <- seq(0.1,10,by=delta)
par(mfrow=c(1,2))
image(log10(alpha),log10(beta),dist1,xaxt="n",yaxt="n",xlab="Alpha",ylab="Beta")
axis(1,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
axis(2,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
points(0,0,col="blue",cex=1,pch=19)
image(log10(alpha),log10(beta),cred1$cred,xaxt="n",yaxt="n",xlab="Alpha",ylab="Beta")
axis(1,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
axis(2,at=c(-2,-1,0,1,2),labels=c("10^-2","10^-1","10^0","10^1","10^2"))
points(0,0,col="blue",cex=1,pch=19)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-2-25 16:40:49
