data(MS, package = "gemR")
gem <- GEM(proteins ~ MS * group, data = MS)
# Extract interaction between 'MS' and 'group
ER.isolated <- knock.in(gem, 'MS:group')
# Remove main effect of 'group'
ER.cleaned <- knock.out(gem, 'group')
Run the code above in your browser using DataLab