Getting started#

Installing and getting started with FPsim is quite simple.

To install, just type pip install fpsim. If it worked, you should be able to import FPsim with import fpsim as fp.

The basic design philosophy of FPsim is: common tasks should be simple. For example:

  • Defining parameters

  • Running a simulation

  • Plotting results

This tutorial walks you through how to define parameters and run the simulation. The next tutorial will show you how to plot the results of a simulation.

An interactive version of this notebook is available on Google Colab or Binder.

Hello world#

To create, run, and plot a sim with default options is just:

[1]:
import fpsim as fp

sim = fp.Sim()
sim.run()
fig = sim.plot()
/home/docs/checkouts/readthedocs.org/user_builds/institute-for-disease-modeling-fpsim/envs/latest/lib/python3.9/site-packages/tqdm/auto.py:21: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html
  from .autonotebook import tqdm as notebook_tqdm
  Running 1960 of 2020 ( 0/721) (0.09 s)  ———————————————————— 0%
  Running 1961 of 2020 (12/721) (0.17 s)  ———————————————————— 2%
  Running 1962 of 2020 (24/721) (0.27 s)  ———————————————————— 3%
  Running 1963 of 2020 (36/721) (0.37 s)  •——————————————————— 5%
  Running 1964 of 2020 (48/721) (0.47 s)  •——————————————————— 7%
  Running 1965 of 2020 (60/721) (0.57 s)  •——————————————————— 8%
  Running 1966 of 2020 (72/721) (0.66 s)  ••—————————————————— 10%
  Running 1967 of 2020 (84/721) (0.76 s)  ••—————————————————— 12%
  Running 1968 of 2020 (96/721) (0.86 s)  ••—————————————————— 13%
  Running 1969 of 2020 (108/721) (0.97 s)  •••————————————————— 15%
  Running 1970 of 2020 (120/721) (1.07 s)  •••————————————————— 17%
  Running 1971 of 2020 (132/721) (1.17 s)  •••————————————————— 18%
  Running 1972 of 2020 (144/721) (1.27 s)  ••••———————————————— 20%
  Running 1973 of 2020 (156/721) (1.37 s)  ••••———————————————— 22%
  Running 1974 of 2020 (168/721) (1.48 s)  ••••———————————————— 23%
  Running 1975 of 2020 (180/721) (1.58 s)  •••••——————————————— 25%
  Running 1976 of 2020 (192/721) (1.68 s)  •••••——————————————— 27%
  Running 1977 of 2020 (204/721) (1.79 s)  •••••——————————————— 28%
  Running 1978 of 2020 (216/721) (1.89 s)  ••••••—————————————— 30%
  Running 1979 of 2020 (228/721) (2.00 s)  ••••••—————————————— 32%
  Running 1980 of 2020 (240/721) (2.10 s)  ••••••—————————————— 33%
  Running 1981 of 2020 (252/721) (2.21 s)  •••••••————————————— 35%
  Running 1982 of 2020 (264/721) (2.31 s)  •••••••————————————— 37%
  Running 1983 of 2020 (276/721) (2.42 s)  •••••••————————————— 38%
  Running 1984 of 2020 (288/721) (2.52 s)  ••••••••———————————— 40%
  Running 1985 of 2020 (300/721) (2.63 s)  ••••••••———————————— 42%
  Running 1986 of 2020 (312/721) (2.74 s)  ••••••••———————————— 43%
  Running 1987 of 2020 (324/721) (2.85 s)  •••••••••——————————— 45%
  Running 1988 of 2020 (336/721) (2.96 s)  •••••••••——————————— 47%
  Running 1989 of 2020 (348/721) (3.07 s)  •••••••••——————————— 48%
  Running 1990 of 2020 (360/721) (3.18 s)  ••••••••••—————————— 50%
  Running 1991 of 2020 (372/721) (3.29 s)  ••••••••••—————————— 52%
  Running 1992 of 2020 (384/721) (3.41 s)  ••••••••••—————————— 53%
  Running 1993 of 2020 (396/721) (3.52 s)  •••••••••••————————— 55%
  Running 1994 of 2020 (408/721) (3.64 s)  •••••••••••————————— 57%
  Running 1995 of 2020 (420/721) (3.75 s)  •••••••••••————————— 58%
  Running 1996 of 2020 (432/721) (3.87 s)  ••••••••••••———————— 60%
  Running 1997 of 2020 (444/721) (3.98 s)  ••••••••••••———————— 62%
  Running 1998 of 2020 (456/721) (4.09 s)  ••••••••••••———————— 63%
  Running 1999 of 2020 (468/721) (4.21 s)  •••••••••••••——————— 65%
  Running 2000 of 2020 (480/721) (4.32 s)  •••••••••••••——————— 67%
  Running 2001 of 2020 (492/721) (4.43 s)  •••••••••••••——————— 68%
  Running 2002 of 2020 (504/721) (4.55 s)  ••••••••••••••—————— 70%
  Running 2003 of 2020 (516/721) (4.66 s)  ••••••••••••••—————— 72%
  Running 2004 of 2020 (528/721) (4.77 s)  ••••••••••••••—————— 73%
  Running 2005 of 2020 (540/721) (4.89 s)  •••••••••••••••————— 75%
  Running 2006 of 2020 (552/721) (5.01 s)  •••••••••••••••————— 77%
  Running 2007 of 2020 (564/721) (5.12 s)  •••••••••••••••————— 78%
  Running 2008 of 2020 (576/721) (5.24 s)  ••••••••••••••••———— 80%
  Running 2009 of 2020 (588/721) (5.35 s)  ••••••••••••••••———— 82%
  Running 2010 of 2020 (600/721) (5.47 s)  ••••••••••••••••———— 83%
  Running 2011 of 2020 (612/721) (5.59 s)  •••••••••••••••••——— 85%
  Running 2012 of 2020 (624/721) (5.70 s)  •••••••••••••••••——— 87%
  Running 2013 of 2020 (636/721) (5.82 s)  •••••••••••••••••——— 88%
  Running 2014 of 2020 (648/721) (5.93 s)  ••••••••••••••••••—— 90%
  Running 2015 of 2020 (660/721) (6.05 s)  ••••••••••••••••••—— 92%
  Running 2016 of 2020 (672/721) (6.17 s)  ••••••••••••••••••—— 93%
  Running 2017 of 2020 (684/721) (6.29 s)  •••••••••••••••••••— 95%
  Running 2018 of 2020 (696/721) (6.41 s)  •••••••••••••••••••— 97%
  Running 2019 of 2020 (708/721) (6.53 s)  •••••••••••••••••••— 98%
  Running 2020 of 2020 (720/721) (6.65 s)  •••••••••••••••••••• 100%

