computes the phylogenetically based transcriptome evolutionary index (TEI)
Examples
## load example sequence data
data("PhyloExpressionSetExample", package="myTAI")
spmat <- as(data.matrix(PhyloExpressionSetExample[,-c(1,2)]), "sparseMatrix")
rownames(spmat) <- PhyloExpressionSetExample$GeneID
ps <- setNames(PhyloExpressionSetExample$Phylostratum, PhyloExpressionSetExample$GeneID)
rcpp_tei_parallel(spmat, ps)
#> $sumx
#> [1] 60413160 59665009 57246260 55866953 55152111 57140339 63197410
#>
#> $teisum
#> [1] 195130989 192456282 177871849 174120165 169537228 181506933 214260302
#>
#> $tei
#> [1] 3.229942 3.225614 3.107135 3.116693 3.073993 3.176511 3.390334
#>