2016-12-28 24 views
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       

^上面已經縮短

回答

0

這些信息是由alias命令給出。這裏重複的例子:

library(vegan) 
data(dune, dune.env) 
m <- rda(dune ~ ., dune.env) 
## the following gives the names of aliased variables 
alias(m, names=TRUE) 
# [1] "Manure^4" 
## the following gives the aliasing equation 
alias(m) 
Complete : 
    A1 Moisture.L Moisture.Q Moisture.C ManagementHF ManagementNM ManagementSF Use.L Use.Q Manure.L Manure.Q Manure.C 
Manure^4            8.366600        5.291503 -4.472136 2.645751 

輸出是寬的,但它是說,當我們有ManagementNMManure一個層次是完全別名。更多細節很容易被辨別交叉表:

with(dune.env, table(Management, Manure)) 
      Manure 
Management 0 1 2 3 4 
     BF 0 2 1 0 0 
     HF 0 1 2 2 0 
     NM 6 0 0 0 0 
     SF 0 0 1 2 3 

我們只有Manure水平0當我們有Management水平NM最後級的`糞便將被別名。

+0

謝謝Jari。這對沙丘的例子有意義。但是我對大數據集有點困惑。別名(ca)輸出爲我提供了一個很長的列表,列出了所有使用數字列出的所有別名變量的變量。我如何辨別別名是什麼?示例: – shu251

+0

**我已經使用示例輸出更新了我的問題 – shu251

相關問題