sleuth_wt: Wald test for a sleuth model

Description Usage Arguments Value See Also Examples

Description

This function computes the Wald test on one specific 'beta' coefficient on every transcript.

Usage

1
sleuth_wt(obj, which_beta, which_model = "full")

Arguments

obj

a sleuth object

which_beta

a character string of denoting which grouping to test. For example, if you have a model fit to 'treatment,' with values of neg_ctl, pos_ctl, and drug, you would need to run sleuth_wt once each for pos_ctl and drug

which_model

a character string of length one denoting which model to use

Value

an updated sleuth object

See Also

models to view which models have been fit and which coefficients can be tested, sleuth_results to get back a data.frame of the results

Examples

1
2
3
4
5
6
7
8
9
# Assume we have a sleuth object with a model fit to both genotype and drug,
models(so)
# formula:  ~genotype + drug
# coefficients:
#   (Intercept)
#   genotypeKO
#   drugDMSO
so <- sleuth_wt(so, 'genotypeKO')
so <- sleuth_wt(so, 'drugDMSO')

pachterlab/zika documentation built on May 24, 2019, 5:58 p.m.