ここからの参照(PDF の 19 ページ) を使用して、直接標準化された発生率の信頼区間を計算しました。引用された章はTable of 95% Confidence Limits for the Number of Cases, for 1-29 Cases, Based on the Poisson Distribution、ここに添付されています。
パッケージ内のci.poisson関数を使用して番号を生成しようとしました:epicalc
> ci.poisson(0,29)
events person.time incidence se exact.lower95ci exact.upper95ci
0 29 0 0 0 0.1272069
出力は PDF のテーブルとは異なります。これは0.0000-である必要があります3.6889。数値を生成するために間違った関数を使用しているのではないかと思います。Rを使用してそのテーブルを構築するにはどうすればよいですか? ありがとう。