- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在从 JAGS 迁移到 LaplacesDemon 并尝试重写我的一些代码。我读过 LaplacesDemon Tutorial和 LaplacesDemon Examples小插图,我对小插图中的一些例子有点困惑。
在LaplacesDemon Tutorial (p.5)中的简单例子中,模型写成:
Model <- function(parm, Data)
{beta <- parm[Data$pos.beta]
sigma <- interval(parm[Data$pos.sigma], 1e-100, Inf)
parm[Data$pos.sigma] <- sigma
beta.prior <- dnormv(beta, 0, 1000, log=TRUE)
sigma.prior <- dhalfcauchy(sigma, 25, log=TRUE)
mu <- tcrossprod(beta, Data$X)
LL <- sum(dnorm(Data$y, mu, sigma, log=TRUE))
LP <- LL + sum(beta.prior) + sigma.prior
Modelout <- list(LP=LP, Dev=-2*LL, Monitor=LP,
yhat=rnorm(length(mu), mu, sigma), parm=parm)
return(Modelout)}
此处,beta.prior
是针对 LP
求和的,因为有多个 beta 参数。
但我发现在 LaplacesDemon Example vignette 中更高级的示例中,它似乎并不总是遵循规则。如示例 87 (p.162):
Model <- function(parm, Data)
{### Log-Prior
beta.prior <- sum(dnormv(beta[,1], 0, 1000, log=TRUE), dnorm(beta[,-1], beta[,-Data$T], matrix(tau, Data$K, Data$T-1), log=TRUE))
zeta.prior <- dmvn(zeta, rep(0,Data$S), Sigma[ , , 1], log=TRUE)
phi.prior <- sum(dhalfnorm(phi[1], sqrt(1000), log=TRUE), dtrunc(phi[-1], "norm", a=0, b=Inf, mean=phi[-Data$T], sd=sigma[2], log=TRUE))
### Log-Posterior
LP <- LL + beta.prior + zeta.prior + sum(phi.prior) + sum(kappa.prior) + sum(lambda.prior) + sigma.prior + tau.prior
Modelout <- list(LP=LP, Dev=-2*LL, Monitor=LP, yhat=rnorm(prod(dim(mu)), mu, sigma[1]), parm=parm)
return(Modelout)}
(由于示例代码较长,只放了部分代码)
这里,zeta
不止一个,但没有在 Log-Prior
或 Log-Posterior
部分求和,beta
不止一个,在 Log-Prior
中求和,phi
也不止一个参数,但在 Log-Prior 中求和Prior
和 Log-Posterior
部分。
而在第 167 页的下一个示例中,它似乎又有所不同。
我想知道在什么情况下我们应该对先验密度求和?非常感谢!
最佳答案
您是否尝试过逐行运行代码?您会了解到,没有什么可以求和的,因为 dmvn
是多元正态分布的密度函数,它返回一个值——观察向量 zeta
的概率密度。所有总和的原因是为了获得一起观察两个独立事件的概率,我们将它们的边际概率相乘(或求和它们的对数)。因此,我们将观察到所有先验的概率相乘以获得它们的联合分布。
关于r - 拉普拉斯恶魔 : when should I sum prior density?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/43115842/
我正在尝试使用R中的“密度”函数查找具有未知分布的连续数据的最佳PDF。现在,给定一个新的数据点,我想基于核密度找到该数据点的概率密度从“密度”函数结果得到的估计量。 我怎样才能做到这一点? 最佳答案
我在尝试绘制 density 对象时遇到问题。考虑例如 require(grDevices) set.seed(43) d0 = density(rexp(1e5,rate=1)) d1 = dens
我正在从 JAGS 迁移到 LaplacesDemon 并尝试重写我的一些代码。我读过 LaplacesDemon Tutorial和 LaplacesDemon Examples小插图,我对小插图中
猜猜看,另一个 Android-Bitmap-OOM 问题! 背景 压力测试our application已经注意到,在持续、大量使用(像猴子奔跑者一样)后,可能会最大化应用程序的进程内存分配,并在随
这两个在数学方面做的一样吗?它们产生相同的结果吗? Java: pixel=10; margin=px * getResources().getDisplayMetrics().density; la
有没有办法让 R 中的 density() 函数使用计数与概率? 例如,在使用直方图函数 hist 检查密度分布时,我有两个选择。 : hist(x,freq=F) #"graphic is a r
今年夏天我有一个项目即将完成,该项目可能会涉及大量的图像数据用于显示。我们讨论的是给定应用程序 session 中的数百个 640x480 左右的图像(显示时缩放到较小的分辨率),以及一次大量的非常大
我正在尝试使用 ggridges 包(基于 ggplot2)创建一个 joyplot。一般的想法是 joyplot 创建很好缩放的堆叠密度图。但是,我似乎无法使用加权密度生成其中之一。在创建 joyp
我正在尝试使用 ggridges 包(基于 ggplot2)创建一个 joyplot。一般的想法是 joyplot 创建很好缩放的堆叠密度图。但是,我似乎无法使用加权密度生成其中之一。在创建 joyp
我为 list 文件中的 Activity 添加了 android:configChanges="density"。但我没有收到 onConfigurationChanged() 回调
我有这个图像映射,它使用 CSS 实现在框中显示信息的悬停效果。基本上,您将鼠标悬停在其中一个 map 标记上,然后它会显示有关该 map 区域的信息。可以在下面的链接上查看。 http://phan
某些设备的设备密度返回 0。以下代码用于计算设备密度: WindowManager wm = (WindowManager) context .getSystemService(Con
float 越大(无论是正数还是负数),我们用于编码十进制数字的位数就越少,这是真的吗? 我们能否在 21 和 22 之间编码比 216 和 232 之间更多的十进制数字? 这两个范围之间的值计数是否
我正在用 python 编码。我有 3 个数组 x、y 和 z,我想用 colorbar 绘制平面 (x,y) 中 z 值的二维密度图。 所以在我的图中,点 x[0] 和 y[0] 的颜色将由 z[0
我正在阅读 The R Graphics Cookbook 并遇到了这组代码: library(gcookbook) library(ggplot2) p <- ggplot(faithful, ae
我想最终做到这一点: library(ggplot2) density=TRUE if (density) {ggplot(diamonds,aes(x=price)) + geom_histogra
在我的 android 中,avd 密度是 120,但是当我在 getResources().getDisplayMetrics().density 的帮助下计算它时,它是 .75。 为什么它们不同?
这是我要运行的代码 import sc2 from sc2 import run_game, maps, Race, Difficulty from sc2. player import Bot, C
我创建了一个依赖于某些屏幕指标来运行的应用程序。它适用于所有设备,除了 Galaxy Note 5。在 Note 5 上,它报告屏幕尺寸为 1080x1920,而屏幕实际上是 1440x2560。我还
我希望构建一个静态 KML(Google 地球标记)文件,该文件以 [lat, lon, density] 元组的形式显示一些给定数据集的热图样式渲染。 我拥有的一个非常简单的数据集是人口密度。 我的
我是一名优秀的程序员,十分优秀!