Schefferville - transitions; Indicator value analysis

Tammy L. Elliott

Date: May 8, 2015

R version 3.1.0

In all analyzes MoeMac, SteLon, PhyCae, and ParKot dropped because of very low frequency affecting SCBD and indval analyzes

# set global chunk options: 
library(knitr)
opts_chunk$set(cache=FALSE, fig.align='center')

Correlations between SCBD for beta diversity and phylogenetic beta diversity in five environmental distance clusters

# Pearson's correlation coefficients

cor.env.dist.bd.vasc<-cor(env.dist.ord,abd.vasc.ord)
cor.env.dist.ph.bd.vasc<-cor(env.dist.ord, ph.abd.vasc.ord)
cor.bd.ph.bd.vasc<-cor(abd.vasc.ord, ph.abd.vasc.ord)

# Correlation between Gower's distance and beta diversity
cor.env.dist.bd.vasc

## [1] 0.86038

# Correlation between Gower's distance and phylogenetic beta diversity
cor.env.dist.ph.bd.vasc

## [1] 0.9147992

# Correlation between beta diversity and phylogenetic beta diversity
cor.bd.ph.bd.vasc

## [1] 0.9234183

#drop values for env.dist.ord to get angiosperms only
env.dist.ord.angio<-env.dist.ord[-match(c("LycAno", "HupApr", "DipCom", "SelSel", "CysMon", "DryExp", "EquSyl", "EquArv", "EquSci", "EquVar",
 "JunCom", "PicMar", "PicGla", "AbiBal", "LarLar"), names(env.dist.ord))]



# Pearson's correlation coefficients for angiosperms
cor.env.dist.bd.angio<-cor(env.dist.ord.angio, abd.angio.ord)
cor.env.dist.ph.bd.angio<-cor(env.dist.ord.angio, ph.abd.angio.ord)
cor.bd.ph.bd.angio<-cor(abd.angio.ord, ph.abd.angio.ord)

# Correlation between Gower's distance and beta diversity
cor.env.dist.bd.angio

## [1] 0.7289247

# Correlation between Gower's distance and phylogenetic beta diversity
cor.env.dist.ph.bd.angio

## [1] 0.8060046

# Correlation between beta diversity and phylogenetic beta diversity
cor.bd.ph.bd.angio

## [1] 0.9391508

SCBD ploted onto vascular(left) and angiosperm trees (right)

#dev.new(fig.width=11.8, fig.height=8)
par(mfrow=c(1,2))
par(mai=c(0.05,0,0.35,0))
 #Botany 2014 tree

  plot(trans.one.tree.short,x.lim=620,show.tip.label = FALSE,edge.width=2)
    nspecies <- length(trans.one.tree.short$tip.label)
    segments(rep(430,nspecies),1:nspecies,rep(430,nspecies)+(400*env.dist.ord), 1:nspecies,lwd=4, col="grey15")
    segments(rep(490,nspecies),1:nspecies,rep(490,nspecies)+(400*abd.vasc.ord), 1:nspecies,lwd=4, col="grey65")
    segments(rep(575,nspecies),1:nspecies,rep(575,nspecies)+(400*ph.abd.vasc.ord), 1:nspecies,lwd=4, col="grey30")

par(mai=c(0.35,0.2,0.95,0.05))
  plot(angio.one.tree,x.lim=330,show.tip.label = FALSE,edge.width=2)

#angiosperms
    nspecies <- length(angio.one.tree$tip.label)
    segments(rep(215,nspecies),1:nspecies,rep(215,nspecies)+(180*env.dist.ord.angio), 1:nspecies,lwd=4, col="grey15")
    segments(rep(240,nspecies),1:nspecies,rep(240,nspecies)+(180*abd.angio.ord), 1:nspecies,lwd=4, col="grey80")
    segments(rep(290,nspecies),1:nspecies,rep(290,nspecies)+(180*ph.abd.angio.ord), 1:nspecies,lwd=4, col="grey50")

