问题标签 [s]

R和S有什么区别？

有人对从 S-Plus 调用 R 的好方法有任何建议吗？理想情况下，我只想将代码传递给 R 并取回数据，而不必编写任何过于复杂的东西来集成它们。

我应该补充一点，我熟悉 Omegahat 上的RinS包，但我没有使用过它。我的印象是 Insightful 在 Tibco 接管之前已经努力整合环境。

编辑：事实证明，RinS 在 Windows 上不起作用。我发现最简单的解决方案就是使用 Rscript。我可以使用system()命令从 S-Plus 调用它。例如，这是一个简单的脚本：

从 S-Plus 调用它：

然后我只是将结果存储到一个文本文件中，并在脚本运行后将其导入 S-Plus。

我有一个 R / S / 非线性回归相关问题，我不是 R 程序员，所以我有点需要帮助。

我有两个数组 - tt 和 td。

我需要找到参数 a、b 和 c，因此对于非线性函数，最小二乘之和是最小的：

我不知道该怎么做。我尝试了nls()函数，nls2()但没有运气......

提前致谢。

编辑：

我的数据：

使用下面答案中的方法，我得到随机数据的 ok 值，但我使用的数据返回缺失值或评估模型消息时产生的无穷大。

很抱歉没有尽快提供数据。

我试图在 R/S-PLUS 中一次性获得按分类列分组的多个汇总统计信息。我发现了几个函数，但它们每次调用都会做一个统计，比如aggregate().

我正在寻找的是在一次调用中获取同一组的多个统计信息，例如平均值、最小值、最大值、标准差等，这可行吗？

在这里改写的问题：

我已经取得了一些进展，但从 R 那里得到了奇怪的行为......

这是我开始的xts

我想生产什么：

我写了一个这样的函数：

这成功地添加了 thisDate （在命令行中逐步运行代码）。但是添加 nextDate 的位似乎覆盖了它！我似乎也得到了一排意想不到的 NA。还在做这个...

我在第一列中添加了“无标题”，以表明它是 xts 日期索引，而不是向量/矩阵的一部分。

关于删除额外行的一点是因为我还没有解决覆盖问题并且正在试验。它不需要在最终答案中出现，但这是我目前要做的。

最后，当我询问这个​​结果并尝试将 nextDate 转换为我得到的日期时......

所以我有点糊涂...

下面的原始问题：

我在 R 中有一个时间序列（我学得很快，但显然不够快！）像这样

我想在每行的新列中使用 NEXT 行中的日期索引创建它的派生词：

对于价值（使用滞后）很容易做到，但对于日期索引本身却不是。

我可能可以使用各种查找等来解决如何做到这一点，但它很混乱。您必须在其他字段上进行匹配，或者摆弄那些感觉不太“符合 R”的行号。

有没有一种漂亮、整洁、优雅的方式来做到这一点？

我很确定我会去“D'OH！” 只要有人给出答案！但到目前为止，我还没有在这个网站上找到滞后日期索引的答案。

我想这样做的原因是我想连续使用每一对日期来询问另一个系列。所以可能有更好的方法来做到这一点。

我需要测试每个因素的“比索”值（见下面的复制代码）。一个因素是否达到“比索”总和的 50%，每个因素的值应粘贴到一个新对象“结果”中，否则，R 应评估哪个因素具有“比索”的最低聚合值，并考虑下一列中的因子再次汇总“比索”。基本上，此过程将最低得分的因素替换为下一个因素。该过程应该重复，直到一个因素超过 50% 的阈值。所以我的问题是，我从哪里开始？

它应该像

因此，因子 8 获得了更大的份额 4/(4+2+3) = 0.4444444，但在 a 轮中没有达到 50% 的阈值。因此，我需要更多的东西：重复聚合，但现在考虑“b”列中的因子 7 而不是“a”列中的因子 9，因为它在第一轮中获得了最低的聚合值。

我对 R 相当陌生。我正在使用 R 中的 ROCR 包来计算 AUC，我可以为一个预测器做这件事。我要做的是为 100 个不同的变量执行许多 AUC 计算。

到目前为止，我所做的如下：

但是，公式是文本格式，而 as.formula 给了我一个错误。任何帮助表示赞赏！提前致谢！

list(...)我对 R 和 S-plus的工作方式感到困惑。对于以下代码

在 S-Plus 中，它可以工作，但在 R 中，它给出“ Error in print(length(list(...))) : argument is missing, with no default”

我对它在 R 中的工作原理以及如何list(...)在 R 函数中获取值更感兴趣。

第一篇文章，所以轻松一点。

在 GLMing 的保险领域，经典的方法是对索赔频率和平均严重程度进行建模。考虑到这一点，我建立了几个模型来为自己进行实验，现在有一个问题。

有人可以解释一下 GLM 如何处理数据集的不同级别的汇总，特别是在错误估计方面？

考虑下面的例子。数据显示两个变量的严重程度趋势： - A 的索赔比 B 高 - 福特 > 起亚 > 沃克斯 > 捷豹

我为数据集的未汇总和汇总版本拟合了一个模型，因此 GLM 在两种情况下都拟合了相同的参数

然而，GLM 为未汇总的数据表明了一个拟合良好的模型。但是当我总结并使用加权平均值时，即平均严重性，模型拟合得很差。也许这正如您所期望的那样，毕竟未汇总的数据有更多可用于建模的点。此外，加权平均值似乎用于表示相对强度，因此在这里，指定加权平均值是没有意义的，因为它们都是相同的权重。

但更根本的是，我不能用 GLM 模拟平均严重性吗？我的意思是，我知道将 GLM 拟合到未汇总的数据集的结果将是平均严重性，但我希望将模型拟合到已经汇总的数据。似乎在聚合数据集上建模并不能真正表明模型拟合。

抱歉，如果这是一个愚蠢的问题，我不是统计学家，所以不完全理解 Hessian 矩阵。

请看下面的代码：

谢谢，

缺口

我使用本教程使用 ggplot2 绘制热图

http://learnr.wordpress.com/2010/01/26/ggplot2-quick-heatmap-plotting/

现在我想使用具有相同长度的另一个因子来订购垂直因子。

有人知道它是如何工作的吗？

谢谢