DotPlot 的宽高自动设置 | 线性拟合

1. 线性模型计算width和height

输入是基因集合

scRNA.markers <- FindAllMarkers(scRNA, only.pos = TRUE, min.pct = 0.25, logfc.threshold = 0.25)
#scRNA@misc[["markers"]]=scRNA.markers
#scRNA.markers %>% group_by(cluster) %>% top_n(n = 2, wt = avg_log2FC)

auto.markers=unique(c(
  scRNA.markers %>% group_by(cluster) %>% top_n(n = 4, wt = avg_log2FC) %>% pull(gene)
))

(1).pdf width

length(auto.markers)
# x, y
# 55,15
# 87,22
# 129, 30
# lm(y~x, data=data.frame(x=c(55,87,129), y=c(15,22,30)))
# y=0.2021*x+4.0787
pdf_width = 0.2021*length(auto.markers) +4.0787

如果基因集合是自定义的list，还需要考虑list中每个key之间的间距：

pdf_width = 0.2021*length(markers |> unlist() ) +4.0787 + length(markers)*0.2 #add space between gene groups