In anthony-aylward/selex: Multinomial Logistic Regrestion for SELEX sequencing data

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

The selex package includes an example of SELEX count data.

library(selex)
example_counts

First, fit a multinomial logit regression model to the counts.

fit <- selex_multinom(example_counts, weights = c(8, 1, 2, 4, 8), ref = "c")
summary(fit)

Then, numerically compute p-values for the coefficients.

selex_pvals(fit)

anthony-aylward/selex documentation built on May 21, 2019, 1:42 a.m.