R语言 SMPracticals包 exp.gibbs()函数中文帮助文档(中英文对照)

loveR

exp.gibbs(SMPracticals)
exp.gibbs()所属R语言包：SMPracticals

                                        Gibbs Sampling for Two Truncated Exponential Variables, Practical 11.3
                                         Gibbs采样截断指数变量，实用11.3

                                         译者：生物统计家园网 机器人LoveR

描述----------Description----------

Performs Gibbs sampling for problem with two truncated exponential variables. See Practical 11.3 of Davison (2003) for details.
执行Gibbs抽样的问题有两个截断指数变量。有关详细信息，看到实际的11.3戴维森（2003年）。


用法----------Usage----------


exp.gibbs(u1 = NULL, u2 = NULL, B, I = 100, S = 100)



参数----------Arguments----------

参数：u1
Initial values for variable 1  
1变量的初始值


参数：u2
Initial values for variable 2
2变量的初始值


参数：B
Value at which exponential distribution is truncated
价值指数分布被截断


参数：I
Number of iterations of sampler
取样器的迭代数


参数：S
Number of replicates of sampler
重复采样


Details

详细信息----------Details----------

This is provided simply so that  readers spend less time typing. It is not intended to be robust and general code.
这是简单地让读者花更少的时间打字。它的目的不是是坚固和一般的代码。


值----------Value----------

A 2 x S x I array containing the values of the variables for the successive iterations
A 2个S x我的连续迭代的数组，包含的变量的值


（作者）----------Author(s)----------


Anthony Davison (<code>anthony.davison@epfl.ch</code>)



参考文献----------References----------



实例----------Examples----------


add.exp.lines <- function( exp.out, i, B=10)
{
  dexp.trunc <- function( u, lambda, B )
     dexp(u, rate=lambda)/(1-exp(-lambda*B))
  S <- dim(exp.out)[2]
  I <- dim(exp.out)[3]
  u <- seq(0.0001,B,length=1000)
  fu <- rep(0,1000)
  for (s in 1:S) fu <- fu + dexp.trunc(u,exp.out[3-i,s,I],B)/S
  lines(u,fu,col="red")
  invisible()
}
par(mfrow=c(3,2))
B <-10; I <- 15; S <- 500
exp.out <- exp.gibbs(B=B,I=I,S=S)
hist(exp.out[1,,I],prob=TRUE,nclass=15,xlab="u1",ylab="PDF",xlim=c(0,B),ylim=c(0,1))
add.exp.lines(exp.out,1)
hist(exp.out[2,,I],prob=TRUE,nclass=15,xlab="u2",ylab="PDF",xlim=c(0,B),ylim=c(0,1))
add.exp.lines(exp.out,2)

转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。


注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。