Ctrl+K
nightly v0.1.2 v0.1.3 v0.1.4 v0.1.5 v0.1.6 v0.1.7
nightly v0.1.2 v0.1.3 v0.1.4 v0.1.5 v0.1.6 v0.1.7

Using Python#

Note

It is possible to run the whole analysis in python without touching the Dash GUI.

Once you prepared your data as described in the prepare your data section, you can use the python API to analyze the data. First you need to load your data in python. Note that the test data has log transformed intensities, which is specified via the logbase parameter in the RAPDORData construction

import pandas as pd
from RAPDOR.datastructures import RAPDORData
intensities = pd.read_csv("https://raw.githubusercontent.com/domonik/RAPDOR/main/RAPDOR/tests/testData/testFile.tsv", sep="\t")
design = pd.read_csv("https://raw.githubusercontent.com/domonik/RAPDOR/main/RAPDOR/tests/testData/testDesign.tsv", sep="\t")
rapdordata = RAPDORData(intensities, design, logbase=2)

Next you want to normalize the fractions and calculate samplewise distances. This is done via a single call. Here we use the recommended Jensen-Shannon-Distance and smooth the initial intensities with an averaging kernel of size 3.

rapdordata.normalize_and_get_distances(method="Jensen-Shannon-Distance", kernel=3)

Since we only have 3 replicates per condition for this dataset, it is not reasonable to calculate a p-value. Instead we will calculate stats like the Jensen-Shannon-Distance of the means and the ANOSIM R value. We will then rank the table based on the R value and the mean distance, which only takes a second for the test data. This is done via:

rapdordata.calc_all_scores()
rapdordata.rank_table(["ANOSIM R", "Mean Distance"], ascending=[False, False])

You can always check your calculated parameters via accessing the rapdordata.extra_df or the rapdordata.df fields. The df contains also the data columns and might be messy. For the test data you can print the head of the table:

print(rapdordata.extra_df.head())

The output should look as follows:

   id  RAPDORid  ...  position strongest shift  Rank
0   0      1118  ...                      10.0    12
1   1      1115  ...                       4.0     7
2   2      1114  ...                       2.0     1
3   3      1721  ...                       7.0     2
4   4      1103  ...                       6.0     8
[5 rows x 12 columns]

You can see that the protein with the RAPDORid 1114 is ranked on place 1.

Finally we want to export the file as a csv table and as a JSON file. The JSON is further needed to display your data in a Dash webserver interface. You can learn more about it in the server setup section.

rapdordata.export_csv(file="test.tsv", sep= "\t")
rapdordata.to_json(file="test.json")
Show Source