Plotting fitted dynamics

[huizizhang949]

Hi,

I have fitted the dynamic model to my dataset and would like to do some extra analysis using R. I try to reproduce the plot of the fitted curve
scv.pl.scatter(adata, basis='Igfbpl1', ncols=5, frameon=True)

Here is how I try to reproduce the above plot in R. I extrace the estimates for the rate parameters alpha, beta and gamma (below), simulate a range of time between [0,20], and assign them to either induction or repression by comparing to 'fit_t_'. Then I use the rate equations to compute the values of u and s.

However, the curve I obtain in R is very different and much more flat. Therefore, I wonder if there is anything that I miss or misunderstand? Could you please briefly explain how the learned dynamics is obtained in your function and the differences to mine? Thank you!