この例では、あるチャレンジ（例えば、宝くじやくじ引きなど）を24万5000回試みた結果、1等が8回当たったという状況が与えられています。これを統計的に評価するために、p値（pバリュー）を計算する方法を考えます。 ### 仮定 1. **帰無仮説 (H₀):**…

#統計 P値関数の解説が載っている教科書は少ないのですが、P値関数とはaに仮説θ=aのP値を対応させる関数のことです。 1つ上の投稿にある演習問題を解けばP値関数についてよく理解できると思います。私もその演習問題を地道に解きました。

#統計 すでにパラメータθの区間推定法が与えられているとき(例えば通常の線形回帰ではt分布を使った区間推定法が与えられている)、仮説θ=aのP値は 信頼水準1-αの信頼区間の端点の1つがaになるαの値 として得られます。これが線形回帰の場合に具体的にどのように書けるかを求めるのは良い演習問題。

#統計 P値が最大の1にならない回帰係数の値達には、データの数値への適合度が最大ではない回帰直線が対応しています。 最尤法でもP値を適切に定義することによって、適合度が最大でない回帰直線達も考えることができるようになります。そのことを利用して区間推定も可能になります。 続く x.com/ds_school_1/st…

#統計 データの数値への回帰係数の適合度(最近の私は「相性の良さの程度」と表現することが多い)の指標の1つはP値と呼ばれており、最尤法で求めた回帰係数はデータの数値への適合度が最大なのでそのP値は1になります。 続く x.com/ds_school_1/st…

統計学を使う学術研究は、基本的にはp値を0.05（5％）に設定するので、極端な話、20回に1回は間違えるものを正しいものとする。であるから、ワクチンなどのヒトに関わる研究の場合、統計手法を変えたり、対象者を変えたり、プラセボを置いたりして、信頼性を高める。しかし、100％安全は無理である。

さてさて、試合終わったので このフレップ神話について3通りの解析をしました。 ①勝てない（負けor引き分け） ②勝たない（負けのみ） ③引き分けを除いた結果 いずれも有意差ありでした。 つまりフレップがいると… p値 ①0.00321 ②0.00777 ③0.00537 x.com/satorippi35/st…

#統計 1つ上の投稿に書いた枠組みを理解していないと、統計学における基本概念である P値と尤度(ゆうど) について全く理解できなくなります。 P値も尤度も ①データの数値 (確率変数ではない) と ②モデル内確率分布 (確率変数が出て来る) の相性の良さ(compatibility)の程度の指標になります。

返信先:@genkuroki質問へのお返事ですが、以下のコードです。Wald法ではstrata.score のところが strata.wald で計算していて、正しい値が得られています。 # p値計算 calculate_p_value <- function(estimate, se) { if (se == 0) { if (estimate == 1) 1 else 0 # オッズ比とリスク比の場合 } else {…

返信先:@PPubmed#統計 EMUYN広報さんに質問 添付画像①でのスコア法のP値はどのようにして計算していますか？ 私はepiR::epi.2by2で計算したのだと解釈して添付画像②のように計算してみましたが、添付画像①のスコア法のP値が見当たりません。 私はepiRに失礼なことを間違って言ってしまったかもしれません。 pic.x.com/hx7ubalgqw x.com/ppubmed/status…

#統計 別スレッドは以下のリンク先から始まります。 EMUYN広報さんが実装したスコア法のP値関数と私が実装したスコア法のP値関数が一致するどうかも知りたいところ。 私の実装の場合 ↓ github.com/genkuroki/publ… それぞれRD=0.1, RR=1.1, OR=1.1という帰無仮説(検定仮説)のP値を計算しています。 ↓ pic.x.com/glbo1sbnhy x.com/ppubmed/status…

返信先:@PPubmed#統計 私の数学的に素直な実装の計算結果(添付画像①②(2×2のデータが互いに転置の関係))と添付画像③を比較すると、適当に数字わせすると、信頼区間の数値は一致しています。 100(1-α)%信頼区間達(αを動かす)から逆にP値関数が得られるので、その意味では「本質的に同じ」だともみなせます。続く pic.x.com/h8lpququ0s

