Example of a MIC Assay Analysis¶
Collecting the inputs:¶
The first step in the analysis of the MIC assay is to collect all variable inputs .. Link MIC_Input.xlsx template in the repo here MIC_Input.xlsx consists of multiple sheets
Substances
Dataset |
External ID |
Internal ID |
Origin Rack |
Origin Position 96 |
MP Barcode 96 |
MP Position 96 |
|---|---|---|---|---|---|---|
Cooperation Partner A |
12345 |
1 |
Organisms
Organism |
Rack |
|---|---|
Organism A ST1234 |
1 |
Organism B ST5678 |
2 |
Organism C ST9101 |
3 |
… |
… |
Dilutions
Concentration |
|---|
50 |
25 |
12.5 |
6.25 |
3,125 |
1.56 |
0.78 |
0.38 |
0.2 |
0.1 |
0.05 |
Controls
Dataset |
Internal ID |
Position |
|---|---|---|
Negative Control |
Bacteria + Medium |
A23 |
Negative Control |
Bacteria + Medium |
C23 |
Negative Control |
Bacteria + Medium |
E23 |
Negative Control |
Bacteria + Medium |
G23 |
Negative Control |
Bacteria + Medium |
I23 |
Negative Control |
Bacteria + Medium |
K23 |
Negative Control |
Bacteria + Medium |
M23 |
Negative Control |
Bacteria + Medium |
O23 |
Blank |
Medium |
A24 |
Blank |
Medium |
B24 |
Blank |
Medium |
C24 |
Blank |
Medium |
D24 |
Blank |
Medium |
E24 |
Blank |
Medium |
F24 |
Blank |
Medium |
G24 |
Blank |
Medium |
H24 |
… |
… |
… |
Reference |
Rifampicin |
N23 |
Reference |
Vancomycin |
P23 |
1import rda_toolbox as rda
2
3input_mapping = rda.mic_process_inputs(
4 # Input specifications excel:
5 "../data/input/MIC_Input.xlsx",
6 # Barcode reader file which shows Motherplate to AsT plate mapping
7 "../data/input/DiS_MP_AsT_2024-10-08.txt",
8 # Barcode reader file which shows Ast plate to AcD plate mapping
9 "../data/input/AmA_AsT_AcD_20241009.txt",
10)
11input_mapping.to_csv("../data/processed/prepared_input_mapping_table.csv", index=False)
Dataset |
External ID |
Original Rack |
Origin Position |
Internal ID |
MP Barcode 96 |
MP Position 96 |
Row_384 |
Col_384 |
AsT Barcode 384 |
Concentration |
Position |
AcD Barcode 384 |
Replicate |
Organism |
Coop Partner A |
123456789 |
2233445566 |
A1 |
987654321 |
3456677889 |
A1 |
||||||||
1rawfiles = rda.parse_readerfiles("../data/raw/")
2rawfiles.to_csv("../data/processed/rawdata.csv", index=False)
Row_384 |
Col_384 |
Raw Optical Density |
AcD Barcode 384 |
1processed_data = rda.preprocess(
2 rawfiles,
3 input_mapping,
4 substance_id="Internal ID",
5 measurement="Optical Density",
6 negative_controls="Bacteria + Medium",
7 blanks="Medium",
8 norm_by_barcode="AcD Barcode 384"
9)