mappingQuality.ssp 4.17 KB
Newer Older
Peter van 't Hof's avatar
Peter van 't Hof committed
1
#import(nl.lumc.sasc.biopet.utils.summary.db.SummaryDb)
2
#import(nl.lumc.sasc.biopet.utils.summary.db.Schema._)
Peter van 't Hof's avatar
Peter van 't Hof committed
3 4
#import(nl.lumc.sasc.biopet.utils.summary.db.SummaryDb._)
#import(nl.lumc.sasc.biopet.utils.summary.db.SummaryDb.Implicts._)
5 6 7
#import(nl.lumc.sasc.biopet.core.report.ReportPage)
#import(nl.lumc.sasc.biopet.pipelines.bammetrics.BammetricsReport)
#import(java.io.File)
Peter van 't Hof's avatar
Peter van 't Hof committed
8 9 10
<%@ var summary: SummaryDb %>
<%@ var sampleId: Option[Int] = None %>
<%@ var libId: Option[Int] = None %>
11
<%@ var sampleLevel: Boolean = false %>
Peter van 't Hof's avatar
Peter van 't Hof committed
12
<%@ var rootPath: String %>
13 14 15 16
<%@ var outputDir: File %>
<%@ var showPlot: Boolean = false %>
<%@ var showTable: Boolean = true %>
<%@ var showIntro: Boolean = true%>
Peter van 't Hof's avatar
Peter van 't Hof committed
17 18
<%@ var runId: Int %>
<%@ var fields: List[String] = List("min", "max", "mean", "median", "modal")%>
19 20
<%@ var allSamples: Seq[Sample] %>
<%@ var allLibraries: Seq[Library] %>
21
<%@ var mappingQualityPlotTables: Array[Map[String, Array[Any]]] %>
22
#{
23 24 25 26
    val samples = sampleId match {
    case Some(id) => allSamples.filter(_.id == id).toList
    case _ => allSamples.toList
    }
27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
}#
#if (showIntro)
    <br/>
    <div class="row">
    <div class="col-md-1"></div>
    <div class="col-md-6">
        <p>
        #if (sampleId.isDefined && libId.isDefined)
            This plot shows the mapping quality distribution for all libraries combined in sample <b>${sampleId}</b>.
        #elseif(sampleId.isDefined)
            This plot shows the mapping quality distribution for the libraries of sample <b>${sampleId}</b>.
        #else
            This plot shows the mapping quality distribution for each of the <b>${samples.size}</b> samples.
        #end
        </p>
    </div>
    </div>
#end

#if (showPlot)
47
    #{ BammetricsReport.mappingQualityPlot(outputDir, "mapping_quality", mappingQualityPlotTables)}#
48 49 50 51 52 53 54 55 56 57 58 59

    <div class="panel-body">
    <img src="mapping_quality.png" class="img-responsive" />
    </div>
    <div class="panel-footer">
        #if (showTable)
            <button type="button" class="btn btn-info" data-toggle="collapse" data-target="#mapping_qualityTable">
                <i class="glyphicon glyphicon-eye-close"></i> Hide table</button>
        #else
            <button type="button" class="btn btn-info" data-toggle="collapse" data-target="#mapping_qualityTable">
            <i class="glyphicon glyphicon-eye-open"></i> Show table</button>
        #end
60
    <a href="mapping_quality.tsv"><button type="button" class="btn btn-info"><i class="glyphicon glyphicon-cloud-download"></i> Data points</button></a>
61 62 63 64 65 66 67 68 69 70 71 72 73 74 75

    </div>
#end

<div class="panel-body collapse #if (showTable)in#end" id="mapping_qualityTable">
<!-- Table -->
<table class="table">
    <thead><tr>
        <th data-sorted="true" data-sorted-direction="ascending">Sample</th>
        #if (!sampleLevel) <th>Library</th> #end
        #for (field <- fields)
            <th>${field.replaceAll("_", " ")}</th>
        #end
    </tr></thead>
    <tbody>
Peter van 't Hof's avatar
Peter van 't Hof committed
76
        #for (sample <- samples.sortBy(_.name))
77
            #{
Peter van 't Hof's avatar
Peter van 't Hof committed
78 79 80
                val libs: List[Option[Int]] = (libId, sampleLevel) match {
                case (_, true) => List(None)
                case (Some(_), _) => List(libId)
81
                case _ => allLibraries.filter(_.sampleId == sample.id).map(x => Some(x.id)).toList
82 83
                }
            }#
84
            <tr><td rowspan="${libs.size}"><a href="${rootPath}Samples/${sample.name}/index.html">${sample.name}</a></td>
85
            #for (libId <- libs)
86
                #{ val libName = libId.map(l => allLibraries.find(_.id == l).get.name) }#
87
                #if (libs.head != libId) <tr> #end
88
                    #if (!sampleLevel) <td><a href="${rootPath}Samples/${sample.name}/Libraries/${libName}/index.html">${libName}</a></td> #end
89
                    #{
Peter van 't Hof's avatar
Peter van 't Hof committed
90
                        val statsPaths = fields.map(x => x -> List("mapping_quality", "general", x)).toMap
Peter van 't Hof's avatar
Peter van 't Hof committed
91
                        val results = summary.getStatKeys(runId, "bammetrics", "bamstats", sample = sample.id, library = libId.map(LibraryId).getOrElse(NoLibrary), keyValues = statsPaths)
92
                    }#
Peter van 't Hof's avatar
Peter van 't Hof committed
93 94
                    #for (field <- fields)
                        <td>${results(field)}</td>
95 96 97 98 99 100 101 102
                    #end
                </tr>
            #end
        #end
    </tbody>
</table>

</div>