Return test statistics from the mgcv package.

getSmootherTestStats(models)

Arguments

models

the GAM models, typically the output from fitGAM. Note that this function only works when models is a list.

Value

a matrix with the wald statistics associated with each lineage's smoother. The matrix has one row per gene where the fitting procedure converged.

Examples

data(gamList, package = "tradeSeq") getSmootherPvalues(gamList)
#> s(t1):l1 s(t2):l2 #> 1700012B09Rik 1.0000000 1.0000000 #> 1700017B05Rik 0.1989893 0.6537023 #> 1700020L24Rik 1.0000000 1.0000000 #> 1810009A15Rik 0.9462776 0.4699110 #> 1810011H11Rik 1.0000000 1.0000000 #> 1810027O10Rik 0.5469890 0.5261559 #> 1810031K17Rik 0.1953402 0.3535604 #> 1810033B17Rik 1.0000000 1.0000000 #> 1810037I17Rik 0.4508156 0.4812300 #> 2010002N04Rik 1.0000000 1.0000000 #> 2010011I20Rik 0.9131269 0.8924305