与群体分散的情节

时间:2018-02-16 02:08:54

标签: r ggplot2 grouping scatter-plot

您好我想做一个不同组的散点图或点图。 我可以让它适用于单个图,但因为ggplot2需要x轴的坐标。多数民众赞成我已经陷入困境。 这是我的结构:

# A tibble: 2 x 33
  gene_id     gene     N1    N2    N3    N4    N5    N6    N7    T1    T2    T3    T4
  <chr>       <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
1 ENSMUSG000… RNS     198   182   206   183   177   194   193   173   191   167   200
2 ENSMUSG000… RNS2    199   198   216   252   273   159   164   159   162   151   199
# ... with 20 more variables: T5 <dbl>, T6 <dbl>, T7 <dbl>, T8 <dbl>, T9 <dbl>,
#   T10 <dbl>, T11 <dbl>, T12 <dbl>, T13 <dbl>, T14 <dbl>, M1 <dbl>, M2 <dbl>,
#   M3 <dbl>, M4 <dbl>, M5 <dbl>, M6 <dbl>, M7 <dbl>, M8 <dbl>, M9 <dbl>, M10 <dbl>

我想首先绘制彼此相邻的所有个体,并将N,T和Ms组合在一起,并用不同的颜色绘制它们。我很确定一旦绘图本身有效,我就可以分配颜色。

ggplot(a, aes(y=a[,3],x=1))+
geom_point()

在这种情况下,我开始在porcess中非常低。;;;

感谢您的帮助!

1 个答案:

答案 0 :(得分:0)

您需要将列名称转换为变量。我使用您发布的数据,但我不确定您是否期望这样的情节

data <- read.csv('tabla.csv')
data

初始数据

 gene_id gene  N1  N2  N3  N4  N5  N6  N7  T1  T2  T3  T4
1 ENSMUSG000…  RNS 198 182 206 183 177 194 193 173 191 167 200
2 ENSMUSG000… RNS2 199 198 216 252 273 159 164 159 162 151 199

融化数据

library(reshape)

newData <- melt(data, id=c("gene_id", "gene"))
newData

你得到一个像这样的数据框

gene_id gene variable value
1  ENSMUSG000…  RNS       N1   198
2  ENSMUSG000… RNS2       N1   199
3  ENSMUSG000…  RNS       N2   182
4  ENSMUSG000… RNS2       N2   198
5  ENSMUSG000…  RNS       N3   206

使用每个组的名称(N,T或M)创建一个新的列

count = 0
for(value in newData$variable){

        # To get the group of each value

        if(startsWith(value, 'N')){

                group <- 'N'

        }else if(startsWith(value, 'T')){

                group <- 'T'

        }else if(startsWith(value, 'M')){

                group <- 'M'
        }

        # Storing the group name into a vector
        if(count == 0){

                vectorGroup <- group
        }else{

                vectorGroup <- c(vectorGroup, group)

        }

        count <- count + 1

}

newData$group <- vectorGroup

最后制作散点图

library(ggplot2)

ggplot(newData, aes(x = gene, y = value, color= group)) + theme_minimal() + geom_point(width = 0.5) + scale_fill_brewer(palette = 'Paired') 

你得到这样的情节 Gene by Group