class: title-slide, middle, left ## Session 5: How to perform quality control of your data ### Dr Shazia Ruybal-Pesántez Presented at the [covidClassifyR Shiny app workshop](https://shaziaruybal.github.io/covidClassifyR-workshop) for researchers from PNGIMR and partner institutions 2022-03-11 --- class: left # Recap -- ####
**[Session 1](https://shaziaruybal.github.io/covidClassifyR-workshop/sessions/2022-03-10-session1/index.html)** gave you an overview of the Luminex technology and the COVID-19 multi-antigen serological assay that was established in PNG -- ####
**Session 2** gave you an overview of the purpose of serosurveillance and insights into the COVID-19 sero-surveys that have been carried out in PNG -- ####
**[Session 3](https://shaziaruybal.github.io/covidClassifyR-workshop/sessions/2022-03-10-session3/index.html)** gave you an overview of the [`covidClassifyR`](https://shaziaruybal.shinyapps.io/covidclassifyr) Shiny app -- ####
**[Session 4](https://shaziaruybal.github.io/covidClassifyR-workshop/sessions/2022-03-11-session4/index.html)** gave you an overview of how to prepare and import your data to the [`covidClassifyR`](https://shaziaruybal.shinyapps.io/covidclassifyr) Shiny app -- .footnote[ For all the workshop materials see the [workshop website
](https://shaziaruybal.github.io/covidClassifyR-workshop/materials.html) ] --- class: center #
# Today we will cover: -- ### How to perform quality control of your Luminex data -- ### How to interpet the QC plots -- ### You can find the app [here
](https://shaziaruybal.shinyapps.io/covidClassifyR) --- class: inverse, middle, center #
## Quality control -- ### Note all of the plots in this interactive section will be downloaded in your QC report! --- ## Check your standard curves .left-column[ The standard curve plots are generated from the antibody data from the standards you indicated in your plate layout (e.g. S1-S10) and median fluorescent intensity (MFI) units are displayed in log10-scale. Remember that there are 11 antigens in the COVID-19 assay panel, all of the COVID-19 antigens start with a "C". **
Tip:** In general your standard curves should look relatively linear (only when the y-axis is on logarithmic scale). ] .right-column[ ![](img/std_curves.png) ] --- # Check your plate bead counts .left-column[ A summary of the bead counts for each plate well are displayed, with blue indicating there are sufficient beads (>15) or red when there are not enough. If any of the wells are red, they should be double-checked manually and re-run on a new plate if required. ] .right-column[ ![](img/plate_qc.png) ] --- # Check your blank samples .left-column[ The median fluorescent intensity (MFI) units for each antigen is displayed for your blank samples. In general, each blank sample should have <50 MFI for each antigen, if they are higher they should be cross-checked manually. **
Tip:** If you have more than one blank sample, make sure you label them differently in your plate layout (e.g. Blank1, Blank2), otherwise the figure will display cumulative MFI values for all your blank samples. ] .right-column[ ![](img/blanks.png) ] --- .pull-left[ ## Automated curve fitting for our Luminex multiplex assay The standard curve on each plate allows us to use the median fluorescence intensity (MFI) values recorded for each standard known concentration/dilution to estimate the antibody concentration for each unknown sample. When plotting the relationship between MFI and antibody concentration for our multiplex serological assay, it results in "S"-shaped curves, which is due to the non-linear relationship between MFI and antibody concentration. To account for the assymetry in this non-linear relationship, we can fit a logistic regression model using 4 parameters (4PL). .footenote[ .bottom[ _[Resource
](https://www.biolegend.com/en-us/blog/curve-fitting-for-immunoassays-legendplex)_ ]] ] .pull-right[ <img src="img/4pl_model.png" width="650" height="550" /> ] --- .pull-left[ ## Automated curve fitting for our Luminex multiplex assay However, for more complex multiplex assays it is more appropriate to also account for the asymmetry to provide a better fit (by including the asymmetry factor). This is called a 5-parameter logistic (5PL) regression model. .footnote[ .bottom[ _[Resource
](https://www.biolegend.com/en-us/blog/curve-fitting-for-immunoassays-legendplex)_ ]] ] .pull-right[ <img src="img/5pl_model.png" width="650" height="350" /> ] --- # Check your model results .left-column[ The automated data processing in this app allows you to convert your median fluorescence intensity (MFI) data into relative antibody units (RAU) by fitting a 5-parameter logistic (5PL) regression model to the standard curve on a per-antigen level. This log-log model is used to obtain a more linear relationship. **
Tip:** In general your model results for each antigen should look relatively linear. ] .right-column[ ![](img/model.png) ] --- # Check the processed data .left-column[ The results from the data processing are displayed in an interactive table within the app. The raw median fluorescent intensity (MFI) units and converted relative antibody units (RAU) are displayed for each antigen and for each sample.
Scroll to the right to explore data for all antigens
Navigate to different pages to view all samples
Sort & filter as needed! ] .right-column[ ![](img/sample_results.png) **
Tip:** Check that the number of entries is as expected baed on the number of unknown samples in your plate layout and that RAU values have been calculated for all of your antigens. ] --- #
Demo <video width="1330" height="410" controls> <source src="screen_records/qc_demo.mp4" type="video/mp4"> </video> --- class: inverse, middle, center #
## Next session: ### Download your processed data and QC report --- # Acknowledgments - Dr Maria Ome-Kaius and Dr Fiona Angrisano - PNGIMR and partner institutions - WEHI & Burnet Institute - All of you for attending! *We are extremely grateful for financial support to develop and host the covidClassifyR Shiny web application, and to host these virtual workshops through the [Regional Collaborations Programme COVID-19 Digital Grant](https://www.science.org.au/news-and-events/news-and-media-releases/regional-research-set-get-digital-boost) from the Australian Academy of Science and Australian Department of Industry, Science, Energy and Resources.* The scripts and functions used in [`covidClassifyR`](https://shaziaruybal.shinyapps.io/covidclassifyr) were developed by Shazia Ruybal-Pesántez, with contributions from the following researchers: Eamon Conway, Connie Li Wan Suen, Narimane Nekkab and Michael White. .footnote[ _These slides were created using the R packages: [xaringan](https://github.com/yihui/xaringan), [xaringanthemer](https://github.com/gadenbuie/xaringanthemer), [xaringanExtra](https://github.com/gadenbuie/xaringanExtra)_ ] --- name: contact class: inverse .pull-left[ .center[ ### Dr Shazia Ruybal-Pesántez <img style="border-radius: 50%;" src="https://shaziaruybal.github.io/covidClassifyR-workshop/images/shazia.png" width="250px"/> #### Contact details [
ruybal.s@wehi.edu.au](mailto:ruybal.s@wehi.edu.au) [
@DrShaziaRuybal](https://twitter.com/DrShaziaRuybal) ]] .pull-right[ ### Session 5 Resources: ### [
Recording](https://youtu.be/bglxoyCM2LY) ### [
`covidClassifyR`](https://shaziaruybal.shinyapps.io/covidClassifyR) ### [
Workshop materials](https://shaziaruybal.github.io/covidClassifyR-workshop/materials.html) ### [
Slides for Session 5](https://shaziaruybal.github.io/covidClassifyR-workshop/slides/session5/session5_slides.html) ] --- class: inverse, middle, center #
Questions?