Skip to contents

Test for reductive hourglass patterns in transcriptomic data by comparing early and late developmental stages to mid developmental stages.

Usage

stat_reductive_hourglass_test(phyex_set, modules, ...)

Arguments

phyex_set

A PhyloExpressionSet object

modules

A named list with elements 'early', 'mid', and 'late' containing stage indices for each developmental module

...

Additional arguments passed to stat_generic_conservation_test

Value

A ConservationTestResult object with reductive hourglass test results

Details

The reductive hourglass test evaluates whether mid developmental stages show lower transcriptomic index values (indicating older genes) compared to both early and late stages. This creates an hourglass-shaped pattern where ancient genes dominate during mid-development. The test computes a score based on the minimum difference between early vs. mid and late vs. mid TXI values.

Examples

# Define developmental modules
# modules <- list(early = 1:3, mid = 4:6, late = 7:9)
# result <- stat_reductive_hourglass_test(phyex_set, modules)