contoh

`

# Load library
library(ggplot2)
library(reshape2)

# Membuat matriks korelasi
korelasi_matrix <- matrix(c( NA,    0.001,  0.001, 0.001, 0.018,
                             0.001,  NA,    0.001, 0.001, 0.009,
                             0.001, 0.001, NA,  0.001, 0.008,
                             0.001, 0.001, 0.001, NA, 0.004,
                             0.018, 0.009, 0.008, 0.004, NA), 
                          nrow = 5, byrow = TRUE)

# Menamai baris dan kolom
rownames(korelasi_matrix) <- colnames(korelasi_matrix) <- c("Cluster 1", "Cluster 2", "Cluster 3",
                                                            "Cluster 4", "Cluster 5")

# Konversi ke data frame untuk ggplot
data <- melt(korelasi_matrix)
colnames(data) <- c("Cluster1", "Cluster2", "Korelasi")

# Menentukan warna berdasarkan signifikansi (asumsi p-value tidak tersedia, hanya nilai korelasi)
data$Warna <- ifelse(data$Korelasi > 0.05|data$Korelasi=="NA", "white", "#99d8c9")

# Plot heatmap
ggplot(data, aes(x = Cluster1, y = Cluster2, fill = Warna)) +
  geom_tile(color = "black") +
  geom_text(aes(label = round(Korelasi, 3)), size = 6) +  # Menampilkan nilai korelasi
  scale_fill_manual(values = c("#99d8c9", "white")) +
  labs(title = "Heatmap p-value Cluster", x = "", y = "") +
  theme_minimal() +
  theme(
    axis.text.x = element_text(size=15),
    axis.text.y = element_text(size=15),
    legend.position = "right"  # Menampilkan legenda di sebelah kanan
    )

## Warning: Removed 5 rows containing missing values or values outside the scale range
## (`geom_text()`).