# #——回声= FALSE,警告= FALSE,消息= FALSE - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -库(HIBAG)库(ggplot2)版本< - c (“v1.24.0 \ n(基线)”、“v1.24.0 \ n (POPCNT)”、“v1.26.1 \ n (SSE2)”、“v1.26.1 \ n (SSE4&POPCNT)”、“v1.26.1 \ n (AVX)”、“v1.26.1 \ n (AVX2)”、“v1.26.1 \ n (AVX512F)”、“v1.26.1 \ n (AVX512BW)”)和< -因子(版本,版本)画< -函数(dat) {ggplot (dat, aes (x = intrinsic, y =加速,填补=基因))+ xlab(“包版本和intrinsic”) + ylab(“加速因素”)+ geom_bar(统计=“身份”,宽度= 0.9,颜色=“白色”,位置= position_dodge ()) + geom_hline (yintercept = 1,线型= 2,颜色=“gray33”) + geom_text (aes(标签=加速),vjust = -0.55, = position_dodge位置(0.9),大小= 2)+ scale_fill_brewer(面板=“Dark2”) + theme_bw()} # # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - #继续而不中断IgnoreError < -函数(cmd) tryCatch (cmd,错误=函数(e){消息(“不支持”);看不见的()})# # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - IgnoreError (hlaSetKernelTarget (“sse4”)) IgnoreError (hlaSetKernelTarget (“avx”)) IgnoreError (hlaSetKernelTarget (“avx2”)) IgnoreError (hlaSetKernelTarget (“avx512f”)) IgnoreError (hlaSetKernelTarget (“avx512bw”)) # #——回声= FALSE, fig.align =‘中心’,fig.width = 9, fig.height = 3.2 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - s < -“抗原HLA-B HLA-C HLA-DRB1 1.0 1.0 1.0 1.0 1.7 1.6 1.6 1.6 1.2 1.1 1.0 1.0 2.3 2.2 2.2 2.2 2.7 2.5 2.8 2.6 3.2 2.7 2.8 2.7 3.3 2.8 3.6 2.9 4.2 3.5 4.6 3.9”b < -阅读。表(文本=年代,头= TRUE) colnames (b) < - gsub (”。”、“——”、colnames (b),固定= TRUE) b1 <——data.frame (intrinsic =代表(vs, ncol (b)),基因= as.factor(代表(名字(b),每个= nrow (b))),加速= unname (unlist (b)))画(b1) + ggtitle(“构建HIBAG模型使用~ 1000样本:“)# #——回声= FALSE, fig.align =‘中心’,fig.width = 9, fig.height = 3.2 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - s < -“抗原HLA-B HLA-C HLA-DRB1 HLA-DQA1 HLA-DQB1 1.0 1.0 1.0 1.0 1.0 1.0 1.6 1.6 1.6 1.6 1.5 1.6 1.1 1.0 1.0 1.0 1.0 1.1 2.2 2.2 2.2 2.2 2.3 2.3 2.6 2.8 2.9 2.8 2.9 2.9 2.7 2.8 2.9 2.9 3.0 3.0 2.9 3.7 3.5 3.4 4.1 3.8 3.5 4.7 4.7 4.6 5.3 5.2”b < -阅读。表(文本=年代,头= TRUE) colnames (b) < - gsub (”。”、“——”、colnames (b),固定= TRUE) b2 <——data.frame (intrinsic =代表(vs, ncol (b)),基因=因素(代表(名字(b),每个= nrow (b))、名称(b)),加速= unname (unlist (b)))画(b2) + ggtitle(“构建HIBAG模型使用~ 5000样本:“)# #——回声= FALSE,警告= FALSE, fig.align =‘中心’,fig.width = 9, fig.height = 3.2 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - s < -“抗原HLA-B HLA-C HLA-DRB1 HLA-DQA1 HLA-DQB1 1.0 1.0 1.0 1.0 1.0 1.0 1.5 1.7 1.7 1.7 1.8 1.7 1.2 1.2 1.2 1.1 1.2 1.2 1.9 2.3 2.3 2.3 2.3 2.4 2.2 2.9 2.8 3.0 3.0 2.9 3.3 3.6 3.6 3.6 3.7 3.7 4.1 4.1 4.4 4.3 4.4 4.5 5.4 6.0 6.4 6.5 6.9 7.0”b < -阅读。表(文本=年代,头= TRUE) colnames (b) < - gsub (”。”、“——”、colnames (b),固定= TRUE) b3 <——data.frame (intrinsic =代表(vs, ncol (b)),基因=因素(代表(名字(b),每个= nrow (b))、名称(b)),加速= unname (unlist (b)))画(b3) + ggtitle(“构建HIBAG模型使用~ 10000样本:“)# #——回声= FALSE,警告= FALSE, fig.align =‘中心’,fig.width = 9, fig.height = 3.2 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - s < -“HLAgenes 1线程8线程16线程1线程8线程线程抗原1.0 7.5 13.8 1.3 9.5 17.5 1.0 HLA-B HLA-C 7.6 14.6 1.6 12.0 22.6 1.0 7.5 13.9 1.6 11.3 20.1 1.0 HLA-DRB1 HLA-DQA1 7.5 14.2 1.6 11.4 21.3 1.0 7.1 12.2 1.9 11.9 19.0 HLA-DQB1 1.0 7.2 12.9 1.7 11.4 19.0”b < -阅读。表(文本=年代,头= TRUE) x < - unname (unlist (b [1])) nt < - c(代表(粘贴(“AVX2: c(1、8、16)),每个= nrow (b)),代表(粘贴(“AVX512BW: c(1、8、16)),每个= nrow (b))) nt < -因子(nt), c (“AVX2: 1”、“AVX2: 8”,“AVX2: 16”,“AVX512BW: 1”、“AVX512BW: 8”,“AVX512BW: 16”)) (nt) ss党卫军< -水平[1 l] < -“AVX2: 1(基线)”水平(nt) < - ss dat < data.frame(加速= x,线程= nt,基因=因素(代表(b HLAgenes美元,6),b HLAgenes美元))ggplot (dat, aes (x =基因,y =加速,填补=线程))+ xlab (“HLA基因/ HIBAG_v1.26”) + ylab(“加速因素”)+ geom_bar(统计=“身份”,宽度= 0.9,颜色=“白色”,位置= position_dodge ()) + geom_hline (yintercept = 1,线型= 2,颜色=“gray33”) + geom_text (aes(标签=加速),vjust = -0.55, = position_dodge位置(0.9),大小= 2)+ scale_fill_brewer(面板=“Dark2”) + theme_bw () + ggtitle(“构建HIBAG模型使用~ 5000样品和AVX2 / AVX512BW intrinsic:“) # # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - sessionInfo ()