Scypy: Wilcoxon测试:比较分布与单个值
Scipy: Wilcoxon测试是一种非参数统计检验方法，用于比较两个相关样本或配对样本的差异。它的目的是判断两个样本是否来自同一总体分布，而不需要对总体分布做出任何假设。
Wilcoxon测试的步骤如下：
- 收集两个相关样本或配对样本的数据。
- 对每对数据进行差异计算。
- 对差异值进行排序，并计算排序的绝对值秩次。
- 计算秩次和，以及秩次和的正负符号。
- 根据秩次和的正负符号，查找Wilcoxon秩次和表格，确定显著性水平下的临界值。
- 比较计算得到的秩次和与临界值，判断两个样本是否存在显著差异。
Wilcoxon测试的优势在于它不需要对数据的分布做出任何假设，适用于非正态分布或小样本量的情况。它也可以用于比较有序分类变量或等级数据。
Wilcoxon测试的应用场景包括但不限于：
- 医学研究中比较两种治疗方法的效果。
- 心理学研究中比较两组受试者的表现。
- 生物学研究中比较不同基因表达的差异。
- 工程领域中比较两种产品的性能。
腾讯云提供了一系列与数据分析和统计相关的产品和服务，可以用于支持Wilcoxon测试的实施和分析，例如：
- 腾讯云数据万象（https://cloud.tencent.com/product/ci）：提供图像和视频处理能力，可用于处理多媒体数据。
- 腾讯云人工智能（https://cloud.tencent.com/product/ai）：提供各种人工智能服务，如图像识别、语音识别等，可用于处理和分析数据。
- 腾讯云数据库（https://cloud.tencent.com/product/cdb）：提供各种数据库服务，如云数据库MySQL、云数据库MongoDB等，可用于存储和管理数据。
- 腾讯云云服务器（https://cloud.tencent.com/product/cvm）：提供云服务器实例，可用于部署和运行应用程序。
请注意，以上仅为腾讯云相关产品的示例，其他云计算品牌商也提供类似的产品和服务，具体选择应根据实际需求和预算进行评估。
相关·内容
1回答
我正在对实验前后每个人的平均心痛数据进行wilcoxon测试。大约有62人参加。
我跑：
statistic, p-value = wilcoxon(dataBefore - dataAfter)
P值为0.0005，而统计量为498.这似乎是一个异常大的数字。我注意到在这里归还的东西有些不清楚。我以为this使用了t检验，但我也看到人们把它作为z统计数据存储。它还回来了什么？
这本手册对我来说没有多大意义：它说：“统计:浮点--是零上或零下差值的总和，以较小的值为准。”但这对我来说没什么意义。这是人们通常报告的价值吗？
浏览 6提问于2016-10-27得票数 0
回答已采纳
1回答
scipy.stats.ranksums，wilcoxon，mannwhitneyu都给出了与Matlab的ranksum不同的结果。
浏览 0提问于2015-07-29得票数 2
回答已采纳
我使用以下函数来计算数据帧(X)中数据的t-stat：
wilcox.test.all.genes<-function(x,s1,s2) {
x1<-x[s1]
x2<-x[s2]
x1<-as.numeric(x1)
x2<-as.numeric(x2)
wilcox.out<-wilcox.test(x1,x2,exact=F,alternative="two.sided",correct=T)
out<-as.numeric(wilcox.out$statisti
浏览 2提问于2011-11-12得票数 1
我运行了两种不同的元启发式算法25次，我想检查哪个算法的结果优于其他算法。我决定使用Wilcoxon测试，但我不能很好地理解给定函数的输出：
from scipy.stats import ranksums
rng = np.random.default_rng()
sample1 = rng.uniform(-1, 1, 200)
sample2 = rng.uniform(-0.5, 1.5, 200)
ranksums(sample1, sample2)
Output: RanksumsResult(statistic=-7.887059, pvalue=3.09390448e-15)
浏览 12提问于2022-07-25得票数 0
我执行"C:\Program Files\QlikView\qv.exe“/R "D:\QlikViewAdm\qlik_file.qvw”来启动加载程序。当没有区段访问权限时，它可以工作。然而，当我运行一个文件的节访问QV要求我登录和密码。我可以链接到Windows用户凭据的访问权限，但它只能在我的终端上工作。
是否有方法将QV指向具有凭据的本地文件？
浏览 0提问于2021-05-28得票数 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
回答已采纳
我想对我的配对数据进行单侧wilcoxon秩检验，因为我感兴趣的是其中一个样本是否明显大于另一个样本。
Scipy提供
scipy.stats.wilcoxon(x,y)
对成对的样本x和y进行双边检验。因为我不能假设正态(对称)分布，所以我不能从双边p值推导出单侧p值。
现在有没有人可以用python的方法来获得单边测试的p值？
谢谢!
浏览 3提问于2013-04-30得票数 14
回答已采纳
我写信是想问你对如何解释这个案例有什么看法。我有两个向量"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
回答已采纳
2回答
我正在使用R中的ggplot2包比较两组长度(不同的个体)和盒图。我想比较两个分布，但到目前为止，我发现使用wilcoxon测试的唯一方法是"ggpubr“包中的stat_compare_means。这是比较分布的正确方式吗？我可以比较分布而不是均值吗？正如您所看到的，我是统计世界中的一个新手。谢谢!
浏览 57提问于2021-03-15得票数 0
1回答
我有一组样本，每个样本都含有元素A的浓度，在下面的柱状图X轴上用对数标度显示。相同浓度的样品数显示在Y轴上.正如您在直方图中看到的那样，分布类似于多模分布。基于我正在使用的实验，左模态，我知道一个事实，那就是仪器噪声(，这是我如何定义噪声:那些数据点大于data+3*standard偏差的平均值)，但右模态是真实的数据。所以，基本上，对于那些在左模态上显示浓度的样本，实际浓度值为零。
我的问题是；
1-我需要一个测量值(p值)，它给我一个概率值，从右边模态中的每个数据点属于(或不)属于左模态(适合或不适合)。
2-如果我在右边只有一个数据点，而不是在右边的另一个分布(正确的模式)，我将如何度量
浏览 2提问于2014-09-19得票数 1
因此，我遇到了一个奇怪的差异，这取决于我如何使用R中的wilcox.test()函数来分析相同的数据，在这个例子中，我正在比较两个组中的值。我可以将它们作为两个独立的向量交给wilcox.test函数，也可以给函数一个data.frame，并使用一个公式来指定我想要进行的比较。奇怪的是，我最终得到了测试统计量(W)的不同值，这取决于我使用的输入法。下面是一个例子(R v3.3.1)：
#Prepare test data
wt_exp = c(0.59, 0.56, 0.45, 0.59, 0.54, 0.13, 0.25, 0.10, 0.15)
kd_exp = c(0.27, 0.27,
浏览 2提问于2017-04-12得票数 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
回答已采纳
2回答
正态得分
、、
具有以下分布(实际分布和预测)，Hist 1到3(左到右)。
📷
我想得到一个分数范围从0-1的实际分布是多么接近正常。我发现了一些统计正态性检验：
夏皮罗-威尔克试验
D‘’Agostino^2检验
我的DataSet很大，因此我决定检查斜度和峰度统计数据，并得到以下结果：
hist-1 Skewness is 0.028386209063816035 and Kurtosis is 2.4224694251429764 <-- Most normal
hist-2 Skewness is 3.7702212103585246 and Kurtosis is 15.2145679750
浏览 0提问于2021-06-28得票数 1
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
1回答
我有一个数据帧，我试图在每一行上应用T-test，但它给了我nan。
代码：
from scipy.stats import ttest_ind, ttest_rel
import pandas as pd
df_stat = df_stat[['day', 'hour', 'CallerObjectId', 'signals_normalized', 'presence_normalized']]
def ttest(a, b):
t = ttest_ind(a, b)
return t
d
浏览 1提问于2021-04-27得票数 0
我有一个30码的分发。
这是用相同的超参数连续训练30次神经网络的结果(例如，ROC-AUC)，但是由于它们是随机初始化的，所以结果总是有点不同。
然后，我使用其他超参数来训练同一个网络，并且只想在较少的运行情况下这样做。让我们说5次。
我的零假设是，较小的运行分布不小于有30次运行的分布(一侧测试)。
哪种统计显着性检验是比较这些小分布最好的方法？
PS:目前我正在使用测试。还有什么更好的吗？
浏览 0提问于2020-09-11得票数 0
1回答
我正在使用scipy.stats.ttest_ind和scipy.stats.wilcoxon进行t检验和Wilcoxon检验，但我不知道如何为它们设定重要水平。在正式文件中，没有提到任何关于这一问题的内容，这两种方法似乎没有任何争议。有人知道怎么设置吗？
正式文件：
浏览 9提问于2022-07-25得票数 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
我在R非常新，可能没有经验排除这个问题是否与以前的帖子密切相关。如果是的话，请接受我的道歉。
我有前后测量(分别为X0variablename和X1variablename )。变量是导入的.CSV数据集中的列。我对每个变量进行了11次观察，并希望检验是否存在真正的差异(meadian是使用describe (心理包)计算的)。
dataset中的某些列不适合此测试，因此我手动编写了所有这些列：
changeipss <- wilcox.test(mydata$X0ipss, mydata$X1ipss)
changeipss
changeqol <- wilcox.test(m
浏览 0提问于2018-09-06得票数 0
2回答
基于目标范围的神经网络回归评价
、、、、
目前，我正在拟合一个神经网络来预测一个从1到10之间的连续目标。然而，样本在整个数据集中的分布并不均匀:目标为1-3的样本代表性很低(仅占数据的5%左右)。然而，它们具有很大的兴趣，因为目标的低范围是某种临界范围。
有没有办法知道我的模型是如何预测这些低范围样本的呢？我知道，在进行多类分类时，我可以检查回忆，以了解模型在某个类上的表现。对于分类用例，我还可以在Keras中设置类权重参数，以考虑类的不平衡，但这显然不可能用于回归。
到目前为止，我使用了典型的指标，如MAE，MSE，RMSE，并取得了令人满意的结果。不过，我想知道模型如何在“关键”示例上执行。
浏览 5提问于2019-12-03得票数 1
回答已采纳
我在比较成对的列数据和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
回答已采纳
我必须对不同的数据挖掘算法进行三种不同的比较。
唯一有问题的比较类型是最基本的比较，在单个数据集上的两种算法----对我来说是有问题的。
我知道Diettrich (1998)的论文提到了McNemar和5x2CV作为选择的选择，并指出重放t检验是不可行的。由于分析是使用子样本( 60:40培训:测试-拆分和总成本)作为性能度量的更大设置的一部分，我不能使用这些。
在这种情况下，还有哪些其他选项来评估性能？
符号测试：只计算情况数，其中两种算法的性能都更好，然后使用二项分布检查p值。就像很弱一样有问题。
Wilcoxon-signed-rank-test:作为t检验的非参数替代，我想到
浏览 7提问于2013-09-10得票数 1
回答已采纳
1回答
最近，我想比较两组点，并找到了R软件包‘k’。函数kde.test (基于内核密度的全局双样本比较测试，而不是KS测试)是我所需要的。例如
library(MASS)
data(crabs)
x1 <- crabs[crabs$sp=="B", c(4,6)]
x2 <- crabs[crabs$sp=="O", c(4,6)]
kt <- kde.test(x1=x1, x2=x2)
但不幸的是，我对R并不熟悉，我的大部分工作都是使用Python完成的，因此很难将我的脚本传输到R。同时，'ks‘的函数(主要是kde.test)对我的工
浏览 2提问于2022-06-01得票数 0
我正在对我的数据进行威尔科克森秩和检验。数据具有不同的单位和很大程度上不同的尺度。在运行测试之前，我需要键入这个命令吗？
## where x is my data frame
scale(x, center = TRUE, scale = TRUE)
或者，Wilcoxon秩和检验本身就是这样做的？
浏览 0提问于2013-02-14得票数 0
回答已采纳
当一个人使用pnorm命令生成一个图时，它会生成一个带有单元的图：
Y轴：Normal F[(Variable Name-m)/s]
X轴：P[i] = i/(N+1)
用手计算X轴似乎是合理的.我搞不懂Y轴的单位是什么意思？
正常的Normal F[(Variable Name-m)/s]是如何分解的？m代表平均值，s代表标准差吗？如果是这样的话，Normal F()代表什么函数？
浏览 4提问于2017-01-26得票数 1
回答已采纳
我正在研究一个数据不平衡的二进制分类问题(阳性类为17%)。
问题如下:当我在测试集上进行预测时，我的三个模型(我没有标签)给出了与列车集非常相似的分布。
但是对这些模型进行建模，在给出较好的结果(F1-评分)的同时，它大大改变了测试集上的分布，从~20%到5%。
我的问题是:
我很困惑于选择一个最好的个体模型，这个模型保持了几乎相同的分布，但失去了一些效率，或者是给出了真正不同的分布的集合模型。
我不知道测试集的分布。
谢谢你的帮助
浏览 0提问于2021-01-31得票数 0
回答已采纳
我正在寻找一种在保留原始组信息的同时使用group_split()或summarise()语法的方法。我见过一些以前的页面，如和使用这些方法，但它们不保留分组信息。有办法这样做吗？当然，我可以加入数据，但我希望避免使用这种方法。
> set.seed(22)
> # Create fake data
> flavor <- data.frame(
+ temperature = sample(x = c('hot','cold'), size = 500, replace = TRUE),
+ color = sample(c(&
浏览 4提问于2022-09-05得票数 1
回答已采纳
我正试图为我的数据获得最好的发行版。拟合完成如下图所示，但我需要一个测量，以选择最好的模式。我将拟合优度与a-平方值进行了比较，并用Kolmogorov-Smirnov (KS)检验了观测值与拟合分布的显着性差异。我搜索了一些潜在的解决方案-- 、、，但没有得到以下结果：
如果p-value大于k-statistic，是否意味着我们可以接受假设或数据很好地符合分布？
Alternatively，是否可以将level of significance(a=0.005)与p-value进行比较并决定是否接受或拒绝假设？如果p-value低于a，那么这两种分布很可能是不同的。
For Kolmogor
浏览 3提问于2020-09-29得票数 3
我正在使用PSPP，并希望比较一个样本的中值与给定的中值值，但我发现的所有选项都是比较两个变量的相互比较。通过定义一个充满引用值的变量，我尝试了一种解决方法：
NPAR测试我的实际数据集的/Wilcoxon变量，变量中填充了我的引用值。
NPAR测试我的实际数据集的/SIGN变量，变量中填充了我的引用值。
但这是一个不可靠的解决办法；( b)结果与我用SPSS得到的Wilcoxon检验结果相去甚远。
什么是正确的语法，还是更好的解决办法？
浏览 2提问于2021-04-22得票数 1
回答已采纳
我试图在R和python的scipy.stats包中运行Wilcoxon测试。但是我得到了不同的结果，有人能解释吗？
R中的我的代码
> des2
[1] 6.2151308 4.7956451 4.7473738 5.4695828 6.3181463 2.8617239
[7] -0.8105824 3.9456856 4.6735000 4.1067193 5.7656002 2.2237666
[13] 1.0354143 4.9547707 5.3156348 4.8163154 3.4024776 4.2876854
[19]
浏览 2提问于2015-11-07得票数 1
回答已采纳
1回答
我正试图计算KS测试，将其指定为一个数组，但是，我遇到了意想不到的结果。通过进一步的评估，我发现了不同的结果，这取决于我是将CDF指定为可调用的字符串还是数组。我的代码如下：
import scipy.stats as st
random_variables = st.norm.rvs(loc=1, scale=1,size=1000000)
cdf_data = st.norm.cdf(random_variables, loc=1,scale=1)
params = st.norm.fit(data=random_variables)
display(params)
print('
浏览 1提问于2021-03-15得票数 1
回答已采纳
我得到了这个数据集，(mydata.csv)，我想做一个非参数配对测试(wilcoxon)，但当我运行代码时，出现了这个错误，我不明白它的含义和如何解决： import pandas as pd
import pingouin as pg
df = pd.read_csv("mydata.csv")
id price review score
0 7949480 99.0 Check-in 10.0
1 6627449 125.0
浏览 136提问于2020-12-10得票数 0
1回答
我使用了Network核心应用程序来获取Cytoscape上的无定向网络的基本参数。对各节点的度、中心度、网络直径等参数进行了满意的测量，给出了每个节点的聚类系数为0.0，计算出网络的总体聚类系数为0.0。接下来我要比较我的网络和随机网络，网络效率是我想比较的一个关键指标，以证明我的网络是一个无标度的网络。可能出了什么问题。在我的网络中有361个节点和695个边。任何想法都会受到赞赏。
浏览 1提问于2017-04-19得票数 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
我怎么才能拉z？
渐近Wilcoxon-Pratt签名-秩检验
data: y by x (pos, neg)
stratified by block
Z = 2.5113, p-value = 0.01203
alternative hypothesis: true mu is not equal to 0
浏览 2提问于2021-04-27得票数 0
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
我不知道如何在Scipy中进行两个样本的KS测试。
在阅读文档之后
我可以看到如何测试一个分布是否与标准正态分布相同。
from scipy.stats import kstest
import numpy as np
x = np.random.normal(0,1,1000)
test_stat = kstest(x, 'norm')
#>>> test_stat
#(0.021080234718821145, 0.76584491300591395)
这意味着在p值为0.76时，我们不能拒绝两个分布相同的零假设。
然而，我想要比较两个分布，看看我是否可
浏览 4提问于2012-06-05得票数 99
回答已采纳
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
回答已采纳
我使用statannot对一些基本数据进行统计测试，但是统计测试的结果似乎不正确。也就是说，我的一些比较得出了"P_val=0.000e+00 U_stat=0.000e+00"，我认为这是不可能的。我的数据框架和/或代码有什么问题吗？
下面是我使用的数据框架：
这是我的密码：
import pandas as pd
from matplotlib import pyplot as plt
import seaborn as sns
from statannot import add_stat_annotation
import scipy.stats as sp
da
浏览 3提问于2021-06-09得票数 0
回答已采纳
我有两个不同长度的数据帧(一个是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
我制作了一个随机数发生器，在0到10007范围内生成了N个数，我想测试这个发生器在R中有多好。
浏览 1提问于2014-03-15得票数 0
回答已采纳
我运行了以下模型
m1 <- lmer ( XX ~ XX_cov
+ time + time:group + ZZ_cov + YY_cov
+ age + (1 | ID), data = dat_long0 )
摘要(M1)
> Linear mixed model fit by REML. t-tests use Satterthwaite's method
> [lmerModLmerTest ] Formula: XX ~ XX_cov + time + time:group + ZZ_cov + YY_cov + age
浏览 6提问于2022-10-11得票数 -1
我在4个矩阵中存储了来自4种不同方法的产生(q)值。4个矩阵中的每个矩阵都包含来自不同方法的Q值，如下所示：
Matrix_1 = 1 row x 20 column
Matrix_2 = 100 rows x 20 columns
Matrix_3 = 100 rows x 20 columns
Matrix_4 = 100 rows x 20 columns
列数表示年数。一行将包含与20年相对应的生产值。矩阵2、3和4的其他99行只是不同的实现(或模拟运行)。因此，矩阵2、3和4的其他99行基本上是重复情况(但由于随机数而不是精确值)。
将Matrix_1作为参考事实(或基
浏览 2提问于2010-09-17得票数 1
回答已采纳
我试图使用wilcox.test()函数在几个数据子集上运行lapply()。这些数据按照文本标识符(网站名称)分组在我的数据框架的第一列中，还有另外两列用于2013年和2017年的数据，我正在对这两列进行wilcox测试。在我的500行数据中，约有10%的数据在2013年或2017年(或两者都)的列中缺失。
当我试图运行如下所示的lapply()函数时，会得到以下错误：
df<-read.csv("myfile.csv",header=T)
split.df<-split(df,df$Site)
lapply(split.df, function(g) wilc
浏览 0提问于2019-05-21得票数 0
回答已采纳
在少数数据集上使用了两种算法(方法)，并取得了一定的结果。现在我想检查一下所得到的结果是否具有可比性？我使用了两个样本的K测试，得到了以下结果，现在如何解释测试结果？结论应该是什么？
KstestResult(statistic=0.11320754716981132, pvalue=0.8906908896753045)
D_alpha= 0.2641897545611759
浏览 4提问于2021-07-11得票数 0
回答已采纳
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
我正在使用一个数据集，其中包含申请者的收入、性别和贷款状况(无论此人是否已被批准贷款)。我根据这些数据创建了下面的图表。直方图图是：
📷
核密度估计(KDE)图是：
📷
KDE地块似乎表明，与妇女相比，在某一特定收入中，男子接受的与被拒绝的比率更高。我想进一步调查一下。注(!)数据集中的男性多于女性，因此任何结论都需要考虑到差异。
一个想法:我最初的想法是把收入存起来，计算每个性别在每个垃圾箱中被接受/拒绝的比率。然后，我们可以绘制比率和方差(使用每个垃圾箱中的男女计数)，看看是否有统计意义的依赖于接受/拒绝的性别。
问:上述观点是否合理？我应该提出一个假设检验问题吗？如果是的话，我该怎么做呢
浏览 0提问于2022-01-23得票数 0
回答已采纳
1回答
我有一个问题，我是否可以在一个循环中对所有生成的表进行Wilcoxon测试。
基本上，我想在每个数据集的2个变量之间进行配对的Wilcoxon测试，并且这2个变量对于每个数据集都在相同的位置(如第x列和第y列)。(对于熟悉生物学的人来说，实际上这是一些重复元素的对照和处理样本之间的RPKM值)，我希望我可以为每个数据集的Wilcoxon检验的p值生成一个表。
我已经使用下面的代码生成了所有的表/数据集/数据帧，我想我想为每个数据集做一个Wilcoxon测试，所以我认为我需要继续这个循环，但我不知道怎么做：
data=sample_vs_norm
filter=unique(data$famil
浏览 30提问于2020-05-27得票数 1
回答已采纳
1回答
用C++实现字符串的二进制搜索
、、、、
我实现了这个二进制搜索算法，目的是找到包含在数组中的所需字符。考虑到ASCII表，我已经尝试用标准方法进行' if‘else’比较：[if( searchElement > arrayCharmid ...等]。这似乎不起作用，所以我意识到这可能与字符串比较有关。我现在使用函数strcmp并使用它的返回值进行比较：
(如果string1 < string2 ->负值。If string1 > string2 ->正值。If string1 == string2 -> 0.)
看起来效果不是很好。
#include <stdlib.h>
#
浏览 0提问于2021-01-10得票数 0
1回答
我想查看特定的图像或数据集的分布，并查看它们是否不同。
只是简单地写了这样的东西：
# mydataset.shape = (50k,32,32,3)
plt.hist(mydataset.reshape(-1))
玩这个把戏吗？还是我该做点别的？例如，在cifar10上这样做给了我这个情节：
📷
然而，它看起来不对，有50K的训练图像，我不知道如何解释这一点，即使这是正确的方法，在第一！
如果我这样做了：
#dataset shape is (50k, 3072)
plt.hist(mydataset.reshape(-1,32*32*3))
#and testset which is (
浏览 0提问于2019-02-17得票数 6
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
