小编Xoc*_*lli的帖子

如果我根据卡方分布绘制我的多元数据的马氏距离的多元 qqplot，我希望伴随的 qqline 是截距为 0 和斜率为 1 的线。但是如果我运行以下代码：

scores<-matrix(rnorm(100*3),nc=3)
mah_scores = mahalanobis(scores,center=colMeans(scores),cov=cov(scores))
chi_scores = qchisq(ppoints(nrow(scores)),df=3)
qqplot(x=chi_scores,y=mah_scores,asp=1)
abline(a=0,b=1)
qqline(mah_scores, distribution = function(p) qchisq(p,df = 3),col="red")

我得到下图：

我希望 qqline（红色）与截距 0 和斜率 1（黑色）的线相同。谁能向我解释为什么这两行不匹配？

（我正在运行 R 版本 2.15.3 (2013-03-01)）

我想知道为什么 Python 的输入模块不支持isinstance(<obj>, Any)并引发TypeError. 我希望它总是会回来True。它不总是返回有什么具体原因吗True？

• 此处引发 TypeError 。
• 类型错误示例：

>>> from typing import Any
>>> isinstance(1, Any)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/home/peter/miniconda3/envs/prefect_utils/lib/python3.9/typing.py", line 338, in __instancecheck__
    raise TypeError(f"{self} cannot be used with isinstance()")
TypeError: typing.Any cannot be used with isinstance()

是否可以使用新值扩展数据类型？

例如:以下编译:

data Axes2D = X | Y
data Axes3D = Axes2D | Z

但是,以下内容:

data Axes2D = X | Y deriving (Show, Eq)
data Axes3D = Axes2D | Z deriving (Show, Eq)

type Point2D = (Int, Int)
type Point3D = (Int, Int, Int)

move_along_axis_2D :: Point2D -> Axes2D -> Int -> Point2D
move_along_axis_2D (x, y) axis move | axis == X = (x + move, y)
                                    | otherwise = (x, y + move)

move_along_axis_3D :: Point3D -> Axes3D -> …

我正在自学 PyMC3，并且正在复制这个变更点检测示例。我想通过随时间对 2 个不同泊松率参数的后验采样来扩展示例。我该怎么做呢？

提供更多细节：如果我有y示例中的时间序列。我想先用观察采样y[0:2]，然后y[0:3], y[0:4], ... , y[0:100]。我的目标是绘制mu0和mu1（泊松率参数）随时间 2..100 的分布。显然，我应该能够通过创建大约 98 个不同的模型并从每个模型中独立采样来做到这一点。但是，这似乎效率很低，我想知道是否有更好的方法来仅使用新观察更新模型并继续跟踪？

我在文档或在线找不到任何相关内容，因此我尝试仅替换变量，但随后出现“ValueError: Variable name changepoint already exists.”。

TLDR；有没有一种有效的方法可以随着时间的推移更新观察结果并使用这些新观察结果继续从先前的跟踪中采样？