#textClick.bold("(A)", cex=1.5)
#textClick.bold("(B)", cex=1.5)
#textClick.bold("(C)", cex=1.5)
#textClick.bold("(D)", cex=1.5)
#textClick.bold("(E)", cex=1.5)
#textClick.bold("(F)", cex=1.5)

Correlations with species frequencies across plots

#Calculate number of occurrences per species
abd.sp.allsp.pa<-abd.sp.allsp
abd.sp.allsp.pa[abd.sp.allsp.pa > 0] <- 1 
freq.sp<-(rowSums(abd.sp.allsp.pa))

angio.sp.abd.pa<-angio.sp.abd
angio.sp.abd.pa[angio.sp.abd.pa > 0] <- 1 
angio.freq.sp<-(colSums(angio.sp.abd.pa))

#Correlations with SCBD
# Environmental distance and frequency
(env.freq.corr<-cor(freq.sp, env.dist.ord))

## [1] -0.1546929

# Vascular abundance with SCBD
(vasc.freq.corr<-cor(freq.sp,abd.vasc.ord))

## [1] -0.07957138

# Vascular phylogenetic abundance with SCBD
(vasc.ph.freq.corr<-cor(freq.sp,ph.abd.vasc.ord))

## [1] -0.1233261

# Angiosperms with SCBD
(angio.freq.corr<-cor(angio.freq.sp,abd.angio.ord))

## [1] 0.07933482

# Angiosperms with angiosperm SCBD (phylogenetic)
(angio.ph.freq.corr<-cor(angio.freq.sp,ph.abd.angio.ord))

## [1] 0.05562471

Indicator species analysis - correlations of species with different habitats

Order for the following plots is top - vasculars; bottom -angiosperms

1 - environmental distance; 2 - vasc.BD; 3 - vasc.phBD; 4 - angioBD; 5 -angioPh.BD

#dev.new(width=11.8, height=11.8)
par(mfrow=c(2,3))
par(mai=c(0.05,0,0.25,0))
  plot(vasc.tree.ind, x.lim=800,show.tip.label = FALSE,edge.width=2)
    nspecies <- length(vasc.tree.ind$tip.label)
    segments(rep(450,nspecies),1:nspecies,rep(450,nspecies)+(100*env.dist.ord.ind.1), 1:nspecies,lwd=3, col="black")
     segments(rep(520,nspecies),1:nspecies,rep(520,nspecies)+(100*env.dist.ord.ind.2), 1:nspecies,lwd=3, col="black")
     segments(rep(600,nspecies),1:nspecies,rep(600,nspecies)+(100*env.dist.ord.ind.3), 1:nspecies,lwd=3, col="black")
     segments(rep(680,nspecies),1:nspecies,rep(680,nspecies)+(100*env.dist.ord.ind.4), 1:nspecies,lwd=3, col="black")
     segments(rep(760,nspecies),1:nspecies,rep(760,nspecies)+(100*env.dist.ord.ind.5), 1:nspecies,lwd=3, col="black")

par(mai=c(0.05,0,0.25,0))
  plot(vasc.tree.ind, x.lim=800,show.tip.label = FALSE,edge.width=2)
    nspecies <- length(vasc.tree.ind$tip.label)
    segments(rep(450,nspecies),1:nspecies,rep(450,nspecies)+(100*vasc.abd1.ord.ind.1), 1:nspecies,lwd=3, col="black")
     segments(rep(540,nspecies),1:nspecies,rep(540,nspecies)+(100*vasc.abd1.ord.ind.2), 1:nspecies,lwd=3, col="black")
     segments(rep(600,nspecies),1:nspecies,rep(600,nspecies)+(100*vasc.abd1.ord.ind.3), 1:nspecies,lwd=3, col="black")
     segments(rep(670,nspecies),1:nspecies,rep(670,nspecies)+(100*vasc.abd1.ord.ind.4), 1:nspecies,lwd=3, col="black")
     segments(rep(740,nspecies),1:nspecies,rep(740,nspecies)+(100*vasc.abd1.ord.ind.5), 1:nspecies,lwd=3, col="black")

 par(mai=c(0.05,0,0.25,0))
  plot(vasc.tree.ind, x.lim=800,show.tip.label = FALSE,edge.width=2)
    nspecies <- length(vasc.tree.ind$tip.label)
    segments(rep(450,nspecies),1:nspecies,rep(450,nspecies)+(100*vasc.ph.abd1.ord.ind.1), 1:nspecies,lwd=3, col="black")
     segments(rep(540,nspecies),1:nspecies,rep(540,nspecies)+(100*vasc.ph.abd1.ord.ind.2), 1:nspecies,lwd=3, col="black")
     segments(rep(600,nspecies),1:nspecies,rep(600,nspecies)+(100*vasc.ph.abd1.ord.ind.3), 1:nspecies,lwd=3, col="black")
     segments(rep(670,nspecies),1:nspecies,rep(670,nspecies)+(100*vasc.ph.abd1.ord.ind.4), 1:nspecies,lwd=3, col="black")
     segments(rep(740,nspecies),1:nspecies,rep(740,nspecies)+(100*vasc.ph.abd1.ord.ind.5), 1:nspecies,lwd=3, col="black")


