Question

好的，这让我绕弯道。

我有一个data.frame df，它从txt文件读入。在其中，我在每个样本标记的列中有一个结果列表（称为df[1]到df[15]。每个结果都按染色体（df$chr）列出。

我想绘制df[i]和df[15]之间的相关性，但是由chr组织。目前，我正在绘制图表上绘制的所有相关性，并生成每条染色体的图。然而，这些图看起来完全一样，所以我怀疑它给出了整个样本的相关性，而不是按染色体分组。你能帮忙吗？

df <- read.delim("~/Desktop/Plots/blabl.txt", header=TRUE)
ind = unique(as.character(df$chr))
for(jj in ind) {
  for(i in 3:14) {
    indic = which(df$chr == jj)
    corrie = (cor(df[,i], df[,15]))
    points(i,corrie)
   }
  plot(i, corrie, xlim=c(0,22), ylim=c(0.7,1), pch=19, main=jj)
}

Answer 1

行。我设法弄清楚这个：＃获取数据显示为df＆lt; - read.delim（＆＃34;＆＃34;〜/ Desktop / Plots / blabl.txt＆＃34;，header = TRUE）查看（df）ind = unique（as.character（df $ chr）） for（jj in ind）{indic = which（df $ chr == jj） plot（i，corrie，xlim = c（0,22），ylim = c（0.7,1），pch = 19，type =＆＃34; s＆＃34;，main = jj） for（i in 5:12）{corrie =（cor（df [，i] [indic]，df [，15] [indic]））points（i，corrie） } }

如何从循环中向点添加点

1 个答案: