0
2つのサブセットの2つの共分散行列の等価性をテストする必要があります。 列車[1:5、c(3,6,8,11)]と列車[6:10、c(3,6,8 、11)]。 boxMテストなどを使用しますか?ありがとうRの2つの共分散行列を比較するにはどうすればよいですか?
2つのサブセットの2つの共分散行列の等価性をテストする必要があります。 列車[1:5、c(3,6,8,11)]と列車[6:10、c(3,6,8 、11)]。 boxMテストなどを使用しますか?ありがとうRの2つの共分散行列を比較するにはどうすればよいですか?
covmat
からcompareCov
を使用できます。あなたがパッケージのビネットでより多くの例と詳細を見つけることができます
library(covmat)
(Cl <- cor(longley))
compareCov(Cl, Cl, labels = c("Robust Croux", "Robust"))
:
例として、私はstats
パッケージからlongely
データセットを使用します。
https://cran.r-project.org/web/packages/covmat/vignettes/CovarianceEstimation.pdf
あなたが最も適切な統計的検定を選択する方法を求めている場合はベンBolkerだけのコメントのように、それが、Cross Validatedのための質問の多くはです。
私は、彼らがCrossValidatedのためにこれをより多くの質問にする何を*統計的なテストに使うべきか質問していると思います... –
ありがとう、ありがとう、私はこれとcompareCovのboxMテストを使用できますか?もしそうなら、コーディングはどのようなものでしょうか?ありがとうございました – JSlocombe95
@ JSlocombe95私はそうは思いませんが、BoxのM-testを使用したい場合、 'heplots'を含むいくつかのパッケージがあります。これは、あなたのためにすべてのダウンボックスMを破ります:https://www.rdocumentation.org/packages/heplots/versions/1.3-1/topics/boxM 'boxM(cbind(Sepal.Length、Sepal.Width、Petal.Length、Petal。幅)〜種、データ=アイリス) '。あなたの質問に答えるには、チェックマークをクリックしてください。 –