我一直在使用这个psych包来比较两个相关矩阵cortest。
现在我想尝试需要类对象的cortest.matandcortest.jennrich函数phych和sim. 我已经尝试转换 mi 相关矩阵，sim.structure从而产生此类类的对象，但是在运行任一函数时都会出错。
这是我尝试使用随机数的方法：
Random<-cor(matrix(rnorm(400, 0, .25), nrow=(20), ncol=(20)))
SimRandom<-sim.structure(Random)
class(SimRandom)
cortest.jennrich(SimRandom,SimRandom,n1=400, n2=400)
产生以下结果：
if (dim(R1)[1] != p) { 中的错误：参数长度为零
我确定我做错了，因为错误消息和 Random 和 SimRandom 中的值不完全相同。
将相关矩阵转换为类型 -phych, sim- 以用作运行 cortest.mat 的输入的正确方法是哪种？
提前致谢。
编辑：关于我想要做什么的简短解释。使用随机数只是一个例子。要比较的实际相关矩阵如下完成。我有一个巨大的文件列表，每个文件都包含针对特定遗传位置的 100 个观察值。这些文件可以根据已知的遗传关系分为 20 个文件，因此我使用这些文件组，将它们作为列加载到矩阵中并计算 cor()。这给出了一个相关矩阵。作为控件，我加载随机文件并以相同的方式处理它们。该矩阵包含真实数据，但分组是随机完成的。最后，我有两个相关矩阵 1-包含预选文件的相关性和 2-包含随机加载文件之间的相关性。两个矩阵大小相同。
我想做的是比较两个相关矩阵，以了解分组是否对观察到的相关值有影响。
很抱歉没有早点解释这一点，我想避免冗长的解释并保持问题简单。