Statistical Methods for High Dimensional ODE Models with Applications to Dynamic Gene Regulatory Network Identification