A**H 发帖数: 4797 | 1 我有近200个样品,每个样品有两组数据,一个是浓度Conc(0.01~50.00),一个是效
果值Read(0~99999)
因为有些效果值超大或者超小,我怀疑是样品的浓度有点问题,所以我想把浓度和效果
值之间做一个paired Wilconxon test
> wilcox.test(t$Read, t$Conc, paired=TRUE)
Wilcoxon signed rank test with continuity correction
data: t$Read and t$Conc
V = 4656, p-value < 2.2e-16
alternative hypothesis: true location shift is not equal to 0
这里显示p value极小
我把两组数据打乱,我还是得到极小的p value
> wilcox.test(sample(t$Read), sample(t$Conc), paired=TRUE)
Wilcoxon signed rank test with continuity correction
data: sample(t$Read) and sample(t$Conc)
V = 4641, p-value < 2.2e-16
alternative hypothesis: true location shift is not equal to 0
我用样品号码,结果也得到极小p value
> wilcox.test(sample(t$SampleNum), rev(sample(t$Conc)), paired=TRUE)
Wilcoxon signed rank test with continuity correction
data: sample(t$Barcode) and rev(sample(t$Conc))
V = 4656, p-value < 2.2e-16
alternative hypothesis: true location shift is not equal to 0
这是怎么回事?求指点 |
s**r 发帖数: 669 | 2 paired Wilconxon test其实和paired T test的null hypothesis都是一样的,就是H0:
mean(Read)=mean(Conc),从data的range来看p-value很小是很正常的。 |
a*******g 发帖数: 80 | 3 paired test 实际是个one sample test
H0:
mean(read-conc)=0
是不是你样本量比较大?样本量大可以造成加阳性。
另外你看看你那几个极值 是不是read conc同大 同小 如果是 这个极值没啥影响
又仔细看下题 是有点奇怪哦 楼主是不是说 随机乱排结果还是一样 你试了很多遍吗
?每次都这样?
H0:
【在 s**r 的大作中提到】 : paired Wilconxon test其实和paired T test的null hypothesis都是一样的,就是H0: : mean(Read)=mean(Conc),从data的range来看p-value很小是很正常的。
|
A**H 发帖数: 4797 | 4 因为担心你所说的问题,我把浓度都乘以了1000,这样就基本上都达到了效果值的水平
了,也还是一样让人困惑的结果,就是说不管怎么样sample, reverse的搞p value都
很小 。。。
为毛?
H0:
【在 s**r 的大作中提到】 : paired Wilconxon test其实和paired T test的null hypothesis都是一样的,就是H0: : mean(Read)=mean(Conc),从data的range来看p-value很小是很正常的。
|
A**H 发帖数: 4797 | 5 200个样本
我试试看你说的
【在 a*******g 的大作中提到】 : paired test 实际是个one sample test : H0: : mean(read-conc)=0 : 是不是你样本量比较大?样本量大可以造成加阳性。 : 另外你看看你那几个极值 是不是read conc同大 同小 如果是 这个极值没啥影响 : 又仔细看下题 是有点奇怪哦 楼主是不是说 随机乱排结果还是一样 你试了很多遍吗 : ?每次都这样? : : H0:
|
a*******g 发帖数: 80 | 6 刚看见你的read 和 conc 差这么大 那当然p很小了 因为paired test 是test paird的
差值是不是0 你这个不可能是0啊
【在 A**H 的大作中提到】 : 200个样本 : 我试试看你说的
|
a*******g 发帖数: 80 | 7 哦 你把乘完以后的conc和read 相减 画个histogram 或box plot 看看location在不在
0
【在 A**H 的大作中提到】 : 因为担心你所说的问题,我把浓度都乘以了1000,这样就基本上都达到了效果值的水平 : 了,也还是一样让人困惑的结果,就是说不管怎么样sample, reverse的搞p value都 : 很小 。。。 : 为毛? : : H0:
|
x****6 发帖数: 27 | 8 我有点蒙了,楼主不是要看浓度对效果有没有影响吗?不是应该把浓度分成两组,浓和
稀,然后比较这两组的效果有没有差别吗?怎么会浓度的值和效果的值比呢?? |
A**H 发帖数: 4797 | 9 谢谢anneeding and sqrr两位
我发现我把浓度乘以1000还是不太够效果值的范围,改成了10000,发现p value变小了
很多,随机话之后就只有0.01了,所以应该还是这个效果值-浓度的差异太大的问题 |
A**H 发帖数: 4797 | 10 我这不是要设计实验,可以搞两个浓度,而是浓度和效果值都已经出来了,现在怀疑出
问题是因为浓度值变化引起的效果值变化
而且也不是算比值,是算差值和rank
【在 x****6 的大作中提到】 : 我有点蒙了,楼主不是要看浓度对效果有没有影响吗?不是应该把浓度分成两组,浓和 : 稀,然后比较这两组的效果有没有差别吗?怎么会浓度的值和效果的值比呢??
|
g******2 发帖数: 234 | |