R z<；-是什么意思排序（x%*%y）怎么办？

我试图理解TA在stats类中发布的问题集的解决方案，他使用以下命令计算问题中的p值：

FisherMeans <- sort(m1%*%Yobs)/NT # FisherMeans contains all possible outcomes.

---

m1

1的每一行中的值正好是一半，而另一半是-1？如果是这样，那么

（m1%*%Yobs）/NT

相当于底部的代码。但是，上面的代码正在计算排序版本，而下面的代码没有排序。这个代码等价的原因实际上是一个数学观察，而不是编码观察。写出几行矩阵积，

m1%*%Yobs

我想你会看到等价性。

mean（Yobs[m1[ii，]==1]）

表示“在

ii

m1第四行中有1的位置，

Yobs

值的平均值。”.如果你不明白是什么东西把它分解了，看看每一步都在做什么。首先键入

？“%*%”

然后查看每个步骤的输出是什么

头（m1%*%Yobs）

头（排序（m1%*%Yobs））

头（排序（m1%*%Yobs）/NT）

FisherMeans <- rep(NA,nrow(m1))
for (ii in 1:nrow(m1)) {
    FisherMeans[ii] <- mean(Yobs[m1[ii,]==1]) - mean(Yobs[m1[ii,]==-1])
}