r - R の Marascuilo プロシージャ

Question

プロポーションの違いを比較するためにマラスキーリオ法を行っています。私は次のコードを使用しています（このチュートリアルからコピーして適応させました：

## Set the proportions of interest.
p = c(0.3481, 0.1730, 0.4788)
N = length(p)
value = critical.range = c()

## Compute critical values.
for (i in 1:(N-1))
{ for (j in (i+1):N)
{
  value = c(value,(abs(p[i]-p[j])))
  critical.range = c(critical.range,
                     sqrt(qchisq(.95,3))*sqrt(p[i]*(1-p[i])/12000 + p[j]*(1-p[j])/12000))
}
}
round(cbind(value,critical.range),3)

出力にカテゴリのラベルも出力する必要があります (たとえば、どのカテゴリが正確に比較されているか)。

したがって、カテゴリが分離されたベクトルでリストされている場合、たとえばcategories <- c("cat1", "cat2", cat"3)、比較はcat1-cat2、cat1-cat3、およびcat2-cat3です。

これらのラベルを出力に追加するにはどうすればよいですか?

    value critical.range
[1,] 0.175          0.016
[2,] 0.131          0.018
[3,] 0.306          0.016