「$SUPER」の編集履歴（バックアップ）一覧はこちら

「$SUPER」(2006/08/23 (水) 00:45:43) の最新版変更点

追加された行は緑色になります。

削除された行は赤色になります。

　経口投与後のデータを発生させて，静脈内投与後のモデルで解析するシミュレーションの例．

-$TABLE で作成したデータを次の $PROB の $DATA で読み込むテクニック
-$SUB での LIBRARY の利用
-$SIM で乱数の種を指定する方法，および，シミュレーションごとに異なる種を使う方法

　この二点に気を配ればよい．
　なお，途中で計算エラーが発生した場合，そこで全てが止まってしまう．悲しい．私は解決方法を知りません．

----

 $PROBLEM DUMMY SIMULATION
 $DATA TEMPLATE-1.CSV IGNORE=@
 $INPUT ID TIME EVID MDV CMT AMT DV
 $SUB ADVAN2 TRANS2 LIBRARY INFN=INFN.FOR
 $PK
  CL=THETA(1)*EXP(ETA(1))
  V =THETA(2)*EXP(ETA(2))
  KA=THETA(3)*EXP(ETA(3))
  SC=V
 $ERROR
  Y=F*EXP(ERR(1))
  IPRED=F
 $THETA
  (0 1)
  (0 10)
  (0 0.7)
 $OMEGA 0.04 0.04 0.25
 $SIGMA 0.01
 $SIM (2601) ONLYSIM
 $TABLE
  NOPRINT NOHEADER FILE=DATA1.TXT
  NOAPPEND
  ID TIME EVID MDV CMT AMT IPRED DV
 ;
 
 $SUPER SCOPE=2 ITERATIONS=2
 $PROBLEM SIMULATION
 $DATA DATA1.TXT (8E12.0) NOOPEN
 $INPUT ID TIME EVID MDV CMT AMT IPRD DV
 $PK
  CL=THETA(1)*EXP(ETA(1))
  V =THETA(2)*EXP(ETA(2))
  KA=THETA(3)*EXP(ETA(3))
  SC=V
  ;
  CMT2=2
  ;
 $ERROR
  Y=F*EXP(ERR(1))
  IPRED=F
 $THETA
  (0 1)
  (0 10)
  (0 0.7)
 $OMEGA 0.04 0.04 0.25
 $SIGMA 0.01
 $SIM (-1) ONLYSIM
 $TABLE
  NOPRINT NOHEADER FILE=DATA2.TXT
  NOAPPEND
  ID TIME EVID MDV CMT2 AMT IPRED DV
 ;
 
 $PROBLEM ESTIMATION
 $DATA DATA2.TXT (8E12.0) NOOPEN
 $INPUT ID TIME EVID MDV CMT AMT IPRD DV
 $THETA
  (0 1)
  (0 10)
  (0 FIXED)
 $OMEGA 0.04 0.04 (0 FIXED)
 $SIGMA 0.01
 $EST NOABORT MAXE=9999 METHOD=ZERO
 $COV
 $TABLE
  NOPRINT ONEHEADER FILE=DATA3.TXT
  ID TIME EVID MDV CMT AMT IPRED
 ;

　経口投与後のデータを発生させて，静脈内投与後のモデルで解析するシミュレーションの例．

-$TABLE で作成したデータを次の $PROB の $DATA で読み込むテクニック
-$SUB での LIBRARY の利用
-$SIM で乱数の種を指定する方法，および，シミュレーションごとに異なる種を使う方法

　この三点に気を配ればよい．
　なお，途中で計算エラーが発生した場合，そこで全てが止まってしまう．悲しい．私は解決方法を知りません．
----

　この手のシミュレーションには $SUPER を使おうと言っているわけではない．$SUPER の使い方を覚えるくらいなら，SAS か S-Plus/R を覚えましょう(^^)．

----

 $PROBLEM DUMMY SIMULATION
 $DATA TEMPLATE-1.CSV IGNORE=@
 $INPUT ID TIME EVID MDV CMT AMT DV
 $SUB ADVAN2 TRANS2 LIBRARY INFN=INFN.FOR
 $PK
  CL=THETA(1)*EXP(ETA(1))
  V =THETA(2)*EXP(ETA(2))
  KA=THETA(3)*EXP(ETA(3))
  SC=V
 $ERROR
  Y=F*EXP(ERR(1))
  IPRED=F
 $THETA
  (0 1)
  (0 10)
  (0 0.7)
 $OMEGA 0.04 0.04 0.25
 $SIGMA 0.01
 $SIM (2601) ONLYSIM
 $TABLE
  NOPRINT NOHEADER FILE=DATA1.TXT
  NOAPPEND
  ID TIME EVID MDV CMT AMT IPRED DV
 ;
 
 $SUPER SCOPE=2 ITERATIONS=2
 $PROBLEM SIMULATION
 $DATA DATA1.TXT (8E12.0) NOOPEN
 $INPUT ID TIME EVID MDV CMT AMT IPRD DV
 $PK
  CL=THETA(1)*EXP(ETA(1))
  V =THETA(2)*EXP(ETA(2))
  KA=THETA(3)*EXP(ETA(3))
  SC=V
  ;
  CMT2=2
  ;
 $ERROR
  Y=F*EXP(ERR(1))
  IPRED=F
 $THETA
  (0 1)
  (0 10)
  (0 0.7)
 $OMEGA 0.04 0.04 0.25
 $SIGMA 0.01
 $SIM (-1) ONLYSIM
 $TABLE
  NOPRINT NOHEADER FILE=DATA2.TXT
  NOAPPEND
  ID TIME EVID MDV CMT2 AMT IPRED DV
 ;
 
 $PROBLEM ESTIMATION
 $DATA DATA2.TXT (8E12.0) NOOPEN
 $INPUT ID TIME EVID MDV CMT AMT IPRD DV
 $THETA
  (0 1)
  (0 10)
  (0 FIXED)
 $OMEGA 0.04 0.04 (0 FIXED)
 $SIGMA 0.01
 $EST NOABORT MAXE=9999 METHOD=ZERO
 $COV
 $TABLE
  NOPRINT ONEHEADER FILE=DATA3.TXT
  ID TIME EVID MDV CMT AMT IPRED
 ;