par(mai=c(0.05,0,0.25,0))
  plot(angio.tree.ind, x.lim=450,show.tip.label = FALSE,edge.width=2)
    nspecies <- length(angio.tree.ind$tip.label)
    segments(rep(220,nspecies),1:nspecies,rep(220,nspecies)+(40*angio.abd1.ord.ind.1), 1:nspecies,lwd=3, col="black")
     segments(rep(265,nspecies),1:nspecies,rep(265,nspecies)+(40*angio.abd1.ord.ind.2), 1:nspecies,lwd=3, col="black")
     segments(rep(300,nspecies),1:nspecies,rep(300,nspecies)+(40*angio.abd1.ord.ind.3), 1:nspecies,lwd=3, col="black")
     segments(rep(350,nspecies),1:nspecies,rep(350,nspecies)+(40*angio.abd1.ord.ind.4), 1:nspecies,lwd=3, col="black")
     #segments(rep(740,nspecies),1:nspecies,rep(740,nspecies)+(100*angio.abd1.ord.ind.5), 1:nspecies,lwd=3, col="black")

 par(mai=c(0.05,0,0.25,0))
  plot(angio.tree.ind, x.lim=450,show.tip.label = FALSE,edge.width=2)
    nspecies <- length(angio.tree.ind$tip.label)
       segments(rep(220,nspecies),1:nspecies,rep(220,nspecies)+(40*angio.ph.abd1.ord.ind.1), 1:nspecies,lwd=3, col="black")
     segments(rep(265,nspecies),1:nspecies,rep(265,nspecies)+(40*angio.ph.abd1.ord.ind.2), 1:nspecies,lwd=3, col="black")
     segments(rep(300,nspecies),1:nspecies,rep(300,nspecies)+(40*angio.ph.abd1.ord.ind.3), 1:nspecies,lwd=3, col="black")
     segments(rep(350,nspecies),1:nspecies,rep(350,nspecies)+(40*angio.ph.abd1.ord.ind.4), 1:nspecies,lwd=3, col="black")
     segments(rep(380,nspecies),1:nspecies,rep(380,nspecies)+(40*angio.ph.abd1.ord.ind.5), 1:nspecies,lwd=3, col="black")

Phylogenetic signal for each five cluster grouping

