### R代码来自Vignette Source'pvca.rnw'############################################################################################(“ Biocmanager”)## Biocmanager :: install(“ PVCA”)############################################################## ###代码块编号2:示例#########################################################“ all.aml”,“ bm.pb”,“ source”)pvcaobj <-pvcabatchassess(golub_merge,batch.fectors.factors,pct_threshold)##########################################################################################barplot(pvcaobj $ dat,xlab =“效果”,## ylab =“加权平均比例变量”,## ylim = c(0,1.1),col = c(“ blue”),las = 2,## main=“ PVCA估计条形图”)## axis(1,at = bp,标签= pvcaobj $ label,xlab =“效果”,cex.axis = 0.5,las = 2)## values = pvcaobj $ dat ## new_values = new_values = round(values,3)## text(bp,pvcaobj $ dat,labels =new_values,pos = 3,cex = 0.8)## ##################################################################################