adaptive_MC_sampling

Core implementation of the adaptive Monte Carlo sampling algorithm. This internal
function explores the parameter space by updating the sampling distribution
based on evaluated likelihoods. It is called by the main <code>run_adaptive_sampling</code>
function.

internal

A robust implementation of Topolow algorithm. It embeds objects into a low-dimensional Euclidean space from a matrix of pairwise dissimilarities, even when the data do not satisfy metric or Euclidean axioms. The package is particularly well-suited for sparse, incomplete, and censored (thresholded) datasets such as antigenic relationships. The core is a physics-inspired, gradient-free optimization framework that models objects as particles in a physical system, where observed dissimilarities define spring rest lengths and unobserved pairs exert repulsive forces. The package also provides functions specific to antigenic mapping to transform cross-reactivity and binding affinity measurements into accurate spatial representations in a phenotype space.
Key features include:
* Robust Embedding from Sparse Data: Effectively creates complete and consistent maps (in optimal dimensions) even with high proportions of missing data (e.g., >95%).
* Physics-Inspired Optimization: Models objects (e.g., antigens, landmarks) as particles connected by springs (for measured dissimilarities) and subject to repulsive forces (for missing dissimilarities), and simulates the physical system using laws of mechanics, reducing the need for complex gradient computations.
* Automatic Dimensionality Detection: Employs a likelihood-based approach to determine the optimal number of dimensions for the embedding/map, avoiding distortions common in methods with fixed low dimensions.
* Noise and Bias Reduction: Naturally mitigates experimental noise and bias through its network-based, error-dampening mechanism.
* Antigenic Velocity Calculation (for antigenic data): Introduces and quantifies "antigenic velocity," a vector that describes the rate and direction of antigenic drift for each pathogen isolate. This can help identify cluster transitions and potential lineage replacements.
* Broad Applicability: Analyzes data from various objects that their dissimilarity may be of interest, ranging from complex biological measurements such as continuous and relational phenotypes, antibody-antigen interactions, and protein folding to abstract concepts, such as customer perception of different brands.
Methods are described in the context of bioinformatics applications in Arhami and Rohani (2025a) <doi:10.1093/bioinformatics/btaf372>, and mathematical proofs and Euclidean embedding details are in Arhami and Rohani (2025b) <doi:10.48550/arXiv.2508.01733>.

Omid Arhami

topolow

Force-Directed Euclidean Embedding of Dissimilarity Data

adaptive_MC_sampling function

<dl><dt>samples_file</dt>
<dd>Path to the CSV file with samples for the current job.</dd>
<dt>dissimilarity_matrix</dt>
<dd>The dissimilarity matrix to be fitted.</dd>
<dt>iterations</dt>
<dd>Number of sampling iterations per job.</dd>
<dt>mapping_max_iter</dt>
<dd>Maximum optimization iterations for the embedding.</dd>
<dt>relative_epsilon</dt>
<dd>Convergence threshold for the optimization.</dd>
<dt>folds</dt>
<dd>Number of cross-validation folds.</dd>
<dt>num_cores</dt>
<dd>Number of cores for parallel processing.</dd>
<dt>scenario_name</dt>
<dd>Name for output files, used for context.</dd>
<dt>verbose</dt>
<dd>Logical. If TRUE, prints progress messages.</dd></dl>

Arguments

Perform Adaptive Monte Carlo Sampling (Internal) — adaptive_MC_sampling

<dl>

<dt>samples_file</dt>
<dd>Path to the CSV file with samples for the current job.</dd>


<dt>dissimilarity_matrix</dt>
<dd>The dissimilarity matrix to be fitted.</dd>


<dt>iterations</dt>
<dd>Number of sampling iterations per job.</dd>


<dt>mapping_max_iter</dt>
<dd>Maximum optimization iterations for the embedding.</dd>


<dt>relative_epsilon</dt>
<dd>Convergence threshold for the optimization.</dd>


<dt>folds</dt>
<dd>Number of cross-validation folds.</dd>


<dt>num_cores</dt>
<dd>Number of cores for parallel processing.</dd>


<dt>scenario_name</dt>
<dd>Name for output files, used for context.</dd>


<dt>verbose</dt>
<dd>Logical. If TRUE, prints progress messages.</dd>

</dl>

adaptive_MC_sampling: Perform Adaptive Monte Carlo Sampling (Internal)

Description

Usage

Value

Arguments