copynumber 是一個功能強(qiáng)大拷貝數(shù)分析工具，這里僅限于自己用的的部分功能，看到可用hg38版本，就嘗試運(yùn)行了一下，識別CNV break points。

remotes::install_github("ShixiangWang/copynumber")
library('sequenza')
data("lymphoma")
sub.lymphoma <- subsetData(lymphoma,sample=1:3)
wins.lymph <- winsorize(sub.lymphoma)

#單個樣本
pcf.segments <- pcf(data=wins.lymph,gamma=12,Y=sub.lymphoma, assembly = "hg38")

pcf.segments <- pcf(data=wins.lymph,gamma=12,Y=sub.lymphoma, assembly = "hg19")


##多個樣本：
pcf.segments11 <- multipcf(data= wins1,gamma=12,Y=dat1, assembly = "hg38")
pcf.segments22 <- multipcf(data= wins1,gamma=12,Y=dat2, assembly = "hg38")