testpackage包做kmeans趋势分析

a <- read.table('kmeans20220919.txt',header = T)
a <- smalltools::col_to_rowname(a,1)
group <- read.table('e:group20220919.txt',header = T)
colnames(a) <- group$Sample
a <- data.frame(a)
data <- a
data <- testpackage::means_calculate(data,group)

testpackage::kmeans_plot(
  data,
  centers = 6,
  row = 2,
  col = 3,
  center_col = "blue",
  las = 2
)

在这里插入图片描述

testpackage::kmeans_plot2(
  data,
  centers = 6,
  row = 2,
  col = 3,
  center_col = "blue",
  las = 2
)

在这里插入图片描述

testpackage::kmeans_plot3(
  data,
  centers = 6,
  row = 2,
  col = 3,
  center_col = "blue",
  las = 2
)

在这里插入图片描述

testpackage::kmeans_plot3(
  data = data,
  centers = 6,
  row = 2,
  col = 3,
  center_col = "blue",
  las = 2,
  box = FALSE,
  point_plot = TRUE

在这里插入图片描述
数据链接：https://pan.baidu.com/s/1LYArBC77M_E9GIG38HK-2Q
提取码：cg83