欠測の母数を推定する
母数の推定 経済統計 鹿野研究室
... 復習問題 出席確認用紙に解答し（用紙裏面を用いても良い） 、退出時に提出せよ。 1. 未知の母数 θ （例えば母分散）に対し、標本 X 1 , X 2 , . . . , X n の使い方次第で二つの推定量 ˆθ 1 、 ˆθ 2 ...
6
高橋将宜 伊藤孝之 : 経済調査における欠測値補定方法について はじめに 1 個人や世帯を対象とする調査と異なり 企業の売上高といった経理項目を対象とする経済調査では 回答ユニットの形態が多岐にわたることもあり データの広がりが大きく かつ 記入漏れや記入誤りが生じやすい その結果 調査から得られる
... Y の欠測パターンを予測することがで きない。X2（サイコロ）の値を見ても、Y のどの値が欠測するのか予測ができない。しかし、 MCAR の場合とは異なり、網掛けとなっている Y ...
68
母数の仮説検定 経済統計 鹿野研究室
... ⊲ 母数（未知）と仮説値の差 θ − θ ∗ は測れない ⇒ 推定量 ˆθ と仮説値の差 ˆθ − θ ∗ に注目。 ⊲ 母数の推定（講義ノート #19 ）とは別の、統計的推測のスタイル。 例：あるビール工場では、品質基準として製品の平均アルコール度数を ...
6
2. 方法対象とする期間と地域対象期間は 2002 年から 2010 年までとした 対象地域は兵庫県本州部とする 用いたデータ推定には以下のデータを使用した 有害捕獲数 ( 年度 )i_yugai[i]:i 年度の有害許可による捕獲数 個体数を反映する指標として用いる 目撃効率 spue[i]:i
... 5. 2009 年の生息個体数の自然対数値 lnNins：2009 年度の個体数に関する事前の情報は ないため、分散を大きめに設定した。 6. 目撃効率、狩猟捕獲数、有害捕獲数の期待値からの誤差分散 v_spue v_ryo v_yugai： ...
12
欠測を含む順序カテゴリカル経時データの解析 -GEE プロシジャの有用性 - 駒嵜弘 1 藤原正和 2 ( 1 マルホ株式会社 2 塩野義製薬株式会社 ) Ordinal longitudinal data analysis with missing data -Usefulness of Proc
... GEEプロシジャ の順序カテゴリカルデータの GEE法は、Heagerty and Zeger(1996)の手法に基づいている。 →相関パラメータ算出の際、alternating logistic regression( ALR )を使用 主効果 の推定方程式は Liang and Zeger (1986) 、 ...
31
臨床試験の欠測データの取り扱いに関する最近の展開と今後の課題について
... サンプルサイズ計算で良く行われてきた方法は，欠測データを考慮せずにサンプルサイズ を一旦計算し，そこで得られた被験者数から想定した中止割合を考慮し増加させる (i.e., 必 要被験者 = 被験者数／ (1 – 想定する中止割合)という単純なものである ∗ ．NAS レポート ...
49
審査資料で欠測データが議論された事例
... れると考えられた．MMRM の結果のみに基づいて治療効果のロバストな推定値が得られた との申請者の主張は受け入れられない．特に，DPM-CF-301 試験では初期の中止割合が大き いため，事前に計画された MMRM による結果は治療効果を過大評価していると考えられた． ...
32
欠測を伴う標本調査に基づく統計的仮説検定について
... 定精度の低下幅は正確に評価できる。しかし、欠測バイアスの問題に関しては、バイアスを 欠測率により一般的に定量評価することは容易ではない。例えば、標本平均による母集団平 ...
21
「欠測値補完に関する調査研究報告書」 欠測データの統計的処理 その１
... った。これに対して、図２－１－２では、消費支出については、支出額の大きい世帯で 欠測が生じ、食料品支出については、支出額の小さい世帯で欠測が生じている。つま り、図２－１－２の利用可能ケース分析では、消費支出の平均値は支出額が相対的に ...
23
空間的自己相関モデルにおける距離行列のべき数の推定と地価決定要因分析への応用
... 一方、空間的自己相関モデルでは空間的な影響の受け方 を表す空間構造を空間重み付け行列で表すが、重み付け行 列の構造によって推定結果が左右されるために適切な重み 付け行列を設定することが重要となる。しかし空間構造は 観測できない未知のものであるから空間構造の恣意性を排 ...
6
スギ 埼 衛 研 所 報 第 47 号 213 年 結 果 及 び 考 察 1 平 成 25 年 の 測 定 結 果 (1) スギ 花 粉 毎 日 のスギ 花 粉 飛 散 数 の 測 定 結 果 を 表 1-1~ 表 1-5 に 示 した.また, 表 2 に 各 観 測 地 点 での 初 観 測 日
... 「さいたま」は同じであり， 「深谷」は 5 日早く， 「秩父」 は 4 日遅かった．一方，飛散終了日は「さいたま」は 20 日， 「深谷」は 3 日， 「秩父」は 15 日程度早かった．特 に「さいたま」はもっとも早く「秩父」は 2 番目に早か った．そのため，飛散期間は比較的短かった (表 3)． 例年，花粉数の推移は最大飛散日を頂点とした山型と なっているが，本年は異なっており，3 ...
7
討を行った.UWB MIMO 伝送方式の開発においては IR (Impulse Radio) 方式を用いた MIMO 伝送によるインプラント通信の伝送実験を行い, その特性評価を行う. UWB 帯パルスを用いた位置推定方式においては, パルスの到来時間を推定することによって送受信間距離を推定する方式
... 測定するためにベクトルネットワークアナライザ の Port1 に送信アンテナ，Port2 に受信アンテナを それぞれ接続した．測定環境を図 5 に示す．図 6 に 製作した送信アンテナと受信アンテナの反射係数 特性を示す．製作した送受信アンテナの反射係数 特性は ...
5
7 魚類の生息数調査カジカ及びサツキマス ( アマゴ ) について 生息数を推定した カジカは 対象事業実施区域を流れる 3 本の河川すべてに生息し 西側流域の推定数が多かった 一方 中央を流れる流域は上流の地点からは確認されず その下流に分布上限があると考えられた なお カジカの再捕獲数が少ないた
... 0 種 調整池工事により、その周辺に 生息する魚類の一部の個体が 移動する可能性がある。 ただし、水域全体に対する改変 区域は限られているとともに、 非改変区域にも同様の種が生 息しており、一定期間を経た後 に回復するものと考えられる ことから、工事による直接的影 響は小さいと予測する。 現段階で対象事業実施区域と ...
45
1. 気温と産業の関係 2. 気温と販売数の関係の分析 過去の気温データをダウンロードする 時系列グラフを描く 気温と販売の関係を調べる 散布図を描く 定量的な関係を求める 気温から販売数を推定する 2 週間先の気温予測を取得し 活用する 気温以外の要素の影響 3. 予報精度 過去の 1 か月予報
... 10月から12月までは気温の下降に伴い販売数が増加する。 1∼2月は気温が低い状態が続くが、販売数は減少する。 • 横軸は平均気温。赤破線は18℃を示す。 • 縦軸は店舗当たり販売数。 ...
53
ケプストラムによる音声のピッチ周波数推定 : 評価信号の提案と分析可能範囲の考察
... まとめ ここでは、ケプストラム法によるピッチ周波 数推定可能範囲について検討した。まず、ピッチ 周波数の推定可能範囲を測定するための評価信 号について考察を行った。正弦波、三角波、方形 波などはケプストラム法のアノレゴリズム上、本 法では正確なピッチが計算できないことを示し た。そして、スベクトノレの幅についての理論的 考察から、鋸歯状波に周波数変調したものを評 価信号[r] ...
11
希少疾病用医薬品の該当性 ( 推定対象患者数 推定方法についても記載する ) 国内の承認内容 ( 適応外薬のみ ) 併用において 通常 成人にはベバシズマブ ( 遺伝子組換え ) として 1 回 15mg/kg( 体重 ) を点滴静脈内注射する 投与間隔は 3 週間以上とする が国内で承認されている
... ２）Pujade-Lauraine E et al. Bevacizumab combined with chemotherapy for platinum-resistant recurrent ovarian cancer: The AURELIA open-label randomized phase III trial. J Clin Oncol. 2014; 32: 1302-8 4) ...
22
以下の内容について説明する 1. VAR モデル推定する 2. VAR モデルを用いて予測する 3. グレンジャーの因果性を検定する 4. インパルス応答関数を描く 1. VAR モデルを推定する ここでは VAR(p) モデル : R による時系列分析の方法 2 y t = c + Φ 1 y t
... vars を使用できる様に なる。 次に使用するデータを読み込む。教科書の著者のホームページからダウンロードした G7 国の 日次の MSCI(Morgan Stanley Capital International)データ (これは MSCI が作成している株式指 数) の(対数階差×100 ...
11
平成24年度知床世界自然遺産地域ヒグマ個体数推定のための解析業務報告書 [PDF]
... 体識別を行い、その情報をもとに個体数推定を行う調査手法である。そのため、外見上 の特徴からヒグマを個体識別できるかどうかは調査の重要なポイントである。 今回の調査では、胸部斑紋の情報に各種情報（構成や体色、体サイズ、耳タグ、子が ...
30
F σ 1 exp (3-1) ここで 形状母数 m と尺度母数 σ0 は材料定数とされ ワイブル応力 σw は Eq.(3-2) で定義される σ V V dv σ V V (3-2) ここで 積分領域 V は破壊が発生する可能性のある領域 体積要素 V0 は 最弱リンク機構 を構成する要素 σi
... 3mm/min の荷重速度で破壊するまで負荷した。負荷中の荷重と変位は、 KYOWA 製 Sensor Interface PCD-320A（以下、測定機と称す）を介して 10Hz で断続的に測定した。破断した試験片は 八光電機製作所製 A 型投込みヒーター（以下、加熱機と称す）で加熱したエタノール溶液 ...
143
2. 方法 対象とする期間と地域対象期間は 2002 年から 2010 年までとした 対象地域は淡路島とする 用いたデータ推定には以下のデータを使用した 有害捕獲数 ( 年度 )i_yugai[i]:i 年度の有害許可による捕獲数 個体数を反映する指標として用いる 目撃効率 spue[i]:i 年度
... 自然増加率や個体数の推定は、上記のデータと時系列的な関係を記述する階層ベイズモデル を構築し、パラメータの推定にはマルコフ連鎖モンテカルロ法によるランダムサンプリング ...
12