Parameter estimation in systems of differential equations based on three-way data arrays and random time changes