问题描述
最近我正在研究用于 PCA 分析的数据，在分析中我必须考虑 3 个组成部分。所以，我使用了 scatterplot3d 库。但是这个库的问题是我不能使组的eclipse和命名变量几乎是不可能的。
使用的数据
dput(NB)
structure(list(Serial = 1:21,Group = c("On-Shore","On-Shore","Off-Shore","On-Shore"),Serial.1 = 32:52,Station = c("NDDZ91","NDDZ92","NDDZ94","NDDZ95","NDDZ96","NDDZ97","NDDZ98","NDDZ99","NDDZ103","NDDZ105","NDDZ106","NDDZ107","NDDZ112","NDDZ113","NDDZ114","NDDZ115","NDDZ116","NDDZ117","NDDZ119","NDDZ122","NDDZ123"),Longitude = c(119.6875,119.6641667,119.7163889,119.7425,119.69,119.7719444,119.6622222,119.7544444,119.6644444,119.6875,119.7369444,119.7672222,119.6583333,119.6327778,119.5991667,119.6769444,119.7019444,119.7419444,119.6216667,119.6722222,119.6408333
),Latitude = c(26.34944444,26.35972222,26.36,26.36027778,26.37583333,26.38277778,26.38638889,26.39277778,26.41055556,26.41972222,26.42833333,26.42861111,26.43361111,26.44416667,26.45,26.44805556,26.46638889,26.47166667,26.4775),Petroleum = c(20.26723856,35.12860786,38.00639611,485.9158071,6.801285888,130.2613429,26.41691724,21.47363409,75.24255785,9.088017875,17.83869051,106.4365614,367.7372615,26.95847112,42.01753212,17.56582294,72.00256098,382.5705508,150.0004825,452.8679126,61.51137016),Sulfide = c(15.16,47.67,30.23,75.05,10.51,38.75,21.64,35.97,122.53,6.95,58.52,75.92,92.23,10.85,72.1,28.46,81.19,195.59,42.54,277.4,81.67),PAH = c(0.0479,0.0638,0.1312,0.1293,0.0589,0.1032,0.2725,0.0712,0.1173,0.0703,0.2097,0.1261,0.0685,0.0796,0.1353,0.0354,0.0323,0.034,0.0269,0.0952,0.0269),OCP = c(0.01071845,0.004247812,0.002020821,0.055799207,0.003690873,0.019063717,0.173323781,0.010429405,0.013486877,0),PES = c(0,44.15139004,17.922209,35.47228577,66.9289217,11.62692114,30.51552678,57.2039725,35.11363161,49.43776022,83.9088046,45.56461801,8.461490717,29.69151766,76.46222766,34.05399385),Acrylic = c(0,41.9544023,PP = c(0,5.263920851,186.2662335,248.4375,4.230745359,Rayon = c(74.23500995,58.86852005,53.766627,70.94457153,26.31960425,69.76152682,61.03105356,114.407945,70.22726323,136.693854,112.339267,36.33232944,16.92298143,59.38303531,305.8489106,102.1619815)),class = "data.frame",row.names = c(NA,-21L))
用于此的代码
NB = read.csv("PCA_Luoyuan.csv")
pca = prcomp(NB[,7:13],scale.=TRUE)
group = NB[,2]
zz = scatterplot3d(pca$rotation[,1:3],xlab = "PC-1",ylab = "PC-2",zlab = "PC-3",pch = 16,color = "blue",grid = TRUE)
我想要这样的图：
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）