Skip to contents

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

Usage

stat_reverse_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 reverse hourglass test results

Details

The reverse hourglass test evaluates whether mid developmental stages show higher transcriptomic index values (indicating younger genes) compared to both early and late stages. This creates a reverse hourglass pattern where recently evolved genes dominate during mid-development. The test computes a score based on the minimum difference between mid vs. early and mid vs. late TXI values.

Examples

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