我的数据已经过两种处理（添加或不添加解决方案）中的一种，并且处于随机区组设计中。我有 8 个块（编号为 1-8），但在 aov 中，df=1。R 不知道有 8 个，所以我给它们写了字母（B1、B2 等）。现在双向 aov 中没有 F 值或 p 值，尽管现在 df=7。

起初：

输出

现在：

输出

我想在 python 中使用 RI 计算 p 值，我正在使用这个包 rpy2。我正在动态生成 count_a 和 count_b，并与它一起计算 p 值。当我运行脚本时，python 意外关闭，并收到以下错误消息：

“错误：'rho' 必须是一个不为 NULL 的环境：在启动期间在 C 级评估中检测到 - 警告消息：

中止陷阱：6"

数据如下：

下面是我的代码：

有什么建议或选项可以在 python 中根据一对值计算 p 值吗？

计算 binom 对象：

数据：c(94L, 107L) 成功次数 = 94，试验次数 = 201，p 值 = 0.3974 备择
假设：成功的真实概率不等于 0.5 95% 置信区间：0.3971286 0.5391627
样本估计值：
成功概率
0.4676617

但是，在提取 p 值时，我收到此错误：

文件“/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/rpy2/robjects/vectors.py”，第 233 行，在getitem res = super(Vector, self) 中。getitem (i) TypeError: 'str' 对象不能被解释为索引

惠，

我有一张包含 nb 汇总数据的表格。个体、平均值和标准差。由于大量 NA（我不能用 0 替换），我无法使用 rawdata，我无法进行成对比较或 Anova。所以我从 BSDA 包中选择 t.test 和 tsum.test。

这给出了一个结果，例如

我的问题是：“如何使用此 tsum.test 应用 Bonferonni 调整？” 如何修改p值？一种相当于 p.adjust.method

我试过了，但这不起作用（我将 n=2 用于 2 次比较，但我在真实数据集中得到 12 次比较）

谢谢大家

我不确定如何从 gofstat 函数中获取 p 值。我确实看到了可能有“未计算”、“未拒绝”和“拒绝”的结论的**测试。任何帮助都会很棒。假设使用的数据是连续的。

使用下面的代码

我想运行逻辑回归的因变量（在我的数据集中它是：dat$admit），其中包含所有可用变量、对和三元组（3 个独立变量），每个回归都有不同的自变量与因变量。我想要返回的结果是一个连续的每个回归摘要的列表：系数、p 值、AUC、CI 95%。使用下面提交的数据集应该有 7 个回归：

这是一个样本数据集（其中 dat$admit 是逻辑回归因变量）：

根据@marek 评论，输出应该是这样的（仅针对女性以及来自女性和 apcalc ）： # Intercept Estimate P-Value (Intercept) P-Value (Estimate) AUC # female 0.000000e+00 0.000000e+00 1 1 0.5

女性+apcalc 0.000000e+00 0.000000e+00 1 1 0.5

@David Arenburg 编写了一个很好的代码，它可以生成统计数据，但没有配对和三重奏的模型创建，所以我想知道如何添加模型创建。这是大卫阿伦伯格的代码吗？

知道如何创建此列表吗？谢谢，罗恩

在 R 中，我可以使用phyper () 函数计算超几何分布的 p 值，其中返回数组中的第一个值是 p 值。

我想知道 Go / Golang 中是否有任何包可以让我完全在 Go 中进行计算？

我正在运行许多回归，并且只对一个特定变量的系数和 p 值的影响感兴趣。因此，在我的脚本中，我希望能够从 glm 摘要中提取 p 值（获取系数本身很容易）。我知道查看 p 值的唯一方法是使用 summary(myReg)。还有其他方法吗？

例如：

我尝试将fit$coefficients其视为矩阵，但仍然无法简单地提取 p 值。

是否有可能做到这一点？

谢谢！

如何确定具有 n 个自由度的 t 分布的 P 值。

对这个主题的研究将我指向这个堆栈交换答案：https ://stackoverflow.com/a/17604216

我假设 np.abs(tt) 是 T 值，但是我如何在自由度上工作，那是 n-1 吗？

提前致谢

我正在尝试使用 和 获得两个光栅砖之间的像素相关性和显着性（p 值cor）cor.test。我的数据在这里：

• 砖 1
• 砖 2

它们都相当小，总共不到 2MB。

我从之前关于 StackOverflow 和 r-sig-geo 的讨论中找到了以下两个代码（均来自 Robert Hijmans）：

两个代码都按预期工作，cor函数产生相关网格。但是，我尝试用corwith代替cor.test以提取 p 值：

我遇到了以下错误（在 RStudio 中有引用）：

在之前的 r-sig-geo 讨论中，用户曾询问过此错误，但没有得到答复。所以我再次询问并收到了对我的询问的回复，其中一位用户指出cor能够输入矩阵并且cor.test不能，但即使在将数据转换为数字向量之后：

我面临以下错误：

我想知道是否有人可以帮助解决这个问题？

我的 sessionInfo() 如下：

谢谢！

我知道有很多关于如何从aov. 但是，我有一个包含数千个样本的列表。我为每个样本做了一个 aov 来比较两种不同的治疗方法，现在我正在寻找一种方法来获取包含所有 p 值的列表，因为我不能一个一个地复制它们。

这甚至可能吗？对于 a 创建的 p 值，我没有问题ttest：

数据是我导入的 .csv 和 [1:3] 表示与列 [4:6] 比较的列

所以这真的不是问题，但似乎不可能为 aov 做类似的事情：

我无法获得包含所有 p 值的列表（称为 Pr(>F).. 这有点令人沮丧.. 希望你明白我想要做什么，