Final population size: 3491.
Run finished for "None" after 6.7 s
../_images/tutorials_T1_intro_3_2.png

Defining parameters and running simulations#

Parameters are defined as a dictionary. In FPsim, we categorize our parameters as:

  • Basic parameters

  • Age limits

  • Durations

  • Pregnancy outcomes

  • Fecundity and exposure

  • MCPR

The most common category of parameters to change in FPsim is the basic category, which includes the location (i.e. Senegal, northern India), the starting population, the starting year, and the initial number of agents. We can define thiese as:

[2]:
pars = dict(
    n_agents   = 10_000,
    location   = 'senegal',
    start_year = 1960,
    end_year   = 2020,
)

Running a simulation is pretty easy. In fact, running a sim with the parameters we defined above is just:

[3]:
sim = fp.Sim(pars)
sim.run()
  Running 1960 of 2020 ( 0/721) (0.06 s)  ———————————————————— 0%
  Running 1961 of 2020 (12/721) (0.19 s)  ———————————————————— 2%
  Running 1962 of 2020 (24/721) (0.34 s)  ———————————————————— 3%
  Running 1963 of 2020 (36/721) (0.61 s)  •——————————————————— 5%
  Running 1964 of 2020 (48/721) (0.76 s)  •——————————————————— 7%
  Running 1965 of 2020 (60/721) (0.92 s)  •——————————————————— 8%
  Running 1966 of 2020 (72/721) (1.08 s)  ••—————————————————— 10%
  Running 1967 of 2020 (84/721) (1.24 s)  ••—————————————————— 12%
  Running 1968 of 2020 (96/721) (1.40 s)  ••—————————————————— 13%
  Running 1969 of 2020 (108/721) (1.57 s)  •••————————————————— 15%
  Running 1970 of 2020 (120/721) (1.74 s)  •••————————————————— 17%
  Running 1971 of 2020 (132/721) (1.91 s)  •••————————————————— 18%
  Running 1972 of 2020 (144/721) (2.08 s)  ••••———————————————— 20%
  Running 1973 of 2020 (156/721) (2.26 s)  ••••———————————————— 22%
  Running 1974 of 2020 (168/721) (2.43 s)  ••••———————————————— 23%
  Running 1975 of 2020 (180/721) (2.62 s)  •••••——————————————— 25%
  Running 1976 of 2020 (192/721) (2.81 s)  •••••——————————————— 27%
  Running 1977 of 2020 (204/721) (2.99 s)  •••••——————————————— 28%
  Running 1978 of 2020 (216/721) (3.18 s)  ••••••—————————————— 30%
  Running 1979 of 2020 (228/721) (3.38 s)  ••••••—————————————— 32%
  Running 1980 of 2020 (240/721) (3.58 s)  ••••••—————————————— 33%
  Running 1981 of 2020 (252/721) (3.78 s)  •••••••————————————— 35%
  Running 1982 of 2020 (264/721) (3.98 s)  •••••••————————————— 37%
  Running 1983 of 2020 (276/721) (4.18 s)  •••••••————————————— 38%
  Running 1984 of 2020 (288/721) (4.39 s)  ••••••••———————————— 40%
  Running 1985 of 2020 (300/721) (4.60 s)  ••••••••———————————— 42%
  Running 1986 of 2020 (312/721) (4.83 s)  ••••••••———————————— 43%
  Running 1987 of 2020 (324/721) (5.05 s)  •••••••••——————————— 45%
  Running 1988 of 2020 (336/721) (5.27 s)  •••••••••——————————— 47%
  Running 1989 of 2020 (348/721) (5.50 s)  •••••••••——————————— 48%
  Running 1990 of 2020 (360/721) (5.73 s)  ••••••••••—————————— 50%
  Running 1991 of 2020 (372/721) (5.96 s)  ••••••••••—————————— 52%
  Running 1992 of 2020 (384/721) (6.20 s)  ••••••••••—————————— 53%
  Running 1993 of 2020 (396/721) (6.54 s)  •••••••••••————————— 55%
  Running 1994 of 2020 (408/721) (6.78 s)  •••••••••••————————— 57%
  Running 1995 of 2020 (420/721) (7.02 s)  •••••••••••————————— 58%
  Running 1996 of 2020 (432/721) (7.27 s)  ••••••••••••———————— 60%
  Running 1997 of 2020 (444/721) (7.52 s)  ••••••••••••———————— 62%
  Running 1998 of 2020 (456/721) (7.77 s)  ••••••••••••———————— 63%
  Running 1999 of 2020 (468/721) (8.01 s)  •••••••••••••——————— 65%
  Running 2000 of 2020 (480/721) (8.27 s)  •••••••••••••——————— 67%
  Running 2001 of 2020 (492/721) (8.52 s)  •••••••••••••——————— 68%
  Running 2002 of 2020 (504/721) (8.79 s)  ••••••••••••••—————— 70%
  Running 2003 of 2020 (516/721) (9.05 s)  ••••••••••••••—————— 72%
  Running 2004 of 2020 (528/721) (9.32 s)  ••••••••••••••—————— 73%
  Running 2005 of 2020 (540/721) (9.58 s)  •••••••••••••••————— 75%
  Running 2006 of 2020 (552/721) (9.85 s)  •••••••••••••••————— 77%
  Running 2007 of 2020 (564/721) (10.11 s)  •••••••••••••••————— 78%
  Running 2008 of 2020 (576/721) (10.39 s)  ••••••••••••••••———— 80%
  Running 2009 of 2020 (588/721) (10.66 s)  ••••••••••••••••———— 82%
  Running 2010 of 2020 (600/721) (10.94 s)  ••••••••••••••••———— 83%
  Running 2011 of 2020 (612/721) (11.22 s)  •••••••••••••••••——— 85%
  Running 2012 of 2020 (624/721) (11.50 s)  •••••••••••••••••——— 87%
  Running 2013 of 2020 (636/721) (11.79 s)  •••••••••••••••••——— 88%
  Running 2014 of 2020 (648/721) (12.08 s)  ••••••••••••••••••—— 90%
  Running 2015 of 2020 (660/721) (12.37 s)  ••••••••••••••••••—— 92%
  Running 2016 of 2020 (672/721) (12.66 s)  ••••••••••••••••••—— 93%
  Running 2017 of 2020 (684/721) (12.96 s)  •••••••••••••••••••— 95%
  Running 2018 of 2020 (696/721) (13.26 s)  •••••••••••••••••••— 97%
  Running 2019 of 2020 (708/721) (13.56 s)  •••••••••••••••••••— 98%
  Running 2020 of 2020 (720/721) (13.86 s)  •••••••••••••••••••• 100%

