User Tools

Site Tools


en:rda_cca_rscript

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
en:rda_cca_rscript [2019/02/10 01:17]
David Zelený
en:rda_cca_rscript [2019/02/10 01:34]
David Zelený
Line 1: Line 1:
 <code rsplus> <code rsplus>
 +# CCA - LC vs WA scores
 +
 library (vegan) library (vegan)
 +library (weimea)
 data (dune) data (dune)
-CA <- cca (dune) +data (dune.env) 
-png ('CA-scaling1-2.png', width = 8, height = 4, units = '​in',​ res = 300) +data (vltava) 
-par (mfrow = c(1,2)) + 
-plot (CAscaling ​1, main = 'CA on dune dataScaling 1') +png ('CCA-LC-WC-scores.png', width = 8, height = 8, units = '​in',​ res = 300) 
-plot (CAscaling ​2, main = 'CA on dune dataScaling ​2')+ 
 +par (mfrow = c(2,2)) 
 + 
 +CCA.dune <- cca (dune ~ Managementdata dune.env) 
 +ordiplot (CCA.dune, display = c('​sp',​ '​wa',​ '​cn'​), main = 'CCA on dune: WA scores'​xlim = c(-3, 2), ylim = c(-2, 2)) 
 +title (main = '(a)', adj = 0
 +ordiplot ​(CCA.dunedisplay ​c('​sp',​ '​lc',​ '​cn'​), main = 'CCA on dune: LC scores',​ xlim = c(-3, 2), ylim = c(-2, 2)) 
 +title (main = '(b)', adj = 0) 
 + 
 +E32 <- round (vltava$env$COVERE32/​10)*10 
 +CCA.vlt <- cca (log1p (vltava$herbs$spe) ~ E32) 
 + 
 +ordiplot (CCA.vlt, display = c('​sp',​ '​wa',​ '​bp'​),​ main = 'CCA on vltava: WA scores'​) 
 +title (main = '​(c)',​ adj = 0) 
 +ordiplot (CCA.vlt, display = c('​sp',​ '​lc',​ '​cn'​),​ main = 'CCA on vltava: WA scores'​) 
 +title (main = '​(d)',​ adj = 0)
 dev.off () dev.off ()
 </​code>​ </​code>​
en/rda_cca_rscript.txt · Last modified: 2019/02/10 01:34 by David Zelený