- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试使用带有 R 的 MH 算法来实现一个简单的 MCMC,问题是我得到了这个错误(我试图计算 alpha 并且它不是 NA 问题)
Error in if (runif(1) <= alpha) { : missing value where TRUE/FALSE needed
PoissonMetropolisHastingRW = function(n=100000,lambda=10,p=0.5,x0=0){
x=rep(0,n); y=0; alpha = 0
x[1]=x0
for(i in 2:n){
if (x[i-1] == 0){
y = sample(c(0,1),1, prob=c(0.5,0.5))
alpha = min(1,((lambda^y)*x[i-1]*p)/((lambda^x[i-1])*y*(1-p)))
#alpha = min(1, ( ((lambda^y)*x[i-1])/( (lambda^x[i-1])*y) )*(p/(1-p)) ))
if(runif(1)<=alpha) {x[i]=y}
else {x[i]= x[i-1]}
}
if (x[i-1] > 0){
y = sample(c(x[i-1]-1,x[i-1]+1), 1, prob=c(1-p,p))
alpha = min(1,((lambda^y)*x[i-1]*p)/((lambda^x[i-1])*y*(1-p)))
#alpha = min(1, (((lambda^y)*x[i-1]/((lambda^x[i-1])*y))*(p/(1-p))))
if(runif(1) <= alpha) {x[i]=y}
else {x[i]= x[i-1]}
}
}
return(x)
}
最佳答案
如 y
恰好为 0(并且每次迭代的概率为 0.5,这肯定会发生),然后 alpha
是 0/0(因为 x[i-1] == 0
)。它给你 NaN
.条件something <= NaN
提供 NA
.
关于r - Metropolis-Hastings MCMC 与 R,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/29698831/
运行 npm run test 时出现以下错误: jest-haste-map: Haste module naming collision: {{name}} The following fil
我已经创建了一个自定义的 npm 模块(将使用 xxx 而不是它的名称)并使用 npm install 手动链接它。 我非常努力地搜索: [Workarounds] Packager unable t
当我安装 react-native-popup-menu 时,我遇到了这样的错误。 这是错误代码 Looking for JS files in E:\react-native\...........
我在运行./gradlew assembleRelease时遇到了这个问题 Error: jest-haste-map: Haste module naming collision: Duplicat
我对 Haskell 这种语言有些熟悉,但对工具链不太熟悉。 (在 cabal 和 stack 存在之前,我曾使用 Haskell。)有人告诉我,stack 是我应该用来管理 Haskell 项目的工
我正在尝试在我的应用程序中集成 LinkedIn。所以,我做了如下: 安装的依赖如下: npm install react-native-linkedin --save; 因此,在启动应用程序时,出现
如果是这样,那么我可以让人们通过更改其(Haskell)配置来玩弄我的项目:在浏览器中重新编译它,然后在浏览器中重新运行它。我不在乎速度。我只想能够在浏览器中演示一个难以设置的桌面项目。因此,可能的新
我正在测试 Haste-compiler,我需要访问事件和/或原始元素,但 Event 类型不提供该信息。 我尝试使用下面的函数进行 ffi 回调,但它返回给回调的“事件”是未定义的。我做错了什么,或
我正在尝试在 C 中实现 Metropolis-Hastings 算法以实现简单的线性回归(不使用其他库(boost、Eigen 等)并且不使用二维数组)*。为了更好地测试代码/跟踪图的评估,我通过保
在运行 ReactNative Android 应用程序时出现以下异常,并且在 Bundle 时出现相同的异常。 react-native bundle --platform android --de
React Native 环境信息: System: OS: macOS High Sierra 10.13.6 CPU: x64 Intel(R) Core(TM) i3-7100 CPU
我正在尝试运行现有的博览会项目。我已经克隆了它并运行了 yarn ,然后博览会开始了。并收到此错误: 无法从 /path-to-project/node_modules/react-native/Li
我正在尝试使用 Metropolis-Hastings 进行贝叶斯回归。测试数据生成如下(Python代码,我没有复制整个代码): trueA = 5 ; trueB = 7 ;trueSD = 10
我一直在阅读 Metropolis-Hastings (MH) 算法。从理论上讲,我了解该算法的工作原理。现在,我正在尝试使用 python 实现 MH 算法。 我遇到了以下 notebook .它完
我是 Node.js 和 react-native 的新手。我按照 send_telemetry.js 上的示例进行操作完全正确,但是当我运行我的 react native 应用程序时,我收到错误:“
我正在尝试创建一个 react-native 包并收到以下错误: Unable to resolve module `expo` from `C:\Users\\Projects\myProject\
在我的 Metro Bundler 控制台输出中看到: Error: Cannot find module 'jest-haste-map' at Function.Module._resolveFi
我正在使用Haskell Stack对于一个项目,我想包括 Haste它编译客户端逻辑。我喜欢这样一个事实:Stack 抽象了环境中不同的构建和安装问题,如果它在我的机器上构建,它也会在其他人的机器上
我正在尝试使用带有 R 的 MH 算法来实现一个简单的 MCMC,问题是我得到了这个错误(我试图计算 alpha 并且它不是 NA 问题) Error in if (runif(1) 0){
我有一个相对简单的函数,具有三个未知输入参数,我只知道它们的上限和下限。我也知道我所有数据的输出 Y 应该是什么。 到目前为止,我已经在 python 中完成了一个简单的网格搜索,遍历所有可能的参数组
我是一名优秀的程序员,十分优秀!