0

我正在运行一个 R 代码(不是我的),旨在读取目录并绘制图表:

library(ggplot2)

dat <- at <- read.table(file="User/alexanderzamani/Documents/results.txt", 
                        sep = "\t", header= TRUE)

#let's just look at a single mutation rate
r <- dat[which(dat$u==0.00010),]

#now let's plot it out
p <- ggplot(r, aes(x=S, y=k, group=pN))
p + geom_line(aes(color=pN), size=2, alpha=0.5) + geom_point() +
    xlab("Overlapping generations (proportion of offspring that survive)") + 
    ylab("substitution rate") + 
    scale_colour_gradient(name="Degree of\npopulationsize\nfluctuation")

安装后ggplot2,我运行了代码,它显示了一个空白的石英窗口和这个错误:

Error in if (nrow(layer_data) == 0) return() : argument is of length zero

我以前没有使用 R 的经验,任何帮助将不胜感激。

- 编辑

这是我正在使用的上传的 .txt 文件的链接 - http://txtup.co/NW6lU

4

1 回答 1

6

当您执行以下操作时,您正在dat创建rwhere 列u等于 0.0001 的子集:

r <- dat[which(dat$u==0.0001),]

您是否可能没有任何值为 0.0001 ....?

我做了0.1并得到:

在此处输入图像描述

于 2012-12-03T23:04:12.727 回答