Skip to content

Instantly share code, notes, and snippets.

@michaelbarton
Created June 19, 2017 20:25
Show Gist options
  • Save michaelbarton/101960391430635fd4ff24c99451ae55 to your computer and use it in GitHub Desktop.
Save michaelbarton/101960391430635fd4ff24c99451ae55 to your computer and use it in GitHub Desktop.
```{r}
library(lme4)
null <- glm.nb(gene_l1_norm ~ 1, data = data)
genome_null <- glm.nb(gene_l1_norm ~ biological_source_name, data = data)
n_misassemblies <- glm.nb(gene_l1_norm ~ n_misassemblies, data = data)
n_misassemblies_genome <- glmer.nb(gene_l1_norm ~ n_misassemblies + (1|biological_source_name), data = data)
nas50 <- glm.nb(gene_l1_norm ~ nas50, data = data)
nas50_genome <- glmer.nb(gene_l1_norm ~ nas50 + (1|biological_source_name), data = data)
nas50_n_misassemblies <- glm.nb(gene_l1_norm ~ nas50 + n_misassemblies, data = data)
nas50_n_misassemblies_genome <- glmer.nb(gene_l1_norm ~ nas50 + n_misassemblies + (1|biological_source_name), data = data)
nas50_n_misassemblies_int <- glm.nb(gene_l1_norm ~ nas50*n_misassemblies, data = data)
nas50_n_misassemblies_int_genome <- glmer.nb(gene_l1_norm ~ nas50*n_misassemblies + (1|biological_source_name), data = data)
BIC(
null,
genome_null,
n_misassemblies,
n_misassemblies_genome,
nas50,
nas50_genome,
nas50_n_misassemblies,
nas50_n_misassemblies_genome,
nas50_n_misassemblies_int,
nas50_n_misassemblies_int_genome)
```
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment