Signaturization

The main feature of the CC is the automatic conversion of virtually any compound-related data into a standard format ready-to-be-used by machine learning algorithms. All CC data is stored as HDF5 files following a folder structure organized by datasets.

The central type of data are the signatures (one numerical vector per molecule), which are of four types:

• sign0 Signatures type 0: A sufficiently-processed version of the raw data, containing TF-IDF weightnings if applicable. They usually show explicit knowledge, which enables connectivity and interpretation.
• sign1 Signatures type 1: A PCA/LSI-projected version of the data, retaining 90% of the variance. They keep most of the complexity of the original data and they can be used for almost-exact similarity calculations.
• sign2 Signatures type 2: Network embedding of the similarity matrix derived from signatures type 1. They have fixed length, which is convenient for machine learning, and capture both explicit and implicit similarity relationships in the data.
• sign3 Signatures type 3: Network embedding of observed and inferred similarity networks. Their added value, compared to signatures type 2, is that they can be derived for virtually any molecule in any dataset. ⚠ These signatures are not calculated yet, and won't be in the near future.

Besides, there are other (auxiliary) types of data that may be of interest. The asterisk * denotes correspondence with signatures type 0-3. Typically, we have the 1 version of the following in the repository.

• sims* Similarity vectors: Full similarities stored as light int8 data. Each molecule receives one such similarity vector per dataset. They may be observed *_obs or predicted *_prd similarities. These signatures are only applicable to exemplary datasets.
• neig* Nearest neighbors: XXXX. Currently, we consider the 1000-nearest neighbors, which is more than sufficient in any realistic scenario. For now, we only keep neig1.
• nprd* Predicted nearest neighbors: XXXX
• clus* Clusters: XXXX
• proj* 2D Projections: XXX

I consider the numbering 0-3 to be conceptually closed. However, further auxiliary data types may be introduced in the future. Note that all names have a 4-character code followed by a digit. Future data should stick to this nomenclature.

Commonalities

The following applies to all data types except the similarity vectors sims*, which are of a very different nature due to their organization in the folder structure.

I suggest that

Signatures type 0

Signatures type 1

Signatures type 2

Signatures type 3

Similarity vectors

Nearest neighbors

Predicted nearest neighbors

Clusters

2D projections

Below, I list an schematic proposal of the classes:

class 

As you know, we have most of the CC data stored as HDF5 files. I think HDF5 format is good and we have to stick to this file format. However, I think that these files should be accessible through some classes. These classes must not load all data into memory.

Signature classes

• Applies to:
• Signatures Type 0
• Signatures Type 1
• Signatures Type 2
• Signatures Type 3
• 2D projections
• Every class must have at least the following attributes:
• V: the values
• keys: the keys, sorted alphabetically
• metric: the distance used
• pvalues: (distance, p-value) array
• PATH: the path where everything is stored
• Every class must have at least the following methods:
• __iter__: smart, batch iteration, if necessary
• __getattr__: returns the vector corresponding to the key. Works fast with bisect, but should return None if the key is not in keys (ideally, keep a set to do this).
• fit: Not sure this is necessary... Maybe we can just do it as part of the pipeline.
• predict: For the new samples, we should be able to produce the corresponding V vectors. This will be, by far, the most tricky part. One should access the models folder and use them correspondingly. To increase speed in this part, probably one should just predict for the ones that are not already in the reference. Sometimes, it will be necessary to learn a mapping functions, for instance via AdaNet; for example, in the case of Signature Type 2, as node2vec does not allow for out-of-sample mapping.
• validate: I'm thinking of a folder where we have validation files (for now, MoA and ATC), and then automatically outputting AUROC and KS metrics, among others.
• background: Not sure this is necessary... Just like fit, this has to be done only at initiation.

Other classes

• We have other data types, such as the nearest neighbors produced by Oriol and the clusters produced by myself.
• These must also have at least the following methods:
• __iter__
• __getattr__
• predict: As always, we want to be able to predict for new molecules using the models stored.
• validate: Here we will not use AUROC and KS, but other statistics, depending on the case.

Here I put a scheme of the first part of the Chemical Checker pipeline: