在R

Question

我正在尝试在R中计算。我有以下值。

nb <- 100
tb <- 25
ns <- 90
ts <- 15
A0 <- 1
S_norm <- 0.4

R <-tb/ts
y_meas <- (ns-nb/R)/A0
sigma_meas = sqrt(ns+(nb+1)/R^2)/A0

我对如何整合L（psi）非常困惑，比如-10到10.因为我正在整合日志A.

Answer 1

您可以替换logA，并且您可以按以下方式整合psi的固定值：

psi <- 5
integrate(function(x) exp(-0.5*(((x/A0)/S_norm)^2 + ((psi-y_meas*A0/exp(x))/sigma_meas)^2)),
                                                                                   -10, 10)
# 0.1775989 with absolute error < 6.6e-05

Answer 2

除了@SandipanDey的优秀答案之外，如果你可以将限制扩展到-Infinity ... + Infinity，有一种更好的方法将函数与e ^{-x 2 集成}内核：Gauss-Hermite quadrature，并且有一个 R 包。

简单示例：

AFRAME.registerComponent('foo', {
  init: function () {
    this.baz = 'hello';
    console.log(this.baz); // hello
  },
  doStuff: function () {
    console.log(this.baz); // hello
    this.baz = 'bonjour';
    console.log(this.baz); // bonjour
  }
});