$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