# NOT RUN {
data(PhyloExpressionSetExample)
# compute partial TAI values for each Phylostratum
partialStrata <- pStrata(PhyloExpressionSetExample)
# show that colSums of pStrata is equavalent to the TAI values
all.equal(colSums(partialStrata),TAI(PhyloExpressionSetExample))
# show that colSums of pStrata is equavalent to colSums of pMatrix(PhyloExpressionSetExample)
all.equal(colSums(partialStrata),colSums(pMatrix(PhyloExpressionSetExample)))
# }
Run the code above in your browser using DataLab