在散点图中显示与变量有关的PCA结果

时间:2019-04-09 13:40:20

标签: r scatter-plot pca

我正在对一个大型数据集(14个变量,2500个观察值)进行主成分分析,目的是找出变量之间的相关性。我使用FactoMineR库。此刻,我仅使用以下代码行:

plot.PCA(X2.PCA,axes=c(1,3),unselect=1,choix="var" ,lim.cos2.var =0.5,  autolab="yes" )

其中X2.PCA是我的PCA分析的结果。为了检查相关性,我手动修改了axis变量中的值,这很麻烦。我想在一个简单的散点图中显示所有PCA图。我尝试了以下方法:

pairs(X2.PCA$var$cor,upper.panel=NULL)

但是,这不会产生与plot.PCA相同的图,并且我在改善它方面有些懈怠。有人可以告诉我一种更好的方式来表示我所有的图。PCA图形吗?

这是我的数据集的前20行:

  

X2 <-structure(list(CO2_Lounge = c(11.2049941666667,12.2042366666667,

     

6.76906、11.03748250025、12.6510383333333、12.6775925、13.0000608333333,   13.6306966666667、8.669745、6.55921333333333、11.5527725、11.6418625,   10.3574658333333、10.4998775、7.97028416666667、11.540345、11.4443291666667,   12.2817125,9.51236,7.06556833333333),RH_Ext = c(8.30166666666667,   8.23433333333333、8.28441666666667、8.55523333333333、7.73623333333333,   7.77083333333333,8.20058333333333、8.4955、8.54416666666667,   8.72533333333333、8.81558333333333、7.64041666666667、7.24775,   7.22408333333333,7.4675,7.96825,7.90991666666667,8.10833333333333,   7.9985,7.90133333333333),Temp_Ext = c(6.315,6.125,6.10333333333333,   6.1125、6.35916666666667、6.11083333333333、5.5025、5.0875、5.02333333333333,   5.01666666666667、5.813333333333333、11.045、11.5866666666667,   11.5366666666667、10.6383333333333、9.3825、9.6775、8.7275、8.73166666666667,   8.86333333333333),MVHR_Sup_RH = c(6.30316666666667,6.43075,   6.07666666666667、6.45391666666667、5.54908333333333、5.60091666666667,   5.27175、5.49183333333333、5.53658333333333、5.81375、6.30241666666667,   6.01791666666667、6.042、6.11508333333333、5.77391666666667,   5.87225、6.103833333333333、6.29066666666667、6.20125、6.38491666666667),MVHR_Sup_Temp = c(22.4733333333333、22.4016666666667,   21.8633333333333,   21.5675、24.3808333333333、24.0341666666667、24.2133333333333,   23.5108333333333,23.1733333333333,22.7741666666667,22.2175,   21.6883333333333、21.625、21.5441666666667、21.2758333333333,   20.9891666666667、20.945、20.9833333333333、20.8458333333333,   20.0966666666667)RH_Bed = c(7.04791666666667,7.187,7.15883333333333,   7.04916666666667、6.91641666666667、6.886725、6.86、6.83366666666667,   7.15516666666667、7.15858333333333、7.16416666666667、6.62466666666667,   6.22441666666667、6.07291666666667、5.83141666666667、6.35733333333333,   6.4435,6.71916666666667,6.965,7.17441666666667),Temp_Bed = c(22.0466666666667,   21.95、21.5733333333333、21.3016666666667、21.795、22.0641666666667,   22.0341666666667、21.8033333333333、22.0308333333333、22.175,   21.8333333333333、20.8741666666667、20.6475、20.4716666666667,   19.8658333333333、19.8841666666667、19.9516666666667、20.165,   20.2558333333333,19.8791666666667),RH_Lounge = c(5.28233333333333,   5.37475,5.00816666666667,5.45675,5.27658333333333,5.2095,   5.18416666666667,5.07325,4.79691666666667,4.689,5.12458333333333,   5.5205,5.516,5.49433333333333,5.27533333333333,5.4145,5.48558333333333,   5.69675、5.56591666666667、5.48875),Temp_Lounge = c(22.5275,   22.0525、20.985、21.4066666666667、25.015、24.8241666666667,   24.5533333333333、24.3766666666667、23.3425、22.5308333333333,   22.1216666666667、22.3766666666667、22.1741666666667、22.1816666666667,   21.445、22.1216666666667、22.025、22.9666666666667、22.3591666666667,   21.145),Delta_Elec_Mains = c(8.77、5.25、4.4、5.2、4.38、4.74,   4.85、4.91、2.94、2.78、6.22、7.44、4.59、4.6、5.55、8.09、15.65,   3.1、2.91、3.42),MVHR_Delta_elec = c(9、8、9、8、9、9、9、9、8、9、8、9、8、8、9、9、9、8、9 ),Delta_PV = c(0.01,0,0,0,0,0,0,   0,0,0,0,2.83,2.91,1.93,0.01,0,0,0,0,0),Delta_Water =   c(2.8,2.4,1,0,0,0,0,0.9,0,0,0,0.1,   0.2,0.4,0.6,1.2,3.1,0.4,0,0),Delta_Gas = c(0.5,0,0.4,0,0,0,0,0,0,0,0,0,0.06,0.24,0.3 ,0,0.5,0,0,0)),row.names =   c(13050L,13051L,13054L,13055L,13058L,13059L,13060L,13062L,   13063L,13064L,13098L,13118L,13119L,13120L,13122L,13125L,   13126L,13132L,13133L,13150L),类=“ data.frame”)

谢谢

0 个答案:

没有答案