python中两个数据帧之间的Wilcoxon秩和检验
在Python中，Wilcoxon秩和检验是一种非参数统计方法，用于比较两个相关样本或配对样本之间的差异。它基于样本的秩次而不是具体数值，因此不需要对数据的分布做出假设。
Wilcoxon秩和检验可以用于以下情况：
- 比较两个相关样本的差异：当我们有两组相关的样本数据时，可以使用Wilcoxon秩和检验来判断这两组样本是否存在显著差异。
- 比较两个配对样本的差异：当我们有一组样本数据在两个不同的条件下进行了测量时，可以使用Wilcoxon秩和检验来判断这两个条件下的测量结果是否存在显著差异。
在Python中，可以使用SciPy库中的wilcoxon函数来进行Wilcoxon秩和检验。该函数的使用方法如下：
from scipy.stats import wilcoxon
# 假设我们有两个数据帧 df1 和 df2
# 假设我们要比较它们的某一列数据 col 的差异
# 提取两个数据帧的 col 列数据
data1 = df1['col']
data2 = df2['col']
# 进行Wilcoxon秩和检验
statistic, p_value = wilcoxon(data1, data2)
# 输出结果
print("Wilcoxon秩和检验统计量：", statistic)
print("p值：", p_value)推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云云服务器（CVM）：https://cloud.tencent.com/product/cvm
- 腾讯云云数据库 MySQL 版：https://cloud.tencent.com/product/cdb_mysql
- 腾讯云人工智能平台（AI Lab）：https://cloud.tencent.com/product/ailab
- 腾讯云物联网平台（IoT Hub）：https://cloud.tencent.com/product/iothub
- 腾讯云移动开发平台（MPS）：https://cloud.tencent.com/product/mps
- 腾讯云对象存储（COS）：https://cloud.tencent.com/product/cos
- 腾讯云区块链服务（BCS）：https://cloud.tencent.com/product/bcs
- 腾讯云元宇宙服务（Tencent XR）：https://cloud.tencent.com/product/xr
相关·内容
1回答
我有两个列表，我想检查一下wilcoxon秩和检验。
我看到有库，但它只显示了双面测试。
如何在python中进行单面测试？
浏览 18提问于2020-02-29得票数 1
回答已采纳
SAS编码：-我对两组(独立但来自同一人群)的差异进行测试。差异量和t-stat匹配的符号(即两组之间的数学差异为负数，tstat为负值)。或者，如果两组之间的数学差异为正，则tstat为正)。
然而，当我进行wilcoxon秩和检验时，我的z-分数的符号与组差的符号(-/+)不匹配。(也就是说，两组之间的数学差异是负的，而z-得分是正的)。如果两组之间的数学差异为正，则z得分为负数)。
我尝试过对数据集进行常规排序和降序排序。
这是我的密码：
*proc sort data = fundawin3t;
by vb_nvb_TTest;
run;
**Wilcoxon rank su
浏览 1提问于2018-03-29得票数 0
回答已采纳
我正在对我的数据进行威尔科克森秩和检验。数据具有不同的单位和很大程度上不同的尺度。在运行测试之前，我需要键入这个命令吗？
## where x is my data frame
scale(x, center = TRUE, scale = TRUE)
或者，Wilcoxon秩和检验本身就是这样做的？
浏览 0提问于2013-02-14得票数 0
回答已采纳
我想为一个小样本(n=32)计算一个秩和测试统计量。scipy.stats.mannwhitneyu对测试统计数据和相关的p值使用正态近似(参见)。对于小样本来说，这可能会产生误导。
在python中是否有一种方法来计算准确的测试统计数据和/或相关的p-值(正如描述的，当然不是在excel中)？
[R有这样的选项，参见，但我更愿意呆在python。]
浏览 2提问于2021-03-15得票数 2
1回答
我有这个数据集。我试图找出连续变量expr与范畴变量WHO_Grade之间的相关性。
> dput(tmp)
structure(list(expr = c(3.72491159808923, 7.8316405937301, 4.1302793124001,
6.81536170645658, 6.68352582647051, 6.0974581720256, 6.81642917136002,
6.52282686468863, 6.95033151442703, 7.40122305409127, 6.734502473652,
4.52338197246748, 5.6
浏览 5提问于2022-03-24得票数 0
python的scipy.stats.ranksums和R的wilcox.test都应该计算用于Wilcoxon秩和检验的双边p值。但是当我对相同的数据运行这两个函数时，我得到的p值相差了几个数量级：
R：
> x=c(57.07168,46.95301,31.86423,38.27486,77.89309,76.78879,33.29809,58.61569,18.26473,62.92256,50.46951,19.14473,22.58552,24.14309)
> y=c(8.319966,2.569211,1.306941,8.450002,1.624244,1.8871
浏览 0提问于2012-10-09得票数 12
回答已采纳
我有一个长数据集(N=499)。在这里，我比较了一个索引在8个不同的处理中的表现(每个处理中有不同的样本数)。
我已经做了一个Kruskal，它是显着的(p值< 2.2e-16)。
现在，对于后特别测试，我正在考虑使用邓恩的测试，但我读到威尔克森也可能有用，有什么建议吗？
非常感谢大家。
浏览 1提问于2015-07-15得票数 1
我想对我的配对数据进行单侧wilcoxon秩检验，因为我感兴趣的是其中一个样本是否明显大于另一个样本。
Scipy提供
scipy.stats.wilcoxon(x,y)
对成对的样本x和y进行双边检验。因为我不能假设正态(对称)分布，所以我不能从双边p值推导出单侧p值。
现在有没有人可以用python的方法来获得单边测试的p值？
谢谢!
浏览 3提问于2013-04-30得票数 14
回答已采纳
1回答
我已经创建了两种使用统一来执行相同功能的模拟，但是使用了两种不同的编程范式(面向对象和面向数据)来演示一种模式相对于另一种模式的性能优势。我收集了两种模拟的CPU速度值(用于测量仿真中每个循环之间的时间)。我已经为模拟的300次迭代编译了两个包含这些CPU速度值的数组，我希望对这两个数组进行统计分析，以确定模拟之间是否存在显着的性能差异。
我目前的问题是，我不知道如何分类我的数据。我最好的猜测是，我的数据是参数化的(我对此不太确定)，我应该使用独立的t检验，但我也读过一些关于Kruskal Wallis和Mann Whitney的东西，这让我怀疑t检验是否是解决这个问题的合适途径。
任何帮助都
浏览 0提问于2022-03-22得票数 1
我深入研究了这项测试的scipy代码Wilcoxon符号秩检验()，我发现在scipy中，他们计算正差值的秩和，以及负差值的秩和。然后，他们选择了较小的一个，并将其用作W。这与维基百科或其他网站()中的测试解释有很大不同，其中W=sum(所有差异)。
这种方法有效吗？
浏览 0提问于2013-06-19得票数 2
回答已采纳
1回答
我正在尝试执行两个数据帧之间的Wilcoxon秩和检验。我想只在行之间执行测试。例如，测试应该只在df1 (A，1，2，3)和df2 (A ,10，12 ,13)中的第1行，df1 (B ,4，5，6)和df2 (B ,14，15，16)中的第2行之间进行，依此类推。 df1=pd.DataFrame(np.array([['A',1, 2, 3], ['B',4, 5, 6], ['C',7, 8, 9]]),
columns=['Details','a', 'b&
浏览 12提问于2020-11-04得票数 0
回答已采纳
1回答
预聚合数据的正态性检验
、、、、
使用spark I聚合每个组(队列)的数据，以仅包含平均值、标准差和方差。 现在，在使用python的第二步中，我想测试正态性(https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.normaltest.html)，然后使用t-测试stats.ttest_ind或stats.wilcoxon等级测试进行显著性测试。 然而，所有这些方法都希望数据以面向记录的原始值的形式输入。如何将它们与预聚合的数据一起使用？
浏览 21提问于2019-09-06得票数 0
回答已采纳
1回答
我正在对实验前后每个人的平均心痛数据进行wilcoxon测试。大约有62人参加。
我跑：
statistic, p-value = wilcoxon(dataBefore - dataAfter)
P值为0.0005，而统计量为498.这似乎是一个异常大的数字。我注意到在这里归还的东西有些不清楚。我以为this使用了t检验，但我也看到人们把它作为z统计数据存储。它还回来了什么？
这本手册对我来说没有多大意义：它说：“统计:浮点--是零上或零下差值的总和，以较小的值为准。”但这对我来说没什么意义。这是人们通常报告的价值吗？
浏览 6提问于2016-10-27得票数 0
回答已采纳
1回答
如何检验以下假设？票价越大，顾客就越有可能独自出海。
利用下面的数据，我们怎样才能检验假设呢？
import seaborn as sns
# dataset
df= sns.load_dataset('titanic')
df[['fare','alone']].head()
fare alone
0 7.2500 False
1 71.2833 False
2 7.9250 True
3 53.1000 False
4 8.0500 True
更新
#subset for alone = True
浏览 0提问于2022-04-13得票数 1
我正在使用PSPP，并希望比较一个样本的中值与给定的中值值，但我发现的所有选项都是比较两个变量的相互比较。通过定义一个充满引用值的变量，我尝试了一种解决方法：
NPAR测试我的实际数据集的/Wilcoxon变量，变量中填充了我的引用值。
NPAR测试我的实际数据集的/SIGN变量，变量中填充了我的引用值。
但这是一个不可靠的解决办法；( b)结果与我用SPSS得到的Wilcoxon检验结果相去甚远。
什么是正确的语法，还是更好的解决办法？
浏览 2提问于2021-04-22得票数 1
回答已采纳
2回答
我正在使用R中的ggplot2包比较两组长度(不同的个体)和盒图。我想比较两个分布，但到目前为止，我发现使用wilcoxon测试的唯一方法是"ggpubr“包中的stat_compare_means。这是比较分布的正确方式吗？我可以比较分布而不是均值吗？正如您所看到的，我是统计世界中的一个新手。谢谢!
浏览 57提问于2021-03-15得票数 0
下面是我的数据集：
stereograms <- structure(list(Time = c(47.2, 22, 20.4, 19.7, 17.4, 14.7, 13.4,
13, 12.3, 12.2, 10.3, 9.7, 9.7, 9.5, 9.1, 8.9, 8.9, 8.4, 8.1,
7.9, 7.8, 6.9, 6.3, 6.1, 5.6, 4.7, 4.7, 4.3, 4.2, 3.9, 3.4, 3.1,
3.1, 2.7, 2.4, 2.3, 2.3, 2.1, 2.1, 2, 1.9, 1.7, 1.7, 19.7, 16.
浏览 6提问于2020-09-11得票数 0
我已经进行了Wilcoxon秩和检验，以查看三个疾病样本和三个对照样本之间的598019个基因的表达是否有任何显着差异。我在R。
当我看到有多少基因的p值<0.05时，我总共得到了41913。我将Wilcoxon的参数设置如下；
wilcox.test(currRow[4:6], currRow[1:3], paired=F, alternative="two.sided", exact=F, correct=F)$p.value
(这是在apply函数中，如果有必要，我可以提供全部代码，我有点不确定alternative="two.sided“是否正确)。
但是
浏览 22提问于2019-11-12得票数 1
回答已采纳
1回答
我有一个csv文件，有744行186列，格式如下：
Label 1 1 0 0 1
TaxID P_ERR161 P_ERR162 P_ERR163 P_ERR164 P_ERR165
333046 0.05 0 22.33 0.06 7.32
1049 0.03 0 0.04 0.01 0.02
337090 0.01 0 9.79 45.88 3.
浏览 23提问于2018-01-31得票数 1
回答已采纳
4回答
我有两个数据帧：
D9 <- as.data.frame(DF$As,DF$Cd,DF$Cu,DF$Cr,DF$Ni,DF$Pb,DF$Zn)
D10 <- as.data.frame(DO$As,DO$Cd,DO$Cu,DO$Cr,DO$Ni,DO$Pb,DO$Zn)
并希望对每个列(DF$As、DO$As)应用wilcox测试，依此类推。我尝试了以下代码：
lapply(ncol(D9), function(i) {wilcox.test((D9[,i]),(D10[,i]))})
输出结果是：
[[1]]
Wilcoxon rank sum test
data: (D9[,
浏览 7提问于2017-10-10得票数 0
我使用以下语句运行以下代码，以检查一系列day1AverageReturnQuartiles中变量day1TotalStake的中位数是否有任何差异：
proc npar1way data = hsb2a wilcoxon;
class day1AverageReturnQuartile;
var day1TotalStake;
run;
但是，SAS将仅输出Kruskal-Wallace统计数据。它确实报告了每个四分位的秩和，但没有整体统计数据。
我应该调用其他函数吗？
谢谢。
浏览 0提问于2013-06-10得票数 1
回答已采纳
1回答
我根据'yr‘值中的值将'cnt’变量分成两组，这样我就可以用Wilcoxon秩和检验来比较这两组变量。我一直收到这样的错误：“样本x和y必须是一维的。”有人能帮我解决这个问题吗？我定义了两个数组cnt_yr0和cnt_yr1，它们是1xn数组。
#Here is some of my code:
from scipy import stats as sc
cnt_yr0 = np.transpose(np.array(data.loc[data['yr']==0,['cnt']]))
cnt_yr1 = np.transpose(np.arr
浏览 0提问于2020-01-24得票数 1
1回答
我有下面的代码来运行Wilcoxon秩和测试
print stats.ranksums(pre_course_scores, during_course_scores)
RanksumsResult(statistic=8.1341352369246582, pvalue=4.1488919597127145e-16)
然而，我感兴趣的是从结果中提取pvalue。我找不到这方面的教程。有人能帮忙吗？
浏览 3提问于2017-01-18得票数 6
回答已采纳
1回答
scipy.stats.ranksums，wilcoxon，mannwhitneyu都给出了与Matlab的ranksum不同的结果。
浏览 0提问于2015-07-29得票数 2
回答已采纳
我正在构建一个预测模型，该模型旨在预测组织内部的自然减员。我正试图决定是否在我的模型中添加某些预测器。我用Kruskal秩和检验了我的几个预测变量和我的反应变量之间的相关性，发现预测因子与我的反应无关。我还应该在我的模型中包括这些变量吗？由于与我的响应变量缺乏相关性，我倾向于否定，但不想放弃决策树中的潜在“拆分”。我的数据集包含8,225条观测数据和173列，我一直使用它们作为预测器。
浏览 0提问于2018-10-17得票数 1
我有两个不同长度的数据帧(一个是16，另一个是28)。我想用scipy.stats.wilcoxon在这两者之间做一个Wilcoxon测试。为此，我创建了一个函数： def wilcoxon_test(df1, df2):
list_col_1 = df1.columns
list_col_2 = df2.columns
for i in range(0, len(list_col_1)):
name = list_col_1[i]
for j in range(0, len(list_col_2)):
name_check = list_col_
浏览 9提问于2021-07-22得票数 0
我在比较成对的列数据和wilcoxon秩和检验，在大多数比较中，我得到的是完全相同的p值。我在想，从数据判断我是犯了错还是一切都还好。以下是一些比较。
这是我用的比较
wtresult<-wilcox.test(datachunk[,i],datachunk[,(i+1)],paired=FALSE)
下面是上面使用的数据的结果。
X1 X2 X3
339.53 354.11 435.56 425.34 434.64 436.08
X1 X2 X3
312.1 282.2 2
浏览 4提问于2015-07-21得票数 1
回答已采纳
1回答
我正在使用scipy.stats.ttest_ind和scipy.stats.wilcoxon进行t检验和Wilcoxon检验，但我不知道如何为它们设定重要水平。在正式文件中，没有提到任何关于这一问题的内容，这两种方法似乎没有任何争议。有人知道怎么设置吗？
正式文件：
浏览 9提问于2022-07-25得票数 0
回答已采纳
2回答
我想对两个样本x和y执行Wilcoxon测试，其中x和y是长度为n的数值向量。
给定一些实验设计，我想给x中的一些条目赋予一些权重。如何将这样的权重包含在Wilcoxon检验中？我找到了各种各样的包："survey"，"Hmisc“
但是，没有权重的简单测试不会返回标准的Wilcoxon结果:例如：
x=rnorm(n=100,mean=0,sd=1)
y=rnorm(n=100,mean=0.1,sd=1)
wilcox.test(x,y)
data: x and y
W = 4389, p-value = 0.1358
使用WWest函数: wwest(x，y)
浏览 1提问于2014-07-09得票数 0
你好，我正在尝试使用一个看起来像这样的列表。
list <- 0.281097 0.245194 0.172255 0.136205 0.435323 0.431997 0.332863 0.371120 0.428497 0.414634 0.369541 0.303205 0.342573 0.328254 0.358234 0.364887 0.120683 0.118715 0.206876 0.387825 0.059059 0.075471 0.173021 0.359584 0.047235 0.083811 0.159509 0.361754 0.085833 0.162
浏览 0提问于2013-07-13得票数 0
我有两个独立样本的非正态分布数据，将患者分为对照组和治疗组。希望验证两组“控制”和“治疗”之间是否存在差异，并测量这一差异，因此我使用以下代码：
wilcox.test(data.to.work$disease ~ data.to.work$group)
对我的测试来说，我的怀疑是:我可以用Cohen的d来测量效果大小吗？
我还用密码做了一个测试：
cohens_d(data.to.work$disease ~ data.to.work$group)
rcompanion::wilcoxonR(data.to.work$disease, g=data.to.work$group, ci=T)
浏览 15提问于2022-07-29得票数 -1
我正在使用scipy.stats.wilcoxon()运行一个wilcoxon测试，它返回一个W和p值。但是，我也想要z统计数据。有办法用W和p来计算吗？
浏览 1提问于2019-06-20得票数 2
我写信是想问你对如何解释这个案例有什么看法。我有两个向量"a“和"b”，我试图比较它们。
wilcoxon检验给出了a/b的p值为5.139217e-303，另一个选择是“更大”。现在，如果我对它们中的每一个做一个总结，我会得到以下内容
> summary(a)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.0000000 0.0001411 0.0002381 0.0002671 0.0003623 0.0012910
> summary(c)
Min. 1st Qu.
浏览 1提问于2012-09-20得票数 0
回答已采纳
我必须对不同的数据挖掘算法进行三种不同的比较。
唯一有问题的比较类型是最基本的比较，在单个数据集上的两种算法----对我来说是有问题的。
我知道Diettrich (1998)的论文提到了McNemar和5x2CV作为选择的选择，并指出重放t检验是不可行的。由于分析是使用子样本( 60:40培训:测试-拆分和总成本)作为性能度量的更大设置的一部分，我不能使用这些。
在这种情况下，还有哪些其他选项来评估性能？
符号测试：只计算情况数，其中两种算法的性能都更好，然后使用二项分布检查p值。就像很弱一样有问题。
Wilcoxon-signed-rank-test:作为t检验的非参数替代，我想到
浏览 7提问于2013-09-10得票数 1
回答已采纳
我试图在来自数据帧的两个数据子集上运行一个wilcox.test()。它们的长度不一样(48对260)。我想看看活橡树和水橡树的胸径(胸径)是否有差别。
Pine_stand <- read.csv("Pine_stand.csv")
live_oaks <- subset(Pine_stand,Species=="live oak",select=c("dbh"));live_oaks
water_oaks <- subset(Pine_stand,Species=="water oak",select=c(&
浏览 2提问于2017-12-10得票数 3
回答已采纳
我有一个30码的分发。
这是用相同的超参数连续训练30次神经网络的结果(例如，ROC-AUC)，但是由于它们是随机初始化的，所以结果总是有点不同。
然后，我使用其他超参数来训练同一个网络，并且只想在较少的运行情况下这样做。让我们说5次。
我的零假设是，较小的运行分布不小于有30次运行的分布(一侧测试)。
哪种统计显着性检验是比较这些小分布最好的方法？
PS:目前我正在使用测试。还有什么更好的吗？
浏览 0提问于2020-09-11得票数 0
我对R程序非常陌生，但我一直在学习:)
我想对测量前后进行配对的Wilcoxon秩和检验。我的受试者接受了4种不同类型的治疗(即作为因子w/ 4水平的治疗)。我的数据是有序的。
我有一个包含3列的数据框架：
在此之前
之后
治疗类型
我是否可以运行wilcox.test，同时对不同类型的治疗进行细分？
我运行wilcox.test的方式如下：
wilcox.test(Before,After,paired=TRUE)
如何将此代码编辑为我的治疗类型的子集？
希望有人能帮我..。:)
BeforeAfterTreatment22SH21SH22SH11SH22OH22OH22
浏览 1提问于2017-05-26得票数 0
回答已采纳
2回答
我有一个信号，我想测量它与一个特定值的离散度。 有没有一个Matlab函数可以做到这一点？ ? 这是我的信号样本，我想知道这个信号与1V的直流值有多远
浏览 20提问于2019-05-18得票数 0
2回答
基于目标范围的神经网络回归评价
、、、、
目前，我正在拟合一个神经网络来预测一个从1到10之间的连续目标。然而，样本在整个数据集中的分布并不均匀:目标为1-3的样本代表性很低(仅占数据的5%左右)。然而，它们具有很大的兴趣，因为目标的低范围是某种临界范围。
有没有办法知道我的模型是如何预测这些低范围样本的呢？我知道，在进行多类分类时，我可以检查回忆，以了解模型在某个类上的表现。对于分类用例，我还可以在Keras中设置类权重参数，以考虑类的不平衡，但这显然不可能用于回归。
到目前为止，我使用了典型的指标，如MAE，MSE，RMSE，并取得了令人满意的结果。不过，我想知道模型如何在“关键”示例上执行。
浏览 5提问于2019-12-03得票数 1
回答已采纳
我的问题是关于SAS程序NPAR1WAY的。具体地说，我想使用Fligner-Policello测试，该测试应该通过"FP“选项调用。该检验是对两样本的标准Wilcoxon秩和检验的一种调整。
SAS帮助网站上记录了此选项；但是，没有给出任何示例。我还找不到任何正确语法的例子来使用这个测试。
下面是我到目前为止尝试过的语法，这似乎是SAS文档建议的语法：
proc npar1way fp data=input;
class stratum;
var varname;
run;
以下是错误消息：
R 22-322: Syntax error, expecting one
浏览 0提问于2014-04-07得票数 0
关于我的威尔克森签名等级测试的结果，我有一个疑问：
我的数据包括两组(配对)的试验，其中使用了一种治疗。结果以%为单位。团体由131人组成。
当我在R中运行测试时，我得到了以下结果：
wilcox.test(no.treatment, with.treatment, paired=T)
# Wilcoxon signed rank test with continuity correction
# data: no.treatment and with.treatment V = 3832, p-value = 0.7958
# alternative hypothesis: true
浏览 0提问于2017-07-21得票数 7
回答已采纳
我正在尝试使用for语句创建一个循环，该语句将由统计测试(wilcoxon.test)分析的变量的名称添加到数据帧单行中所有相关测试统计信息的左侧。 我发现除了包含'i‘引用的第二个' for’语句之外的代码工作正常-所以统计数据是正确的，但是当我添加代码(包括引用'i‘的for语句，并在'rbind’下指定的向量下列出'i‘时，所有以前正确的统计信息都会被一组来自单个测试的统计结果覆盖。 dta <- c()
max_col_A1.1_OR_df_transp <- ncol(A1.1_OR_df_transp)
for(j in
浏览 19提问于2019-04-17得票数 0
我真的是新手，助教让我来这里为RStudio寻求帮助。这周我们将学习威尔科克森秩和精确检验。我知道(我认为)测试返回的P值是，如果P< 0.05，那么数据来自不同的分布。但是，作业问题是这样问的： 对于因杀婴而失去一只幼仔的雌性，下一只幼仔的存活时间是(更长/更短)的_。 我被卡住了，我不能解决这个问题，我需要一些帮助！ 这是一些背景下的完整家庭作业问题。 “当入侵的雄狮接管一群雌性狮子时，它们通常会杀死群体中的大部分或所有幼狮。在接管之后，这种骄傲的稳定性是不可预测的，科学家们假设，雌性狮子可能会推迟排卵，直到不确定性过去。为了验证这一点，一项针对塞伦盖蒂地区狮子的研究测量了雌性狮子
浏览 51提问于2021-03-18得票数 0
我使用来自https://en.wikipedia.org/wiki/Wilcoxon_签名等级_测试的数据，W值是9。但是，对于下面的代码W值是18，原因是什么？
# Wilcoxon signed-rank test
from numpy.random import seed
from numpy.random import randn
from scipy.stats import wilcoxon
# seed the random number generator
seed(1)
# generate two independent samples
a=np.array([125
浏览 0提问于2018-11-06得票数 0
回答已采纳
1回答
我试图用Python做一些统计。我有几个缺失值的数据，填充了np.nan，我不确定是手动删除它，还是can能够处理它。所以我两个都试过了：
import scipy.stats, numpy as np
a = [0.75, np.nan, 0.58337, 0.75, 0.75, 0.91663, 1.0, np.nan, 0.663, 0.837, 0.837, 1.0, 0.663, 1.0, 1.0, 0.91663, 0.75, 0.41669, 0.58337, 0.663, 0.75, 0.58337]
b = [0.837, np.nan, 0.663, 0.583
浏览 3提问于2015-05-10得票数 3
回答已采纳
1回答
我想使用ggpubr包来引用全局环境中字符串中列出的数据帧列名，但是ggpubr似乎不接受变量，只接受硬编码的列名。有什么办法我可以做任何改变，这样它就能做到这一点吗？
vars = c('var1', 'var2')
controls = c('a', 'w')
df = data.frame(subject = 1:100,
value = rnorm(100, 100, 10),
var1 = rep(c('a', 'b'),
浏览 2提问于2019-08-08得票数 1
回答已采纳
1回答
：不等N
、、
Rs wilcox.test可以采取不同的长度向量，但是scipy.stats的wilcoxon不能:我得到一个unequal N错误消息。
from scipy.stats import wilcoxon
wilcoxon(range(10), range(12))
有没有一种方法可以在Python中获取Rs行为？
浏览 0提问于2015-11-24得票数 8
回答已采纳
我有一个暗淡的15x555的数据帧，其中行方向的样本号与列方向的蛋白质名称相对应。该数据帧的最后3列保存映射信息，即处理、Treatment_Time和月份，并且被如此标记。
在以列的方式遍历数据帧时，我希望使用wilcoxon.test执行wilcoxon测试，并根据映射信息直接从df获取信息。
粗略的例子：
pre_post <- vector()
for(i in names(df[,1:552])){
pre_post <- append(pre_post, wilcox.test(df[df$Treatment_Time %in% "Pre",
浏览 0提问于2018-03-10得票数 1
我正在努力计算R中数据帧中两组之间的t-test。下面的示例代码产生了一个有两列的数据帧:变量和值。有两个变量："M“和"F”。
data <- data.frame(variable = c("M", "F", "F"), value = c(10,5,6))
我需要证明M和F的值在统计上是不同的。换句话说，10在统计上与5和6的平均值不同。我需要在此数据框中添加另一列来显示p值。当我运行下面的代码时，它给出了以下错误：
result <- data %>% mutate(newcolumn = t.test(
浏览 5提问于2020-02-27得票数 0
1回答
我是这方面的新手。在使用统计测试对算法进行比较时，我面临一些问题。本文给出了几种分类算法的Gmean结果。鲍鱼、平衡秤、Car、国际象棋是这里的数据集，ROS、RUS、RFS、NoS是算法.
如果我想用一些统计检验(如t检验、Friedmen检验、Wilcoxon检验等)来比较哪个数据集比其他数据集更好，那么我能用下表比较算法吗？
ROS RUS RFS NoS
Abalone 0.003 0.0036 0.0039 0
Balance-scale 0.8858 0.8065 0.8
浏览 0提问于2020-01-10得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
