Utilizei a seguinte linha de comando do ggplot2
em R para gerar o gráfico em anexo:
g1 <- ggplot(count, aes(x=cluster, y=Reads)) +
geom_point(shape=1) +
facet_grid(. ~ Spp)
ggplot(count, aes(x=cluster, y=Reads, fill=Spp, col=Spp)) +
geom_point(shape=1)
Esse gráfico diz respeito a quantidade de sequências de DNA em cada cluster de sequencias similares de duas espécies (Ca e Cp). Os clusters localizados na faixa de 0 reads de uma espécie são clusters enriquecidos na outra espécie. Para tentar deixar isso mais visual gostaria de plotar o gráfico da espécie Ca (os resultados que estão na cor alaranjada) na região superior e manter o da espécie Cp na região inferior. Literalmente inverter o gráfico de Ca junto, inclusive, com uma escala de Y no outro oposto, também invertida. A ideia é fazer como no esquema, mas claro, com os dados dispersos e não em uma reta.
structure(list(cluster = 1:20, Spp = structure(c(2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label =
c("Ca", "Cp"), class = "factor"), Reads = c(10808L, 16118L, 7928L,
4171L, 6525L, 2480L, 4151L, 7870L, 2508L, 2141L, 1318L, 6072L, 6459L,
10270L, 2068L, 55L, 1108L, 2116L, 4219L, 3939L)), row.names = c(NA,
20L), class = "data.frame")
Alguém sabe como posso fazer isso utilizando o pacote ggplot2
?
dput(count)
ou, se a base for muito grande, dedput(head(count, 20))
?g1
.print(g1)