Final population size: 30670.
Run finished for "None" after 14.0 s
[3]:
Sim(<no label>; n=10000; 1960-2020; results: b=37491 ☠=13893 pop=30670)

This will generate a results dictionary sim.results. For example, the number of pregnancies in the sim can be found using sim.results['pregnancies'].

Rather than creating a parameter dictionary, any valid parameter can also be passed to the sim directly. For example, exactly equivalent to the above is:

[4]:
sim = fp.Sim(n_agents=10e3, location='senegal', start_year=1960, end_year=2020)
sim.run()
  Running 1960 of 2020 ( 0/721) (0.20 s)  ———————————————————— 0%
  Running 1961 of 2020 (12/721) (0.33 s)  ———————————————————— 2%
  Running 1962 of 2020 (24/721) (0.48 s)  ———————————————————— 3%
  Running 1963 of 2020 (36/721) (0.63 s)  •——————————————————— 5%
  Running 1964 of 2020 (48/721) (0.79 s)  •——————————————————— 7%
  Running 1965 of 2020 (60/721) (0.95 s)  •——————————————————— 8%
  Running 1966 of 2020 (72/721) (1.11 s)  ••—————————————————— 10%
  Running 1967 of 2020 (84/721) (1.27 s)  ••—————————————————— 12%
  Running 1968 of 2020 (96/721) (1.43 s)  ••—————————————————— 13%
  Running 1969 of 2020 (108/721) (1.60 s)  •••————————————————— 15%
  Running 1970 of 2020 (120/721) (1.77 s)  •••————————————————— 17%
  Running 1971 of 2020 (132/721) (1.94 s)  •••————————————————— 18%
  Running 1972 of 2020 (144/721) (2.11 s)  ••••———————————————— 20%
  Running 1973 of 2020 (156/721) (2.29 s)  ••••———————————————— 22%
  Running 1974 of 2020 (168/721) (2.47 s)  ••••———————————————— 23%
  Running 1975 of 2020 (180/721) (2.65 s)  •••••——————————————— 25%
  Running 1976 of 2020 (192/721) (2.83 s)  •••••——————————————— 27%
  Running 1977 of 2020 (204/721) (3.02 s)  •••••——————————————— 28%
  Running 1978 of 2020 (216/721) (3.21 s)  ••••••—————————————— 30%
  Running 1979 of 2020 (228/721) (3.40 s)  ••••••—————————————— 32%
  Running 1980 of 2020 (240/721) (3.60 s)  ••••••—————————————— 33%
  Running 1981 of 2020 (252/721) (3.80 s)  •••••••————————————— 35%
  Running 1982 of 2020 (264/721) (4.01 s)  •••••••————————————— 37%
  Running 1983 of 2020 (276/721) (4.21 s)  •••••••————————————— 38%
  Running 1984 of 2020 (288/721) (4.42 s)  ••••••••———————————— 40%
  Running 1985 of 2020 (300/721) (4.64 s)  ••••••••———————————— 42%
  Running 1986 of 2020 (312/721) (4.85 s)  ••••••••———————————— 43%
  Running 1987 of 2020 (324/721) (5.07 s)  •••••••••——————————— 45%
  Running 1988 of 2020 (336/721) (5.30 s)  •••••••••——————————— 47%
  Running 1989 of 2020 (348/721) (5.52 s)  •••••••••——————————— 48%
  Running 1990 of 2020 (360/721) (5.75 s)  ••••••••••—————————— 50%
  Running 1991 of 2020 (372/721) (5.98 s)  ••••••••••—————————— 52%
  Running 1992 of 2020 (384/721) (6.36 s)  ••••••••••—————————— 53%
  Running 1993 of 2020 (396/721) (6.59 s)  •••••••••••————————— 55%
  Running 1994 of 2020 (408/721) (6.84 s)  •••••••••••————————— 57%
  Running 1995 of 2020 (420/721) (7.08 s)  •••••••••••————————— 58%
  Running 1996 of 2020 (432/721) (7.33 s)  ••••••••••••———————— 60%
  Running 1997 of 2020 (444/721) (7.58 s)  ••••••••••••———————— 62%
  Running 1998 of 2020 (456/721) (7.83 s)  ••••••••••••———————— 63%
  Running 1999 of 2020 (468/721) (8.08 s)  •••••••••••••——————— 65%
  Running 2000 of 2020 (480/721) (8.34 s)  •••••••••••••——————— 67%
  Running 2001 of 2020 (492/721) (8.59 s)  •••••••••••••——————— 68%
  Running 2002 of 2020 (504/721) (8.86 s)  ••••••••••••••—————— 70%
  Running 2003 of 2020 (516/721) (9.12 s)  ••••••••••••••—————— 72%
  Running 2004 of 2020 (528/721) (9.38 s)  ••••••••••••••—————— 73%
  Running 2005 of 2020 (540/721) (9.65 s)  •••••••••••••••————— 75%
  Running 2006 of 2020 (552/721) (9.92 s)  •••••••••••••••————— 77%
  Running 2007 of 2020 (564/721) (10.19 s)  •••••••••••••••————— 78%
  Running 2008 of 2020 (576/721) (10.46 s)  ••••••••••••••••———— 80%
  Running 2009 of 2020 (588/721) (10.74 s)  ••••••••••••••••———— 82%
  Running 2010 of 2020 (600/721) (11.02 s)  ••••••••••••••••———— 83%
  Running 2011 of 2020 (612/721) (11.31 s)  •••••••••••••••••——— 85%
  Running 2012 of 2020 (624/721) (11.60 s)  •••••••••••••••••——— 87%
  Running 2013 of 2020 (636/721) (11.89 s)  •••••••••••••••••——— 88%
  Running 2014 of 2020 (648/721) (12.18 s)  ••••••••••••••••••—— 90%
  Running 2015 of 2020 (660/721) (12.48 s)  ••••••••••••••••••—— 92%
  Running 2016 of 2020 (672/721) (12.77 s)  ••••••••••••••••••—— 93%
  Running 2017 of 2020 (684/721) (13.07 s)  •••••••••••••••••••— 95%
  Running 2018 of 2020 (696/721) (13.36 s)  •••••••••••••••••••— 97%
  Running 2019 of 2020 (708/721) (13.67 s)  •••••••••••••••••••— 98%
  Running 2020 of 2020 (720/721) (13.97 s)  •••••••••••••••••••• 100%

