Question

有关于R语言中t.test()函数在两样本t检验应用时的参数设置问题。

各位大侠好，小弟用R语言写的一个程序中用到了两样本t检验，目的在于比对两组小样本之间的均值是否相同。由于有很多对对应的组，迫于无奈以及懒惰，在R中写了循环，当循环进行到其中一对比对的组时，由于两组的数据全都是0或者NA所以检验得到的p值就变成了NA循环报错，我检查过报错的那一对数据，连接是我截下来的图：http://xiangce.baidu.com/picture/album/list/80d4e3d2f164c466e033557085593b22f3287667

求问：如何设置t.test()参数使两组数据分布情况与上述情况类似时仍然可以求得p值。
ps：全部为0的情况视为两组数据的分布不存在差异，p值应当大于0.05.
我试图用正太性检验排除数据中包括全都是0在内的数值等于恒量的数据，但是出现了一下问题：

http://xiangce.baidu.com/picture/detail/68793c5fc47fbf3831368738ad1d07e01ce9ee1b?from=dialog

ps：请原谅我一次问多个问题。

Answer 1

问题关键在于数据全为0，这样在算t检验统计里时，其分母为0，故得到NaN。
因此在t.test设置任何参数都是没有用的，你应该在t.test之有加条件判断出现这种情况 应该 先令
p>0.05才行

追问

是的，因此我昨天晚上在循环程序之前附加了这种全是0的情况排除，但很可惜这批数据中还有全是1的情况，这让我想到数据中有一些全是恒量的数据并不是正太分布的，实际上这些数据也需要排除掉，于是我想到了在循环之前先进行正态性检验的问题。但是似乎这也不容易。

我已经将出现的错误截图，大侠您有空的话请看问题补充，同时我也在寻找答案。
ps：我是学医的，统计程序只学到些皮毛，但只要有线索我就能继续挖下去的。

追答

你就干脆先用t.test出现p=nan再判断

追问

对 我接着以您给的算法做了程序，但是如果数据是其他恒量比如说全是1的时候它t.test()是可以的到p.value但是由于数据都是恒量，依然会报错，导致循环无法继续。具体情况我已经截图放在问题补充之中了，还请大侠不吝赐教。

追答

数据全一样的道理是一样的，方差都为零，导致t检验统计量没有意义。
你用is.na只能判断缺失值，而这个实际上t.test自己会判断。
你应该判断是不是恒量，是恒量要就不要调用t.test，不是恒量才调用