问题标签 [mass]

运行这个glmer.nb，我收到了错误信息

边界（奇异）拟合：参见 ?isSingular

警告消息：在 theta.ml(Y, mu, weights = object@resp$weights, limit = limit, : 达到迭代限制

我跑了singular(model)，它又回来了TRUE。然后我跑了summary(model)，它显示了警告信息

警告消息：1：在 vcov.merMod(object, use.hessian = use.hessian) 中：从有限差分 Hessian 计算的方差-协方差矩阵不是正定的或包含 NA 值：回退到从 RX 2 估计的 var-cov ：在 vcov.merMod(object, correlation = correlation, sigm = sig) 中：从有限差分 Hessian 计算的方差-协方差矩阵不是正定的或包含 NA 值：回退到从 RX 估计的 var-cov

我抬头一看，似乎警告消息表明模型不准确。但我不确定我能做些什么来解决这个警告。我将非常感谢任何帮助！

使用包分析序数数据survey，我在尝试使用 raked 数据时遇到了一些问题。没有耙，svyolr()工作没有任何问题，但是当我在耙后尝试分析时，svyolr 遇到错误Error in if (any(y < 0 | y > 1)) stop("y values must be 0 <= y <= 1") : missing value where TRUE/FALSE needed。

可以使用示例数据集 api 重现该问题：

难道我做错了什么？无法svyolr处理倾斜或后分层（相同的错误消息）数据？我在帮助文件中没有看到任何提及，甚至在survey:::svyolr.survey.design2.

是否可以使用survey或任何其他包来执行此操作R？

我正在尝试使用具有 1270 个维度/变量的数据集绘制一个平行坐标图，并且由于背景网格，该图似乎都是一个灰色区域（如果我绘制具有 800 个维度的相同数据集，我仍然得到网格占据了大部分空间，但你可以或多或少地看到情节）。

有谁知道我怎样才能删除这个背景网格？

我正在尝试建立一个泊松混合模型，但残差有些异方差。所以，我想为模型定义一个方差结构。据我所知，我可以通过使用包glmmPQL中的命令来做到这一点MASS。但这由于某种原因不起作用。这是代码：

我认为一切都很酷，但我收到了以下错误消息：

有任何想法吗？谢谢

这是向客户发送带有文件（图片）的消息。您好，请帮忙。现在我只能发送到一个号码。如何向多个号码发送时事通讯。

我尝试了不同的选项，但没有帮助，也许有兴趣的人可以帮助解决方案（添加）邮寄到超过 1 个号码

我正在尝试根据癌症 Recurrence(R) 或 Non-Recurrence(NR) 训练我的基因表达数据。我正在处理的数据集：

87个样本和10个基因（一列称为“条件”，表明每个样本是R或NR。我使用了这些命令：

但我得到这个错误：

FUN(x, aperm(array(STATS, dims[perm]), order(perm)), ...) 中的错误：二进制运算符的非数字参数

请你帮帮我好吗？

在 R 中，使用库中的'geyser'数据集MASS，我们正在构建一个非参数公差区间 cont。90% 的间歇泉喷发 w/ 99% CI

请在下面找到我正在使用的方程式：

这是我正在使用的代码，其中我遗漏了一些东西......：

寻找两者

当我试图在其中一个参数中进行一些计算时，我发生了 mvrnorm() 函数的一个奇怪行为。

这确实可以完成工作：

然而，用 n = (1-0.8)*1000 替换 n = 200 并不能完成这项工作：

为什么会这样？

我有一个数据集，其中包含从 2019 年到 2021 年的每周计数。我想做的是将 2020 年给定一周的每周计数与 2019 年同一周的计数进行比较，并将 2021 年的计数与2019 年同一周。数据如下：

在我的真实数据中，存在显着的过度离散，因此我认为负二项式模型可能最适合。我已经运行了以下内容：

参考类别factor(year)是 2019 年，这就是我想要的。但是，我正在努力解释系数（以及 IRR），week因为第 1 周是参考类别。

有没有更好的方法来做到这一点，以实现周/年比较？我的主要目标是绘制 2020 年和 2021 年相对于 2019 年的每周 IRR。