问题标签 [multivariate-testing]

我想对网站的整个部分进行 A/B 测试或实验。例如，在我的 /blog/ 页面上，其中一种变体将具有时事通讯表格，而另一种变体将具有免费的电子书下载按钮。

问题是我必须使用完整的 URL 路径进行实验，例如 /blog/2013/article/1?var=1 和 /blog/2013/article/1?var=2 使用这种方法我需要创建每篇博文都有一个新的实验。这是不可能的。

关于如何解决这个问题的任何提示？

你好 Stack Overflow 的人，我花了一段时间寻找解决我的问题的方法，但没有找到任何东西，所以我想我会发布。

基本上我有一个按字母顺序列出的 196 个国家的数据集。其中一个变量根据国家所在的地区分配一个 1-10 的数字。例如东欧 = 1、西欧 = 2、中东 = 3、南美 = 4 等等。

这是数据集的可视化表示：

国家名称-----国家地区------婴儿死亡率

阿富汗------------3------------180

阿根廷---------------4------------65

法国-------------------2------------------------12

德国---------------2-------------10

波兰-----------------1------------------------16

我需要做的是将这 10 个区域分成各自的虚拟变量，以便我通过多元回归对它们进行运行，以确定它们对婴儿死亡率的个体影响。

我想知道创建虚拟变量（1 = 东欧，0 = 其他等）的必要代码是什么，然后如何在单独和多元回归中测试它们的效果。

抱歉，如果这看起来很简单或愚蠢的问题，我对使用 R 相当陌生。

我在这里先向您的帮助表示感谢。

编辑：这是要求的 dput 输出：

有谁知道 R 中执行 MANOVA 同时控制系统发育非独立性的包或方法？

谢谢！

第一次在这里提出问题，我会尽力明确 - 但如果我应该提供更多信息，请告诉我！其次，这是一个很长的问题......希望为某人解决简单;）！所以使用“R”，我根据一些论文（Manera et al. 2012）对多元 GARCH 模型进行建模。

我在均值方程中使用外部回归器对恒定条件相关 (CCC) 和动态条件相关 (DCC) 模型进行建模；使用“R”版本 3.0.1 和包“rugarch”版本 1.2-2 用于带有外部回归器的单变量 GARCH，以及用于 CCC/DCC 模型的“ccgarch”包（版本 0.2.0-2）。（我目前正在研究“rmgarch”包——但它似乎只适用于 DCC，我也需要 CCC 模型。）

我的模型的平均方程有问题。在我上面提到的论文中，CCC 和 DCC 模型之间的平均方程的参数估计发生了变化！而且我不知道我将如何在 R 中做到这一点......（目前，查看谷歌并查看 Tsay 的书“金融时间序列分析”和 Engle 的书“预测相关性”以找到我的错误）

我所说的“我的平均方程不会在 CCC 和 DCC 模型之间改变”的意思是：我使用包 rugarch 为我的 n=5 时间序列指定单变量 GARCH。然后，我使用 GARCH 的估计参数（ARCH + GARCH 项）并将它们用于 CCC 和 DCC 函数“eccc.sim()”和“dcc.sim()”。然后，从 eccc.estimation() 和 dcc.estimation() 函数，我可以检索方差方程的估计值以及相关矩阵。但不适用于平均方程。

我只发布了单变量模型和 CCC 模型的 R 代码（可重现的和我的原始代码）。已经谢谢你阅读我的帖子了！！！！！！

注意：在下面的代码中，“data.repl”是dim 843x22 的“zoo”对象（9 个每日商品返回系列和解释变量系列）。多元 GARCH 仅适用于 5 个系列。

可重现的代码：

我的原始代码：

我正在使用 FactoMine R 绘制 MCA 图。我的数据表如下所示：

这是我的第一个红色标志......之后，我执行了 MCA 只是为了看看我会得到什么，这就是代码：

然后我在控制台中得到以下信息：

我是使用 R 的新手（如 1 周），但我有使用 SAS 的经验......我不知道我做错了什么以及为什么 R 将我的数据修复到上述结构中（3L、3L、3L.. .) 有没有人知道如何进行？

