import grandpy as gp
import os


sars = gp.read_grand("https://github.com/maus2310/GrandPy/blob/main/tests/data/sars.tsv", design=("Condition", "dur.4sU", "Replicate"))
sars = sars.filter_genes()
print(sars)
# sars = sars.fit_kinetics(steady_state={"Mock": True, "SARS": False})
# gp.plot_scatter(sars, x="kinetics_Mock_Half-life", y="kinetics_SARS_Half-life",
#                 diagonal=True, log=True, limit=(0,1.5), log_labels=True)
# gp.plot_gene_progressive_timecourse(sars, gene="ATF3", steady_state={"Mock": True, "SARS": False})



