inst/unitTests/pseudoautoIntensityPlot_test.R

In smgogarten/GWASTools: Tools for Genome Wide Association Studies

test_pseudoautoIntensityPlot <- function() {
  data(illuminaScanADF)
  scanAnnot <- illuminaScanADF
  
  blfile <- system.file("extdata", "illumina_bl.nc", package="GWASdata")
  blnc <- NcdfIntensityReader(blfile)
  blData <-  IntensityData(blnc, scanAnnot=scanAnnot)

  scanID <- getScanID(scanAnnot, index=1)
  pseudoautoIntensityPlot(scan.ids=scanID, intenData=blData)

  pseudoautoIntensityPlot(scan.ids=scanID, intenData=blData, plotY=TRUE)
  close(blData)
}