使用R繪製三元相圖(Ternary plot)結合等高線圖(topographic map)

## Version 1.2.1335
## coding: utf-8 
## 繪製三元相圖(Ternary plot)結合等高線圖(topographic map)
# 安裝相關套件
install.packages("AlgDesign")
install.packages('ggtern')
install.packages('ggplot2')
install.packages('RColorBrewer')
#載入相關套件
library(AlgDesign)
library(ggplot2)
library(RColorBrewer)
library(ggtern)

data_ggtern <- read.csv(file.choose()) #以視窗方式取得檔案位址
head(data_ggtern) #查看數據前6筆數據
## 有時候數據抬頭為中文名稱時讀取後會是亂碼
##可以利用colnames()重新命名
colnames(data_ggtern) <- c("A", "B", "C", "D")

p <- ggtern(data = data_ggtern,
         aes(
           x = A,
           y = B,
           z = C,
           value = data_ggtern$D
         )) + theme_bw() +
  geom_interpolate_tern(
    data = data_ggtern,
    mapping = aes(value = data_ggtern$D, color = ..level..),
    method = lm,
    #等高線趨勢資料計算
    formula = expand.formula(value ~ cubicS(x, y) + quad(x, y)),
    base = "identity",
  ) +
  #添加資料點
  geom_point(size = 2, aes(color = data_ggtern$D)) +
  #設定顏色漸層
  scale_color_gradient(low = 'green', high = 'red') +
  labs(fill = "D")
p