diff --git a/README.md b/README.md
index 2012f813c..a24038fa3 100644
--- a/README.md
+++ b/README.md
@@ -1,6 +1,6 @@
 # A dynamic benchmark for gene regulatory network (GRN) inference
 
-[![Documentation Status](https://readthedocs.org/projects/grn-inference-benchmarking/badge/?version=latest)](https://grn-inference-benchmarking.readthedocs.io/en/latest/?badge=latest)
+The full documentation is hosted on [ReadTheDocs](https://openproblems-grn-task.readthedocs.io/en/latest/index.html). [![Documentation Status](https://readthedocs.org/projects/grn-inference-benchmarking/badge/?version=latest)](https://grn-inference-benchmarking.readthedocs.io/en/latest/?badge=latest)
 
 <!--
 This file is automatically generated from the tasks's api/*.yaml files.
diff --git a/docs/source/add_stuff.rst b/docs/source/add_stuff.rst
deleted file mode 100644
index 61f33c9b4..000000000
--- a/docs/source/add_stuff.rst
+++ /dev/null
@@ -1,4 +0,0 @@
-Adding datasets or algorithms
-===
-
-The general documentation that includes instructions on how to add new methods for the task, as well as new datasets is hosted `here <https://openproblems.bio/documentation/>`_.
diff --git a/docs/source/evaluation.rst b/docs/source/evaluation.rst
index 1dd40abc3..f202b0437 100644
--- a/docs/source/evaluation.rst
+++ b/docs/source/evaluation.rst
@@ -1,16 +1,30 @@
 Evaluation
-=====
+==========
 
 The evaluation is done with the help of pertubation data, using two different approaches:
-1. Regression from GRN regulations to target expression
-2. Regression from TF expression of predicted regulators to target expression
 
-Regression from GRN regulations to target expression
-----------------
+#. Regression from GRN regulations to target expression
+#. Regression from TF expression of predicted regulators to target expression
 
+|
 
-Regression from TF expression of predicted regulators to target expression
-----------------
+.. image:: images/regressions.png
+   :width: 100 %
+   :alt: overview of the two regression evaluation approaches
+   :align: center
+
+|
+|
+
+
+Evaluation 1: Regression from GRN regulations to target expression
+------------------------------------------------------------------
+The first approach we used is similar to GRaNPA and the multivariate decision tree in Decoupler, where regulatory weights from the GRN form the feature space to predict perturbation data. In this method, we train one model per sample. The feature space matrix has dimensions of genes by transcription factors (TFs), with values being the regulatory weights from the GRN or 0 if the link is absent. The target space matrix represents the perturbation data for each sample. We evaluate the model's predictive performance using a 5-fold cross-validation scheme and the coefficient of determination (R²) as the metric. LightGBM is used for computational efficiency.
+
+
+Evaluation 2: Regression from TF expression of predicted regulators to target expression
+----------------------------------------------------------------------------------------
+In the second approach, instead of using regulatory weights, we utilized the expression of putative regulators (TFs) from the perturbation data to construct the feature space. We fit one model per gene, selecting regulators based on the regulatory weights suggested by the GRNs. This method is similar to many modern GRN inference techniques.
 
 
 
diff --git a/docs/source/extending.rst b/docs/source/extending.rst
new file mode 100644
index 000000000..ec351153d
--- /dev/null
+++ b/docs/source/extending.rst
@@ -0,0 +1,21 @@
+Extending the pipeline
+======================
+
+Currently the perturbation dataset from the Open Problems: Single Cell Perturbation 2023 data science competition is used.
+It provides single-cell perturbation gene expression for peripheral blood mononuclear cells (PBMCs), along with integrated multi-omics data of scRNA-seq and scATAC-seq for the baseline compound from the same experiment.
+It includes 146 perturbations, making it the largest drug perturbation study on primary human tissue with donor replicates. 
+
+Currently, the following six enhancer aware GRN inference methods (eGRN methods) are implemented in the pipeline:
+
+#. Scenic+ (`Paper <https://doi.org/10.1038/s41592-023-01938-4>`_)
+#. CellOracle (`Paper <https://doi.org/10.1038/s41586-022-05688-9>`_)
+#. FigR (`Paper <https://doi.org/10.1016/j.xgen.2022.100166>`_)
+#. scGLUE (`Paper <https://doi.org/10.1038/s41587-022-01284-4>`_)
+#. GRaNIE (`Paper <https://doi.org/10.15252/msb.202311627>`_)
+#. ANANSE (`Paper <https://doi.org/10.1093/nar/gkab598>`_)
+
+To add a method to the repository, follow the instructions in the ``scripts/add_a_method.sh`` script.
+
+
+
+
diff --git a/docs/source/images/overview.png b/docs/source/images/overview.png
new file mode 100644
index 000000000..3a8d9bd6b
Binary files /dev/null and b/docs/source/images/overview.png differ
diff --git a/docs/source/images/regressions.png b/docs/source/images/regressions.png
new file mode 100644
index 000000000..5ed729bd8
Binary files /dev/null and b/docs/source/images/regressions.png differ
diff --git a/docs/source/index.rst b/docs/source/index.rst
index cd763b3c8..b46b8fa81 100644
--- a/docs/source/index.rst
+++ b/docs/source/index.rst
@@ -1,14 +1,30 @@
-Welcome to the documentation for the GRN inference task on OpenProblems Benchmarking platform!
-===================================
+Welcome to the documentation for the GRN inference task on the OpenProblems Benchmarking platform!
+==================================================================================================
 
 This task is one of many other tasks hosted on the `OpenProblems benchmarking platform <https://openproblems.bio/>`_.
 The interactive benchmarking results are hosted `here <https://openproblems.bio/results/>`_.
 
-The *GRN inference* task focuses on the inference of gene regulatory networks (GRN) from RNA-Seq expression or chromatin accessibility data (ATAC-Seq) or both. 
+The **GRN inference** task focuses on the inference of gene regulatory networks (GRN) from RNA-Seq expression or chromatin accessibility data (ATAC-Seq) or both. 
 The pipeline evaluates inferred GRNs against pertubation data, by training two types of regression models. This type of evaluation is closer to evaluating the biological knowledge that a GRN should represent, instead of evaluating the presence of edges in a statistical way only, as commonly done by using metrics, such as AUPRC or AUROC.
 
-Jump to the :doc:`evaluation` section to get a deeper explanation on how the benchmarking task is setup.
-If you want to add your own datasets or algorithms to the benchmark, check our the :doc:`add_stuff` section.
+Jump to the :doc:`overview` section to get a first summary of the pipeline.
+If you want to add your own datasets or algorithms to the benchmark, check our the :doc:`extending` section.
+
+
+.. list-table:: Authors & contributors
+   :widths: 25 25
+   :header-rows: 1
+
+   * - name
+     - roles
+   * - Jalil Nourisa
+     - author
+   * - Antoine Passemiers
+     - author
+   * - Robrecht Cannoodt
+     - author
+   * - Marco Stock
+     - author
 
 .. note::
 
@@ -19,5 +35,8 @@ Contents
 
 .. toctree::
 
+   overview
+   objects
    evaluation
-   add_stuff
\ No newline at end of file
+   extending
+   
\ No newline at end of file
diff --git a/docs/source/objects.rst b/docs/source/objects.rst
new file mode 100644
index 000000000..031ae6065
--- /dev/null
+++ b/docs/source/objects.rst
@@ -0,0 +1,327 @@
+Objects
+=======
+
+This provides an overview of the objects used in the GRN inference task.
+
+.. toctree::
+
+   objects
+
+Perturbation
+------------
+
+Perturbation dataset for benchmarking.
+
+Example file: ``resources/grn-benchmark/perturbation_data.h5ad``
+
+Format:
+
+.. list-table:: AnnData object
+   :header-rows: 1
+
+   * - Attribute
+     - Description
+   * - ``obs``
+     - 'cell_type', 'sm_name', 'donor_id', 'plate_name', 'row', 'well', 'cell_count'
+   * - ``layers``
+     - 'n_counts', 'pearson', 'lognorm'
+
+Slot description:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Slot
+     - Type
+     - Description
+   * - ``obs["cell_type"]``
+     - ``string``
+     - The annotated cell type of each cell based on RNA expression.
+   * - ``obs["sm_name"]``
+     - ``string``
+     - The primary name for the (parent) compound (in a standardized representation) as chosen by LINCS. This is provided to map the data in this experiment to the LINCS Connectivity Map data.
+   * - ``obs["donor_id"]``
+     - ``string``
+     - Donor id.
+   * - ``obs["plate_name"]``
+     - ``string``
+     - Plate name 6 levels.
+   * - ``obs["row"]``
+     - ``string``
+     - Row name on the plate.
+   * - ``obs["well"]``
+     - ``string``
+     - Well name on the plate.
+   * - ``obs["cell_count"]``
+     - ``string``
+     - Number of single cells pseudobulked.
+   * - ``layers["n_counts"]``
+     - ``double``
+     - Pseudobulked values using mean approach.
+   * - ``layers["pearson"]``
+     - ``double``
+     - Normalized values using pearson residuals.
+   * - ``layers["lognorm"]``
+     - ``double``
+     - Normalized values using shifted logarithm.
+
+Control Method
+--------------
+
+Path:
+``src/control_methods``
+
+A control method.
+
+Arguments:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Name
+     - Type
+     - Description
+   * - ``--perturbation_data``
+     - ``file``
+     - Perturbation dataset for benchmarking.
+   * - ``--layer``
+     - ``string``
+     - Which layer of perturbation data to use to find tf-gene relationships. Default: ``lognorm``.
+   * - ``--prior_data``
+     - ``file``
+     - Prior data used for GRN benchmark.
+   * - ``--prediction``
+     - ``file``
+     - (*Output*) GRN prediction.
+
+Label
+-----
+
+Path:
+``src/metrics``
+
+A metric to evaluate the performance of the inferred GRN
+
+Arguments:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Name
+     - Type
+     - Description
+   * - ``--perturbation_data``
+     - ``file``
+     - Perturbation dataset for benchmarking.
+   * - ``--prediction``
+     - ``file``
+     - GRN prediction.
+   * - ``--score``
+     - ``file``
+     - (*Output*) File indicating the score of a metric.
+   * - ``--reg_type``
+     - ``string``
+     - (*Optional*) name of regression to use. Default: ``ridge``.
+   * - ``--subsample``
+     - ``integer``
+     - (*Optional*) number of samples randomly drawn from perturbation data. Default: ``-1``.
+
+
+GRN
+---
+
+GRN prediction
+
+Example file: ``resources/grn-benchmark/collectri.csv``
+
+Format:
+
+.. list-table:: Tabular data
+   :header-rows: 1
+
+   * - Attribute
+     - Description
+   * - ``Columns``
+     - 'source', 'target', 'weight'
+
+Slot description:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Column
+     - Type
+     - Description
+   * - ``source``
+     - ``string``
+     - Source of regulation.
+   * - ``target``
+     - ``string``
+     - Target of regulation.
+   * - ``weight``
+     - ``float``
+     - Weight of regulation.
+
+
+Score
+-----
+
+File indicating the score of a metric.
+
+Example file: ``resources/grn-benchmark/score.csv``
+
+Format:
+
+.. list-table:: Tabular data
+   :header-rows: 1
+
+   * - Attribute
+     - Description
+   * - ``Columns``
+     - 'accuracy', 'completeness'
+
+Slot description:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Column
+     - Type
+     - Description
+   * - ``accuracy``
+     - ``string``
+     - (*Optional*) some explanation.
+   * - ``completeness``
+     - ``double``
+     - (*Optional*) some explanation.
+
+
+Prior data
+----------
+
+Prior data used for grn benchmark
+
+Example file: ``resources/grn-benchmark/prior_data.h5ad``
+
+Format:
+
+.. list-table:: AnnData object
+   :header-rows: 1
+
+   * - Attribute
+     - Description
+   * - ``uns``
+     - 'tf_list'
+
+Slot description:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Slot
+     - Type
+     - Description
+   * - ``uns["tf_list"]``
+     - ``list``
+     - List of known TFs obtained from https://resources.aertslab.org/cistarget.
+
+
+Multiomics RNA
+--------------
+
+RNA expression for multiomics data.
+
+Example file: ``resources/grn-benchmark/multiomics_rna.h5ad``
+
+Format:
+
+.. list-table:: AnnData object
+   :header-rows: 1
+
+   * - Attribute
+     - Description
+   * - ``obs``
+     - 'cell_type', 'donor_id'
+
+Slot description:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Slot
+     - Type
+     - Description
+   * - ``obs["cell_type"]``
+     - ``string``
+     - The annotated cell type of each cell based on RNA expression.
+   * - ``obs["donor_id"]``
+     - ``string``
+     - Donor id.
+
+Method
+------
+
+Path:
+``src/methods``
+
+A GRN inference method
+
+Arguments:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Name
+     - Type
+     - Description
+   * - ``--multiomics_rna``
+     - ``file``
+     - RNA expression for multiomics data.
+   * - ``--multiomics_atac``
+     - ``file``
+     - (*Optional*) Peak data for multiomics data.
+   * - ``--prediction``
+     - ``file``
+     - (*Output*) GRN prediction.
+
+
+Multiomics ATAC
+---------------
+
+Peak data for multiomics data.
+
+Example file: ``resources/grn-benchmark/multiomics_atac.h5ad``
+
+Format:
+
+.. list-table:: AnnData object
+   :header-rows: 1
+
+   * - Attribute
+     - Description
+   * - ``obs``
+     - 'cell_type', 'donor_id'
+
+Slot description:
+
+.. list-table:: 
+   :header-rows: 1
+
+   * - Slot
+     - Type
+     - Description
+   * - ``obs["cell_type"]``
+     - ``string``
+     - The annotated cell type of each cell based on RNA expression.
+   * - ``obs["donor_id"]``
+     - ``string``
+     - Donor id.
+
+
+
+
+
+
+
+
+
diff --git a/docs/source/overview.rst b/docs/source/overview.rst
new file mode 100644
index 000000000..7928b39a4
--- /dev/null
+++ b/docs/source/overview.rst
@@ -0,0 +1,27 @@
+Overview
+========
+
+The pipeline evaluates inferred GRNs against perturbation data. The evaluation is done by training two types of regression models. 
+This type of evaluation is closer to evaluating the biological knowledge that a GRN should represent, instead of evaluating the presence of edges in a statistical way only, as commonly done by using metrics, such as AUPRC or AUROC. 
+The general setup is shown in the figure below:
+
+.. image:: images/overview.png
+   :width: 100 %
+   :alt: overview of pipeline
+   :align: center
+
+The pipeline can evaluate algorithms that leverage only one of the multi-omic data types (RNA-Seq or ATAC-Seq) or both.
+It also evaluates the performance of two controls:
+
+#. As a *negative control*, the pipeline evaluates the performance of a random network.
+#. As a *positive control*, the pipeline evaluates the performance of a network derived from correlation of genes in the perturbation dataset used for evaluation.
+
+The two types of regression models are:
+
+#. Regression from GRN regulations to target expression
+#. Regression from TF expression of predicted regulators to target expression
+
+More details can be found in section :doc:`evaluation`.
+
+In the future, other prior knowledge can be incorporated in the pipeline to allow for evaluation with binary classification metrics and as an additional control method. 
+