12/04/08 15:12:11.67 UM3uS8YD
> library(qcc)
x <- c(33, 29, 36, 53, 39, 39, 37, 37, 34, 40, 42, 34, 24, 41, 44, 30, 44, 41, 28, 29, 32, 33, 31, 33, 42, 45,
40, 38, 36, 37, 53, 96, 37, 34, 39, 40, 41, 48, 33, 37, 34, 38, 32, 39, 36, 32, 22, 41, 54, 35, 39, 35, 35, 38,
49, 37, 31, 28, 25, 33, 41, 32, 47, 33, 35, 34, 23, 27, 38, 41, 36, 33, 29, 32, 21, 28, 38, 30, 37, 28, 36, 32,
35, 29, 32, 44, 33, 37, 34, 31, 39, 31, 50, 37, 30, 39, 43, 30, 40, 30, 36, 37, 32, 38, 31, 37, 38, 39, 37, 33,
39, 39, 40, 31, 31, 37, 37, 26, 34, 36, 36, 34, 40, 31, 38, 33, 44, 41, 37, 26, 36, 41, 32, 33, 27, 34, 39, 41,
36, 28, 18, 35, 35, 40, 33, 32, 23, 28, 40, 30, 34, 45, 43, 31, 41, 53, 26, 34, 35, 39, 32, 33, 33, 36, 44, 28,
48, 40, 29, 34, 41, 35, 23, 40, 33, 32, 36, 31, 35, 29, 40, 36, 38, 44, 36, 24, 36, 33, 35, 32, 38, 43, 28, 24,
31, 27, 24, 41, 24, 35, 43, 47, 41, 41, 28, 28, 37, 50, 38, 35, 26, 32, 30, 35, 35, 35, 32, 24, 37, 32, 37, 41,
41, 37, 33, 35, 36, 41, 34, 29, 35, 29, 37, 42, 29, 33, 17, 31, 28, 32, 32, 28, 49, 34, 30, 37, 40, 47, 37, 28,
38, 23, 43, 38, 26, 29, 33, 36, 25, 36, 33, 39, 31, 37, 36, 30, 31, 27, 31, 30, 43, 40, 28, 34, 35, 39, 42, 39,
36, 42, 43, 31, 28, 38, 37, 35, 33, 39, 34, 40, 34, 38, 30, 30, 36, 38, 43, 36, 38, 37, 26, 34, 31, 34, 36, 43,
40, 40, 38, 26, 27, 26, 39, 30, 26, 32, 30, 36, 36, 45, 39, 41, 44, 33, 42, 47, 37, 49, 38, 37, 39, 31, 27, 36,
40, 25, 34, 31, 28, 36, 38, 39, 51, 34, 23, 44, 39, 39, 41, 44, 32, 42, 36, 31, 37, 37, 26, 35, 47, 33, 31, 32,
37, 42, 43, 65, 45, 35, 44, 32, 35, 50, 39, 22, 28, 42, 31, 27, 27, 37, 29, 33, 37, 27, 36, 39, 30, 30, 30, 36,
42, 17, 35, 40, 53, 30, 44, 47, 37, 52, 38, 32, 42, 25, 41, 41, 30, 38, 34, 41, 34, 38, 31, 29, 38, 28, 22, 34,
26, 26, 39, 33, 28, 38, 23, 39, 40, 40, 28, 41, 43, 35, 32, 29, 29, 32, 28, 43, 37, 24, 30, 53, 24, 33, 42, 32,
29, 40, 31, 36, 36, 36, 39, 43, 34, 41, 40, 40, 25, 48, 36, 38, 33, 49, 39, 44, 30, 41, 33, 24, 32, 37, 38, 36,
29, 35, 41, 26, 42, 33, 47, 40, 42, 31, 30, 25, 36, 32, 32, 42, 43, 41, 41, 40, 28, 37, 46, 32, 38, 38, 44, 32,
33, 37, 37, 36, 26, 39, 29, 45, 28, 19, 35, 34, 34, 36, 29, 37, 32, 45, 36, 31, 25, 27, 34, 29, 33, 28, 33, 30,
34, 40, 29, 34, 28, 38, 37, 33, 42, 36, 38, 34, 31, 29, 45, 54, 33, 33, 42, 38, 26, 36, 24, 22, 30, 56, 31, 30,
29, 24, 35, 37, 33, 24, 37, 26, 25, 27, 43, 35, 46, 40, 33, 31, 41, 32, 44, 28, 40, 33, 35, 38, 33, 36, 38, 40,
40, 42, 32, 31, 47, 37, 36, 38, 37, 37, 31, 28, 37, 31, 40, 30, 36, 34, 36, 32, 39, 35, 20, 32, 38, 38, 28, 40,
29, 31, 41, 37, 41, 39, 29, 32, 38, 35, 33, 39, 26, 37, 27, 26, 31, 31, 46, 36, 39, 32, 25, 28, 38, 30, 32, 44,
32, 36, 39, 33, 49, 43, 38, 41, 30, 42, 41, 39, 35, 39, 33, 36, 36, 42, 32, 25, 33, 35, 34, 44, 41, 37, 39, 36,
30, 38, 34, 33, 34, 40, 28, 34, 48, 30, 40, 40, 27, 37, 33, 42, 40, 35, 37, 28, 26, 27, 35, 27, 44, 37, 34, 37,
41, 37, 33, 33, 36, 34, 41, 37, 37, 52, 31, 39, 46, 40, 27, 24, 33, 37, 29, 40, 46, 40, 19, 37, 21, 27, 28, 28,
34, 35, 41, 38, 30, 32, 36, 42, 56, 30, 24, 42, 36, 37, 46, 36, 41, 35, 48, 33, 45, 34, 37, 25, 29, 43, 34, 34,
29, 28, 26, 31, 37, 38, 40, 38, 40, 45, 34, 41, 24, 37, 29, 42, 23, 30, 43, 34, 38, 32, 36, 28, 32, 37, 37, 37,
29, 34, 38, 33, 39, 35, 31, 52, 48, 34, 43, 43, 39, 32, 40, 30, 56, 87)
15:名無しに影響はない(栃木県)
12/04/08 15:16:23.53 UM3uS8YD
weight <- matrix(x,nrow=160)
sample <- rep(1:160, each=5)
d <- qcc.groups(weight, sample)
colnames(d) <- c("1st","2nd", "3th", "4th", "5th")
> qcc(data=d, type="xbar.one")
Call:
qcc(data = d, type = "xbar.one")
xbar.one chart for d
Summary of group statistics:
Min. 1st Qu. Median Mean 3rd Qu. Max.
17.00 31.00 36.00 35.47 39.00 96.00
Group sample size: 5
Number of groups: 800
Center of group statistics: 35.47125
Standard deviation: 6.584028
Control limits:
LCL UCL
15.71917 55.22333
>
334 CPM が 1 uSv/hですので、測定値(単位CPM)からuSV/hへの換算はできるでしょう。
16:名無しに影響はない(やわらか銀行)
12/04/08 22:08:53.34 39gPaL8L
>>15
たぶんほとんどの人が意味がわからないと思います。
１．測定対象
２．測定結果
３.補足コメント
４. データ貼り付け
のようにしたほうがわかりやすいと思います。
17:名無しに影響はない(栃木県)
12/04/09 20:07:17.87 WOmGBIk+
>>16
１．測定対象
バソコンを置いてある部屋の大気(2012.03.26日夕方から翌日にかけて)を測定したもの。単位CPM。インスペクターは 334 CPM が 1 uSv/h なので、シーベルト単位の値が必要な方は換算してください。
測定場所は栃木県内のどこか。官公庁の発表は 0.05-0.12 uSV/h 程度、インスペクターのユーザーが極端に少ないのでこれ以上の情報はご勘弁ください。
２．測定結果
>>14
３.補足コメント
問題なのが、５箇所ぐらい、１σ内に５点ぐらい連続して集まっている場所があります。
これは、おそらく、何かしらかの原因で外気が室内に入って、線量が下がったものと思われます。
何点が上側管理限界を超えている点がありますが、これは猫の騒音でしょう。
一般に、上下３σ内から外れる確率は 約2/1000で、管理限界を超えた時点で即異常値と判断できます。
しかし、そのようなわかりきった内容で管理図を眺めると肝心な内容を見落としてしまう。
２σ線(X管理図の場合にはUCLとLCLの間を6分割して、下から 3σ, 2σ, 1σ, CL, 1σ, 2σ, 3σ線)内(上下2σ線の間)に全データ数の95%が入っているか、
１σ線内に66%が入っているか、上下が対象か、に注意してください。
２σ線と３σ線の間に集まっていて、１σ線内の数が少ない場合には、2つの群に分けなければならない測定値をひとまとめにして測定していることがわかります。
上下が対照でなく特定の範囲に固まって存在する場合にも同様なことがいえます。
高濃度な空気と低濃度な空気が混在していて、変な偏りがあり、微量な分析には適さないのです。
日本国内で使われている「Ｘ管理図」と呼ばれている物は３種類あります。そのうちのひとつが qcc で描かれる管理図で、連続する何個か(ここでは５個を使用)をひとつの群(group)として取り扱い、
合計800個のデータを、160群に分けて使用し、中心線(CL)・管理限界(UCL, LCL)を計算する方法のようです。
計算結果を見ると、CL=35.47125, UCL=55.22933, LCL=15.71917 です。σ=6.584028 という計算結果がありますので、これを使うと、
１σ線が 35.47125+6.584028 = 42.055278 と 35.47125-6.584028 = 28.887222 の２本になります。同様に２σ線が 35.47125+6.584028*2 = 48.639306 と 35.47125-6.584028*2 = 22.303194 に、
３σ線は既にある UCL, LCLと同じ 55.22933 と 15.71917 になります。
４. データ貼り付け
> library(qcc)
>x <- c(33, (中略) , 56, 87) # >>14 生データ。
>weight <- matrix(x,nrow=160)
>sample <- rep(1:160, each=5)
>d <- qcc.groups(weight, sample)
>colnames(d) <- c("1st","2nd", "3th", "4th", "5th")
> qcc(data=d, type="xbar.one")
Call:
qcc(data = d, type = "xbar.one")
xbar.one chart for d
Summary of group statistics:
Min. 1st Qu. Median Mean 3rd Qu. Max.
17.00 31.00 36.00 35.47 39.00 96.00
Group sample size: 5
Number of groups: 800
Center of group statistics: 35.47125
Standard deviation: 6.584028
Control limits:
LCL UCL
15.71917 55.22333
>
18:名無しに影響はない(栃木県)
12/04/10 19:19:15.20 rYU0Fh//
１．測定対象
2012.04.08 採取, 裏の畑ふきのとう, 40度１日乾燥 (11*9*2cm) 19g の測定値。
２．測定結果
まずは生データ。
1. n=31 BG
x <- c(42, 42, 26, 41, 36, 28, 51, 45, 42, 46, 42, 41, 40, 44, 30, 29, 41, 38, 39, 44, 34, 38, 35, 40, 52, 35, 43, 26, 34, 45, 47)
> mean(x) [1] 39.22581
> var(x) [1] 44.84731
2. n=32 試料
x <- c(65, 65, 70, 77, 67, 51, 68, 62, 72, 57, 58, 74, 76, 71, 71, 67, 70, 85, 50, 64, 89, 56, 57, 84, 68, 64, 55, 61, 69, 66, 63, 62)
> mean(x) [1] 66.6875
> var(x) [1] 84.67339
3. n=34 BG
x <- c(39, 49, 45, 32, 35, 39, 35, 33, 37, 40, 45, 35, 42, 33, 38, 32, 36, 35, 38, 34, 36, 30, 47, 35, 33, 34, 33, 39, 42, 25, 41, 31, 38, 31)
> mean(x) [1] 36.67647
> var(x) [1] 26.58913
4. n=33 試料
x <- c(63, 61, 70, 68, 76, 86, 66, 67, 66, 59, 86, 80, 61, 68, 65, 55, 84, 64, 65, 64, 62, 75, 62, 63, 66, 65, 79, 55, 66, 63, 61, 59, 63)
> mean(x) [1] 67.0606
> var(x) [1] 67.24621
5. n=33 BG
x <- c(40, 37, 37, 48, 42, 48, 39, 44, 37, 33, 45, 46, 31, 30, 25, 31, 36, 33, 39, 36, 45, 39, 40, 40, 40, 47, 42, 31, 48, 46, 55, 35, 42)
> mean(x) [1] 39.60606
> var(x) [1] 42.62121
19:名無しに影響はない(栃木県)
12/04/10 19:19:51.52 rYU0Fh//
４. データ貼り付け
全体の分析
x <- c(42, 42, 26, 41, 36, 28, 51, 45, 42, 46, 42, 41, 40, 44, 30, 29, 41, 38, 39, 44, 34, 38, 35, 40, 52, 35, 43, 26, 34, 45, 47 ,
65, 65, 70, 77, 67, 51, 68, 62, 72, 57, 58, 74, 76, 71, 71, 67, 70, 85, 50, 64, 89, 56, 57, 84, 68, 64, 55, 61, 69, 66, 63, 62 ,
39, 49, 45, 32, 35, 39, 35, 33, 37, 40, 45, 35, 42, 33, 38, 32, 36, 35, 38, 34, 36, 30, 47, 35, 33, 34, 33, 39, 42, 25, 41, 31, 38, 31 ,
63, 61, 70, 68, 76, 86, 66, 67, 66, 59, 86, 80, 61, 68, 65, 55, 84, 64, 65, 64, 62, 75, 62, 63, 66, 65, 79, 55, 66, 63, 61, 59, 63 ,
40, 37, 37, 48, 42, 48, 39, 44, 37, 33, 45, 46, 31, 30, 25, 31, 36, 33, 39, 36, 45, 39, 40, 40, 40, 47, 42, 31, 48, 46, 55, 35, 42)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 , 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5)
> mean(x)
[1] 49.79755
> var(x)
[1] 247.36
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 149.7578, 第1自由度 = 4, 第2自由度 = 158, P値 < 2.2e-16
有意差あり、群による差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 12.2524, 自由度 = 4, P値 = 0.01557
有意差あり。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 134.2289, 第1自由度 = 4.000, 第2自由度 = 77.717, P値 < 2.2e-16
有意差あり、群による差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 6 3.680982 3.680982
30 20 12.269939 15.950920
35 29 17.791411 33.742331
40 24 14.723926 48.466258
45 16 9.815951 58.282209
50 4 2.453988 60.736196
55 10 6.134969 66.871166
60 17 10.429448 77.300613
65 18 11.042945 88.343558
70 7 4.294479 92.638037
75 5 3.067485 95.705521
80 3 1.840491 97.546012
85 4 2.453988 100.000000
>
20:名無しに影響はない(栃木県)
12/04/10 19:20:19.75 rYU0Fh//
試料とバックグラウンドの比較。
x <- c(42, 42, 26, 41, 36, 28, 51, 45, 42, 46, 42, 41, 40, 44, 30, 29, 41, 38, 39, 44, 34, 38, 35, 40, 52, 35, 43, 26, 34, 45, 47 ,
65, 65, 70, 77, 67, 51, 68, 62, 72, 57, 58, 74, 76, 71, 71, 67, 70, 85, 50, 64, 89, 56, 57, 84, 68, 64, 55, 61, 69, 66, 63, 62 ,
39, 49, 45, 32, 35, 39, 35, 33, 37, 40, 45, 35, 42, 33, 38, 32, 36, 35, 38, 34, 36, 30, 47, 35, 33, 34, 33, 39, 42, 25, 41, 31, 38, 31 ,
63, 61, 70, 68, 76, 86, 66, 67, 66, 59, 86, 80, 61, 68, 65, 55, 84, 64, 65, 64, 62, 75, 62, 63, 66, 65, 79, 55, 66, 63, 61, 59, 63 ,
40, 37, 37, 48, 42, 48, 39, 44, 37, 33, 45, 46, 31, 30, 25, 31, 36, 33, 39, 36, 45, 39, 40, 40, 40, 47, 42, 31, 48, 46, 55, 35, 42)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 594.8668, 第1自由度 = 1, 第2自由度 = 161, P値 < 2.2e-16
有意差あり、試料とバックグラウンドによる差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 8.4869, 自由度 = 1, P値 = 0.003577
有意差あり。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 522.673, 第1自由度 = 1.00, 第2自由度 = 107.23, P値 < 2.2e-16
有意差あり、試料とバックグラウンドによる差異がある。
>
21:名無しに影響はない(栃木県)
12/04/10 19:21:41.71 rYU0Fh//
繰り返しによる差異
1. BG
x <- c(42, 42, 26, 41, 36, 28, 51, 45, 42, 46, 42, 41, 40, 44, 30, 29, 41, 38, 39, 44, 34, 38, 35, 40, 52, 35, 43, 26, 34, 45, 47 ,
39, 49, 45, 32, 35, 39, 35, 33, 37, 40, 45, 35, 42, 33, 38, 32, 36, 35, 38, 34, 36, 30, 47, 35, 33, 34, 33, 39, 42, 25, 41, 31, 38, 31 ,
40, 37, 37, 48, 42, 48, 39, 44, 37, 33, 45, 46, 31, 30, 25, 31, 36, 33, 39, 36, 45, 39, 40, 40, 40, 47, 42, 31, 48, 46, 55, 35, 42)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x)
[1] 38.46939
> var(x)
[1] 38.72586
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.247, 第1自由度 = 2, 第2自由度 = 95, P値 = 0.1113
有意差あり、繰り返しによる差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 2.4918, 自由度 = 2, P値 = 0.2877
有意差あり。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.581, 第1自由度 = 2.000, 第2自由度 = 61.442, P値 = 0.08389
有意差あり、繰り返しによる差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 6 6.122449 6.122449
30 20 20.408163 26.530612
35 29 29.591837 56.122449
40 24 24.489796 80.612245
45 16 16.326531 96.938776
50 2 2.040816 98.979592
55 1 1.020408 100.000000
>
22:名無しに影響はない(栃木県)
12/04/10 19:22:25.06 rYU0Fh//
2. 試料
x <- c(65, 65, 70, 77, 67, 51, 68, 62, 72, 57, 58, 74, 76, 71, 71, 67, 70, 85, 50, 64, 89, 56, 57, 84, 68, 64, 55, 61, 69, 66, 63, 62 ,
63, 61, 70, 68, 76, 86, 66, 67, 66, 59, 86, 80, 61, 68, 65, 55, 84, 64, 65, 64, 62, 75, 62, 63, 66, 65, 79, 55, 66, 63, 61, 59, 63)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 66.87692
> var(x)
[1] 74.67212
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.0298, 第1自由度 = 1, 第2自由度 = 63, P値 = 0.8634
有意差無し。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.4111, 自由度 = 1, P値 = 0.5214
有意差無し。先の分散分析は有効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.0297, 第1自由度 = 1.000, 第2自由度 = 61.689, P値 = 0.8637
有意差無し。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
50 2 3.076923 3.076923
55 9 13.846154 16.923077
60 17 26.153846 43.076923
65 18 27.692308 70.769231
70 7 10.769231 81.538462
75 5 7.692308 89.230769
80 3 4.615385 93.846154
85 4 6.153846 100.000000
>
23:名無しに影響はない(栃木県)
12/04/10 19:22:48.67 rYU0Fh//
バックグラウンドと試料の比較
1. BG
gr1 <- c(42, 42, 26, 41, 36, 28, 51, 45, 42, 46, 42, 41, 40, 44, 30, 29, 41, 38, 39, 44, 34, 38, 35, 40, 52, 35, 43, 26, 34, 45, 47 ,
39, 49, 45, 32, 35, 39, 35, 33, 37, 40, 45, 35, 42, 33, 38, 32, 36, 35, 38, 34, 36, 30, 47, 35, 33, 34, 33, 39, 42, 25, 41, 31, 38, 31 ,
40, 37, 37, 48, 42, 48, 39, 44, 37, 33, 45, 46, 31, 30, 25, 31, 36, 33, 39, 36, 45, 39, 40, 40, 40, 47, 42, 31, 48, 46, 55, 35, 42)
> mean(x)
[1] 38.46939
> var(x)
[1] 38.72586
gr2 <- c(65, 65, 70, 77, 67, 51, 68, 62, 72, 57, 58, 74, 76, 71, 71, 67, 70, 85, 50, 64, 89, 56, 57, 84, 68, 64, 55, 61, 69, 66, 63, 62 ,
63, 61, 70, 68, 76, 86, 66, 67, 66, 59, 86, 80, 61, 68, 65, 55, 84, 64, 65, 64, 62, 75, 62, 63, 66, 65, 79, 55, 66, 63, 61, 59, 63)
> mean(x)
[1] 66.87692
> var(x)
[1] 74.67212
t.test(gr1, gr2, v=T)
t.test(gr1, gr2)
var.test(gr1, gr2)
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = -24.3899, 自由度 = 161, P値 < 2.2e-16
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -30.70765 -26.10743
標本推定値:
平均値x 平均値y
38.46939 66.87692
有意。平均値に差がある。
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = -22.862, 自由度 = 107.23, P値 < 2.2e-16
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -30.87071 -25.94436
標本推定値:
平均値x 平均値y
38.46939 66.87692
有意。平均値に差がある。
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.5186, 第1自由度 = 97, 第2自由度 = 64, P値 = 0.003399
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.3273010 0.8046984
標本推定値:
分散比
0.518612
有意差あり。分散は等しくない。平均値の差の検定は（Welchの方法）を使用する。
>
24:名無しに影響はない(栃木県)
12/04/10 19:24:10.18 rYU0Fh//
３.補足コメント
本来は、重み付演算を使うのですが、面倒なので普通の算術演算で
BG: 38.46939, 試料: 66.87692 CPM
よって、66.87692 - 38.46939 = 28.40753
本来の統計処理だと、試料のデータ数が100個以下の65個だから有効桁が１つ増えて、小数点以下第一位まで求める。
BGも100個以下の98だから有効桁が1つ増えて、小数点以下第１位まで求めて丸める。
BG: 38.5, 試料: 66.9 CPM
66.9-38.5=28.4
インスペクターによるやさしおの測定では、1.819個 / (分・100mg)なので、インスペクターの計数値からの換算計数は 0.47 Bq/CPM となる。
よって、測定した試料の線量は
28.4 * 0.47 = 13.3 Bq
となる。
試料の寸法は、 縦: 11cm, 横: 9cm, 高さ: 2cm で、風袋を温めた重さが 19g (ここではポリ袋の重さを0とする)より、比重は約 0.096 g/cm3 である。
単位面積あたりの試料の重さが 1g以下なので、試料によるβ線の吸収を0とする。
測定範囲が底面の半径 7.4cm、高さ2cm、上面の半径を4.1cm とする円錐台と考えれば、試料全体がほぼ測定可能範囲内にはいるので、19g全量を計ったとする。
よって、
13.3 * 1000 / 19 = 700
より、700Bq/kgとなる。
ふきの水分量が 96%(URLﾘﾝｸ(food.longseller.org))なので、ふきのとうの水分量を同量とすると、700 * (100-96.1) / 96.1 = 28.4 より
栃木産ふきのとう(自家栽培)の線量は 28 Bq/kg となる。
インスペクターはβ線測定を特長とする測定器ですから、γ線を測定を基準としている正確な測定とは異なります。
25:名無しに影響はない(栃木県)
12/04/10 19:39:24.31 rYU0Fh//
2022.13.25 マイクの集音能力の測定。
2-3分に１回の３連音を受信できず
猫等の生活音を受信する
パソコン値, インスペクター値, 差,備考
08:10 317, 276, 41
08:20 620, 584, 36
08:31 1044, 1010, 34
08:40 1289, 1260, 29
08:50 1644, 1613, 31, 猫音
09:00 2037, 2006, 31
これは、マイクロホンで音を拾って計数するというソフトの性能試験結果です。
26:名無しに影響はない(やわらか銀行)
12/04/11 02:20:45.27 8YEj24yd
１．測定対象
乾燥わかめ １５０ｇ 中国産 直付け ビニル袋で保護
２．測定結果
5分間バックグラウンド 198
5分間乾燥わかめ 251
３.補足コメント
バックグラウンドとの差があった。
カリウムを検出したと思う。三陸産わかめがあれば測ってみて比較したい。
27:名無しに影響はない(やわらか銀行)
12/04/11 02:30:14.57 8YEj24yd
>>24
ベクレルまでわかるんですか。
カルシウムの多い東日本の食品にはストロンチウムが入っているように思います。
干しえびや魚肉ソーセージや煮干しなど。
海に流出した原発の汚染水に含まれている放射性物質がわかればいいですが。
プルトニウムは豆類に蓄積するようです。
28:名無しに影響はない(やわらか銀行)
12/04/11 02:38:08.57 8YEj24yd
>>26はインスペクター＋のトータルタイマーでの5分間のカウント数です。
29:名無しに影響はない(やわらか銀行)
12/04/11 03:16:37.49 8YEj24yd
インスペクター＋ トータルタイマーモード 5分間カウント数
１．測定対象
丸大食品（大阪） フィッシュソーセージ７０ｇ 購入時の袋のまま直付け ビニル袋で保護
２．測定結果 5分間を2回測定
バックグラウンド 195、198
フィッシュソーセージ 201、198
３.補足コメント
数値に差がないようだ。次回はフィッシュソーセージを袋から出して測りたい。
30:名無しに影響はない(やわらか銀行)
12/04/11 03:22:10.59 8YEj24yd
インスペクター＋ トータルタイマーモード 5分間カウント数
１．測定対象
（株）ヒガシマル（鹿児島） そうめん７００ｇ 袋のまま直付け ビニル袋で保護
２．測定結果 5分間を2回測定
バックグラウンド 195、198
そうめん 219、208
３.補足コメント
数値に差があるが誤差の範囲か。 原材料は小麦粉と食塩のみ
31:名無しに影響はない(やわらか銀行)
12/04/11 03:36:54.96 8YEj24yd
インスペクター＋ トータルタイマーモード 5分間カウント数
１．測定対象
（株）ミタカ（熊本） 佐賀県産大麦 麦茶ティーバッグ400ｇ 袋のまま直付け ビニル袋で保護
２．測定結果 5分間
バックグラウンド 195、198
麦茶 230
32:名無しに影響はない(やわらか銀行)
12/04/11 03:41:44.72 8YEj24yd
インスペクター＋ トータルタイマーモード 5分間カウント数
１．測定対象
二豊フーズ（株）（大分県）納豆40ｇ 直付け ビニル袋で保護
２．測定結果 5分間
バックグラウンド 195、198
納豆 229
33:名無しに影響はない(栃木県)
12/04/11 14:51:46.67 3KPqPkME
>>12
ｽﾚﾘﾝｸ(radiation板:357-359番)
＞GJ!
＞やっとこのスレもスッキリしたね！
ってあるから、エーヨッ(共通語訳 いいとも)。
>>27
まず、単位の定義を読んでくれ。
URLﾘﾝｸ(ja.wikipedia.org)
URLﾘﾝｸ(ja.wikipedia.org)
URLﾘﾝｸ(ja.wikipedia.org)
ンデ、測定器のほうで何をやっているのかッテート、
何分間又は何秒間又は何時間か、放射線の数(ガイガーのカリカリ音)を数えて
適切な演算をして、表示しているだけ。
つまり、検出器とカリカリ音を数えることができれば、
空中線量も、表面線量も、物質中の放射線量も数えられるわけ。
だから(精度は別にして)カウンターがついている放射線量計ならば測れることになります。
実際にどういう値になるかは、測定者の力量次第。
ここはにちゃんねるだから、嘘を嘘だとわかる人が利用しているわけで、嘘を見抜けないならば利用するな。
＞プルトニウムは豆類に蓄積するようです。
そうなんだよな。
除去方法がわかれば、大豆でもインゲンでも小豆でもをいくらでも作付けするんだが
わかんないんで、作付けできないでいる。
34:名無しに影響はない(やわらか銀行)
12/04/11 15:32:28.04 8YEj24yd
福島、宮城、茨城が３大汚染地域です。
この地域の大豆や小豆などの豆類からプルトニウムが検出できるかもしれません。
汚染水にも含まれているので海藻や海底魚からも検出されそうです。
プルトニウムを検出する場合は、マイカ窓と対象物の間にビニルなど何もあってはいけません。
マイカ窓を３センチ程度まで近づけた数値と、間に紙一枚挟んだ数値に差があればプルトニウムです。
福島の自動車フィルターからはすでにプルトニウムが検出されています。
URLﾘﾝｸ(www.youtube.com)
35:名無しに影響はない(やわらか銀行)
12/04/11 15:42:00.47 8YEj24yd
納豆菌はほとんどが東日本産。
納豆菌は放射性物質で汚染するのであろうか？
九州で作った外国産大豆納豆はどこの納豆菌を使用しているか不明。
東京都、山形県、宮城県のメーカーが国内３大メーカー。
URLﾘﾝｸ(www.nattou.com)
36:名無しに影響はない(栃木県)
12/04/11 19:26:06.13 3KPqPkME
>>35
たしか、1gで100kgだったかな、使用量が。
だから、菌体が汚染されていたとしても薄まってしまうので、気にする必要なし。
>>34
α線放出核種は、重大な吸入障害を起こすので、呼吸器の防護を必要とする
放射線緊急事態時の評価および対応のための一般的手順
現場監督者の対応 手順書 Ci 2/4ページ内コラム(アクロバットではp40)より。
URLﾘﾝｸ(www.nirs.go.jp)
より。
口をあけて寝ていると、口内から出血をしたりしていますから、ある程度のα線核種が存在します。
また、集落内在住老人が「こんな刺すような鋭い痛みは経験していない」ということから、
2011年6月頃にストロンチウムよりも強いβ線を出す核種、おそらく、ネプツニウムが存在していたと思われます。
つまり、プルトニウムがそれなりの分量で家屋内に存在することが予想されます。
α線核種の場合には、ちょっとした薄膜程度で遮蔽されますので、どうしても、灰化しないと測定できません。
現時点では、ドラフトの入手ができず灰化はできません。
URLﾘﾝｸ(www.n-hakko.com) 灰化
URLﾘﾝｸ(ci.nii.ac.jp) 高温灰化による消失
URLﾘﾝｸ(www.jrias.or.jp) 低温灰化による消失
乾燥の場合にも、核種の希散と思われる現象が観察されています。机を改造した乾燥機を使用していますが
しいたけの乾燥中に、室内の線量が 0.05-0.12(最頻値 0.08) → 0.05-0.15(最頻値 0.12)と上昇しました。
37:名無しに影響はない(栃木県)
12/04/11 19:34:07.62 3KPqPkME
家庭内騒音による誤差。
測定値は以下のとおり。
x <- c( 0, 0, 2, 3, 2, 7, 3, 6, 7, 9, 7, 11, 7, 10, 15, 7, 9, 6, 9, 11, 12, 7, 8, 8, 13, 7, 9, 9, 15, 12, 20, 16, 17, 19, 14, 9, 12, 17, 26, 13, 16, 17, 11, 16, 19)
戸の開け閉め(2-9と10-15)、咳(26)、その他、思いつく限りの音を出しています。
１分間に１回、で画面を見ながら繰り返しを取っています。
38:名無しに影響はない(やわらか銀行)
12/04/11 22:44:22.19 8YEj24yd
インスペクター＋ トータルタイマーモード 5分間カウント数
１．測定対象
長崎県産 小魚（きびなご）１００ｇ 直付け ビニル袋で保護
２．測定結果 5分間
バックグラウンド 195、198
きびなご 210
39:名無しに影響はない(栃木県)
12/04/14 07:52:56.32 PzUMpoQ1
１．測定対象
「JヌードルコーポレーションBK こだわりの熟成 2014.01 B3941204」の測定。20g
１回目の「茹で」は、説明書どおり。水200ccに麺20gをいれて茹でたもの
２回目の「茹で」は、沸騰した200ccのお湯に、茹であがった麺をいれて放置。
１５分後麺が透き通ってきたら加熱を止めて、「茹であがった」として説明書3の通りにした。
「茹であがった麺」は、合金標準バット６号( URLﾘﾝｸ(taniguchi-metal.com) )に標準バット用金網を入れて、
キッチンペーパー( URLﾘﾝｸ(item.rakuten.co.jp) )１組をすいた上に広げた。
キッチンペーパーを丸めて先端を金網とバットの隙間に押し込んで、机を改造した乾燥機(40度)に４日間放置した。
乾燥後、麺は外側のキッチンペーパーを捨て内側のキッチンペーパーごと、
ユニパック F-4( URLﾘﾝｸ(www.seinichi.co.jp) )に入れて測定試料とした。
「袋」はF-4にキッチンペーパー１枚を入れたもの。「未茹で」は、麺20gを6cm-10cm程度に折り、キッチンペーパー１枚で包んだもの。
いずれも、直径10cm以内に入る大きさに折りたたんで、机の上に置き、上空3cmにインスペクターを置いて、３０分程度測定。音をパソコンで集計したもの。
40:名無しに影響はない(栃木県)
12/04/14 07:54:58.04 PzUMpoQ1
２．測定結果
まずは生データ。
1. n=30 BG
x <- c(26, 46, 30, 40, 47, 41, 32, 42, 30, 40, 44, 39, 42, 51, 44, 38, 38, 34, 48, 54, 44, 34, 31, 29, 45, 36, 53, 31, 41, 38)
> mean(x)
[1] 39.6
> var(x)
[1] 53.69655
>
2. n=31 袋
x <- c(34, 42, 40, 35, 38, 39, 29, 54, 39, 40, 35, 50, 35, 35, 40, 41, 33, 38, 33, 36, 47, 40, 33, 38, 48, 40, 42, 33, 25, 29, 39)
> mean(x)
[1] 38.06452
> var(x)
[1] 37.72903
>
3. n=35 未茹で
x <- c(43, 51, 38, 37, 45, 38, 41, 42, 35, 35, 37, 45, 39, 37, 41, 43, 46, 38, 49, 33, 35, 35, 43, 49, 50, 33, 43, 45, 48, 50, 39, 21, 41, 41, 46)
> mean(x)
[1] 40.91429
> var(x)
[1] 38.49244
>
4. n=30 茹で１回
x <- c(37, 38, 40, 51, 37, 38, 54, 36, 25, 36, 38, 45, 53, 32, 42, 50, 37, 36, 46, 60, 34, 50, 48, 34, 30, 37, 34, 31, 34, 39)
> mean(x)
[1] 40.06667
> var(x)
[1] 66.4092
>
5. n=32 茹で２回
x <- c(37, 31, 37, 47, 38, 39, 42, 38, 40, 33, 41, 42, 35, 36, 43, 33, 42, 49, 41, 38, 39, 36, 48, 39, 33, 35, 32, 25, 35, 38, 32, 45)
> mean(x)
[1] 38.09375
> var(x)
[1] 27.44254
>
6. n=31 BG
x <- c(37, 32, 46, 42, 33, 40, 38, 41, 35, 39, 34, 31, 43, 37, 43, 42, 48, 43, 26, 36, 29, 45, 32, 36, 43, 35, 50, 44, 44, 41, 47)
> mean(x)
[1] 39.09677
> var(x)
[1] 34.75699
>
41:名無しに影響はない(栃木県)
12/04/14 07:55:22.73 PzUMpoQ1
7. n=35 袋
x <- c(50, 27, 34, 42, 45, 38, 37, 39, 26, 51, 27, 37, 40, 40, 30, 43, 30, 47, 35, 34, 54, 38, 36, 41, 38, 26, 40, 32, 37, 32, 40, 40, 50, 37, 37)
> mean(x)
[1] 38
> var(x)
[1] 49.82353
>
8. n=31 茹で１回
x <- c(39, 27, 37, 44, 49, 49, 38, 48, 56, 35, 42, 41, 38, 34, 40, 27, 36, 40, 31, 38, 40, 29, 37, 37, 37, 43, 34, 39, 41, 41, 37)
> mean(x)
[1] 38.83871
> var(x)
[1] 39.47312
>
9. n=31 未茹で
x <- c(46, 40, 38, 33, 55, 44, 53, 39, 47, 41, 39, 31, 47, 44, 36, 38, 47, 33, 34, 36, 39, 44, 40, 55, 34, 36, 40, 31, 50, 53, 50)
> mean(x)
[1] 41.70968
> var(x)
[1] 50.67957
>
10. n=30 茹で２回
x <- c(37, 44, 49, 30, 37, 35, 41, 44, 32, 36, 29, 33, 45, 43, 33, 33, 46, 43, 48, 39, 31, 36, 41, 40, 32, 34, 43, 38, 39, 44)
> mean(x)
[1] 38.5
> var(x)
[1] 31.36207
>
11. n=30 BG
x <- c(38, 40, 40, 41, 37, 41, 47, 33, 33, 42, 37, 45, 32, 38, 35, 46, 36, 38, 36, 44, 26, 40, 42, 53, 41, 44, 37, 32, 40, 35)
> mean(x)
[1] 38.96667
> var(x)
[1] 28.72299
>
12. n=31 袋
x <- c(39, 32, 42, 41, 38, 46, 28, 37, 46, 35, 34, 40, 35, 41, 42, 43, 41, 31, 47, 46, 31, 42, 27, 51, 37, 41, 41, 34, 29, 37, 45)
> mean(x)
[1] 38.67742
> var(x)
[1] 36.29247
>
42:名無しに影響はない(栃木県)
12/04/14 07:56:14.68 PzUMpoQ1
４. データ貼り付け
全体の分析
x <- c(26, 46, 30, 40, 47, 41, 32, 42, 30, 40, 44, 39, 42, 51, 44, 38, 38, 34, 48, 54, 44, 34, 31, 29, 45, 36, 53, 31, 41, 38 ,
34, 42, 40, 35, 38, 39, 29, 54, 39, 40, 35, 50, 35, 35, 40, 41, 33, 38, 33, 36, 47, 40, 33, 38, 48, 40, 42, 33, 25, 29, 39 ,
43, 51, 38, 37, 45, 38, 41, 42, 35, 35, 37, 45, 39, 37, 41, 43, 46, 38, 49, 33, 35, 35, 43, 49, 50, 33, 43, 45, 48, 50, 39, 21, 41, 41, 46 ,
37, 38, 40, 51, 37, 38, 54, 36, 25, 36, 38, 45, 53, 32, 42, 50, 37, 36, 46, 60, 34, 50, 48, 34, 30, 37, 34, 31, 34, 39 ,
37, 31, 37, 47, 38, 39, 42, 38, 40, 33, 41, 42, 35, 36, 43, 33, 42, 49, 41, 38, 39, 36, 48, 39, 33, 35, 32, 25, 35, 38, 32, 45 ,
37, 32, 46, 42, 33, 40, 38, 41, 35, 39, 34, 31, 43, 37, 43, 42, 48, 43, 26, 36, 29, 45, 32, 36, 43, 35, 50, 44, 44, 41, 47 ,
50, 27, 34, 42, 45, 38, 37, 39, 26, 51, 27, 37, 40, 40, 30, 43, 30, 47, 35, 34, 54, 38, 36, 41, 38, 26, 40, 32, 37, 32, 40, 40, 50, 37, 37 ,
39, 27, 37, 44, 49, 49, 38, 48, 56, 35, 42, 41, 38, 34, 40, 27, 36, 40, 31, 38, 40, 29, 37, 37, 37, 43, 34, 39, 41, 41, 37 ,
46, 40, 38, 33, 55, 44, 53, 39, 47, 41, 39, 31, 47, 44, 36, 38, 47, 33, 34, 36, 39, 44, 40, 55, 34, 36, 40, 31, 50, 53, 50 ,
37, 44, 49, 30, 37, 35, 41, 44, 32, 36, 29, 33, 45, 43, 33, 33, 46, 43, 48, 39, 31, 36, 41, 40, 32, 34, 43, 38, 39, 44 ,
38, 40, 40, 41, 37, 41, 47, 33, 33, 42, 37, 45, 32, 38, 35, 46, 36, 38, 36, 44, 26, 40, 42, 53, 41, 44, 37, 32, 40, 35 ,
39, 32, 42, 41, 38, 46, 28, 37, 46, 35, 34, 40, 35, 41, 42, 43, 41, 31, 47, 46, 31, 42, 27, 51, 37, 41, 41, 34, 29, 37, 45)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 , 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5 , 6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6 ,
7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7 , 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8 ,
9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9 , 10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10 ,
11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11,11 ,
12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12)
> mean(x)
[1] 39.2122
> var(x)
[1] 41.30591
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 1.0675, 第1自由度 = 11, 第2自由度 = 365, P値 = 0.3866
有意差あり。群による差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 11.9859, 自由度 = 11, P値 = 0.3647
有意。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.9729, 第1自由度 = 11.000, 第2自由度 = 143.248, P値 = 0.4738
有意差あり。群による差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 1 0.265252 0.265252
25 21 5.570292 5.835544
30 64 16.976127 22.811671
35 117 31.034483 53.846154
40 100 26.525199 80.371353
45 46 12.201592 92.572944
50 24 6.366048 98.938992
55 3 0.795756 99.734748
60 1 0.265252 100.000000
>
43:名無しに影響はない(栃木県)
12/04/14 07:56:50.56 PzUMpoQ1
x <- c(26, 46, 30, 40, 47, 41, 32, 42, 30, 40, 44, 39, 42, 51, 44, 38, 38, 34, 48, 54, 44, 34, 31, 29, 45, 36, 53, 31, 41, 38 ,
34, 42, 40, 35, 38, 39, 29, 54, 39, 40, 35, 50, 35, 35, 40, 41, 33, 38, 33, 36, 47, 40, 33, 38, 48, 40, 42, 33, 25, 29, 39 ,
43, 51, 38, 37, 45, 38, 41, 42, 35, 35, 37, 45, 39, 37, 41, 43, 46, 38, 49, 33, 35, 35, 43, 49, 50, 33, 43, 45, 48, 50, 39, 21, 41, 41, 46 ,
37, 38, 40, 51, 37, 38, 54, 36, 25, 36, 38, 45, 53, 32, 42, 50, 37, 36, 46, 60, 34, 50, 48, 34, 30, 37, 34, 31, 34, 39 ,
37, 31, 37, 47, 38, 39, 42, 38, 40, 33, 41, 42, 35, 36, 43, 33, 42, 49, 41, 38, 39, 36, 48, 39, 33, 35, 32, 25, 35, 38, 32, 45 ,
37, 32, 46, 42, 33, 40, 38, 41, 35, 39, 34, 31, 43, 37, 43, 42, 48, 43, 26, 36, 29, 45, 32, 36, 43, 35, 50, 44, 44, 41, 47 ,
50, 27, 34, 42, 45, 38, 37, 39, 26, 51, 27, 37, 40, 40, 30, 43, 30, 47, 35, 34, 54, 38, 36, 41, 38, 26, 40, 32, 37, 32, 40, 40, 50, 37, 37 ,
39, 27, 37, 44, 49, 49, 38, 48, 56, 35, 42, 41, 38, 34, 40, 27, 36, 40, 31, 38, 40, 29, 37, 37, 37, 43, 34, 39, 41, 41, 37 ,
46, 40, 38, 33, 55, 44, 53, 39, 47, 41, 39, 31, 47, 44, 36, 38, 47, 33, 34, 36, 39, 44, 40, 55, 34, 36, 40, 31, 50, 53, 50 ,
37, 44, 49, 30, 37, 35, 41, 44, 32, 36, 29, 33, 45, 43, 33, 33, 46, 43, 48, 39, 31, 36, 41, 40, 32, 34, 43, 38, 39, 44 ,
38, 40, 40, 41, 37, 41, 47, 33, 33, 42, 37, 45, 32, 38, 35, 46, 36, 38, 36, 44, 26, 40, 42, 53, 41, 44, 37, 32, 40, 35 ,
39, 32, 42, 41, 38, 46, 28, 37, 46, 35, 34, 40, 35, 41, 42, 43, 41, 31, 47, 46, 31, 42, 27, 51, 37, 41, 41, 34, 29, 37, 45)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 , 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5 , 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,
2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 , 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 , 5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.6641, 第1自由度 = 4, 第2自由度 = 372, P値 = 0.03232
有意差あり。群による差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 5.5387, 自由度 = 4, P値 = 0.2364
有意。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.6132, 第1自由度 = 4.000, 第2自由度 = 174.657, P値 = 0.03704
有意差あり。群による差異がある。
44:名無しに影響はない(栃木県)
12/04/14 07:57:54.42 PzUMpoQ1
繰り返しによる差異
1. n=91 BG
x <- c(26, 46, 30, 40, 47, 41, 32, 42, 30, 40, 44, 39, 42, 51, 44, 38, 38, 34, 48, 54, 44, 34, 31, 29, 45, 36, 53, 31, 41, 38 ,
37, 32, 46, 42, 33, 40, 38, 41, 35, 39, 34, 31, 43, 37, 43, 42, 48, 43, 26, 36, 29, 45, 32, 36, 43, 35, 50, 44, 44, 41, 47 ,
38, 40, 40, 41, 37, 41, 47, 33, 33, 42, 37, 45, 32, 38, 35, 46, 36, 38, 36, 44, 26, 40, 42, 53, 41, 44, 37, 32, 40, 35)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x)
[1] 39.21978
> var(x)
[1] 38.21783
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.0862, 第1自由度 = 2, 第2自由度 = 88, P値 = 0.9175
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 3.0282, 自由度 = 2, P値 = 0.22
有意。先の分散分析は無効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.0745, 第1自由度 = 2.000, 第2自由度 = 57.737, P値 = 0.9283
有意ではない。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 5 5.494505 5.494505
30 16 17.582418 23.076923
35 23 25.274725 48.351648
40 31 34.065934 82.417582
45 11 12.087912 94.505495
50 5 5.494505 100.000000
>
45:名無しに影響はない(栃木県)
12/04/14 07:58:22.34 PzUMpoQ1
2. n=97 袋
x <- c(34, 42, 40, 35, 38, 39, 29, 54, 39, 40, 35, 50, 35, 35, 40, 41, 33, 38, 33, 36, 47, 40, 33, 38, 48, 40, 42, 33, 25, 29, 39 ,
50, 27, 34, 42, 45, 38, 37, 39, 26, 51, 27, 37, 40, 40, 30, 43, 30, 47, 35, 34, 54, 38, 36, 41, 38, 26, 40, 32, 37, 32, 40, 40, 50, 37, 37 ,
39, 32, 42, 41, 38, 46, 28, 37, 46, 35, 34, 40, 35, 41, 42, 43, 41, 31, 47, 46, 31, 42, 27, 51, 37, 41, 41, 34, 29, 37, 45)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x)
[1] 38.23711
> var(x)
[1] 40.87027
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.1069, 第1自由度 = 2, 第2自由度 = 94, P値 = 0.8988
有意。繰り返しによる差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.9797, 自由度 = 2, P値 = 0.6127
有意。先の分散分析は無効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.1126, 第1自由度 = 2.000, 第2自由度 = 62.638, P値 = 0.8937
有意。繰り返しによる差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 10 10.309278 10.30928
30 16 16.494845 26.80412
35 29 29.896907 56.70103
40 26 26.804124 83.50515
45 9 9.278351 92.78351
50 7 7.216495 100.00000
>
46:名無しに影響はない(栃木県)
12/04/14 07:58:51.52 PzUMpoQ1
3. n=66 未茹で
x <- c(43, 51, 38, 37, 45, 38, 41, 42, 35, 35, 37, 45, 39, 37, 41, 43, 46, 38, 49, 33, 35, 35, 43, 49, 50, 33, 43, 45, 48, 50, 39, 21, 41, 41, 46 ,
46, 40, 38, 33, 55, 44, 53, 39, 47, 41, 39, 31, 47, 44, 36, 38, 47, 33, 34, 36, 39, 44, 40, 55, 34, 36, 40, 31, 50, 53, 50)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 41.28788
> var(x)
[1] 43.68508
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.2353, 第1自由度 = 1, 第2自由度 = 64, P値 = 0.6293
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.5951, 自由度 = 1, P値 = 0.4404
有意。先の分散分析は無効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.2313, 第1自由度 = 1.000, 第2自由度 = 59.987, P値 = 0.6323
有意ではない。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 1 1.515152 1.515152
25 0 0.000000 1.515152
30 8 12.121212 13.636364
35 20 30.303030 43.939394
40 16 24.242424 68.181818
45 12 18.181818 86.363636
50 7 10.606061 96.969697
55 2 3.030303 100.000000
>
47:名無しに影響はない(栃木県)
12/04/14 07:59:20.24 PzUMpoQ1
4. n=61 茹で１回
x <- c(37, 38, 40, 51, 37, 38, 54, 36, 25, 36, 38, 45, 53, 32, 42, 50, 37, 36, 46, 60, 34, 50, 48, 34, 30, 37, 34, 31, 34, 39 ,
39, 27, 37, 44, 49, 49, 38, 48, 56, 35, 42, 41, 38, 34, 40, 27, 36, 40, 31, 38, 40, 29, 37, 37, 37, 43, 34, 39, 41, 41, 37)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 39.44262
> var(x)
[1] 52.21749
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.4361, 第1自由度 = 1, 第2自由度 = 59, P値 = 0.5116
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 1.9459, 自由度 = 1, P値 = 0.1630
有意。先の分散分析は無効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.4324, 第1自由度 = 1.00, 第2自由度 = 54.52, P値 = 0.5136
有意ではない。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 4 6.557377 6.557377
30 10 16.393443 22.950820
35 23 37.704918 60.655738
40 11 18.032787 78.688525
45 6 9.836066 88.524590
50 5 8.196721 96.721311
55 1 1.639344 98.360656
60 1 1.639344 100.000000
>
48:名無しに影響はない(栃木県)
12/04/14 07:59:44.72 PzUMpoQ1
5. n=62 茹で２回
x <- c(37, 31, 37, 47, 38, 39, 42, 38, 40, 33, 41, 42, 35, 36, 43, 33, 42, 49, 41, 38, 39, 36, 48, 39, 33, 35, 32, 25, 35, 38, 32, 45 ,
37, 44, 49, 30, 37, 35, 41, 44, 32, 36, 29, 33, 45, 43, 33, 33, 46, 43, 48, 39, 31, 36, 41, 40, 32, 34, 43, 38, 39, 44)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 38.29032
> var(x)
[1] 28.89794
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.0871, 第1自由度 = 1, 第2自由度 = 60, P値 = 0.7689
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.1314, 自由度 = 1, P値 = 0.717
有意ではない。先の分散分析は有効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.0867, 第1自由度 = 1.000, 第2自由度 = 58.971, P値 = 0.7694
有意ではない。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 2 3.225806 3.225806
30 14 22.580645 25.806452
35 22 35.483871 61.290323
40 16 25.806452 87.096774
45 8 12.903226 100.000000
>
49:名無しに影響はない(栃木県)
12/04/14 08:00:46.27 PzUMpoQ1
バックグラウンドとの比較
1. n=91 BG
gr1 <- c(26, 46, 30, 40, 47, 41, 32, 42, 30, 40, 44, 39, 42, 51, 44, 38, 38, 34, 48, 54, 44, 34, 31, 29, 45, 36, 53, 31, 41, 38 ,
37, 32, 46, 42, 33, 40, 38, 41, 35, 39, 34, 31, 43, 37, 43, 42, 48, 43, 26, 36, 29, 45, 32, 36, 43, 35, 50, 44, 44, 41, 47 ,
38, 40, 40, 41, 37, 41, 47, 33, 33, 42, 37, 45, 32, 38, 35, 46, 36, 38, 36, 44, 26, 40, 42, 53, 41, 44, 37, 32, 40, 35)
> mean(x)
[1] 39.21978
> var(x)
[1] 38.21783
2. n=97 袋
gr2 <- c(34, 42, 40, 35, 38, 39, 29, 54, 39, 40, 35, 50, 35, 35, 40, 41, 33, 38, 33, 36, 47, 40, 33, 38, 48, 40, 42, 33, 25, 29, 39 ,
50, 27, 34, 42, 45, 38, 37, 39, 26, 51, 27, 37, 40, 40, 30, 43, 30, 47, 35, 34, 54, 38, 36, 41, 38, 26, 40, 32, 37, 32, 40, 40, 50, 37, 37 ,
39, 32, 42, 41, 38, 46, 28, 37, 46, 35, 34, 40, 35, 41, 42, 43, 41, 31, 47, 46, 31, 42, 27, 51, 37, 41, 41, 34, 29, 37, 45)
> mean(x)
[1] 38.23711
> var(x)
[1] 40.87027
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = 1.0702, 自由度 = 186, P値 = 0.2859
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -0.8288018 2.7941355
標本推定値:
平均値x 平均値y
39.21978 38.23711
有意。平均値に差がある。
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = 1.0713, 自由度 = 185.825, P値 = 0.2854
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -0.8268647 2.7921983
標本推定値:
平均値x 平均値y
39.21978 38.23711
有意。平均値に差がある。
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.9351, 第1自由度 = 90, 第2自由度 = 96, P値 = 0.7489
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.6219275 1.4098211
標本推定値:
分散比
0.9351008
有意。（Welchの方法）を使用する。
>
39.2-38.2 = 1
1 * 0.47 = 0.47 Bq/１袋(袋による遮蔽)
50:名無しに影響はない(栃木県)
12/04/14 08:01:16.53 PzUMpoQ1
3. n=66 未茹で
gr2 <- c(43, 51, 38, 37, 45, 38, 41, 42, 35, 35, 37, 45, 39, 37, 41, 43, 46, 38, 49, 33, 35, 35, 43, 49, 50, 33, 43, 45, 48, 50, 39, 21, 41, 41, 46 ,
46, 40, 38, 33, 55, 44, 53, 39, 47, 41, 39, 31, 47, 44, 36, 38, 47, 33, 34, 36, 39, 44, 40, 55, 34, 36, 40, 31, 50, 53, 50)
> mean(x)
[1] 41.28788
> var(x)
[1] 43.68508
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = -2.0097, 自由度 = 155, P値 = 0.0462
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -4.10089286 -0.03530427
標本推定値:
平均値x 平均値y
39.21978 41.28788
有意。平均値に差がある。
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = -1.9883, 自由度 = 134.536, P値 = 0.04881
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -4.12522061 -0.01097653
標本推定値:
平均値x 平均値y
39.21978 41.28788
有意。平均値に差がある。
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.8748, 第1自由度 = 90, 第2自由度 = 65, P値 = 0.553
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.5505676 1.3668892
標本推定値:
分散比
0.8748485
有意。（Welchの方法）を使用する。
>
41.3 - 39.2 = 2.1
2.1 * 0.47 = 0.99
51:名無しに影響はない(栃木県)
12/04/14 09:46:06.53 PzUMpoQ1
連続投稿規制に引っかかったので中断したものを再開。
4. n=61 茹で１回
gr2 <- c(37, 38, 40, 51, 37, 38, 54, 36, 25, 36, 38, 45, 53, 32, 42, 50, 37, 36, 46, 60, 34, 50, 48, 34, 30, 37, 34, 31, 34, 39 ,
39, 27, 37, 44, 49, 49, 38, 48, 56, 35, 42, 41, 38, 34, 40, 27, 36, 40, 31, 38, 40, 29, 37, 37, 37, 43, 34, 39, 41, 41, 37)
> mean(x)
[1] 39.44262
> var(x)
[1] 52.21749
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = -0.2034, 自由度 = 150, P値 = 0.839
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -2.387191 1.941505
標本推定値:
平均値x 平均値y
39.21978 39.44262
有意ではない。平均値に差があるがどうかわからない
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = -0.1973, 自由度 = 114.881, P値 = 0.844
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -2.460391 2.014706
標本推定値:
平均値x 平均値y
39.21978 39.44262
有意ではない。平均値に差があるがどうかわからない
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.7319, 第1自由度 = 90, 第2自由度 = 60, P値 = 0.1781
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.4543745 1.1533524
標本推定値:
分散比
0.7318971
有意。（Welchの方法）を使用する。
>
差が有意ではないので、Bq値を求められない。
39.4-39.2 = 0.2
0.2 * 0.47 = 0.09
52:名無しに影響はない(栃木県)
12/04/14 09:46:43.62 PzUMpoQ1
5. n=62 茹で２回
gr2 <- c(37, 31, 37, 47, 38, 39, 42, 38, 40, 33, 41, 42, 35, 36, 43, 33, 42, 49, 41, 38, 39, 36, 48, 39, 33, 35, 32, 25, 35, 38, 32, 45 ,
37, 44, 49, 30, 37, 35, 41, 44, 32, 36, 29, 33, 45, 43, 33, 33, 46, 43, 48, 39, 31, 36, 41, 40, 32, 34, 43, 38, 39, 44)
> mean(x)
[1] 38.29032
> var(x)
[1] 28.89794
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = 0.9616, 自由度 = 151, P値 = 0.3378
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -0.9803271 2.8392424
標本推定値:
平均値x 平均値y
39.21978 38.29032
有意。平均値に差がある。
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = 0.9874, 自由度 = 142.202, P値 = 0.3251
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -0.9313193 2.7902346
標本推定値:
平均値x 平均値y
39.21978 38.29032
有意。平均値に差がある。
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 1.3225, 第1自由度 = 90, 第2自由度 = 61, P値 = 0.2459
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.8234076 2.0803142
標本推定値:
分散比
1.322511
有意。（Welchの方法）を使用する。
>
39.2-38.2 = 1
1 * 0.47 = 0.47
バックグラウンドとの比較。試料20g
袋*, 未茹で*,茹で１回 ,茹で２回*
===============================================
差(CPS) -1.0 , 2.1 , 0.2 , -1.0
量(Bq ) -0.47, 0.99, 0.09, -0.47
(*印;有意, 無印;有意ではない)
53:名無しに影響はない(栃木県)
12/04/14 09:47:09.40 PzUMpoQ1
回帰分析
y <- c(43, 51, 38, 37, 45, 38, 41, 42, 35, 35, 37, 45, 39, 37, 41, 43, 46, 38, 49, 33, 35, 35, 43, 49, 50, 33, 43, 45, 48, 50, 39, 21, 41, 41, 46 ,
37, 38, 40, 51, 37, 38, 54, 36, 25, 36, 38, 45, 53, 32, 42, 50, 37, 36, 46, 60, 34, 50, 48, 34, 30, 37, 34, 31, 34, 39 ,
37, 31, 37, 47, 38, 39, 42, 38, 40, 33, 41, 42, 35, 36, 43, 33, 42, 49, 41, 38, 39, 36, 48, 39, 33, 35, 32, 25, 35, 38, 32, 45 ,
39, 27, 37, 44, 49, 49, 38, 48, 56, 35, 42, 41, 38, 34, 40, 27, 36, 40, 31, 38, 40, 29, 37, 37, 37, 43, 34, 39, 41, 41, 37 ,
46, 40, 38, 33, 55, 44, 53, 39, 47, 41, 39, 31, 47, 44, 36, 38, 47, 33, 34, 36, 39, 44, 40, 55, 34, 36, 40, 31, 50, 53, 50 ,
37, 44, 49, 30, 37, 35, 41, 44, 32, 36, 29, 33, 45, 43, 33, 33, 46, 43, 48, 39, 31, 36, 41, 40, 32, 34, 43, 38, 39, 44)
x <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 , 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,
2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 , 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> result<-lm(y~x)
> result
Call:
lm(formula = y ~ x)
Coefficients:
(Intercept) x
41.179 -1.502
> plot(x,y,main="Regression line")
> abline(result,col="red")
>
３.補足コメント
インスペクターによるやさしおの測定では、1.819個 / (分・100mg)なので、インスペクターの計数値からの換算計数は 0.47 Bq/CPM となる。
1.502 * 0.47 = 0.70
0.70 Bq/20g・回 減少する。
なお、べつに同様な事を行ったものがある。これは、食べる予定のうどんからキッチンペーパーに入る程度の分量をサンプリングした場合で、
２回茹でたほう(41.8-38.4=3.4 CPM, 12.4g)が、１回(41.7-38.4=1.9 CPM 33.0g)だけよりも高濃度という結果になった。
茹で方に問題があるのか、水道水が汚染されていたか、乾燥中に汚染されたか、等の原因が考えられる。
一番可能性のあるのは、攪拌不備による濃度分布(一部分高濃度なところがあり、高濃度の部分の試料を採取)で、
少量(50g以下)の実験室レベルの実験からテストプラント(100kg以下)に移行するときに現れる問題です。
バックグラウントと袋の比較で、線量の減少が観察された。キッチンペーパー6枚(3回重ね折り)+ポリエチレン 0.04mm３枚で遮蔽可能な放射線はα線だけである。
ポリエチレンはゴムを除くと比較的水素が通りやすい膜である。つまり、α線が透過できることを示している。
54:名無しに影響はない(栃木県)
12/04/14 13:49:48.14 PzUMpoQ1
１．測定対象
「アルミホイル(厚さ 11um, 中国製)」の測定。
インスペクターによる食品測定を目的に色々やっていたが、使用している機材が汚染されているらしく
α線らしきものを検出した。>>39-53 α線かどうかの確認作業を行った。
使用している機材は、以下のとおり。
台、床面より45cmの高さ、35*40cm。この上にダンボール製箱(35*40cm)を置き、上に鉄板(厚さ1mm)をすいて机様に加工。
ダンボール内に水入りポリタンク(URLﾘﾝｸ(item.rakuten.co.jp))を置いてある。
マンテンの棚用の鋼材を井桁に組み、インスペクターを置くと下に約3cmの空間がとれるようにした。
インスペクターは汚染を防ぐために常に食品にも使えるポリ袋(ユニパック I-4, ポリエチレン 0.04mm厚)の中に保管した。
ポリエチレンフィルムのガスの透過性については下記参照。
URLﾘﾝｸ(www.agc.com)
URLﾘﾝｸ(www.as-1.co.jp)
アルミホイルは、インスペクターを包んだポリ袋の外側、井桁の上にすいた。
３.補足コメント
有意差が得られず、α線検出とはならなかった。
紙６枚でβ線を遮蔽できるのか、ちょっとわけがわからない状態。
55:名無しに影響はない(栃木県)
12/04/14 13:53:05.50 PzUMpoQ1
２．測定結果
生データは以下のとおり。連続投稿禁止の制限に引っかかって途中でちぎれるかも。
1. n=30 BG
x <- c(36, 45, 46, 49, 34, 58, 37, 40, 39, 47, 44, 47, 45, 54, 66, 29, 34, 47, 39, 59, 43, 37, 44, 43, 44, 39, 54, 44, 33, 41)
> mean(x)
[1] 43.9
> var(x)
[1] 68.02414
>
2. n=35 試料
x <- c(42, 42, 49, 45, 42, 41, 39, 43, 42, 37, 35, 44, 42, 44, 41, 31, 47, 45, 40, 33, 35, 35, 45, 49, 38, 40, 50, 47, 36, 34, 39, 42, 32, 33, 47)
> mean(x)
[1] 40.74286
> var(x)
[1] 26.96134
>
3. n=33 BG
x <- c(38, 32, 44, 51, 43, 40, 38, 22, 39, 40, 38, 42, 37, 32, 24, 35, 37, 44, 35, 51, 37, 46, 39, 40, 42, 44, 47, 43, 47, 50, 43, 32, 38)
> mean(x)
[1] 39.69697
> var(x)
[1] 44.4678
>
4. n=33 試料
x <- c(40, 38, 41, 40, 50, 50, 42, 39, 41, 45, 40, 44, 40, 36, 42, 36, 40, 27, 35, 40, 36, 51, 34, 37, 25, 47, 43, 36, 36, 40, 40, 41, 43)
> mean(x)
[1] 39.84848
> var(x)
[1] 30.88258
>
5. n=33 BG
x <- c(34, 41, 31, 44, 38, 43, 39, 42, 44, 46, 41, 48, 33, 29, 36, 33, 41, 27, 36, 46, 36, 26, 38, 33, 51, 30, 46, 42, 35, 41, 39, 35, 27)
> mean(x)
[1] 37.90909
> var(x)
[1] 41.52273
>
56:名無しに影響はない(栃木県)
12/04/14 13:54:55.77 PzUMpoQ1
４.データ貼り付け
全体の分析
x <- c(36, 45, 46, 49, 34, 58, 37, 40, 39, 47, 44, 47, 45, 54, 66, 29, 34, 47, 39, 59, 43, 37, 44, 43, 44, 39, 54, 44, 33, 41 ,
42, 42, 49, 45, 42, 41, 39, 43, 42, 37, 35, 44, 42, 44, 41, 31, 47, 45, 40, 33, 35, 35, 45, 49, 38, 40, 50, 47, 36, 34, 39, 42, 32, 33, 47 ,
38, 32, 44, 51, 43, 40, 38, 22, 39, 40, 38, 42, 37, 32, 24, 35, 37, 44, 35, 51, 37, 46, 39, 40, 42, 44, 47, 43, 47, 50, 43, 32, 38 ,
40, 38, 41, 40, 50, 50, 42, 39, 41, 45, 40, 44, 40, 36, 42, 36, 40, 27, 35, 40, 36, 51, 34, 37, 25, 47, 43, 36, 36, 40, 40, 41, 43 ,
34, 41, 31, 44, 38, 43, 39, 42, 44, 46, 41, 48, 33, 29, 36, 33, 41, 27, 36, 46, 36, 26, 38, 33, 51, 30, 46, 42, 35, 41, 39, 35, 27)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 , 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5)
> mean(x)
[1] 40.35976
> var(x)
[1] 44.36671
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 3.6123, 第1自由度 = 4, 第2自由度 = 159, P値 = 0.007568
有意。群別の差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 8.1988, 自由度 = 4, P値 = 0.08456
有意。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.6566, 第1自由度 = 4.000, 第2自由度 = 78.113, P値 = 0.03896
有意。群別の差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 2 1.2195122 1.219512
25 7 4.2682927 5.487805
30 18 10.9756098 16.463415
35 43 26.2195122 42.682927
40 57 34.7560976 77.439024
45 24 14.6341463 92.073171
50 10 6.0975610 98.170732
55 2 1.2195122 99.390244
60 0 0.0000000 99.390244
65 1 0.6097561 100.000000
>
57:名無しに影響はない(栃木県)
12/04/14 13:56:03.19 PzUMpoQ1
バックグラウンドと試料の比較
x <- c(36, 45, 46, 49, 34, 58, 37, 40, 39, 47, 44, 47, 45, 54, 66, 29, 34, 47, 39, 59, 43, 37, 44, 43, 44, 39, 54, 44, 33, 41 ,
42, 42, 49, 45, 42, 41, 39, 43, 42, 37, 35, 44, 42, 44, 41, 31, 47, 45, 40, 33, 35, 35, 45, 49, 38, 40, 50, 47, 36, 34, 39, 42, 32, 33, 47 ,
38, 32, 44, 51, 43, 40, 38, 22, 39, 40, 38, 42, 37, 32, 24, 35, 37, 44, 35, 51, 37, 46, 39, 40, 42, 44, 47, 43, 47, 50, 43, 32, 38 ,
40, 38, 41, 40, 50, 50, 42, 39, 41, 45, 40, 44, 40, 36, 42, 36, 40, 27, 35, 40, 36, 51, 34, 37, 25, 47, 43, 36, 36, 40, 40, 41, 43 ,
34, 41, 31, 44, 38, 43, 39, 42, 44, 46, 41, 48, 33, 29, 36, 33, 41, 27, 36, 46, 36, 26, 38, 33, 51, 30, 46, 42, 35, 41, 39, 35, 27)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.0068, 第1自由度 = 1, 第2自由度 = 162, P値 = 0.9346
有意ではない。バックグラウンドと試料の差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 8.2786, 自由度 = 1, P値 = 0.004012
有意。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.0075, 第1自由度 = 1.000, 第2自由度 = 161.976, P値 = 0.9309
有意ではない。バックグラウンドと試料の差異は不明。
>
58:名無しに影響はない(栃木県)
12/04/14 13:56:49.19 PzUMpoQ1
繰り返しの比較
1. n=96 BG
x <- c(36, 45, 46, 49, 34, 58, 37, 40, 39, 47, 44, 47, 45, 54, 66, 29, 34, 47, 39, 59, 43, 37, 44, 43, 44, 39, 54, 44, 33, 41 ,
38, 32, 44, 51, 43, 40, 38, 22, 39, 40, 38, 42, 37, 32, 24, 35, 37, 44, 35, 51, 37, 46, 39, 40, 42, 44, 47, 43, 47, 50, 43, 32, 38 ,
34, 41, 31, 44, 38, 43, 39, 42, 44, 46, 41, 48, 33, 29, 36, 33, 41, 27, 36, 46, 36, 26, 38, 33, 51, 30, 46, 42, 35, 41, 39, 35, 27)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x)
[1] 40.39583
> var(x)
[1] 55.92588
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 5.7929, 第1自由度 = 2, 第2自由度 = 93, P値 = 0.004256
有意。繰り返しによる差異ある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 2.2143, 自由度 = 2, P値 = 0.3305
先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 5.0623, 第1自由度 = 2.000, 第2自由度 = 60.331, P値 = 0.00928
有意。繰り返しによる差異ある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 2 2.083333 2.083333
25 5 5.208333 7.291667
30 12 12.500000 19.791667
35 26 27.083333 46.875000
40 28 29.166667 76.041667
45 14 14.583333 90.625000
50 6 6.250000 96.875000
55 2 2.083333 98.958333
60 0 0.000000 98.958333
65 1 1.041667 100.000000
>
59:名無しに影響はない(栃木県)
12/04/14 13:57:34.87 PzUMpoQ1
2. n=68 試料
x <- c(42, 42, 49, 45, 42, 41, 39, 43, 42, 37, 35, 44, 42, 44, 41, 31, 47, 45, 40, 33, 35, 35, 45, 49, 38, 40, 50, 47, 36, 34, 39, 42, 32, 33, 47 ,
40, 38, 41, 40, 50, 50, 42, 39, 41, 45, 40, 44, 40, 36, 42, 36, 40, 27, 35, 40, 36, 51, 34, 37, 25, 47, 43, 36, 36, 40, 40, 41, 43)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 40.30882
> var(x)
[1] 28.63455
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.4707, 第1自由度 = 1, 第2自由度 = 66, P値 = 0.4951
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.1498, 自由度 = 1, P値 = 0.6987
有意。先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.4688, 第1自由度 = 1.000, 第2自由度 = 64.946, P値 = 0.496
有意ではない。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 2 2.941176 2.941176
30 6 8.823529 11.764706
35 17 25.000000 36.764706
40 29 42.647059 79.411765
45 10 14.705882 94.117647
50 4 5.882353 100.000000
>
60:名無しに影響はない(栃木県)
12/04/14 13:58:10.33 PzUMpoQ1
バックグラウンドと試料の比較
1. n=96 BG
gr1 <- c(36, 45, 46, 49, 34, 58, 37, 40, 39, 47, 44, 47, 45, 54, 66, 29, 34, 47, 39, 59, 43, 37, 44, 43, 44, 39, 54, 44, 33, 41 ,
38, 32, 44, 51, 43, 40, 38, 22, 39, 40, 38, 42, 37, 32, 24, 35, 37, 44, 35, 51, 37, 46, 39, 40, 42, 44, 47, 43, 47, 50, 43, 32, 38 ,
34, 41, 31, 44, 38, 43, 39, 42, 44, 46, 41, 48, 33, 29, 36, 33, 41, 27, 36, 46, 36, 26, 38, 33, 51, 30, 46, 42, 35, 41, 39, 35, 27)
> mean(x)
[1] 40.39583
> var(x)
[1] 55.92588
2. n=68 試料
gr2 <- c(42, 42, 49, 45, 42, 41, 39, 43, 42, 37, 35, 44, 42, 44, 41, 31, 47, 45, 40, 33, 35, 35, 45, 49, 38, 40, 50, 47, 36, 34, 39, 42, 32, 33, 47 ,
40, 38, 41, 40, 50, 50, 42, 39, 41, 45, 40, 44, 40, 36, 42, 36, 40, 27, 35, 40, 36, 51, 34, 37, 25, 47, 43, 36, 36, 40, 40, 41, 43)
> mean(x)
[1] 40.30882
> var(x)
[1] 28.63455
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = 0.0822, 自由度 = 162, P値 = 0.9346
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -2.004175 2.178194
標本推定値:
平均値x 平均値y
40.39583 40.30882
有意ではない。差があるかどうか不明
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = 0.0869, 自由度 = 161.976, P値 = 0.9309
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -1.891316 2.065336
標本推定値:
平均値x 平均値y
40.39583 40.30882
有意ではない。差があるかどうか不明
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 1.9531, 第1自由度 = 95, 第2自由度 = 67, P値 = 0.004175
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 1.240164 3.022964
標本推定値:
分散比
1.953091
有意。分散が一様ではないので（Welchの方法）を使用する
>
61:名無しに影響はない(やわらか銀行)
12/04/14 15:50:01.69 RDG7UOvf
>>54
３.補足コメント
有意差が得られず、α線検出とはならなかった。
紙６枚でβ線を遮蔽できるのか、ちょっとわけがわからない状態。
--
お疲れ様です。
紙一枚でα線を遮断できます。
紙を６枚以上重ねると低エネルギーの一部のβ線が遮断されるはずです。
62:名無しに影響はない(栃木県)
12/04/14 20:14:16.88 PzUMpoQ1
>>61
汚染されていない紙があればよいのですが。
スーパーで輸入物の紙を探してみます。
さて今回は、
１．測定対象
「しいたけ(2012.04.09 自宅裏より収穫、４日間乾燥)」 の分析結果
２．測定結果
まずは生データ。
1. n=32 BG
x <- c(39, 38, 36, 39, 39, 35, 30, 39, 40, 27, 41, 38, 38, 49, 46, 40, 38, 33, 44, 41, 40, 49, 52, 28, 34, 40, 42, 42, 43, 30, 38, 41)
2. n=41 試料
x <- c(67, 58, 48, 47, 71, 62, 50, 52, 55, 58, 57, 44, 58, 53, 56, 62, 32, 60, 45, 67, 48, 49, 53, 47, 36, 66, 53, 54, 54, 55, 50, 61, 60, 59, 42, 50, 56, 50, 53, 31, 47)
3. n=39 BG
x <- c(44, 42, 50, 32, 46, 48, 44, 29, 35, 52, 43, 46, 33, 54, 35, 48, 35, 43, 43, 42, 42, 37, 30, 40, 51, 30, 50, 36, 41, 46, 52, 45, 31, 45, 33, 41, 34, 36, 31)
4. n=31 試料
x <- c(56, 59, 47, 49, 51, 67, 51, 43, 49, 49, 51, 50, 54, 58, 47, 58, 58, 63, 50, 52, 50, 54, 56, 60, 46, 49, 60, 57, 64, 67, 50)
5. n=31 BG
x <- c(47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45)
63:名無しに影響はない(栃木県)
12/04/14 20:15:26.65 PzUMpoQ1
４. データ貼り付け
全体の分析
x <- c(39, 38, 36, 39, 39, 35, 30, 39, 40, 27, 41, 38, 38, 49, 46, 40, 38, 33, 44, 41, 40, 49, 52, 28, 34, 40, 42, 42, 43, 30, 38, 41 ,
67, 58, 48, 47, 71, 62, 50, 52, 55, 58, 57, 44, 58, 53, 56, 62, 32, 60, 45, 67, 48, 49, 53, 47, 36, 66, 53, 54, 54, 55, 50, 61, 60, 59, 42, 50, 56, 50, 53, 31, 47 ,
44, 42, 50, 32, 46, 48, 44, 29, 35, 52, 43, 46, 33, 54, 35, 48, 35, 43, 43, 42, 42, 37, 30, 40, 51, 30, 50, 36, 41, 46, 52, 45, 31, 45, 33, 41, 34, 36, 31 ,
56, 59, 47, 49, 51, 67, 51, 43, 49, 49, 51, 50, 54, 58, 47, 58, 58, 63, 50, 52, 50, 54, 56, 60, 46, 49, 60, 57, 64, 67, 50 ,
47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 ,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5)
> mean(x)
[1] 45.62644
> var(x)
[1] 90.80184
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 40.26, 第1自由度 = 4, 第2自由度 = 169, P値 < 2.2e-16
有意。群別の差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 10.6791, 自由度 = 4, P値 = 0.03042
有意。先の分散分析が無効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 41.9037, 第1自由度 = 4.000, 第2自由度 = 83.865, P値 < 2.2e-16
有意。群別の差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 3 1.7241379 1.724138
30 20 11.4942529 13.218391
35 26 14.9425287 28.160920
40 37 21.2643678 49.425287
45 28 16.0919540 65.517241
50 28 16.0919540 81.609195
55 17 9.7701149 91.379310
60 9 5.1724138 96.551724
65 5 2.8735632 99.425287
70 1 0.5747126 100.000000
>
64:名無しに影響はない(栃木県)
12/04/14 20:15:59.17 PzUMpoQ1
BGと試料の比較
x <- c(39, 38, 36, 39, 39, 35, 30, 39, 40, 27, 41, 38, 38, 49, 46, 40, 38, 33, 44, 41, 40, 49, 52, 28, 34, 40, 42, 42, 43, 30, 38, 41 ,
67, 58, 48, 47, 71, 62, 50, 52, 55, 58, 57, 44, 58, 53, 56, 62, 32, 60, 45, 67, 48, 49, 53, 47, 36, 66, 53, 54, 54, 55, 50, 61, 60, 59, 42, 50, 56, 50, 53, 31, 47 ,
44, 42, 50, 32, 46, 48, 44, 29, 35, 52, 43, 46, 33, 54, 35, 48, 35, 43, 43, 42, 42, 37, 30, 40, 51, 30, 50, 36, 41, 46, 52, 45, 31, 45, 33, 41, 34, 36, 31 ,
56, 59, 47, 49, 51, 67, 51, 43, 49, 49, 51, 50, 54, 58, 47, 58, 58, 63, 50, 52, 50, 54, 56, 60, 46, 49, 60, 57, 64, 67, 50 ,
47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 160.6906, 第1自由度 = 1, 第2自由度 = 172, P値 < 2.2e-16
有意。バックグラウンドと試料との差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 4.0066, 自由度 = 1, P値 = 0.04532
有意。先の分散分析が無効
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 149.1595, 第1自由度 = 1.000, 第2自由度 = 131.373, P値 < 2.2e-16
有意。バックグラウンドと試料との差異がある。
>
65:名無しに影響はない(栃木県)
12/04/14 20:16:37.61 PzUMpoQ1
繰り返しによる影響
1. n = 102 BG 繰り返し数 =3
x <- c( 39, 38, 36, 39, 39, 35, 30, 39, 40, 27, 41, 38, 38, 49, 46, 40, 38, 33, 44, 41, 40, 49, 52, 28, 34, 40, 42, 42, 43, 30, 38, 41 ,
44, 42, 50, 32, 46, 48, 44, 29, 35, 52, 43, 46, 33, 54, 35, 48, 35, 43, 43, 42, 42, 37, 30, 40, 51, 30, 50, 36, 41, 46, 52, 45, 31, 45, 33, 41, 34, 36, 31 ,
47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x)
[1] 40.07843
> var(x)
[1] 38.52844
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.7926, 第1自由度 = 2, 第2自由度 = 99, P値 = 0.4555
有意。繰り返しによる差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 2.9036, 自由度 = 2, P値 = 0.2341
有意。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.7678, 第1自由度 = 2.000, 第2自由度 = 65.874, P値 = 0.4681
有意。繰り返しによる差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 3 2.9411765 2.941176
30 18 17.6470588 20.588235
35 25 24.5098039 45.098039
40 34 33.3333333 78.431373
45 14 13.7254902 92.156863
50 7 6.8627451 99.019608
55 1 0.9803922 100.000000
>
66:名無しに影響はない(栃木県)
12/04/14 20:17:12.22 PzUMpoQ1
2. n = 72 試料 13.4g 繰り返し数 =2
x <- c( 67, 58, 48, 47, 71, 62, 50, 52, 55, 58, 57, 44, 58, 53, 56, 62, 32, 60, 45, 67, 48, 49, 53, 47, 36, 66, 53, 54, 54, 55, 50, 61, 60, 59, 42, 50, 56, 50, 53, 31, 47 ,
56, 59, 47, 49, 51, 67, 51, 43, 49, 49, 51, 50, 54, 58, 47, 58, 58, 63, 50, 52, 50, 54, 56, 60, 46, 49, 60, 57, 64, 67, 50)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 53.48611
> var(x)
[1] 59.57727
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.2697, 第1自由度 = 1, 第2自由度 = 70, P値 = 0.6051
有意ではない。繰り返しによる差異が不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 3.7321, 自由度 = 1, P値 = 0.05338
有意。先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.2959, 第1自由度 = 1.000, 第2自由度 = 69.781, P値 = 0.5882
有意ではない。繰り返しによる差異が不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
30 2 2.777778 2.777778
35 1 1.388889 4.166667
40 3 4.166667 8.333333
45 14 19.444444 27.777778
50 21 29.166667 56.944444
55 16 22.222222 79.166667
60 9 12.500000 91.666667
65 5 6.944444 98.611111
70 1 1.388889 100.000000
>
67:名無しに影響はない(栃木県)
12/04/14 20:18:27.54 PzUMpoQ1
BGと試料の比較
1. n = 102 BG 繰り返し数 =3
gr1 <- c( 39, 38, 36, 39, 39, 35, 30, 39, 40, 27, 41, 38, 38, 49, 46, 40, 38, 33, 44, 41, 40, 49, 52, 28, 34, 40, 42, 42, 43, 30, 38, 41 ,
44, 42, 50, 32, 46, 48, 44, 29, 35, 52, 43, 46, 33, 54, 35, 48, 35, 43, 43, 42, 42, 37, 30, 40, 51, 30, 50, 36, 41, 46, 52, 45, 31, 45, 33, 41, 34, 36, 31 ,
47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45)
> mean(x)
[1] 40.07843
> var(x)
[1] 38.52844
2. n = 72 試料 13.4g 繰り返し数 =2
gr2 <- c( 67, 58, 48, 47, 71, 62, 50, 52, 55, 58, 57, 44, 58, 53, 56, 62, 32, 60, 45, 67, 48, 49, 53, 47, 36, 66, 53, 54, 54, 55, 50, 61, 60, 59, 42, 50, 56, 50, 53, 31, 47 ,
56, 59, 47, 49, 51, 67, 51, 43, 49, 49, 51, 50, 54, 58, 47, 58, 58, 63, 50, 52, 50, 54, 56, 60, 46, 49, 60, 57, 64, 67, 50)
> mean(x)
[1] 53.48611
> var(x)
[1] 59.57727
68:名無しに影響はない(栃木県)
12/04/14 20:18:58.39 PzUMpoQ1
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = -12.6764, 自由度 = 172, P値 < 2.2e-16
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -15.49540 -11.31996
標本推定値:
平均値x 平均値y
40.07843 53.48611
有意。平均値に差異がある。
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = -12.2131, 自由度 = 131.373, P値 < 2.2e-16
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -15.57936 -11.23600
標本推定値:
平均値x 平均値y
40.07843 53.48611
有意。平均値に差異がある。
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.6467, 第1自由度 = 101, 第2自由度 = 71, P値 = 0.04395
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.4162745 0.9882415
標本推定値:
分散比
0.646697
有意。分散に差異がある。（Welchの方法）を使用する
>
53.5 - 40.1 = 13.4 CPM
13.4 * 0.47 = 6.3 Bq
試料の重さが 13.4g だから、試料による吸収がない、すべての光が受光可能であると近似して
13.4 * 1000 / 13.4 = 1000 Bq/kg
栃木県産(家庭内生産)乾ししいたけの線量は 1000 Bq/kg となった。
しいたけは、91%が水分なので( URLﾘﾝｸ(www.kanesada.co.jp) )
1000 * 9 / 100 = 90 Bq/kg が生しいたけの線量となる。
３.補足コメント
注意、乾燥温度が40度と低い(食品乾燥温度 １３５℃±２℃, ２時間)ので、乾燥不十分により、線量が増えます。
69:名無しに影響はない(栃木県)
12/04/15 12:14:28.89 T3LGXiu9
１．測定対象
「カットわかめ 2012.12.17 中国原産 KKジャパンスパイス」 の測定。
(15 * 10 * 2, 40g) 0.26 g/cm2
２．測定結果
1. n=31 BG
x <- c(25, 37, 46, 41, 38, 38, 46, 36, 43, 36, 47, 40, 34, 28, 44, 52, 38, 37, 35, 37, 25, 30, 32, 35, 38, 37, 40, 52, 42, 40, 33)
> mean(x) [1] 38.12903
> var(x) [1] 43.5828
2. n=31 試料
x <- c(33, 35, 48, 40, 33, 48, 39, 27, 49, 43, 41, 39, 31, 49, 44, 47, 46, 48, 25, 47, 37, 41, 33, 58, 64, 60, 51, 42, 49, 36, 45)
> mean(x) [1] 42.83871
> var(x) [1] 82.33978
3. n=42 BG
x <- c(53, 42, 46, 38, 41, 33, 56, 37, 45, 39, 40, 37, 47, 42, 35, 36, 37, 43, 44, 36, 49, 41, 41, 39, 42, 25, 44, 37, 49, 40, 54, 31, 48, 38, 37, 50, 41, 32, 36, 38, 36, 63)
> mean(x) [1] 41.38095
> var(x) [1] 51.9489
4. n=32 試料
x <- c(35, 37, 42, 45, 34, 35, 53, 42, 41, 43, 38, 48, 46, 45, 49, 36, 38, 43, 40, 47, 38, 38, 41, 41, 38, 40, 35, 37, 50, 44, 54, 51)
> mean(x) [1] 42
> var(x) [1] 30.58065
5. n=37 BG
x <- c(39, 37, 36, 48, 45, 34, 30, 41, 31, 35, 35, 34, 43, 28, 36, 43, 44, 34, 32, 48, 43, 47, 49, 38, 41, 39, 40, 44, 38, 37, 32, 38, 39, 39, 34, 42, 52)
> mean(x) [1] 39.05405
> var(x) [1] 32.55255
70:名無しに影響はない(栃木県)
12/04/15 12:15:55.09 T3LGXiu9
４. データ貼り付け
全体の比較
x <- c(25, 37, 46, 41, 38, 38, 46, 36, 43, 36, 47, 40, 34, 28, 44, 52, 38, 37, 35, 37, 25, 30, 32, 35, 38, 37, 40, 52, 42, 40, 33 ,
33, 35, 48, 40, 33, 48, 39, 27, 49, 43, 41, 39, 31, 49, 44, 47, 46, 48, 25, 47, 37, 41, 33, 58, 64, 60, 51, 42, 49, 36, 45 ,
53, 42, 46, 38, 41, 33, 56, 37, 45, 39, 40, 37, 47, 42, 35, 36, 37, 43, 44, 36, 49, 41, 41, 39, 42, 25, 44, 37, 49, 40, 54, 31, 48, 38, 37, 50, 41, 32, 36, 38, 36, 63 ,
35, 37, 42, 45, 34, 35, 53, 42, 41, 43, 38, 48, 46, 45, 49, 36, 38, 43, 40, 47, 38, 38, 41, 41, 38, 40, 35, 37, 50, 44, 54, 51 ,
39, 37, 36, 48, 45, 34, 30, 41, 31, 35, 35, 34, 43, 28, 36, 43, 44, 34, 32, 48, 43, 47, 49, 38, 41, 39, 40, 44, 38, 37, 32, 38, 39, 39, 34, 42, 52)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 , 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5)
> mean(x) [1] 40.6763
> var(x) [1] 49.69694
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.7227, 第1自由度 = 4, 第2自由度 = 168, P値 = 0.03125
有意差あり。群別の差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 10.5313, 自由度 = 4, P値 = 0.03237
有意。分散の差異があり先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.784, 第1自由度 = 4.000, 第2自由度 = 81.505, P値 = 0.03193
有意差あり。群別の差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 7 4.046243 4.046243
30 20 11.560694 15.606936
35 56 32.369942 47.976879
40 44 25.433526 73.410405
45 30 17.341040 90.751445
50 11 6.358382 97.109827
55 2 1.156069 98.265896
60 3 1.734104 100.000000
>
71:名無しに影響はない(栃木県)
12/04/15 12:16:28.92 T3LGXiu9
BGと試料の比較
x <- c(25, 37, 46, 41, 38, 38, 46, 36, 43, 36, 47, 40, 34, 28, 44, 52, 38, 37, 35, 37, 25, 30, 32, 35, 38, 37, 40, 52, 42, 40, 33 ,
33, 35, 48, 40, 33, 48, 39, 27, 49, 43, 41, 39, 31, 49, 44, 47, 46, 48, 25, 47, 37, 41, 33, 58, 64, 60, 51, 42, 49, 36, 45 ,
53, 42, 46, 38, 41, 33, 56, 37, 45, 39, 40, 37, 47, 42, 35, 36, 37, 43, 44, 36, 49, 41, 41, 39, 42, 25, 44, 37, 49, 40, 54, 31, 48, 38, 37, 50, 41, 32, 36, 38, 36, 63 ,
35, 37, 42, 45, 34, 35, 53, 42, 41, 43, 38, 48, 46, 45, 49, 36, 38, 43, 40, 47, 38, 38, 41, 41, 38, 40, 35, 37, 50, 44, 54, 51 ,
39, 37, 36, 48, 45, 34, 30, 41, 31, 35, 35, 34, 43, 28, 36, 43, 44, 34, 32, 48, 43, 47, 49, 38, 41, 39, 40, 44, 38, 37, 32, 38, 39, 39, 34, 42, 52)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 6.1927, 第1自由度 = 1, 第2自由度 = 171, P値 = 0.01378
有意差あり。試料とバックグラウンドに差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 1.002, 自由度 = 1, P値 = 0.3168
有意。試料とバックグラウンドに分散の差異があり、先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 5.8266, 第1自由度 = 1.000, 第2自由度 = 117.735, P値 = 0.01733
有意差あり。試料とバックグラウンドに差異がある。
72:名無しに影響はない(栃木県)
12/04/15 12:16:52.26 T3LGXiu9
繰り返しによる差異
1. BG
x <- c(25, 37, 46, 41, 38, 38, 46, 36, 43, 36, 47, 40, 34, 28, 44, 52, 38, 37, 35, 37, 25, 30, 32, 35, 38, 37, 40, 52, 42, 40, 33 ,
53, 42, 46, 38, 41, 33, 56, 37, 45, 39, 40, 37, 47, 42, 35, 36, 37, 43, 44, 36, 49, 41, 41, 39, 42, 25, 44, 37, 49, 40, 54, 31, 48, 38, 37, 50, 41, 32, 36, 38, 36, 63 ,
39, 37, 36, 48, 45, 34, 30, 41, 31, 35, 35, 34, 43, 28, 36, 43, 44, 34, 32, 48, 43, 47, 49, 38, 41, 39, 40, 44, 38, 37, 32, 38, 39, 39, 34, 42, 52)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x) [1] 39.68182
> var(x) [1] 44.21893
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.4442, 第1自由度 = 2, 第2自由度 = 107, P値 = 0.09162
有意差あり。繰り返しによる差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 2.0317, 自由度 = 2, P値 = 0.3621
有意差あり。分散の差異があり、先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.1902, 第1自由度 = 2.000, 第2自由度 = 68.786, P値 = 0.1196
有意差あり。繰り返しによる差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 5 4.5454545 4.545455
30 15 13.6363636 18.181818
35 40 36.3636364 54.545455
40 28 25.4545455 80.000000
45 14 12.7272727 92.727273
50 6 5.4545455 98.181818
55 1 0.9090909 99.090909
60 1 0.9090909 100.000000
>
73:名無しに影響はない(栃木県)
12/04/15 12:17:18.74 T3LGXiu9
2. 試料
x <- c(33, 35, 48, 40, 33, 48, 39, 27, 49, 43, 41, 39, 31, 49, 44, 47, 46, 48, 25, 47, 37, 41, 33, 58, 64, 60, 51, 42, 49, 36, 45 ,
35, 37, 42, 45, 34, 35, 53, 42, 41, 43, 38, 48, 46, 45, 49, 36, 38, 43, 40, 47, 38, 38, 41, 41, 38, 40, 35, 37, 50, 44, 54, 51)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 , 2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x) [1] 42.4127
> var(x) [1] 55.3108
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.1977, 第1自由度 = 1, 第2自由度 = 61, P値 = 0.6582
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 7.112, 自由度 = 1, P値 = 0.007657
有意。繰り返しによる分散の差異があり先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.1948, 第1自由度 = 1.000, 第2自由度 = 49.295, P値 = 0.6609
有意。繰り返しによる差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 2 3.174603 3.174603
30 5 7.936508 11.111111
35 16 25.396825 36.507937
40 16 25.396825 61.904762
45 16 25.396825 87.301587
50 5 7.936508 95.238095
55 1 1.587302 96.825397
60 2 3.174603 100.000000
>
74:名無しに影響はない(栃木県)
12/04/15 12:17:42.26 T3LGXiu9
試料とバックグラウンドの比較
1. BG
gr1 <- c(25, 37, 46, 41, 38, 38, 46, 36, 43, 36, 47, 40, 34, 28, 44, 52, 38, 37, 35, 37, 25, 30, 32, 35, 38, 37, 40, 52, 42, 40, 33 ,
53, 42, 46, 38, 41, 33, 56, 37, 45, 39, 40, 37, 47, 42, 35, 36, 37, 43, 44, 36, 49, 41, 41, 39, 42, 25, 44, 37, 49, 40, 54, 31, 48, 38, 37, 50, 41, 32, 36, 38, 36, 63 ,
39, 37, 36, 48, 45, 34, 30, 41, 31, 35, 35, 34, 43, 28, 36, 43, 44, 34, 32, 48, 43, 47, 49, 38, 41, 39, 40, 44, 38, 37, 32, 38, 39, 39, 34, 42, 52)
> mean(x) [1] 39.68182
> var(x) [1] 44.21893
2. 試料
gr2 <- c(33, 35, 48, 40, 33, 48, 39, 27, 49, 43, 41, 39, 31, 49, 44, 47, 46, 48, 25, 47, 37, 41, 33, 58, 64, 60, 51, 42, 49, 36, 45 ,
35, 37, 42, 45, 34, 35, 53, 42, 41, 43, 38, 48, 46, 45, 49, 36, 38, 43, 40, 47, 38, 38, 41, 41, 38, 40, 35, 37, 50, 44, 54, 51)
> mean(x) [1] 42.4127
> var(x) [1] 55.3108
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = -2.4885, 自由度 = 171, P値 = 0.01378
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -4.8970633 -0.5646971
標本推定値:
平均値x 平均値y
39.68182 42.41270
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = -2.4138, 自由度 = 117.735, P値 = 0.01733
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -4.9713019 -0.4904585
標本推定値:
平均値x 平均値y
39.68182 42.41270
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.7995, 第1自由度 = 109, 第2自由度 = 62, P値 = 0.3069
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.5057011 1.2296866
標本推定値:
分散比
0.7994628
試料とBGとの差異がある。
>
３. 補足コメント
.4 - 39.7 = 2.7 CPM
2.7 * 0.47 = 1.3 Bq
容積が、15 * 10 * 2 cm , 面積あたりの重さが 0.28 g/cm なので、底面の半径 7.42 cm、上面の半径 4.14cm 高さ 2cm 円錐台全体(216cm3)が受光範囲で、試料による吸収がないとする。
比重が 40 / (15*10*2) = 0.13 g/cm3 だから、円錐台の質量は 28.8g。よって1kgあたりの線量は
1.3 * 1000 / 29 = 45 Bq/kg
75:名無しに影響はない(栃木県)
12/04/15 16:31:16.62 FM+psz0R
１．測定対象
「富良野産玉葱(輪切り下8mmの乾燥 12.0g)」の分析
２．測定結果
まずは生データ
1. n=31 BG
x <- c(47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45)
> mean(x)
[1] 40.12903
> var(x)
[1] 29.64946
2. n=34 試料
x <- c(54, 36, 60, 46, 40, 61, 63, 41, 50, 40, 52, 52, 45, 55, 56, 45, 50, 50, 60, 47, 56, 50, 51, 49, 51, 51, 62, 37, 43, 45, 38, 45, 49, 49)
> mean(x)
[1] 49.38235
> var(x)
[1] 51.87968
3. n=31 BG
x <- c(41, 53, 43, 47, 41, 42, 49, 50, 37, 40, 36, 32, 34, 45, 44, 38, 37, 40, 38, 46, 39, 33, 46, 36, 37, 40, 39, 44, 47, 48, 37)
> mean(x)
[1] 41.25806
> var(x)
[1] 27.79785
4. n=32 試料
x <- c(40, 54, 40, 55, 45, 49, 40, 53, 36, 57, 40, 39, 60, 50, 43, 38, 52, 59, 59, 43, 63, 38, 34, 47, 47, 37, 38, 35, 54, 55, 40, 46)
> mean(x)
[1] 46.4375
> var(x)
[1] 71.80242
5. n=31 BG
x <- c(40, 34, 42, 39, 41, 35, 42, 43, 50, 40, 50, 38, 47, 38, 36, 38, 40, 42, 34, 47, 44, 46, 39, 37, 29, 41, 35, 34, 45, 46, 40)
> mean(x)
[1] 40.3871
> var(x)
[1] 24.91183
76:名無しに影響はない(栃木県)
12/04/15 16:33:00.19 FM+psz0R
４. データ貼り付け
全体の分析
x <- c(47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45 ,
54, 36, 60, 46, 40, 61, 63, 41, 50, 40, 52, 52, 45, 55, 56, 45, 50, 50, 60, 47, 56, 50, 51, 49, 51, 51, 62, 37, 43, 45, 38, 45, 49, 49 ,
41, 53, 43, 47, 41, 42, 49, 50, 37, 40, 36, 32, 34, 45, 44, 38, 37, 40, 38, 46, 39, 33, 46, 36, 37, 40, 39, 44, 47, 48, 37 ,
40, 54, 40, 55, 45, 49, 40, 53, 36, 57, 40, 39, 60, 50, 43, 38, 52, 59, 59, 43, 63, 38, 34, 47, 47, 37, 38, 35, 54, 55, 40, 46 ,
40, 34, 42, 39, 41, 35, 42, 43, 50, 40, 50, 38, 47, 38, 36, 38, 40, 42, 34, 47, 44, 46, 39, 37, 29, 41, 35, 34, 45, 46, 40)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 ,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4 ,
5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5)
> mean(x)
[1] 43.6478
> var(x)
[1] 54.84985
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 13.5625, 第1自由度 = 4, 第2自由度 = 154, P値 = 1.704e-09
有意。群別の差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 13.3107, 自由度 = 4, P値 = 0.009854
有意。先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 12.4487, 第1自由度 = 4.000, 第2自由度 = 76.685, P値 = 7.49e-08
有意。群別の差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 1 0.6289308 0.6289308
30 13 8.1761006 8.8050314
35 37 23.2704403 32.0754717
40 43 27.0440252 59.1194969
45 30 18.8679245 77.9874214
50 19 11.9496855 89.9371069
55 9 5.6603774 95.5974843
60 7 4.4025157 100.0000000
>
77:名無しに影響はない(栃木県)
12/04/15 16:33:24.13 FM+psz0R
BGと試料の比較
x <- c(47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45 ,
54, 36, 60, 46, 40, 61, 63, 41, 50, 40, 52, 52, 45, 55, 56, 45, 50, 50, 60, 47, 56, 50, 51, 49, 51, 51, 62, 37, 43, 45, 38, 45, 49, 49 ,
41, 53, 43, 47, 41, 42, 49, 50, 37, 40, 36, 32, 34, 45, 44, 38, 37, 40, 38, 46, 39, 33, 46, 36, 37, 40, 39, 44, 47, 48, 37 ,
40, 54, 40, 55, 45, 49, 40, 53, 36, 57, 40, 39, 60, 50, 43, 38, 52, 59, 59, 43, 63, 38, 34, 47, 47, 37, 38, 35, 54, 55, 40, 46 ,
40, 34, 42, 39, 41, 35, 42, 43, 50, 40, 50, 38, 47, 38, 36, 38, 40, 42, 34, 47, 44, 46, 39, 37, 29, 41, 35, 34, 45, 46, 40)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 49.9886, 第1自由度 = 1, 第2自由度 = 157, P値 = 4.786e-11
有意。バックグラウンドと試料との差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 13.6218, 自由度 = 1, P値 = 0.0002236
有意。先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 43.6323, 第1自由度 = 1.000, 第2自由度 = 104.015, P値 = 1.721e-09
有意。バックグラウンドと試料との差異がある。
>
78:名無しに影響はない(栃木県)
12/04/15 16:33:47.68 FM+psz0R
繰り返しによる影響
1. n = 93 BG 繰り返し数 =3
x <- c( 47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45 ,
41, 53, 43, 47, 41, 42, 49, 50, 37, 40, 36, 32, 34, 45, 44, 38, 37, 40, 38, 46, 39, 33, 46, 36, 37, 40, 39, 44, 47, 48, 37 ,
40, 34, 42, 39, 41, 35, 42, 43, 50, 40, 50, 38, 47, 38, 36, 38, 40, 42, 34, 47, 44, 46, 39, 37, 29, 41, 35, 34, 45, 46, 40)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2 ,
3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> mean(x)
[1] 40.5914
> var(x)
[1] 27.0921
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 0.3952, 第1自由度 = 2, 第2自由度 = 90, P値 = 0.6747
有意ではない。繰り返しによる差異は不明。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.2272, 自由度 = 2, P値 = 0.8926
有意ではない。先の分散分析は有効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 0.3813, 第1自由度 = 2.000, 第2自由度 = 59.921, P値 = 0.6846
有意ではない。繰り返しによる差異は不明。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 1 1.075269 1.075269
30 12 12.903226 13.978495
35 27 29.032258 43.010753
40 32 34.408602 77.419355
45 16 17.204301 94.623656
50 4 4.301075 98.924731
55 1 1.075269 100.000000
>
79:名無しに影響はない(栃木県)
12/04/15 16:34:12.23 FM+psz0R
2. n = 66 試料 繰り返し数 =2
x <- c( 54, 36, 60, 46, 40, 61, 63, 41, 50, 40, 52, 52, 45, 55, 56, 45, 50, 50, 60, 47, 56, 50, 51, 49, 51, 51, 62, 37, 43, 45, 38, 45, 49, 49 ,
40, 54, 40, 55, 45, 49, 40, 53, 36, 57, 40, 39, 60, 50, 43, 38, 52, 59, 59, 43, 63, 38, 34, 47, 47, 37, 38, 35, 54, 55, 40, 46)
g <- c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 ,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2)
> mean(x)
[1] 47.95455
> var(x)
[1] 62.78252
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.3234, 第1自由度 = 1, 第2自由度 = 64, P値 = 0.1324
有意。繰り返しによる差異がある。
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.8303, 自由度 = 1, P値 = 0.3622
有意。先の分散分析は無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.3005, 第1自由度 = 1.000, 第2自由度 = 60.999, P値 = 0.1345
有意。繰り返しによる差異がある。
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
30 1 1.515152 1.515152
35 10 15.151515 16.666667
40 11 16.666667 33.333333
45 14 21.212121 54.545455
50 15 22.727273 77.272727
55 8 12.121212 89.393939
60 7 10.606061 100.000000
>
80:名無しに影響はない(栃木県)
12/04/15 16:35:15.38 FM+psz0R
BGと試料の比較
1. n = 93 BG 繰り返し数 =3
gr1 <- c( 47, 33, 47, 37, 35, 42, 44, 42, 42, 35, 55, 43, 44, 34, 33, 45, 37, 36, 37, 44, 43, 41, 32, 42, 32, 39, 33, 44, 43, 38, 45 ,
41, 53, 43, 47, 41, 42, 49, 50, 37, 40, 36, 32, 34, 45, 44, 38, 37, 40, 38, 46, 39, 33, 46, 36, 37, 40, 39, 44, 47, 48, 37 ,
40, 34, 42, 39, 41, 35, 42, 43, 50, 40, 50, 38, 47, 38, 36, 38, 40, 42, 34, 47, 44, 46, 39, 37, 29, 41, 35, 34, 45, 46, 40)
> mean(x)
[1] 40.5914
> var(x)
[1] 27.0921
2. n = 66 試料 繰り返し数 =2
gr2 <- c( 54, 36, 60, 46, 40, 61, 63, 41, 50, 40, 52, 52, 45, 55, 56, 45, 50, 50, 60, 47, 56, 50, 51, 49, 51, 51, 62, 37, 43, 45, 38, 45, 49, 49 ,
40, 54, 40, 55, 45, 49, 40, 53, 36, 57, 40, 39, 60, 50, 43, 38, 52, 59, 59, 43, 63, 38, 34, 47, 47, 37, 38, 35, 54, 55, 40, 46)
> mean(x)
[1] 47.95455
> var(x)
[1] 62.78252
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）
データ: gr1 と gr2
t値 = -7.0703, 自由度 = 157, P値 = 4.786e-11
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -9.420160 -5.306135
標本推定値:
平均値x 平均値y
40.59140 47.95455
有意。試料とバックグラウンドに差異がある
> t.test(gr1, gr2)
二標本t検定（Welchの方法）
データ: gr1 と gr2
t値 = -6.6055, 自由度 = 104.015, P値 = 1.721e-09
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -9.573642 -5.152653
標本推定値:
平均値x 平均値y
40.59140 47.95455
有意。試料とバックグラウンドに差異がある
> var.test(gr1, gr2)
二群の等分散性の検定
データ: gr1 と gr2
F = 0.4315, 第1自由度 = 92, 第2自由度 = 65, P値 = 0.000212
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.2720116 0.6724900
標本推定値:
分散比
0.431523
有意。分散に違いがあるので、（Welchの方法）を使用する。
>