福島県産米の全袋検査：平成24年と25年の比較
-
- いいね！0
今井長兵衛
@medanjin
[米全袋検査 1] 米の全袋検査の結果を両対数グラフにプロットした図がTLに流れてきた。東大の早野さんの「参考出品」。→ @ hayano さんより 福島県の米の全量全袋検査の結果 http://t.co/YU5PKXBLMz
2013-11-11 20:30:59
今井長兵衛
@medanjin
[米全袋検査 2] 変なグラフだと文句を言ったら、「これはヒストグラム」という非公式＠ツイートが帰ってきた。これには、ビックリした。マジかー？？ 早野さんのこのツイは多くの人に「御意」とのツイとともにRTされた。なんだこれは、宗教団体かと二度ビックリ。
2013-11-11 20:32:27
今井長兵衛
@medanjin
[米全袋検査 3] ヒストグラムなら、それなりの表現法があると思いますよと返したら、「御指摘考慮します」。これまで300本ほどの論文を査読し、2000本以上の論文を読んできたが、こんな「ヒストグラム」にお目にかかったことはない。まともな学会誌なら、リジェクトが当然だろう。
2013-11-11 20:34:03
今井長兵衛
@medanjin
[米全袋検査 4] ある人から平成24年と25年の米の全袋検査の結果をまとめたサイトを教えて頂いた。この2つだ。 https://t.co/PmZ2z1OOoG https://t.co/L55WYEVFRA
2013-11-11 20:35:52
今井長兵衛
@medanjin
[米全袋検査 5] サイトによると、ベルトコンベア式全袋検査器でスクリーニング検査を行い、検査器の性能にもよるのだろうが、およそ50 Bq/kg 超となった検体を１kg秤量してゲルマで詳細検査を行う。それでも100 Bq/kg超の米は流通させないというシステムになっている。
2013-11-11 20:37:15
今井長兵衛
@medanjin
[米全袋検査 6] 平成24年は、10,343,757件のうち、スクリーニングで引っかかったのは867件（76－100 Bq/kg の72件と恐らく51～75 Bq/kgの中の値が高かった795件）だった。これらが詳細検査にまわされ、最終的に100 Bg/kg超は71件だった。
2013-11-11 20:38:25
今井長兵衛
@medanjin
[米全袋検査 7] 平成25年は、11月10日までの8,141,290件のうち、スクリーニングで引っかかったのは426件（51－75 Bq/kg の156件と恐らく25～50 Bq/kgの中の値が高かった270件）。詳細検査にまわされ、最終的に100 Bq/kg超は15件だった。
2013-11-11 20:39:29
今井長兵衛
@medanjin
[米全袋検査 8] サイトからも分かるように、スクリーニングで引っかかっても、詳細検査で25 Bq/kg未満というケースもある。30 kg のうちの1 kgだけ詳細検査するので、値にバラツキが出るのだろう。少なくとも、測定値に大きな開きがある場合は、複数回の詳細測定が必要だろう。
2013-11-11 20:40:29
今井長兵衛
@medanjin
[米全袋検査 9] とりあえず、公表されているスクリーニングと詳細検査のデータをまとめると、表のようになる。 http://t.co/mw78Sk2dlq
2013-11-11 20:42:06
拡大
今井長兵衛
@medanjin
[米全袋検査 10] さて、どのように解析するかだが、僕は統計の専門家ではないので、簡単な方法である２行５列のカイ２乗検定を試してみた。結果はきわめて大きなカイ２乗値が得られ、平成24年と25年の結果には統計的に有意な差が検出された。
2013-11-11 20:43:23
今井長兵衛
@medanjin
[米全袋検査 11] 次にどう解析するかだが、ここでも簡単な方法しか思いつかなかった。まず、平成25年にも24年と同数の検体が測定されたと仮定し、平成25年の各濃度階級の件数の割合に応じて総件数10,343,756を按分した。
2013-11-11 20:44:48
今井長兵衛
@medanjin
[米全袋検査 12] そして、各濃度階級の両年の件数を用い、１行２列のカイ２乗検定を行った。その結果、25 Bq/kg未満では平成25年＞平成24年、その他の４つの濃度階級では全て平成25年＜平成24年という有意差が有意水準１％あるいはもっと小さな確率で認められた。
2013-11-11 20:45:45
今井長兵衛
@medanjin
[米全袋検査 13] つまり、平成25年には平成24年より福島県産玄米の放射性セシウム汚染度が有意に改善されていると言うことである。最終的な結論は本年の全検体の検査が結果を待たなければならないが、少なくとも現在までの検査結果からは、そのように結論できる。大変、喜ばしいことである。
2013-11-11 20:46:51
今井長兵衛
@medanjin
[米全袋検査 14] 最後に、早野さんのグラフにケチを付けた落とし前として、正しいヒストグラムを図示しておく。縦軸は対数目盛なので、ほぼ同等に見える25 Bq/kg未満でも、統計的に有意な差があるので、誤解無きようお願いする。統計学の専門家による詳細な解析をお願いしたい。
2013-11-11 20:48:12
今井長兵衛
@medanjin
[米全袋検査 15] 字数オーバーのため、グラフとデータをこのツイに貼っておきます。連投、失礼しました。 http://t.co/31TOu8VdYq
2013-11-11 20:51:34
拡大