compare_models: Perform likelihood ratio test

Description Usage Arguments Value Examples

View source: R/embeddr.R

Description

Perform likelihood ratio test

Usage

1
compare_models(model, null_model)

Arguments

model

The full model y ~ pseudotime returned by AER::tobit

null_model

The null model y ~ 1 returned by AER::tobit

Value

The p-value

Examples

1
2
3
4
5
6
7
library(scater)
data('sc_example_counts') ; sce <- newSCESet(countData = sc_example_counts)
sce <- embeddr(sce)
sce <- fit_pseudotime(sce)
full_model <- fit_pseudotime_model(sce, 1) # fit for first gene
null_model <- fit_null_model(sce, 1) # fit for first gene
pval <- compare_models(full_model, null_model)

kieranrcampbell/embeddr documentation built on May 20, 2019, 9:24 a.m.