返信先:@PPubmed#統計 P値と信頼区間の整合性のためには パラメータθの100(1-α)%信頼区間 と P値≥αとなるパラメータθの値の範囲 を一致させる必要があります。 ユーザー側にとっては整合性のないP値と信頼区間を報告するのは苦しいので、不整合なアウトプットを出すソフトを使い難くなると思います。

#統計 有罪→統計的有意性 有罪→NHST 有罪→ポストP値時代などと言うこと 無罪→P値 ただしP値のcompatibilityによる解釈(ASA声明の原則1の明晰版)によって、P値は推定のための道具であることを理解する必要がある。その解釈は事後分布まで拡張される。

#統計 例えばこれは酷いと思いました。 ckt-hanbai.co.jp/wp/wp-content/… 練習問題1の「出た回数の数値例」が非現実的でかつ(2)の「小数第3位までで答えなさい」が酷い。 練習問題2(2)も酷い。試しにnull P値を計算したら巨大な値に。 pic.x.com/rsymar2ynr

最低用量でのP値1.01なので、（高容量の結果から結果が疑われているとはいえ）この用量でのP3が失敗する確率は低いでしょう

#統計 パラメータθのP値関数から、θの点推定値と信頼区間が (θの点推定値) = (P値関数を最大化するθの値) (θの100(1-α)%信頼区間) = (P値関数≥αとなるθの値の範囲) によって定義される。さらにθ=θ₀が「効果無し」を意味するならば、 (null P値) = (θ=θ₀でのP値関数の値) も定義される。続く

Bigqueryで手軽にウェルチのT検定できるUDFを考えてる。 select t_test(data1_query,data2_query,有意水準) で、 統計量T、p値、有意水準、棄却判定結果 のテーブルを返す

【名目上のp値と検証的なp値】 たくさん検定したらP<0.05の有意差出たものが1つあったー！ この1つの結果だけ出せばいいよね？？ ・・・ダメです！… pic.x.com/fuo33kf9rp

効果量のP値の信頼区間について、EXCELの計算式を追記しました。 ①平均値の差の検定の信頼区間を、プールした標準偏差で割ったものが、効果量の信頼区間 ②効果量のP値の信頼区間は、効果量の信頼区間の上側と下側に対して求めたP値（普通のP値は区間の平均値） data-science.tokyo/ed/edj1-1-3-4-… pic.x.com/tt3fuy3s5n

#統計 モデルのパラメータの値に関する仮説のP値 は モデルのパラメータの値に関する仮説とデータの数値の相性の良さの程度を表す指標の1つ に過ぎません。このことはASA声明の原則1に書いてあるとみなせるのですが、適切に引用している人を見たことがない。続く

P値は、0から1までの範囲の値を取ります。

返信先:@surinobaasan他2人二つの試験の有意差を計算しましたか？ 4価単独では、ギリギリ有意差が出ないのですが、2価は単独で有意差あり、 この二つをメタ解析すると オッズ比5.4 95%信頼区間2.91〜 10 p値1.1× 10⁻⁷ I²検定0% 有意差ありですね。計算が間違えていたら教えてください。

１人反省会 ぴーち→피치 (違和感に気づかない) 피 값→P値!!!!!(間違えた😭) 마안합니다 ㅠㅠ 値も치でギリギリわかってくれたのかな😭😭

【新刊】2024年9月号『グルカゴン・GLP-1・GIPの創薬革命』が本日発行！ 様々な代謝性疾患への有効性がみえてきたグルカゴン，GLP-1，GIP．急速に進む創薬研究について特集！ ☆注目記事　ビッグデータ時代の医学系研究におけるp値 ↓目次・購入はコチラから yodosha.co.jp/jikkenigaku/bo…

#統計 資料 ①二項分布モデルの尤度函数と3種のP値関数 ②正規分布モデルの尤度函数と「母平均はμである」という仮説のP値関数 pic.x.com/9limiccxq3

#統計 尤度に類似のP値は データの数値以上に極端な値がモデルの確率分布で生成される確率またはその近似値 として構成されることが多い。 P値もモデルのパラメータの関数になりますが、尤度と違ってモデルの一部のパラメータの関数になることが多く、その値の範囲は常に0から1になります。

