我一直在使用 psych 包来比较使用函数 cortest 的两个相关矩阵 .
现在我想尝试 cortest.mat 和 cortest.jennrich 函数，它需要类 phych 和 sim 的对象 . 我已经尝试使用 sim.structure 转换mi相关矩阵，这导致了这样的类的对象，但是在运行任一函数时出现错误 .
这是我尝试使用随机数：
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-包含随机加载文件之间的相关性 . 两个矩阵的大小相同 .
我想要做的是比较两个相关矩阵，以了解分组是否对观察到的相关值有影响 .
很抱歉没有解释这个，我想避免冗长的解释并保持简单的问题 .