Simulation PK 1 cpt, perfusion

$PROB PK 1 cpt perfusion, simulation
$INPUT ID TIME DV AMT RATE
$DATA ..\data.csv IGNORE=#
$SUBROUTINE ADVAN1 TRANS2
$PK
CL =THETA(1) * EXP(ETA(1))
V =THETA(2) * EXP(ETA(2))

K=CL/V ; REPARAMETRISATION

$ERROR
IPRED=F
Y=F+F*ERR(1)
IRES=DV-IPRED

$THETA (0, 12.1) ; CL
$THETA (0, 38.7) ; VC

$OMEGA 0.13 ;OM_CL
$OMEGA 0.101 ;OM_VC

$SIGMA 0.08 ;SIM_P

$SIM (123456) ONLYSIM SUB=1000
$TABLE ID TIME DV MDV ONEHEADER NOPRINT FILE=base_sim.fit