Ca/P値によると思います 当院の方針としては ①Calcimimetics(ウパシタ,パーサビブ,オルケディア)でしっかりとi-PTHを下げる ② ①でCa/Pも連動で下がるので、低Caに対して、CaCO3やVDRAで補充 ③低PTHで高P血症はテナパノル塩酸塩などを活用 x.com/ryu04196650ryu…

#統計 ASA声明の原則1の解釈の仕方については以下のリンク先スレッドを参照。 ASA声明の原則1にはP値の正しい解釈の仕方が書いてあり、非常に重要です。ASA声明を引用する人は、以下のリンク先の注意に従って、原則1を引用するべきだと思います。 x.com/genkuroki/stat…

#統計 補足 ベイズ信用区間の場合に限らず、任意の区間推定法には aが閾値100(1-α)%での区間推定の結果の端点のとき (仮説θ=aのP値) = α によってP値を定義できます。例えば、ETI版信用区間には (仮説θ=aのP値) = (事後分布でのθ≤aの確率とθ≥aの確率の小さい方の2倍) が対応している。続く x.com/genkuroki/stat…

相関係数の検定を、表に加えました。 「従来からある相関係数の検定は、計算した相関係数の確からしさを表す」という位置付けにしています。 相関性自体についての検定は、P値として、 「1 - 相関係数の２乗」 を提案しています。 data-science.tokyo/ed/edj1-1-3-5.… pic.x.com/wknfevplgi

数学Bなら正規分布表が使えるので棄却域ですね。 数学1の方では（実質的に）p値です。 peing.net/ja/qs/11879778… #Peing #質問箱

#統計 以上では、リスク比RR=p/qを母比率p,qの違いの指標として採用した場合のP値関数の作り方について説明しましたが、リスク差RD=p-qを採用したり、オッズ比OR=(p/(1-p))/(q/(1-q))を採用した場合についても考えることができます。 そういう複雑な話は後回しにして、まずが二項分布モデルの(文字数)

#統計 その場合には添付画像のような方法を使います。 ルーレット1,2をそれぞれm,n回回したときに出た当たりの回数a,cから、ルーレット1,2での当たりが出る確率p,qの比p/qの推定をP値を使って行う方法の1つ。 対数リスク比に関するWald型P値関数。 pic.x.com/kjsj8ixgjc

動画中の話 ・世界で一番論文捏造してる研究者はぶっちぎり1位で日本人 ・凡ミスでスプレッドシート選択漏れ、偏った結果が出たが気付かない論文 ・白いマウスを黒いマジックで塗って皮膚移植成功論文 ・STAP細胞左右反転画像 ・p値ハッキング ・実は再現実験実はあんまりされない

#統計 二項分布の正規分布近似を使って作ったP値関数の構成は、そういう一般的かつ標準的な方法の1つである中心極限定理を使ったP値関数の構成法(スコア法)の特殊な場合になっています。 一般にはスコア法、Wald法、対数最大尤度比法を使う。

#統計 一般の場合のP値の構成では、モデルが調べたいパラメータ以外の余計なパラメータを含んでいることが問題になります。(例えば、平均μだけに興味があるのに、正規分布モデルは分散σ²もパラメータに持つ。) その問題の一般的かつ標準的な解決法の1つは中心極限定理を上手に使うことです。続く

#統計 P値関数は別に新しい概念ではない。当たる確率pの値に関する仮説を固定せずに動かして考えるだけ。 P値関数については、最初に1つ上の投稿で定義した最も易しいP値関数のグラフ(横軸はp、縦軸はP値)のグラフをコンピュータで描いてみるべき。 グラフのプロットの仕方の学習コストは結構高い。 pic.x.com/5ktvwxqlnh

#統計 P値関数の易しい例は、二項分布の正規分布近似で構成されたP値関数。定義: データの数値「n回中x回当たり」と検定仮説「当たる確率はpである」のP値は 標準正規分布に従う確率変数の絶対値が|x-np|/√(np(1-p))以上になる確率。 与えられたx,nについて、pを動かしてP値のグラフを描けば良い。 pic.x.com/9ldcvhmtxj x.com/genkuroki/stat…