Automatic model building for theophylline PK data