R code of a gene regulatory network simulation of hybrid vigour