маркировка конкретных точек данных на графике R без ggplot - PullRequest
Новая
       42

маркировка конкретных точек данных на графике R без ggplot

0 голосов
/ 06 ноября 2019

Я пытаюсь обозначить точки данных, которые заштрихованы на графике. enter image description here

Вот мои примеры данных: 

    genes     logFC       PValue
1 Arhgap8 -5.492152 2.479473e-99
2    Asns -2.519970 2.731718e-93
3    Bmp4 -1.663583 4.767201e-72
4   Casp1 -1.650139 2.212689e-25
5    Ctgf -1.272772 1.000103e-61
6    Eya4 -2.328052 2.077364e-68

мой графиккод до сих пор: 

plot(sample$logFC,-log10(as.numeric(sample$PValue)),pch = 20,xlab = 'Log2 FoldChange',ylab = '-Log10 p-value',col = 'blue',xlim = c(-10,8),ylim = c(0,300),cex.lab=1.5, cex.axis=1.5)
points(sample$logFC,-log10(as.numeric(sample$PValue)),col = "dark green")
with(subset(sample,genes=='Arhgap8'),points(logFC,-log10(as.numeric(PValue)),pch = 20, col="orange"))

Я пытался использовать приведенную ниже команду, включая текст, но она не показывает мне метку. 

with(subset(sample,genes=='Arhgap8'),points(logFC,-log10(as.numeric(PValue)),pch = 20, col="violet"),text(sample,labels = sample$genes,cex = 0.9,pos = 4))

1 Ответ

1 голос
/ 06 ноября 2019

Правильный способ использования with для выполнения команд: 

with(subset(sample, genes=='Arhgap8'), { 
  points(logFC, -log10(as.numeric(PValue)), pch = 20, col="violet")
  text(logFC, -log10(as.numeric(PValue)), labels = genes, cex = 0.9, pos = 4)
})

Когда вы передаете больше аргументов с помощью with(), они игнорируются. Например 

with(iris, mean(Sepal.Length), stop("not run"))
...