基于另一个数据帧的数据帧子集,在R中具有相同的行名称和列名
我使用Hmisc的rcorr()进行了多次关联(在大约500 x 40的大数据帧上)。此命令返回一个列表,但我只将列表中需要的内容保存为数据帧。因此,我有两个数据帧,一个是Pearson rho,另一个是p值。数据帧之间的所有内容都相同,但单元格值不同基于另一个数据帧的数据帧子集,在R中具有相同的行名称和列名,r,list,dataframe,subset,R,List,Dataframe,Subset,我使用Hmisc的rcorr()进行了多次关联(在大约500 x 40的大数据帧上)。此命令返回一个列表,但我只将列表中需要的内容保存为数据帧。因此,我有两个数据帧,一个是Pearson rho,另一个是p值。数据帧之间的所有内容都相同,但单元格值不同 AdiASigFatCorr <- rcorr(as.matrix(AdiASigFat), type=c("pearson")) AdiASigFatCorrP <- as.data.frame(AdiASigFatCorr$P)
AdiASigFatCorr <- rcorr(as.matrix(AdiASigFat), type=c("pearson"))
AdiASigFatCorrP <- as.data.frame(AdiASigFatCorr$P)
AdiASigFatCorrR <- as.data.frame(AdiASigFatCorr$r)
ADIASIGFATCORR(rho值):
如果相应的p值<0.05,我想将rho值子集
我读了很多书,尝试了很多不同的方法,比如使用subset()等,但我完全被卡住了
如果有其他方法/包用于此目的,请告知我
谢谢大家! 调用rho矩阵
corr.rho
和p值矩阵corr.p
,您可以使用简单的向量运算来查找corr.rho
的所有值,从而corr.p<0.05
。如果数据尚未采用矩阵形式,则可以将其强制转换为具有as.matrix
> corr.rho[corr.p < 0.05]
[1] 0.6838388 0.4349413 0.6459200 -0.4148415
>corr.rho[corr.p<0.05]
[1] 0.6838388 0.4349413 0.6459200 -0.4148415
pre_ISI.kg_ pre_ISI.mmol pre_Weight
pre_BMP3 0.6838388 0.6459200 -0.4148415
pre_CES1 -0.3887246 -0.3653443 0.3685230
pre_CETP 0.4349413 0.3837841 -0.1948498
> corr.rho[corr.p < 0.05]
[1] 0.6838388 0.4349413 0.6459200 -0.4148415