r - 将一个统计程序的输出复制到 R

Question

我组织的一名前员工在与 R 不同的统计程序上执行了所有分析（没有文档），目前没有人知道使用了哪个程序。查看模型输出表并将其与 Google 搜索结果进行比较，我认为他们使用了 Statistica。为了对与我们合作的其他组织保持透明，我正在尝试复制他们的工作并可能重新评估它。

模型：他们模拟了三个变量之间的关系，我称之为A、B、C。变量的选择基于探索性分析（即相关矩阵和 GLM 建模）。参数估计用于预测目的。据我所知，他们使用具有日志链接功能的 GLM 将C建模为A和B的函数。

资料：

A <- c(0.937918714, 1.277501774, 34.46428571, 3.843879361, 5.135520685, 0.324675325, 1.038421599, 0.333333333, 0.058139535, 0.09009009, 0.080515298, 5.174234424, 10.625, 21.9047619, 0.162074554, 2.372881356, 1.084430674, 18.53658537, 6.438631791, 0.172413793, 0.291120815, 9.090909091, 5.882352941)

B <- c(0.416666667, 0.555555556, 0.833333333, 0.4, 0.833333333, 0.4, 0.625, 0.625, 0.294117647, 0.37037037, 0.285714286, 1.111111111, 0.588235294, 0.476190476, 0.555555556, 0.833333333, 0.666666667, 0.476190476, 0.208333333, 0.163934426, 0.163934426, 0.3125, 0.454545455)

C <- c(0.009533367, 0.020812183, 0.056208054, 0.015002587, 0.042735043, 0.013661202, 0.004377736, 0.00635324, 0.001345895, 0.001940492, 0.00446144, 0.043768997, 0.021134594, 0.004471772, 0.023488256, 0.029441118, 0.052287582, 0.003526093, 0.030984508, 0.010891089, 0.020812686, 0.016032064, 0.018145161)

我的方法：

我将每个向量组合成一个数据框（dat）并使用以下模型进行建模：

glm(formula = C ~ A + B, family = binomial(link = logit), data = dat)

问题：

我注意到我们有不同的参数估计；事实上，他们的分析包括“规模”作为一个因素，以及相关的参数估计和标准误差（见下文）。我还没有弄清楚如何包含单独的“比例”因子。我的参数估计值接近这些，但包含一个新变量显然不同。

任何熟悉此 [Statistica] 输出以及如何在 R 中复制它的人？首先，我如何将比例因子合并到我的分析中？

旁注 我也将此发布到 Reddit（r/rstats - 复制在不同软件中执行的分析）。

非常感激！

标签： rstatistica