SASグラフを今の自分が出来る限界までいじってみた。 dm 'log; clear; output; clear'; proc datasets lib = work kill; run; options linesize=120 pagesize = 300 mprint; %let execpath = " "; %let Path = " "; %macro setexecpath; %let execpath = %sy…

同じネタで 5 回も引っ張ってしまった。 何はともあれ、シミュレーション用のプログラムができました。 dm 'log; clear; output; clear'; proc datasets kill; run; option linesize = 130 pagesize = 500; %let execpath = " "; %let Path = " "; %macro se…

問題点 血中濃度のデータを想定しているくせに、前の時点よりも濃度が高くなってしまうシミュレーションデータが生成されてしまう 根本的な解決になるかどうかは分からないが、シミュレーションデータを作ってみた。 dm 'log; clear; output; clear'; proc d…

今回は点滴静注投与の場合を解析してみます。 コンパートメントモデルの部分を変更しただけなので、モデル式さえ書ければいろいろな場合を解析できますね。 dm 'log; clear; output; clear'; proc datasets kill; run; option linesize = 130 pagesize = 500…

今回は誤差分布に対数正規分布を仮定したモデルによる解析を行ってみようと思います。依然シミュレーションデータの問題は時間の都合上修正していませんが、アイデアは浮かんだので近いうちに直す予定です。 シミュレーションデータは誤差分布に対数正規分布…

emf ファイルは Power Point 等で、ベクトルデータのオブジェクトに変換することが出来ます。 SAS で描いた emf ファイルのグラフを使って、きれいなベクトルデータの曲線なんかを使いたいときにたまに使っています。 Univariate Procedure を使ってヒストグ…

ノンパラメトリック回帰は面白いなぁ。 SAS の Loess Procedure を使って解析してみたので、ついでに勉強しておきたい。 しかし、大学には教科書が無く、授業で取り扱うことが無かった。 今回は、竹澤邦夫先生のHP [1] を見てなにやらすごいことが分かり、そ…

分散安定化変換は、確率変数 $Y$ の期待値が $\mu$ で分散が $g(\mu)$、つまり $Y$ の分散が期待値の関数であるとし、$h(Y)$ の分散 $\mathrm{Var}\{h(Y)\}$ が漸近的に一定の値になるような変換 $h$ を考える方法です。 分散安定化変換 デルタ法を用いると…

Harvard School of Public Health PROGRAM ON CAUSAL INFERENCE in Epidemiology and Allied Sciences Harvard School of Public HelthのEpidemiologyとCausal InferenceのHP。 Publicationから、Robins先生の論文がかなり沢山ダウンロードできるようになっ…