Gladstone-Bioinformatics-Wo.../intro-statistics-experimental-design-hypothesis-testing/DoExperiments.R
2022-10-26 19:23:33 -07:00

16 lines
605 B
R

rm(list = ls())
setwd("~/Dropbox (Gladstone)/scripts/Bioinformatics-Workshops/intro-experimental-design")
##load assay functions
source("SImulateData.R")
##do experiment comparing expression of gene between E9.5 and E11.5
PerformAssayGeneExpression(NsamplesPerTime = 100)
##Two models for expression of gene at E9_5: Gamma and Normal
AssayExpression_E9_5_two_models(NsamplesPerTime = 100)
##Repeat the experiment Nsim times to estimate mean gene expression at E9_5 under the two models of
##data generating distribution
Repeat_AssayExpression_E9_5_two_models(NsamplesPerTime = 100, Nrepeat = 10000)