R/CorDdPlot.R

In HierDpart: Partitioning Hierarchical Diversity and Differentiation Across Metrics and Scales, from Genes to Ecosystems

CorDdPlot = function(x, d, ncode) {
    data = COR_DeltaDd(x, d, ncode,nrepet=999)
   # requireNamespace("ggplot2")
     a= plot(data$Dgeo, data$PairwiseDeltaD, xlab = "Geographic Distance", ylab = "Genetic differentiation (DeltaD)")
    abline(lm(data$PairwiseDeltaD ~ data$Dgeo))
    lm = lm(data$PairwiseDeltaD ~ data$Dgeo)
    return(summary(lm))
}