Final population size: 30670.
Run finished for "None" after 14.1 s
[4]:
Sim(<no label>; n=10000; 1960-2020; results: b=37491 ☠=13893 pop=30670)

You can mix and match too – pass in a parameter dictionary with default options, and then include other parameters as keywords (including overrides; keyword arguments take precedence). For example:

[5]:
sim = fp.Sim(pars, n_agents=100) # Use parameters defined above, except start with 100 agents instead of 10,000
sim.run()
  Running 1960 of 2020 ( 0/721) (0.05 s)  ———————————————————— 0%
  Running 1961 of 2020 (12/721) (0.12 s)  ———————————————————— 2%
  Running 1962 of 2020 (24/721) (0.20 s)  ———————————————————— 3%
  Running 1963 of 2020 (36/721) (0.28 s)  •——————————————————— 5%
  Running 1964 of 2020 (48/721) (0.37 s)  •——————————————————— 7%
  Running 1965 of 2020 (60/721) (0.45 s)  •——————————————————— 8%
  Running 1966 of 2020 (72/721) (0.53 s)  ••—————————————————— 10%
  Running 1967 of 2020 (84/721) (0.61 s)  ••—————————————————— 12%
  Running 1968 of 2020 (96/721) (0.69 s)  ••—————————————————— 13%
  Running 1969 of 2020 (108/721) (0.78 s)  •••————————————————— 15%
  Running 1970 of 2020 (120/721) (0.87 s)  •••————————————————— 17%
  Running 1971 of 2020 (132/721) (0.95 s)  •••————————————————— 18%
  Running 1972 of 2020 (144/721) (1.03 s)  ••••———————————————— 20%
  Running 1973 of 2020 (156/721) (1.12 s)  ••••———————————————— 22%
  Running 1974 of 2020 (168/721) (1.20 s)  ••••———————————————— 23%
  Running 1975 of 2020 (180/721) (1.29 s)  •••••——————————————— 25%
  Running 1976 of 2020 (192/721) (1.38 s)  •••••——————————————— 27%
  Running 1977 of 2020 (204/721) (1.46 s)  •••••——————————————— 28%
  Running 1978 of 2020 (216/721) (1.55 s)  ••••••—————————————— 30%
  Running 1979 of 2020 (228/721) (1.64 s)  ••••••—————————————— 32%
  Running 1980 of 2020 (240/721) (1.72 s)  ••••••—————————————— 33%
  Running 1981 of 2020 (252/721) (1.81 s)  •••••••————————————— 35%
  Running 1982 of 2020 (264/721) (1.89 s)  •••••••————————————— 37%
  Running 1983 of 2020 (276/721) (1.98 s)  •••••••————————————— 38%
  Running 1984 of 2020 (288/721) (2.07 s)  ••••••••———————————— 40%
  Running 1985 of 2020 (300/721) (2.15 s)  ••••••••———————————— 42%
  Running 1986 of 2020 (312/721) (2.24 s)  ••••••••———————————— 43%
  Running 1987 of 2020 (324/721) (2.33 s)  •••••••••——————————— 45%
  Running 1988 of 2020 (336/721) (2.42 s)  •••••••••——————————— 47%
  Running 1989 of 2020 (348/721) (2.50 s)  •••••••••——————————— 48%
  Running 1990 of 2020 (360/721) (2.59 s)  ••••••••••—————————— 50%
  Running 1991 of 2020 (372/721) (2.68 s)  ••••••••••—————————— 52%
  Running 1992 of 2020 (384/721) (2.76 s)  ••••••••••—————————— 53%
  Running 1993 of 2020 (396/721) (2.85 s)  •••••••••••————————— 55%
  Running 1994 of 2020 (408/721) (2.94 s)  •••••••••••————————— 57%
  Running 1995 of 2020 (420/721) (3.02 s)  •••••••••••————————— 58%
  Running 1996 of 2020 (432/721) (3.11 s)  ••••••••••••———————— 60%
  Running 1997 of 2020 (444/721) (3.20 s)  ••••••••••••———————— 62%
  Running 1998 of 2020 (456/721) (3.29 s)  ••••••••••••———————— 63%
  Running 1999 of 2020 (468/721) (3.37 s)  •••••••••••••——————— 65%
  Running 2000 of 2020 (480/721) (3.46 s)  •••••••••••••——————— 67%
  Running 2001 of 2020 (492/721) (3.55 s)  •••••••••••••——————— 68%
  Running 2002 of 2020 (504/721) (3.63 s)  ••••••••••••••—————— 70%
  Running 2003 of 2020 (516/721) (3.72 s)  ••••••••••••••—————— 72%
  Running 2004 of 2020 (528/721) (3.80 s)  ••••••••••••••—————— 73%
  Running 2005 of 2020 (540/721) (3.89 s)  •••••••••••••••————— 75%
  Running 2006 of 2020 (552/721) (3.98 s)  •••••••••••••••————— 77%
  Running 2007 of 2020 (564/721) (4.07 s)  •••••••••••••••————— 78%
  Running 2008 of 2020 (576/721) (4.15 s)  ••••••••••••••••———— 80%
  Running 2009 of 2020 (588/721) (4.24 s)  ••••••••••••••••———— 82%
  Running 2010 of 2020 (600/721) (4.33 s)  ••••••••••••••••———— 83%
  Running 2011 of 2020 (612/721) (4.42 s)  •••••••••••••••••——— 85%
  Running 2012 of 2020 (624/721) (4.50 s)  •••••••••••••••••——— 87%
  Running 2013 of 2020 (636/721) (4.59 s)  •••••••••••••••••——— 88%
  Running 2014 of 2020 (648/721) (4.68 s)  ••••••••••••••••••—— 90%
  Running 2015 of 2020 (660/721) (4.77 s)  ••••••••••••••••••—— 92%
  Running 2016 of 2020 (672/721) (4.85 s)  ••••••••••••••••••—— 93%
  Running 2017 of 2020 (684/721) (4.94 s)  •••••••••••••••••••— 95%
  Running 2018 of 2020 (696/721) (5.02 s)  •••••••••••••••••••— 97%
  Running 2019 of 2020 (708/721) (5.11 s)  •••••••••••••••••••— 98%
  Running 2020 of 2020 (720/721) (5.19 s)  •••••••••••••••••••• 100%

