[英]Plot correlation between two columns
有人可以尝试调试以下脚本吗? ggplot找不到我的x和y值。 怎么了?
数据:
data <- read.table(file="./MatureMiRNA_1M_2M_Trim18.csv", sep=",", header=T, check.names=F, row.names=1)
head(data)
1M 2M
hsa-let-7a-2-3p 0 0
hsa-let-7a-3p 0 12
hsa-let-7a-5p 2770 1344
hsa-let-7b-3p 26 27
hsa-let-7b-5p 627 492
hsa-let-7c-3p 0 0
脚本
keep <- rowSums(cpm(data)>3) >=2
data <- data[keep, ]
table(keep)
keep
FALSE TRUE
2381 200
data <- log(data,2)
head(data)
1M 2M
hsa-let-7a-5p 11.435670 10.392317
hsa-let-7b-3p 4.700440 4.754888
hsa-let-7b-5p 9.292322 8.942515
hsa-let-7c-5p 5.906891 4.754888
hsa-let-7d-3p 9.552669 7.882643
hsa-let-7d-5p 8.573647 7.851749
library(ggplot2)
ggplot(data, aes(x=1M,y=2M)) + geom_point()
Error: unexpected symbol in "ggplot(data, aes(x=1M"
R不允许以数字开头的变量名( 例如,参见here )。 但是read.table
命令中的check_names=F
阻止R检查有效名称。 R将check.names
保留为默认值(TRUE),R会通过在变量名称后添加“ X”来使名称合法。 如果要将原始名称用作ggplot中的轴标签,请使用labs
更改轴标签。
data <- read.table(text = " 1M 2M
hsa-let-7a-5p 11.435670 10.392317
hsa-let-7b-3p 4.700440 4.754888
hsa-let-7b-5p 9.292322 8.942515
hsa-let-7c-5p 5.906891 4.754888
hsa-let-7d-3p 9.552669 7.882643
hsa-let-7d-5p 8.573647 7.851749",
header = TRUE,
check.names = TRUE,
row.names = 1)
names(data)
library(ggplot2)
ggplot(data, aes(X1M, X2M)) +
geom_point() +
labs(x = "1M", y = "2M")
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.