我想在 R 中制作一个我认为非常简单的图表——它将显示属于两个类别之一的基因数量。
圆圈的面积必须彼此相关,并显示我的两个类别中的计数数量之间的巨大差异。一类是 15000,另一类是 15。因此,一个圆的面积应该是另一个圆的 1000 倍。是否有一个简单的 R 脚本可用于执行此操作?(画两个圆,其中一个的面积比另一个小X倍)
plotrix
您可以使用包和draw.circle
函数绘制圆圈。所以要回答你的问题,我们只需要计算每个圆的半径。为了进行比较,使第一个圆具有单位面积更容易。所以,
## Calculate radius for given area
get_radius = function(area = 1) sqrt(area/pi)
##Load package and draw blank graph
library(plotrix)
plot(-10:10,seq(-10,10,length=21),type="n",xlab="",ylab="")
## Unit area
draw.circle(0, 0, get_radius())
## 10 times larger
draw.circle(0, 0, get_radius(10))
如本文所示,您可以使用例如shape
包并使用plotcircle
可以选择半径的功能。例子:
require("shape")
emptyplot(c(0, 1))
plotcircle(mid = c(0.2, 0.5), r = 0.1)
plotcircle(mid = c(0.6, 0.5), r = 0.01)