我编写了一个代码来运行多个证券的多个时间序列滚动回归。由于证券数量超过 10,000 个，并且每个证券有超过 200 个滚动窗口，因此顺序设置（使用 foreach %do%）的运行时间约为 30 分钟。

我想使用“doParralel”后端为并行计算实现 foreach %dopar% 。简单地在代码中用 %dopar% 更改 %do% 并不能解决问题。我对这种并行计算方法很陌生，希望能得到一些帮助。

这是 foreach %do% 代码：

为了使用 %dopar% 启用并行计算，我调用并注册了后端“doParallel”。

非常感谢！

更新

这是我的第一次尝试：

这在循环结束之前完美运行。但是，beta.df = do.call('combine', result)给出以下错误：Error in do.call("combine", result) : could not find function "combine"。

如何组合结果的输出。现在它是一个列表而不是一个数据框。

谢谢，

我在 FactoMineR 中遇到了 MFA 问题。我正在处理一个数据集，其中包含在番茄植物中测量的物理、化学和微生物连续变量，这些变量取自 2 种不同的处理方法和 3 个时间点。我已经容纳了这样的数据：

我将变量分类（前 2 个），然后其他 16 个是连续的。但是，我想分别处理这 2 个分类变量。所以我写了以下代码：

但是，它似乎不起作用。因此，我尝试了以下方法：

还有这个：

但我一直遇到同样的问题（“不方便的组定义”）。我可以做些什么来分别保留前 2 个分类组？我非常感谢您对如何正确运行模型的建议！

最好的祝愿，

艾玛

我一直在尝试使用 R 作为 MVSP 的替代品进行聚类分析和 PCA。但是，R 使用我发现的所有函数（包括 dist、bcdist、hclust 和 daisy 函数）给出了与 MVSP 截然不同的输出。我还使用了 MVSP 的距离表输出作为 R 中的距离矩阵，它产生了进一步不同的输出。它需要使用欧几里得距离、索伦森骰子系数和平均/UPGMA 聚类的聚类/树状图。第二个人/输入以及我自己在多台计算机和两个版本上都重复了这些问题。SAS 给出了与 MVSP 相同的结果。

是否有另一个我可以使用的包（替代 dist 或 hclust）或在 R 中查看/测试算法的方法？有什么东西会导致这种情况，并且可以将 R 的版本回溯到更早的作品吗？

编辑; 发现了问题，我尝试的所有函数中使用的算法都没有使用 Sorenson 系数，所以我使用了 Proxy 包，带有 dist 函数（以及内置方法 =“Dice”）。

我正在尝试运行多元线性混合模型，并且需要使用远程工作站来减少计算时间。

当我在我的个人计算机（R 版本 2.15.1、lme4 版本 0.999999-0、64 位 Unix）上从 lme4 运行 glmer() 时，我的模型运行正常。

当我切换到删除工作站（R 版本 3.0.2、lme4 版本 1.0-6、64 位 Linux）时，我的模型无法运行，并且收到错误和警告消息：

有趣的是，如果我对任一因变量运行 lmer()，模型将运行而不会出现错误或警告。

这是示例代码：

这是我个人电脑的输出：

以下是使用较新版本 lme4 的远程工作站的输出：

我想知道模型是否形成不正确，或者一个版本的 lme4 是否存在并发症。非常感谢任何想法或建议。谢谢！！

我正在使用 vegan 包中的 adonis 函数来确定不同环境因素对不同地区森林植物群落组成的影响。我想首先使用 adonis 来消除区域效应，这是为了适合像这样的模型

其中社区是属于不同区域的森林斑块中物种的存在/不存在数据矩阵。

然后我想使用这个模型的残差作为其他一些分析的响应变量。我知道我可以使用地层将区域建模为块变量，但这对我不感兴趣，因为之后我想执行另一种分析，我希望消除区域效应。

我的问题是我无法弄清楚如何从 adonis 模型中获取残差值。

任何可以解决此问题的其他类型的多变量分析的提示都非常受欢迎。