SURVEYSELECT Procedureについて。 いろんな使い道が考えられるが、複雑なサンプリングもできちゃう。 使いそうなのは、ケース・コントロール研究のマッチングとかかな。 単純にリサンプリングを使ったり、ブートストラップを使ったP値を求めるだけならMULTT…

NLMIXED Procedureで一般化線型モデルの最尤推定を行ってみよう。 今回はロジスティック回帰分析。 説明変数xのオッズ比が大体2になるようにダミーデータを発生させ、NLMIXED Procedureでロジスティック回帰分析を行うプログラムを以下に示した。 data dummy…

最近使わなくなってしまったNLMIXED Procedureについて少しだけ書いてみようと思う。 尤度 (対数尤度) を明示的に書けて、ヘンテコなモデル (部分的に変量効果が欠けていたり…) を想定しない場合には簡単に最尤推定して、$H_0: \theta_i = 0$ (パラメータ $\…

SAS Excel一括変換するマクロ (2) 昔(2006年頃)作ったマクロを改良してみました。 SASファイルがあるディレクトリ内のSASデータセットを一括でExcelファイルへと変換するマクロです。 SASファイルの実行パスを取得する方法がSAS Faqに掲載されていたので参考…

ファイル名を桁を揃えた連番にして、複数のグラフをファイル出力をする。 1.xxx, 11.xxx, 111.xxx, 2.xxxみたいな感じにファイル名ソートされるのが気に食わなかったので作った。 %macro graph_plot(data, num); %do i = 1 %to # goptions reset = a…