Limiting a linear model to 3rd level interactions in R

I have a dataset with 14 binary variables. I've already tested for significant single variables, but I'd like to also check for significant interactions. However, I know that higher level interactions are unlikely to be significant and just muddle the model. Is there anyway to run a linear model in R, but tell it to only test for interaction between a maximum of 3 variables?