http://clinical-pk-pd.cocolog-nifty.com/blog/ Pharmacometrics ja-JP 2009-12-29T23:06:05+09:00 http://clinical-pk-pd.cocolog-nifty.com/blog/2009/12/glimmix.html 今日はモデルについての復習。PROC GENMOD, GLIMMIX, NLMI... <p>今日はモデルについての復習。PROC GENMOD, GLIMMIX, NLMIXEDといったものを整理したPPT資料を作った。<br />PROC GLIMMIX data=infection2&nbsp; order=data OR;<br /> class trt clinic;<br /> model cure=trt / dist=bin link = logit solution cl;<br /> random intercept /subject=clinic;<br /> estimate &quot;treatment&quot; trt 1 -1 /exp;<br />RUN;<br />比較のためにNLMIXEDを使った。NLMIXEDは小さいデータなら実行速度が気にならない。QPOINTS=20 TECH=NEWRAPとか、重い設定と聞いていたのに。QPOINTSはGauss積分の分位点の数を決めるので当然多いと実行速度が遅くなるのだが。考えてみれば、Iterationの少ない数で収束できる問題については、気にならない速度になるのかもしれない。<br />あるWebで、「全てをNLMIXEDを使って解析してみよう。そうすれば各統計手法が何をしているのかがわかる」という記述があった。まぁ最近そうかなと感じる。というのもNLMIXEDは結局のところ、尤度を解析者がコードするようなものだし、カテゴリー変数のためのダミー変数も自分で作らないとだめだから、理解は深まるだろうなと。</p> SAS NONMEM user 2009-12-29T23:06:05+09:00 http://clinical-pk-pd.cocolog-nifty.com/blog/2009/12/2009-ead0.html 全て終了。明日から年末年始休みに入る。休みの間はちょっと頭を切り替えていろいろと... <p>全て終了。明日から年末年始休みに入る。休みの間はちょっと頭を切り替えていろいろと勉強したいこともあるのだが、明日はかなり寝ると思う。</p> 雑 NONMEM user 2009-12-28T23:00:46+09:00 http://clinical-pk-pd.cocolog-nifty.com/blog/2009/12/post-ea43.html 久しぶりの休日で特に頭を使わずに統計ソフトをちらりちらりいじっていた。正直大した... <p>久しぶりの休日で特に頭を使わずに統計ソフトをちらりちらりいじっていた。正直大したことをしていない。仕事も一段落し、年末年始は別のことの準備とか勉強とかにあてようかと考えている。来年はBUGSもやりたいのだけど、数理統計の勉強がまずは最優先。</p> 雑 NONMEM user 2009-12-26T19:33:16+09:00 http://clinical-pk-pd.cocolog-nifty.com/blog/2009/12/ess.html EmacsというGNUのテキストエディタから、Rを起動させることができる。そのた... <p>EmacsというGNUのテキストエディタから、Rを起動させることができる。そのためにはESS(Emacs Speaks Statistics)というどこかで聞いたような名前のパッケージを用いる。UNIX環境でなくてもインストールが楽なESSプレインストールドﾊﾞｰｼﾞｮﾝがあって、<a href="http://vgoulet.act.ulaval.ca/en/ressources/emacs/windows">http://vgoulet.act.ulaval.ca/en/ressources/emacs/windows</a>からダウンロードできる。</p> <p>まぁいまさらなのだが、R Booksを読んで使ってみたいかと思って。画面が分割できるようなのだが、いまさっきインストールしたからいまいちわからない。</p> R/S-PLUS NONMEM user 2009-12-25T20:15:56+09:00