繁体   English   中英

两列之间的图相关

[英]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.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM