Объяснение построения байесовских априорных и задних распределений на одной панели с использованием R

Question

Может ли кто-нибудь подробно объяснить, как работает этот код?

require(lattice)
?lattice  # essential reading

 data <- dgamma(seq(from=0.00001,to=0.01,by=0.00001),shape = .1, scale = .01)

dfrm <- data.frame(dgam = data, param="s.1.01")
dfrm <- rbind(dfrm, data.frame(dgam =
                                 dgamma( seq(from=0.00001,to=0.01,by=0.00001), 
                                 shape = .2, scale = .01), 
                               param="s.2.01") )
dfrm <- cbind( dfrm, X.val=seq(from=0.00001,to=0.01,by=0.00001) )
str(dfrm)
#'data.frame':  2000 obs. of  3 variables:
# $ dgam : num  5263 2817 1954 1507 1231 ...
# $ param: Factor w/ 2 levels "s.1.01","s.2.01": 1 1 1 1 1 1 1 1 1 1 ...
# $ X.val: num  1e-05 2e-05 3e-05 4e-05 5e-05 6e-05 7e-05 8e-05 9e-05 1e-04 ...
xyplot( dgam ~ X.val , 
        group=param, 
        data=dfrm, type="l")

Например, в какой части специфицируется предшествующий и задний?

Я нашел код вответ дан здесь Как построить байесовские априорные и задние распределения на одной панели, используя R?

Answer 1

Этот ответ только иллюстрирует, как построить две кривые на одном графике, но вовсе не включает вычисление истинного апостериорного значения с учетом априора и данных.А именно, он просто отображает два гамма-распределения с shape=0.1 и shape=0.2, и они различаются во фрейме данных путем включения столбца фактора param, где "s.1.01" и "s.2.01" указывают соответствующие распределения.