看看 $FAI 文件里面的 染色体ID和你的 all.chr.xpclr 染色体ID是否一致;
还有就是 19226500 这个值表示小于这个长度的染色体会被过滤掉;
xpclr_smooth_line_plot.r -i all.chr.xpclr -F $FAI -f 19226500 -n all.chr.xpclr.smoothline --> Q&A for bioinformatics, please visit the website: https://www.omicsclass.com/ --> R beginners ? I suggest your learning R language: https://study.omicsclass.com/index Loading required package: ggplot2 chr_name = 1 chr_name = 8 chr_name = 15 chr_name = 2 chr_name = 9 chr_name = 16 chr_name = 3 chr_name = 10 chr_name = 17 chr_name = 4 chr_name = 11 chr_name = 18 chr_name = 5 chr_name = 12 chr_name = 19 chr_name = 6 chr_name = 13 chr_name = 20 chr_name = 7 chr_name = 14 chr_name = 21 chr_name = 0 1% 5% 80% 90% 95% 99% -0.2429295 -0.2373629 -0.1200045 -0.1187288 0.4832382 4.2705026 Warning message: Removed 1406463 row(s) containing missing values (geom_path). null device 1 Warning message: Removed 1406463 row(s) containing missing values (geom_path). null device 1