我正在尝试生成相关矩阵 (对称psd)具有预先指定的稀疏结构(由节点上的图形指定)。图中连接的节点具有相关性,其余全部为0,对角线全部为1。
我曾尝试多次生成此矩阵,但很少能获得有效的相关矩阵。
有什么方法可以确保相关矩阵whp吗?请注意,我只能具有正相关,因此等不是一个选择。
任何帮助是极大的赞赏!
也许R中的Matrix软件包的nearPD函数可以提供帮助。
—
niandra82
您为您确定的稀疏度是什么?您的数据应该是二进制数据还是非负连续数据?
—
ttnphns 2015年
@ niandra82:NearPD不好,因为它将破坏矩阵的稀疏性。
—
银翼杀手》 2015年
通常,不存在此问题中描述的矩阵分布。例如,考虑 具有三个系数的情况 。如果 和 , 然后 当且仅当矩阵是正定的。但是,那么你不能两者兼得 和 。
—
ub
那为什么不先生成相关矩阵呢?然后,为该矩阵创建一个对称索引,在其中将索引元素强制为0。稀疏度将由索引的大小指定,并且可以通过r中的sample等函数合并randommess。无论您将多少对角线元素强制设为0,该Matix仍将是pd
—
Zachary Blumenfeld