Search in sources :

Example 11 with DenseDoubleMatrix

use of ubic.basecode.dataStructure.matrix.DenseDoubleMatrix in project Gemma by PavlidisLab.

the class ExpressionExperimentQCController method visualizeCorrMat.

/**
 * @param id         of experiment
 * @param size       Multiplier on the cell size. 1 or null for standard small size.
 * @param text       if true, output a tabbed file instead of a png
 * @param showLabels if the row and column labels of the matrix should be shown.
 * @param os         response output stream
 */
@RequestMapping("/expressionExperiment/visualizeCorrMat.html")
public ModelAndView visualizeCorrMat(Long id, Double size, String contrVal, Boolean text, Boolean showLabels, Boolean forceShowLabels, OutputStream os) throws Exception {
    if (id == null) {
        log.warn("No id!");
        return null;
    }
    ExpressionExperiment ee = expressionExperimentService.load(id);
    if (ee == null) {
        log.warn("Could not load experiment with id " + id);
        return null;
    }
    DoubleMatrix<BioAssay, BioAssay> omatrix = sampleCoexpressionMatrixService.findOrCreate(ee);
    List<String> stringNames = new ArrayList<>();
    for (BioAssay ba : omatrix.getRowNames()) {
        stringNames.add(ba.getName() + " ID=" + ba.getId());
    }
    DoubleMatrix<String, String> matrix = new DenseDoubleMatrix<>(omatrix.getRawMatrix());
    matrix.setRowNames(stringNames);
    matrix.setColumnNames(stringNames);
    if (text != null && text) {
        StringWriter s = new StringWriter();
        MatrixWriter<String, String> mw = new MatrixWriter<>(s, new DecimalFormat("#.##"));
        mw.writeMatrix(matrix, true);
        ModelAndView mav = new ModelAndView(new TextView());
        mav.addObject(TextView.TEXT_PARAM, s.toString());
        return mav;
    }
    /*
         * Blank out the diagonal so it doesn't affect the colour scale.
         */
    for (int i = 0; i < matrix.rows(); i++) {
        matrix.set(i, i, Double.NaN);
    }
    ColorMatrix<String, String> cm = new ColorMatrix<>(matrix);
    this.cleanNames(matrix);
    int row = matrix.rows();
    int cellsize = (int) Math.min(ExpressionExperimentQCController.MAX_HEATMAP_CELLSIZE, Math.max(1, size * ExpressionExperimentQCController.DEFAULT_QC_IMAGE_SIZE_PX / row));
    MatrixDisplay<String, String> writer = new MatrixDisplay<>(cm);
    boolean reallyShowLabels;
    int minimumCellSizeForText = 9;
    if (forceShowLabels != null && forceShowLabels) {
        cellsize = Math.min(ExpressionExperimentQCController.MAX_HEATMAP_CELLSIZE, minimumCellSizeForText);
        reallyShowLabels = true;
    } else {
        reallyShowLabels = showLabels != null && (showLabels && cellsize >= minimumCellSizeForText);
    }
    writer.setCellSize(new Dimension(cellsize, cellsize));
    boolean showScalebar = size > 2;
    writer.writeToPng(cm, os, reallyShowLabels, showScalebar);
    // nothing to return;
    return null;
}
Also used : DecimalFormat(java.text.DecimalFormat) DoubleArrayList(cern.colt.list.DoubleArrayList) ModelAndView(org.springframework.web.servlet.ModelAndView) ColorMatrix(ubic.basecode.graphics.ColorMatrix) MatrixDisplay(ubic.basecode.graphics.MatrixDisplay) ExpressionExperiment(ubic.gemma.model.expression.experiment.ExpressionExperiment) DenseDoubleMatrix(ubic.basecode.dataStructure.matrix.DenseDoubleMatrix) TextView(ubic.gemma.web.view.TextView) BioAssay(ubic.gemma.model.expression.bioAssay.BioAssay) MatrixWriter(ubic.basecode.io.writer.MatrixWriter) RequestMapping(org.springframework.web.bind.annotation.RequestMapping)

Aggregations

DenseDoubleMatrix (ubic.basecode.dataStructure.matrix.DenseDoubleMatrix)11 CompositeSequence (ubic.gemma.model.expression.designElement.CompositeSequence)8 BioMaterial (ubic.gemma.model.expression.biomaterial.BioMaterial)7 ExpressionDataDoubleMatrix (ubic.gemma.core.datastructure.matrix.ExpressionDataDoubleMatrix)5 BioAssay (ubic.gemma.model.expression.bioAssay.BioAssay)5 DoubleMatrix2D (cern.colt.matrix.DoubleMatrix2D)4 ExpressionExperiment (ubic.gemma.model.expression.experiment.ExpressionExperiment)3 DoubleArrayList (cern.colt.list.DoubleArrayList)2 DenseDoubleMatrix2D (cern.colt.matrix.impl.DenseDoubleMatrix2D)2 Algebra (cern.colt.matrix.linalg.Algebra)2 QuantitationType (ubic.gemma.model.common.quantitationtype.QuantitationType)2 DoubleMatrix1D (cern.colt.matrix.DoubleMatrix1D)1 DecimalFormat (java.text.DecimalFormat)1 ArrayList (java.util.ArrayList)1 Map (java.util.Map)1 Test (org.junit.Test)1 RequestMapping (org.springframework.web.bind.annotation.RequestMapping)1 ModelAndView (org.springframework.web.servlet.ModelAndView)1 ColorMatrix (ubic.basecode.graphics.ColorMatrix)1 MatrixDisplay (ubic.basecode.graphics.MatrixDisplay)1