morの推定勉強

小n影響が少ない ”差率”による colSums計算～層限定なし

■ colSumsにより、計算が簡単になる 非曝露試行は、効果を調べる有益な手段かもしれない 式を工夫する r*；曝露発生率差とは違った式を思いついた 差率とよぶ 差率は、cr発生率、RDほどは値が暴れない ‐‐全体；限定なし-- osw事例 yall<- s...

2025/06/21 09:42

差率のみかけを調べ因子を分類する

■ 差率は全体の率に及ぼす因子の影響が際立っているように見えた 差率にもみかけがありうるだろう ■ 独自の生起性 ・vice以外の生起因子がどれなのか cakeが独自生起因子であって他はみかけと考えたのだった colSums(osd[,1]*(1-osd[,13])...

2025/06/14 22:14

ﾒﾓ 非曝露試行と重複を考えあわせ：生起因子以外の”生起因子”を仕分ける vice RD；率差 でなく 発生数差/分母；率

□ 生起因子以外にy1の因子が複数あった 複数の因子はRD colSums(osd[,1]*osd[] )/colSums( osd[]) - colSums( osd[,1]*0^osd[,])/colSums(0^osd[] ) y ...

2025/06/08 22:11

r1* 曝露発生率差” で 非曝露試行し 因子分類

〇 各因子非曝露とする、あらたなrisk値を使って調べる risk値：r1* は、j 因子なしの発生率-j含む発生率で、0を基準に分類できる それから、N個×N行のmatrixを作り、vecで条件付けし、表示し観察する ☆ 記事 ” 非曝露で 発生率mxから,...

2025/05/31 15:08

生起因子らしいか.-2 非FBDを模擬 colSumsで2×2表を記述

■ 過去記事 ”生起因子らしくもない？ 感染症など外の要因” ■ colSumsによる2×2表 acc <- colSums( osd[,1]*osd[,]) bcc<- colSums( osd[,1]*0^osd[,]) ccc <- colSums( 0^osd...

2025/05/22 22:43

メモ 1因子を非曝露として調べる；”非曝露試行”

■ 各因子について”非曝露”としたときの各因子のrisk値を計算しベクトルをつないだmxから調べるということをした 因子の曝露から調べるのとは違い、ある1つが非曝露になることで、変化する他の因子たちのrisk値を調べ、非曝露とした因子の性質をみる方法 ■ 非曝露のパターン...

2025/05/20 20:57

曝露発生率mxを 非限定層1因子非曝露で・・ colSums描画 ＆ cake chocoを読み分ける

〇 層を限定せず、因子性質を調べる r1 で限定のない層をみる すべての因子を1回づつ非曝露にしてr1を計算し、線でつなぐ ベクトルで描画を条件付けし、表示因子を選び、並び順を変える r1 のベクトルから mx kもついでに計算しておく 〇 2つの有力な副生起因子ca...

2025/05/16 22:37

曝露発生率 mxとし、選択・順を変えてplotする記述

■ 曝露群での発生率 r1に対して及ぼされる、因子効果を名目因子について調べる ・各因子削減したときの率mxをつくり、表示の際に選択・順を変える；ベクトルで指定 ・生起因子曝露の下で調べる例 ■ 準備 # k1の小n はじくため k1ベクトルをmxにする記述 ...

2025/05/09 19:27

メモ colSums によるrisk値の式

■risk値・度数 ・RDの表現 colSums(hune[,1]* hune[]) /colSums( hune[] ) # r1 colSums(hune[,1]*(1- hune[]) )/colSums(1- hu...

2025/04/21 23:00

生起性を非曝露層でリストから数え上げる～ 船事例 risk値の信頼区間から 率でIDが切り捨てられ

・みかけを考える時、判明している因子からのみかけを計算である程度みることができる ・MH調整risk値で信頼区間；確率的な水準に基づいてみかけの影響を評価する資料があった 信頼区間により判断することで有意水準に許される範囲の発生IDは 有意性なし と切り捨てる2） 有...

2025/04/05 22:39

ﾒﾓ 信頼区間の幅で切り捨てられる発生～ときに曝露数により広がる

▼ kによる信頼区間の違い ・kが異なるとき、信頼区間の度数は変化する 0.95でみることに代えて、変曲点を度数幅でみる 平均 0.25 において k 小の範囲では plot(xlim=c(1,60),ylim=c(0,0.1),dhyper(1:100,20...

2025/04/01 22:23

近年の ”報告” ･･気になる点、優れた調査を参考に 対処メモ

■流行曲線 一峰であることと 単一曝露を結びつける ・・複数の山も短期に合わされば一峰になるよ なお、単一曝露だと一峰が尤もらしい ■欠席*などの曝露の違いを観察に加えない 2×2表がない 疑わしいイベント；例えば、1日目弁当と2日目弁当 が挙げられたとき、それぞれ...

2025/03/27 19:18

信頼下限 を示せば、 χ2値・P値は いらない･･と思うこのごろ

・納豆オクラ事例の解説資料には、度数、OR信頼下限ほかいろいろ計算した表がある ・信頼下限があれば、他の検定は不要であり、余計なものだと主張したい ・2013 H25 3月審議会 議論あり ・R記述の練習になるので考えてみる ▼ 数メニューのchi2 OR信頼下限、p値をみ...

2025/03/22 22:59

複数menuに高risk値が現れる意味･･発生数の”みかけ”と”実効”を計算する 交叉汚染は一つの可能性

・事例で粗な risk値の大きな、ときに”有意”となる因子が複数見つかることがしばしばある 交叉汚染とみなされる例があるが解明していない ”メモrisk値がうまく働かない”をもう少しきっちり調べる ・2因子の関係のみ考える時、risk値が明らかに大きい因子x1が、も...

2025/03/18 19:14

MHOR mantelhaen.test と 手計算；再

・Rの関数と定義式による手計算 定義式が分かりやすいように再掲載；過去記事削除 ■ 例としてデータ yts_ts[,,] , , sake = 0 tam y 0 1 0 25 18 1 ...

2025/03/17 23:14

メモ risk値がうまく働かない件

・有力な生起因子を探すためにrisk値を計算するのだが、複数；多数のmenuが高risk値になる例がある ・複数menuの汚染も1つの理由のようだ 調整・層化してみかけを数値上処理できる場合もあるが、万能ではない 不調を処理しきれない例もある ・不調になる例についてメモする...

2025/03/14 22:41

”症例研究” 発生度数差 a-b は どのような事例で主役を指し示すか･･ 7事例で調べる

・症例だけから原因に迫れて対策につながった例がある ・CCSにおいては、曝露の偏りを用いて因子の性質を探っている caseだけをみても、曝露の偏りが著しいものがある ・case a-b とは、どんなものか コホートについて実例の数値を調べてみる ・5事例/7事例で原因とされ...

2025/03/11 19:14

”症例研究”！驚くべきアプローチ 散発事例 サルモネラ、鶏卵関連

▼ 症例研究とでもいうべき 対照を置かずに対策が進んだ注目すべき中間報告 平成29年 2017 山形 散発症例 サルモネラ 15名把握、その後対象をあつめ 40 共通的な外食等なしだったが 曝露調査で卵が 多かった N 40 うち30が卵 喫食と. 1IDは”...

2025/03/06 21:06

生起因子との曝露率・RD,MHRD

RDは生起因子の重複の影響は残すが、素直な値になる MHRD計算の記述 RDとMHとの比較 ▼ osw事例 粗なRD mxrd0 MHRD oswmh ** RDによるMHRD plot(xlim=c(-0.3,0.4),ylim=c(-0....

2025/03/04 00:38

「粗なrisk値信頼下限を並べるのは避けよう」事例と 学べる資料の歴史

・粗なrisk値信頼下限、χ2値、p値の羅列でとどまる報告・解説がこのごろ多くみられる 食い合わせによるみかけに触れていないものがある 事例の解析や例示は学習の教材なのでながめる際に、注意したい点を記事にして残したい ・ほか気づいたことも記す ▼ まず目に留まったある事例 ...

2025/03/01 21:10