Analysis: 20250616
!pip install nucleus-cdk==0.5.0rc2 | tail -n2Requirement already satisfied: asttokens in /opt/homebrew/anaconda3/lib/python3.12/site-packages (from stack-data->ipython>=6.1.0->ipywidgets==8.*->jupyter-bokeh<5.0.0,>=4.0.5->nucleus-cdk==0.5.0rc2) (2.0.5)
Requirement already satisfied: pure-eval in /opt/homebrew/anaconda3/lib/python3.12/site-packages (from stack-data->ipython>=6.1.0->ipywidgets==8.*->jupyter-bokeh<5.0.0,>=4.0.5->nucleus-cdk==0.5.0rc2) (0.2.2)
from cdk.analysis.cytosol import platereader as pr
import matplotlib.pyplot as plt
import warnings
# Filter warnings
warnings.filterwarnings('ignore')
# Initialize plotting
pr.plot_setup()Load the data¶
Provide a CSV file containing the data, and a platemap. This function returns both the data with the plate map mapped to it, and the platemap by itself, which is useful for certain tasks.
data, platemap = pr.load_platereader_data("./data/20250616-160504-cytation5-pure-timecourse-gfp-ppk-cp-biotek-cdk.txt", "20250616-PPK CP-CK platemap.csv")Basic Plots¶
Kinetics¶
Kinetic time traces of every well on the plate
p = pr.plot_curves(data, palette="colorblind")
p.savefig("plot1-timeseries")
Steady state¶
Bar graph of steady-state endpoint of each sample. Steady state is calculated as the maximum fluorescence value over a 3-sample rolling average on the data.
p1= pr.plot_steadystate(data, palette="colorblind")
plt.xlabel("Condition")
p1.savefig("plot2-endpoint")
Kinetics Analysis¶
These functions calculate key kinetic parameters of the time series.
pr.plot_kinetics(data)PROVIDING AVERAGED KINETICS
(<seaborn.axisgrid.FacetGrid at 0x14b8af9d0>,
Velocity \
Time
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 0 days 01:37:02.985893445
PURE w/CP 8mM Mg 2+ GFP-Gext 0 days 01:30:11.183020599
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 0 days 01:16:15.056125978
PURExpress NC (-DNA) GFP-Gext 0 days 08:22:11.208282580
PURExpress PC GFP-Gext 0 days 01:35:03.233243742
\
Data Max
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 1142988.56 1380549.46
PURE w/CP 8mM Mg 2+ GFP-Gext 584893.79 793465.68
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 135845.79 143510.00
PURExpress NC (-DNA) GFP-Gext 1363.52 167.92
PURExpress PC GFP-Gext 850639.37 967516.28
Lag \
Time
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 0 days 00:47:22.392950327
PURE w/CP 8mM Mg 2+ GFP-Gext 0 days 00:45:55.787595782
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 0 days 00:19:21.788669307
PURExpress NC (-DNA) GFP-Gext 0 days 00:13:49.447245369
PURExpress PC GFP-Gext 0 days 00:42:16.995848500
\
Data
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 270166.30
PURE w/CP 8mM Mg 2+ GFP-Gext 137313.14
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 32170.94
PURExpress NC (-DNA) GFP-Gext 322.44
PURExpress PC GFP-Gext 200602.45
Steady State \
Time
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 0 days 02:50:11.072913668
PURE w/CP 8mM Mg 2+ GFP-Gext 0 days 02:35:20.507916682
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 0 days 02:40:00.134997899
PURExpress NC (-DNA) GFP-Gext 0 days 20:21:09.831959821
PURExpress PC GFP-Gext 0 days 02:52:44.629645062
\
Data
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 2171678.26
PURE w/CP 8mM Mg 2+ GFP-Gext 1111298.20
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 258106.99
PURExpress NC (-DNA) GFP-Gext 2590.69
PURExpress PC GFP-Gext 1616214.80
Fit \
params
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext [2285977.111810936, 2.4160489944570442, 1.6174...
PURE w/CP 8mM Mg 2+ GFP-Gext [1169787.5764115045, 2.712290228560942, 1.5031...
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext [271691.5708766907, 2.1120079231950206, 1.2708...
PURExpress NC (-DNA) GFP-Gext [2727.047352489795, 0.2459395018310663, 8.3697...
PURExpress PC GFP-Gext [1701278.733888923, 2.274204732121089, 1.58423...
R^2 drift
Name Read
PURE w/CP + PPK + Polyp 18mM Mg 2+ GFP-Gext 1.00 -2951.38
PURE w/CP 8mM Mg 2+ GFP-Gext 0.99 -2785.05
PURE w/PPK + Polyp 18mM Mg 2+ GFP-Gext 0.97 -671.76
PURExpress NC (-DNA) GFP-Gext 1.00 -8.07
PURExpress PC GFP-Gext 0.99 -3113.84 )
We can also calculate the kinetics and display the parameters as a table.
pr.kinetic_analysis(data)Loading...


