0
我已經執行了一個使用素食主義者的CCA並繪製了雙指箭頭。 「ca」輸出告訴我我有「一些約束是別名,因爲它們是共線的(冗餘)」。我知道這意味着這些限制因爲多餘而被放棄。R素食主義者CCA,確定哪些限制是別名?
ca<-cca(CCA_data[3:4567]~.,data=Meta_data_2)
ca
我還可以列出別名沒有問題:
alias(ca, names=TRUE,complete=TRUE)
但是,我怎麼確定哪些約束條件別名或爲共線與其他因素。我需要知道哪些因素與什麼因素是多餘的,所以我可以重新思考我的方法。
這又給了我列出的別名,但沒有雙標點值。而在這裏第二行給我雙標圖值,但別名-ED的人失蹤...
ca$CCA$alias
ca$CCA$biplot
如何得到這些信息的建議「CA」將是有益的。
**更新:
Model :
CCA_data[3:4567] ~ `Assimilatory nitrate reduction, nitrate => ammonia` +
`beta-Oxidation` + `beta-Oxidation, acyl-CoA synthesis` +
`Citrate cycle (TCA cycle, Krebs cycle)` + `Dissimilatory nitrate reduction, nitrate => ammonia` +
`Fatty acid biosynthesis, elongation` + `Fatty acid biosynthesis, elongation, endoplasmic reticulum` +
`Fatty acid biosynthesis, elongation, mitochondria` + `Fatty acid biosynthesis, initiation` +
`Gluconeogenesis, oxaloacetate => fructose-6P` + `Glycolysis (Embden-Meyerhof pathway), glucose => pyruvate` +
`Glyoxylate cycle` + `Nitrate assimilation` + Photorespiration +
`Photosystem I` + `Photosystem II` + `Pyruvate oxidation, pyruvate => acetyl-CoA` +
`Reductive citrate cycle (Arnon-Buchanan cycle)` + `Reductive pentose phosphate cycle (Calvin cycle)` +
`Urea cycle`
Complete :
`Assimilatory nitrate reduction, nitrate => ammonia`
`Photosystem I` 1.23566
`Photosystem II` -0.64830
`Pyruvate oxidation, pyruvate => acetyl-CoA` -0.40986
`Reductive citrate cycle (Arnon-Buchanan cycle)` 22.01804
`Reductive pentose phosphate cycle (Calvin cycle)` 16.10217
`Urea cycle` 1.55773
`beta-Oxidation`
`Photosystem I` 0.24919
`Photosystem II` 1.20864
`Pyruvate oxidation, pyruvate => acetyl-CoA` -0.24875
`Reductive citrate cycle (Arnon-Buchanan cycle)` -1.84576
`Reductive pentose phosphate cycle (Calvin cycle)` -0.74552
`Urea cycle` 5.03302
`beta-Oxidation, acyl-CoA synthesis`
`Photosystem I` 0.08856
`Photosystem II` 1.19711
`Pyruvate oxidation, pyruvate => acetyl-CoA` 0.39169
`Reductive citrate cycle (Arnon-Buchanan cycle)` 0.23713
`Reductive pentose phosphate cycle (Calvin cycle)` -6.34640
`Urea cycle` -1.88670
^上面已經縮短
謝謝Jari。這對沙丘的例子有意義。但是我對大數據集有點困惑。別名(ca)輸出爲我提供了一個很長的列表,列出了所有使用數字列出的所有別名變量的變量。我如何辨別別名是什麼?示例: – shu251
**我已經使用示例輸出更新了我的問題 – shu251