admixture可视化R语言代码

以下是Admixture结果可视化的R语言代码示例：

# 导入Admixture的结果文件
data <- read.table("result.q")

# 绘制堆叠条形图
barplot(t(data), col = rainbow(ncol(data)), border = NA, space = 0, ylab = "Ancestry", xlab = "Sample ID", main = "Admixture Results")

# 添加图例
legend("topright", c(paste0("Ancestry ", 1:ncol(data))), fill = rainbow(ncol(data)), border = NA)

解释：上述代码中，先导入Admixture结果文件，然后使用R语言的`barplot`函数绘制堆叠条形图，其中每一个样本对应一个条形，每个条形被分成不同颜色的部分，代表该样本来自不同祖先的贡献。最后使用`legend`函数添加图例，以便更好地理解图形。