医薬品製造販売後調査に必要な例数を求める

有害事象の発現割合 p%で、95%の確率で少なくとも1例検出するに必要な例数を求めたいとします。

例数を算出する例を示します。

``` # 医薬品販売後調査の例数設定 # 少なくとも1例を見つける場合 power = 0.95 # 検出力 p = 0.001 # 発現割合0.1% = 0.001 n_bin = np.log(1-power)/np.log(1-p) print(n_bin) n_poi = -1*np.log(1-power)/p print(n_poi) ```

例数が決まっている場合の、検出力を算出する例を示します。

``` # 2例以上 from scipy.stats import binom_test,binom p=0.05 expected_num = 9 #ddof=1 # 311例で95を超えるので、この方法では311例が答え。 # 関数によって多少変動します for n in range(300,320,1): # 利用不可 # p_power = binom_test(expected_num, n, p, alternative='greater') # alternative : {'two-sided', 'greater', 'less'}, optional # print(n,p_power) p_power = binom.sf(expected_num, n, p) print(n,p_power) ```

Visionary Imaging Services, Inc.

Imaging CROサービスをご提供させていただきます。

お問い合わせお待ちしております。

https://www.vis-ionary.com/