#Phylogenetic signal in environmental distance - vaculars
(env.dist.1.k<-Kcalc(env.dist.ord.ind.1,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01930172

(env.dist.2.k<-Kcalc(env.dist.ord.ind.2,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.02706511

(env.dist.3.k<-Kcalc(env.dist.ord.ind.3,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01605215

(env.dist.4.k<-Kcalc(env.dist.ord.ind.4,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01943013

(env.dist.5.k<-Kcalc(env.dist.ord.ind.5,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01897711

#Phylogenetic signal in vascular beta diversity
(vasc.abd.1.k<-Kcalc(vasc.abd1.ord.ind.1,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01361344

(vasc.abd.2.k<-Kcalc(vasc.abd1.ord.ind.2,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01764653

(vasc.abd.3.k<-Kcalc(vasc.abd1.ord.ind.3,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01750493

(vasc.abd.4.k<-Kcalc(vasc.abd1.ord.ind.4,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01312226

(vasc.abd.5.k<-Kcalc(vasc.abd1.ord.ind.5,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.02572936

#Phylogenetic signal in vascular phylogenetic beta diversity clusters
(vasc.ph.abd.1.k<-Kcalc(vasc.ph.abd1.ord.ind.1,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.03437222

(vasc.ph.abd.2.k<-Kcalc(vasc.ph.abd1.ord.ind.2,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01399836

(vasc.ph.abd.3.k<-Kcalc(vasc.ph.abd1.ord.ind.3,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01530427

(vasc.ph.abd.4.k<-Kcalc(vasc.ph.abd1.ord.ind.4,vasc.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.02401423

(vasc.ph.abd.5.k<-Kcalc(vasc.ph.abd1.ord.ind.5,vasc.tree.ind,checkdata=TRUE))

##           [,1]
## [1,] 0.0218506

#Phylogenetic signal in angiosperm beta diversity clusters
(angio.abd.1.k<-Kcalc(angio.abd1.ord.ind.1,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.02604287

(angio.abd.2.k<-Kcalc(angio.abd1.ord.ind.2,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01529402

(angio.abd.3.k<-Kcalc(angio.abd1.ord.ind.3,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.00924189

(angio.abd.4.k<-Kcalc(angio.abd1.ord.ind.4,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.04440957

(angio.abd.5.k<-Kcalc(angio.abd1.ord.ind.5,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01339638

#Phylogenetic signal in angiosperm phylogenetic beta diversity clusters
(angio.ph.abd.1.k<-Kcalc(angio.ph.abd1.ord.ind.1,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.01689364

(angio.ph.abd.2.k<-Kcalc(angio.ph.abd1.ord.ind.2,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.03779782

(angio.ph.abd.3.k<-Kcalc(angio.ph.abd1.ord.ind.3,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.05576866

(angio.ph.abd.4.k<-Kcalc(angio.ph.abd1.ord.ind.4,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.02380039

(angio.ph.abd.5.k<-Kcalc(angio.ph.abd1.ord.ind.5,angio.tree.ind,checkdata=TRUE))

##            [,1]
## [1,] 0.02522516

Pearson’s correlation coefficients within and among same number clusters

# Correlations between environmental distance and vascular abundance clusters
(env.vasc.cor.1<-cor(env.dist.ord.ind.1, vasc.abd1.ord.ind.1))

## [1] 0.2104629

(env.vasc.cor.2<-cor(env.dist.ord.ind.2, vasc.abd1.ord.ind.2))

## [1] -0.3552442

(env.vasc.cor.3<-cor(env.dist.ord.ind.3, vasc.abd1.ord.ind.3))

## [1] 0.2711327

(env.vasc.cor.4<-cor(env.dist.ord.ind.4, vasc.abd1.ord.ind.4))

## [1] -0.1466499

(env.vasc.cor.5<-cor(env.dist.ord.ind.5, vasc.abd1.ord.ind.5))

## [1] 0.04541109

# Cross Correlations between environmental distance and vascular abundance clusters
(env.vasc.cor.1.2<-cor(env.dist.ord.ind.1, vasc.abd1.ord.ind.2))

## [1] -0.1603619

(env.vasc.cor.1.3<-cor(env.dist.ord.ind.1, vasc.abd1.ord.ind.3))

## [1] 0.1897604

(env.vasc.cor.1.4<-cor(env.dist.ord.ind.1, vasc.abd1.ord.ind.4))

## [1] 0.2369927

(env.vasc.cor.1.5<-cor(env.dist.ord.ind.1, vasc.abd1.ord.ind.5))

## [1] -0.4691192

(env.vasc.cor.2.3<-cor(env.dist.ord.ind.2, vasc.abd1.ord.ind.3))

## [1] -0.05167944

(env.vasc.cor.2.4<-cor(env.dist.ord.ind.2, vasc.abd1.ord.ind.4))

## [1] 0.07377345

(env.vasc.cor.2.5<-cor(env.dist.ord.ind.2, vasc.abd1.ord.ind.5))

## [1] -0.3578274

(env.vasc.cor.3.4<-cor(env.dist.ord.ind.3, vasc.abd1.ord.ind.4))

## [1] 0.1284278

(env.vasc.cor.4.5<-cor(env.dist.ord.ind.4, vasc.abd1.ord.ind.5))

## [1] 0.8351188

# Correlations between environmental distance and vascular phylogenetic abundance clusters
(env.vasc.ph.cor.1<-cor(env.dist.ord.ind.1, vasc.ph.abd1.ord.ind.1))

## [1] 0.05886758

(env.vasc.ph.cor.2<-cor(env.dist.ord.ind.2, vasc.ph.abd1.ord.ind.2))

## [1] -0.3626696

(env.vasc.ph.cor.3<-cor(env.dist.ord.ind.3, vasc.ph.abd1.ord.ind.3))

## [1] -0.3022219

(env.vasc.ph.cor.4<-cor(env.dist.ord.ind.4, vasc.ph.abd1.ord.ind.4))

## [1] 0.2214922

(env.vasc.ph.cor.5<-cor(env.dist.ord.ind.5, vasc.ph.abd1.ord.ind.5))

## [1] -0.02335573

# Cross correlations between environmental distance and vascular phylogenetic abundance clusters

(env.vasc.ph.cor.1.2<-cor(env.dist.ord.ind.1, vasc.ph.abd1.ord.ind.2))

## [1] 0.1040027

(env.vasc.ph.cor.1.3<-cor(env.dist.ord.ind.1, vasc.ph.abd1.ord.ind.3))

## [1] 0.784136

(env.vasc.ph.cor.1.4<-cor(env.dist.ord.ind.1, vasc.ph.abd1.ord.ind.4))

## [1] -0.5035221

(env.vasc.ph.cor.1.5<-cor(env.dist.ord.ind.1, vasc.ph.abd1.ord.ind.5))

## [1] -0.3882596

(env.vasc.ph.cor.2.3<-cor(env.dist.ord.ind.2, vasc.ph.abd1.ord.ind.3))

## [1] 0.2426948

(env.vasc.ph.cor.2.4<-cor(env.dist.ord.ind.2, vasc.ph.abd1.ord.ind.4))

## [1] -0.1292322

(env.vasc.ph.cor.2.5<-cor(env.dist.ord.ind.2, vasc.ph.abd1.ord.ind.5))

## [1] -0.3738739

(env.vasc.ph.cor.3.4<-cor(env.dist.ord.ind.3, vasc.ph.abd1.ord.ind.4))

## [1] 0.5450795

(env.vasc.ph.cor.4.5<-cor(env.dist.ord.ind.4, vasc.ph.abd1.ord.ind.5))

## [1] 0.7770791

# Correlations between vascular abundance clusters and vascular phylogenetic abundance clusters
(vasc.vasc.ph.cor.1<-cor(vasc.abd1.ord.ind.1, vasc.ph.abd1.ord.ind.1))

## [1] 0.7571733

(vasc.vasc.ph.cor.2<-cor(vasc.abd1.ord.ind.2, vasc.ph.abd1.ord.ind.2))

## [1] 0.7783666

(vasc.vasc.ph.cor.3<-cor(vasc.abd1.ord.ind.3, vasc.ph.abd1.ord.ind.3))

## [1] 0.1817989

(vasc.vasc.ph.cor.4<-cor(vasc.abd1.ord.ind.4, vasc.ph.abd1.ord.ind.4))

## [1] 0.2115801

(vasc.vasc.ph.cor.5<-cor(vasc.abd1.ord.ind.5, vasc.ph.abd1.ord.ind.5))

## [1] 0.8169962

# Cross correlatons between vascular abundance clusters and vascular phylogenetic abundance clusters
(vasc.vasc.ph.cor.1.2<-cor(vasc.abd1.ord.ind.1, vasc.ph.abd1.ord.ind.2))

## [1] -0.3831161

(vasc.vasc.ph.cor.1.3<-cor(vasc.abd1.ord.ind.1, vasc.ph.abd1.ord.ind.3))

## [1] 0.2693246

(vasc.vasc.ph.cor.1.4<-cor(vasc.abd1.ord.ind.1, vasc.ph.abd1.ord.ind.4))

## [1] -0.08701183

(vasc.vasc.ph.cor.1.5<-cor(vasc.abd1.ord.ind.1, vasc.ph.abd1.ord.ind.5))

## [1] -0.4214

(vasc.vasc.ph.cor.2.3<-cor(vasc.abd1.ord.ind.2, vasc.ph.abd1.ord.ind.3))

## [1] -0.317559

(vasc.vasc.ph.cor.2.4<-cor(vasc.abd1.ord.ind.2, vasc.ph.abd1.ord.ind.4))

## [1] -0.2654193

(vasc.vasc.ph.cor.2.5<-cor(vasc.abd1.ord.ind.2, vasc.ph.abd1.ord.ind.5))

## [1] -0.02401277

(vasc.vasc.ph.cor.3.4<-cor(vasc.abd1.ord.ind.3, vasc.ph.abd1.ord.ind.4))

## [1] -0.02184757

(vasc.vasc.ph.cor.4.5<-cor(vasc.abd1.ord.ind.4, vasc.ph.abd1.ord.ind.5))

## [1] -0.2035186

# Correlations between angiosperm abundance clusters and angiosperm phylogenetic abundance clusters
(angio.angio.ph.cor.1<-cor(angio.abd1.ord.ind.1, angio.ph.abd1.ord.ind.1))

## [1] -0.2945521

(angio.angio.ph.cor.2<-cor(angio.abd1.ord.ind.2, angio.ph.abd1.ord.ind.2))

## [1] -0.2901857

(angio.angio.ph.cor.3<-cor(angio.abd1.ord.ind.3, angio.ph.abd1.ord.ind.3))

## [1] -0.5133346

(angio.angio.ph.cor.4<-cor(angio.abd1.ord.ind.4, angio.ph.abd1.ord.ind.4))

## [1] 0.4199509

(angio.angio.ph.cor.5<-cor(angio.abd1.ord.ind.5, angio.ph.abd1.ord.ind.5))

## [1] -0.392256

# Cross correlatons between angiosperm abundance clusters and angiosperm phylogenetic abundance clusters
(angio.angio.ph.cor.1.2<-cor(angio.abd1.ord.ind.1, angio.ph.abd1.ord.ind.2))

## [1] 0.4172365

(angio.angio.ph.cor.1.3<-cor(angio.abd1.ord.ind.1, angio.ph.abd1.ord.ind.3))

## [1] 0.3794714

(angio.angio.ph.cor.1.4<-cor(angio.abd1.ord.ind.1, angio.ph.abd1.ord.ind.4))

## [1] 0.06506476

(angio.angio.ph.cor.1.5<-cor(angio.abd1.ord.ind.1, angio.ph.abd1.ord.ind.5))

## [1] -0.4092137

(angio.angio.ph.cor.2.3<-cor(angio.abd1.ord.ind.2, angio.ph.abd1.ord.ind.3))

## [1] 0.3784006

(angio.angio.ph.cor.2.4<-cor(angio.abd1.ord.ind.2, angio.ph.abd1.ord.ind.4))

## [1] 0.04954705

(angio.angio.ph.cor.2.5<-cor(angio.abd1.ord.ind.2, angio.ph.abd1.ord.ind.5))

## [1] -0.01725278

(angio.angio.ph.cor.3.4<-cor(angio.abd1.ord.ind.3, angio.ph.abd1.ord.ind.4))

## [1] -0.4757189

(angio.angio.ph.cor.4.5<-cor(angio.abd1.ord.ind.4, angio.ph.abd1.ord.ind.5))

## [1] 0.7247129