我想在R中创建一个函数，使用双尾t检验计算成对零假设的所有（k选择2）p值。
该函数将采用两个参数x和y，其中x = k个不同级别的因子变量，y =实数值的向量。 （x和y具有相同的长度）。
每个不同的x值至少有2（x，y）对。 我的想法是我应该先计算t值，计算（m1-m2）/ se，其中se = sqrt（sd1 * sd1 / n1 + sd2 * sd2 / n2）..然后用pt计算p值（ -abs（T），DF = PMIN（N1，N2）-1 ..
我应该使用某种for循环来计算每对中的所有p值吗？
由于这是一个“pval = function（x，y）{＃mycodehere}”函数，我不知道如何为每对进行这样的计算。
任何人都可以帮我这个吗？！ 非常感谢！:)