data(tamoxifen_peaks)
# peakcaller scores PCA
dba.plotPCA(tamoxifen)
# raw count correlation PCA
data(tamoxifen_analysis)
dba.plotPCA(tamoxifen)
#PCA based on normalized data for all sites
dba.plotPCA(tamoxifen,contrast=1,th=1)
#PCA based on DB sites only
p <- dba.plotPCA(tamoxifen,contrast=1)
p <- dba.plotPCA(tamoxifen,contrast=1,attributes=DBA_TISSUE)
p <- dba.plotPCA(tamoxifen,contrast=1,attributes=DBA_TISSUE,label=DBA_CONDITION)
p <- dba.plotPCA(tamoxifen,contrast=1,attributes=DBA_CONDITION,label=DBA_TISSUE)
p <- dba.plotPCA(tamoxifen,contrast=1,attributes=c(DBA_TISSUE,DBA_CONDITION),
label=DBA_REPLICATE)
Run the code above in your browser using DataLab