メモ 船事例データ 関連データ名
■ mxは d という名のmxである MH hunemh MH船事例 RD hurdという、cRD 計算;観光船 1+26因子 名義因子名ごと hurd <- NULL # for(i in 1:27) { # hur<- sum( d[,1]*d[,i]) ... 続きをみる
解析疫学、リスクにまつわるメモや計算
「推定」のまわりをさぐる.率の重ね合わせをどう計算するか、因子の性質をどう分けるか、教科書では「解析はMHにより行う、因子が多ければ重回帰を用いる」という風で詳しい例は少ない.独自(のつもり)な思いつきで具体に試行.
数理を用いるべきアセスメントにも切り込む.
■ mxは d という名のmxである MH hunemh MH船事例 RD hurdという、cRD 計算;観光船 1+26因子 名義因子名ごと hurd <- NULL # for(i in 1:27) { # hur<- sum( d[,1]*d[,i]) ... 続きをみる
・ osw 事例データは因子重複時の効果が生起、抑制とも追加的でない viceが発生に支配的なせいだろう ・生起性因子をすべて仕分ければ、他の因子の生起性はBGとなる 生起性因子を挙げ切れる事例 ・・ これらのことから、モデルが単純になる →有効な層は、v単独と... 続きをみる
・小nな、性質の同じ因子をまとめて1因子とするアイデア ・データをほぼ説明しつくしたモデルの改善 ■【モデル】 ・watまたはmilk がv0cakeを抑制するかを調べる これら2抑制因子は n小のため、あまり取り上げていない 先に作っていたモデルで cake単独の生起性に対する作用を試す ... 続きをみる
■ 結果 ・因子の影響関係を調べるためにモデリングや層化を繰り返し、有効な因子を絞る試みを経て、分割することを考えた ・曝露ベクトルでR記述するのは便利だった ・得られた1つのモデルは 有効とした因子が、生起因子2つ:独立としたもの ただし従な生起因子と主な生起因子との重複部分を主... 続きをみる
・生起因子曝露ベクトルを”抑制されたもの”と、”されないもの”に分けるアイデアは改善に役立つか ■ 因子の分割 1つのa因子をb因子曝露したa因子と、しないa因子に分ける それぞれ独立となる 曝露ベクトルについて処理するのがポイントになる b因子... 続きをみる
■ 独立性 ~発生度数以前に ・生起、抑制因子に対する曝露から考えて効果のある因子との独立性が薄いものをあぶりだす ・単独 n<10 因子はどっちみちはっきりした効果を言い難くなるので最初に切り捨てておこうとするもの 1 vice 重複 vice単独 i単独 ... 続きをみる