Final population size: 272.
Run finished for "None" after 5.2 s
[5]:
Sim(<no label>; n=100; 1960-2020; results: b=350 ☠=157 pop=272)

Now you know how to run a basic simulation in FPsim and change the parameters. Now let’s take a look at the output of the sim.

Explore plotting options for a single sim#

Let’s take a look at the basic suite of plotting options, once we’ve run our initial simulation.

The basic plot function will plot births, deaths, and mcpr over the entire simulation.

There are also pre-defined options that combine similar types of output. For instance, ‘apo’ stands for adverse pregnancy outcomes, and will plot infant deaths, stillbirths, and abortions.

plot() will take any of the following options:

  • ‘cpr’ will plot three different ways to define contraceptive prevalence - mCPR, CPR (includes traditional), and aCPR (includes traditional and restricts denominator to sexually active non-pregnant women)

  • ‘apo’ will plot adverse pregnancy outcomes, including abortion and miscarriage

  • ‘mortality’ will plot mortality-related outcomes

  • ‘method’ plots the method mix over time

[6]:
sim.plot() # the default
sim.plot('cpr')
../_images/tutorials_T1_intro_14_0.png
../_images/tutorials_T1_intro_14_1.png
[6]:
../_images/tutorials_T1_intro_14_2.png

In the next tutorial, you’ll learn how to build intervention scenarios using multiple sims and plot those.