【ガイガー】インスペクター＋統計スレ２

４.　データ貼り付け
全体の分析
x <- c( 49,34,38,36,35,53,26,43,42,30,39,35,31,33,38,32,35,43,52,40,41,36,42,43,41,42,33,34,35,41,31,27,36,35,42 ,
40,35,46,43,41,50,37,55,47,39,43,43,54,50,40,59,45,32,44,43,51,43,45,48,38,38,36,41,40,52,46,40,43,40,50,38,37,43,42,43,43,38,44,31,49,41 ,
35,45,35,37,42,39,32,45,45,39,31,52,34,42,39,35,32,35,32,39,34,28,34,34,35,36,34,44,37,35,42,42,35,32,26,38,30,
39,36,33,41,41,41,35,48,48,43,52,44,58,33,37,36,52,32,32,25,30,50,37,30,40,35,36,43,43,39,53,38,39,29,40,35,41 ,
39,53,44,36,50,42,59,49,30,44,35,49,45,61,37,36,53,35,45,56,56,44,45,53,45,48,51,35,42,34,47,47,37 ,
33,28,28,36,42,43,40,28,38,35,42,45,36,40,34,30,42,43,39,37,34,44,42,32,47,46,55,42,34,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,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,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,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)
> mean(x)
[1] 40.25688
> var(x)
[1] 50.21021
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 9.9357, 第1自由度 = 4, 第2自由度 = 213, P値 = 2.169e-07
有意。群間に差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 3.7347, 自由度 = 4, P値 = 0.4431
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 9.1963, 第1自由度 = 4.000, 第2自由度 = 89.627, P値 = 2.848e-06
有意。群間に差異がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 9 4.1284404 4.12844
30 35 16.0550459 20.18349
35 61 27.9816514 48.16514
40 61 27.9816514 76.14679
45 26 11.9266055 88.07339
50 18 8.2568807 96.33028
55 7 3.2110092 99.54128
60 1 0.4587156 100.00000
>

BGと試料の比較
x <- c( 49,34,38,36,35,53,26,43,42,30,39,35,31,33,38,32,35,43,52,40,41,36,42,43,41,42,33,34,35,41,31,27,36,35,42 ,
40,35,46,43,41,50,37,55,47,39,43,43,54,50,40,59,45,32,44,43,51,43,45,48,38,38,36,41,40,52,46,40,43,40,50,38,37,43,42,43,43,38,44,31,49,41 ,
35,45,35,37,42,39,32,45,45,39,31,52,34,42,39,35,32,35,32,39,34,28,34,34,35,36,34,44,37,35,42,42,35,32,26,38,30,39,36,33,
41,41,41,35,48,48,43,52,44,58,33,37,36,52,32,32,25,30,50,37,30,40,35,36,43,43,39,53,38,39,29,40,35,41 ,
39,53,44,36,50,42,59,49,30,44,35,49,45,61,37,36,53,35,45,56,56,44,45,53,45,48,51,35,42,34,47,47,37 ,
33,28,28,36,42,43,40,28,38,35,42,45,36,40,34,30,42,43,39,37,34,44,42,32,47,46,55,42,34,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,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,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,1,1,1,1,1,1,1,1,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)
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 38.3633, 第1自由度 = 1, 第2自由度 = 216, P値 = 2.926e-09
有意。試料とバックグラウンドに差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 0.323, 自由度 = 1, P値 = 0.5698
有意ではない。先の分散分析は有効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 37.1783, 第1自由度 = 1.000, 第2自由度 = 154.711, P値 = 8.26e-09
有意。試料とバックグラウンドに差異がある。

>

繰り返しによる影響

1. n = 139 BG 繰り返し数 =3
x <- c( 49,34,38,36,35,53,26,43,42,30,39,35,31,33,38,32,35,43,52,40,41,36,42,43,41,42,33,34,35,41,31,27,36,35,42 ,
35,45,35,37,42,39,32,45,45,39,31,52,34,42,39,35,32,35,32,39,34,28,34,34,35,36,34,44,37,35,42,42,35,32,26,38,30,39,
36,33,41,41,41,35,48,48,43,52,44,58,33,37,36,52,32,32,25,30,50,37,30,40,35,36,43,43,39,53,38,39,29,40,35,41 ,
33,28,28,36,42,43,40,28,38,35,42,45,36,40,34,30,42,43,39,37,34,44,42,32,47,46,55,42,34,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,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,2,2,2,2,2,2,2,2,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] 38.18705
> var(x)
[1] 41.05172
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 0.1462, 第1自由度 = 2, 第2自由度 = 136, P値 = 0.8641
有意ではない。繰り返しによる差異は不明。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 0.1473, 自由度 = 2, P値 = 0.929
有意ではない。先の分散分析は有効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 0.1509, 第1自由度 = 2.000, 第2自由度 = 66.709, P値 = 0.8602
有意ではない。繰り返しによる差異は不明。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 9 6.474820 6.47482
30 31 22.302158 28.77698
35 44 31.654676 60.43165
40 36 25.899281 86.33094
45 10 7.194245 93.52518
50 7 5.035971 98.56115
55 2 1.438849 100.00000
>

2. n = 79 栃木県産ジャガイモ 30.5wg 9.8dg 2012.07.06 ６日間乾燥 繰り返し数 =2
x <- c( 40,35,46,43,41,50,37,55,47,39,43,43,54,50,40,59,45,32,44,43,51,43,45,48,38,38,36,41,40,52,46,40,43,40,50,38,37,43,42,43,43,38,44,31,49,41 ,
39,53,44,36,50,42,59,49,30,44,35,49,45,61,37,36,53,35,45,56,56,44,45,53,45,48,51,35,42,34,47,47,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,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] 43.89873
> var(x)
[1] 45.98961
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 1.2622, 第1自由度 = 1, 第2自由度 = 77, P値 = 0.2647
有意。繰り返しによる差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 3.3429, 自由度 = 1, P値 = 0.0675
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 1.1459, 第1自由度 = 1.000, 第2自由度 = 56.064, P値 = 0.289
有意。繰り返しによる差異がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
30 4 5.063291 5.063291
35 17 21.518987 26.582278
40 25 31.645570 58.227848
45 16 20.253165 78.481013
50 11 13.924051 92.405063
55 5 6.329114 98.734177
60 1 1.265823 100.000000
>

BGと試料の比較

1. n = 139 BG 繰り返し数 =3
gr1 <- c( 49,34,38,36,35,53,26,43,42,30,39,35,31,33,38,32,35,43,52,40,41,36,42,43,41,42,33,34,35,41,31,27,36,35,42 ,
35,45,35,37,42,39,32,45,45,39,31,52,34,42,39,35,32,35,32,39,34,28,34,34,35,36,34,44,37,35,42,42,35,32,26,38,30,39,36,
33,41,41,41,35,48,48,43,52,44,58,33,37,36,52,32,32,25,30,50,37,30,40,35,36,43,43,39,53,38,39,29,40,35,41 ,
33,28,28,36,42,43,40,28,38,35,42,45,36,40,34,30,42,43,39,37,34,44,42,32,47,46,55,42,34,45)

2. n = 79 栃木県産ジャガイモ 30.5wg 9.8dg 2012.07.06 ６日間乾燥 繰り返し数 =2
gr2 <- c( 40,35,46,43,41,50,37,55,47,39,43,43,54,50,40,59,45,32,44,43,51,43,45,48,38,38,36,41,40,52,46,40,43,40,50,38,37,43,42,43,43,38,44,31,49,41 ,
39,53,44,36,50,42,59,49,30,44,35,49,45,61,37,36,53,35,45,56,56,44,45,53,45,48,51,35,42,34,47,47,37)
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）

データ: gr1 と gr2
t値 = -6.1938, 自由度 = 216, P値 = 2.926e-09
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -7.529266 -3.894101
標本推定値:
平均値x 平均値y
38.18705 43.89873
有意。試料とバックグラウンドに差異がある。

> t.test(gr1, gr2)
二標本t検定（Welchの方法）

データ: gr1 と gr2
t値 = -6.0974, 自由度 = 154.711, P値 = 8.26e-09
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -7.562136 -3.861231
標本推定値:
平均値x 平均値y
38.18705 43.89873
有意。試料とバックグラウンドに差異がある。

> var.test(gr1, gr2)
二群の等分散性の検定

データ: gr1 と gr2
F = 0.8926, 第1自由度 = 138, 第2自由度 = 78, P値 = 0.5572
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.5947746 1.3106013
標本推定値:
分散比
0.8926301
有意。分散が異なるので（Welchの方法）を使用する。

>

43.89873 - 38.18705 = 5.711681 CPM (2.68Bq)
2.69 * 1000 / 30.5 = 88.0 Bq/wkg
2.69 * 1000 / 9.8 = 273 Bq/dkg

１．測定対象
「栃木県産ホウノキの葉 3.5wg 2.0dg 2012.07.20取得」の分析

２．測定結果
1. n=30 糊台
x <- c( 34,39,42,34,27,43,26,36,45,33,44,48,31,36,33,41,34,34,28,36,45,41,32,37,37,40,34,39,41,44)

2. n=30 栃木県産ホウノキの葉 3.5wg 2.0dg 2012.07.20取得
x <- c( 31,29,50,56,33,53,36,51,45,40,47,43,38,34,37,46,37,37,37,35,43,41,47,47,43,34,30,31,36,51)

3. n=30 糊台
x <- c( 34,47,35,39,34,38,41,24,39,52,47,33,38,34,32,29,40,39,39,24,38,35,51,31,40,24,45,31,46,30)

4. n=31 栃木県産ホウノキの葉 3.5wg 2.0dg 2012.07.20取得
x <- c( 40,34,32,47,43,37,53,48,55,43,40,37,52,43,49,47,42,44,54,51,37,41,39,39,39,41,42,43,47,34,42)

5. n=30 糊台
x <- c( 39,29,47,48,33,44,43,34,35,27,27,31,35,36,45,36,39,43,44,28,43,29,39,38,39,44,33,37,42,36)

４.　データ貼り付け
全体の分析
x <- c( 34,39,42,34,27,43,26,36,45,33,44,48,31,36,33,41,34,34,28,36,45,41,32,37,37,40,34,39,41,44 ,
31,29,50,56,33,53,36,51,45,40,47,43,38,34,37,46,37,37,37,35,43,41,47,47,43,34,30,31,36,51 ,
34,47,35,39,34,38,41,24,39,52,47,33,38,34,32,29,40,39,39,24,38,35,51,31,40,24,45,31,46,30 ,
40,34,32,47,43,37,53,48,55,43,40,37,52,43,49,47,42,44,54,51,37,41,39,39,39,41,42,43,47,34,42 ,
39,29,47,48,33,44,43,34,35,27,27,31,35,36,45,36,39,43,44,28,43,29,39,38,39,44,33,37,42,36)
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 ,
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)
> mean(x)
[1] 39.06623
> var(x)
[1] 47.64892
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 5.1946, 第1自由度 = 4, 第2自由度 = 146, P値 = 0.0006104
有意。群間に差がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 3.7831, 自由度 = 4, P値 = 0.4362
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 5.6152, 第1自由度 = 4.000, 第2自由度 = 72.723, P値 = 0.0005383
有意。群間に差がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 3 1.986755 1.986755
25 10 6.622517 8.609272
30 30 19.867550 28.476821
35 41 27.152318 55.629139
40 35 23.178808 78.807947
45 20 13.245033 92.052980
50 10 6.622517 98.675497
55 2 1.324503 100.000000
>

BGと試料の比較
x <- c( 34,39,42,34,27,43,26,36,45,33,44,48,31,36,33,41,34,34,28,36,45,41,32,37,37,40,34,39,41,44 ,
31,29,50,56,33,53,36,51,45,40,47,43,38,34,37,46,37,37,37,35,43,41,47,47,43,34,30,31,36,51 ,
34,47,35,39,34,38,41,24,39,52,47,33,38,34,32,29,40,39,39,24,38,35,51,31,40,24,45,31,46,30 ,
40,34,32,47,43,37,53,48,55,43,40,37,52,43,49,47,42,44,54,51,37,41,39,39,39,41,42,43,47,34,42 ,
39,29,47,48,33,44,43,34,35,27,27,31,35,36,45,36,39,43,44,28,43,29,39,38,39,44,33,37,42,36)
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 ,
1,1,1,1,1,1,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)
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 18.6345, 第1自由度 = 1, 第2自由度 = 149, P値 = 2.873e-05
有意。試料とバックグラウンドに差がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 0.3483, 自由度 = 1, P値 = 0.5551
有意ではない。先の分散分析は有効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 18.1399, 第1自由度 = 1.000, 第2自由度 = 122.827, P値 = 4.043e-05
有意。試料とバックグラウンドに差がある。

>

繰り返しによる影響

1. n = 90 BG 繰り返し数 =3
x <- c( 34,39,42,34,27,43,26,36,45,33,44,48,31,36,33,41,34,34,28,36,45,41,32,37,37,40,34,39,41,44 ,
34,47,35,39,34,38,41,24,39,52,47,33,38,34,32,29,40,39,39,24,38,35,51,31,40,24,45,31,46,30 ,
39,29,47,48,33,44,43,34,35,27,27,31,35,36,45,36,39,43,44,28,43,29,39,38,39,44,33,37,42,36)
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 ,
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] 37.17778
> var(x)
[1] 40.21523
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 0.0408, 第1自由度 = 2, 第2自由度 = 87, P値 = 0.96
有意ではない。繰り返しによる差異は不明。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 2.3543, 自由度 = 2, P値 = 0.3082
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 0.0386, 第1自由度 = 2.000, 第2自由度 = 57.328, P値 = 0.9621
有意ではない。繰り返しによる差異は不明。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 3 3.333333 3.333333
25 9 10.000000 13.333333
30 21 23.333333 36.666667
35 27 30.000000 66.666667
40 18 20.000000 86.666667
45 10 11.111111 97.777778
50 2 2.222222 100.000000
>

2. n = 61 栃木県産ホウノキの葉 3.5wg 2.0dg 2012.07.20取得 繰り返し数 =2
x <- c( 31,29,50,56,33,53,36,51,45,40,47,43,38,34,37,46,37,37,37,35,43,41,47,47,43,34,30,31,36,51 ,
40,34,32,47,43,37,53,48,55,43,40,37,52,43,49,47,42,44,54,51,37,41,39,39,39,41,42,43,47,34,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 ,2,2,2,2,2,2,2,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.85246
> var(x)
[1] 46.22787
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 2.0378, 第1自由度 = 1, 第2自由度 = 59, P値 = 0.1587
有意。繰り返しによる差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 1.2452, 自由度 = 1, P値 = 0.2645
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 2.0239, 第1自由度 = 1.000, 第2自由度 = 55.852, P値 = 0.1604
有意。繰り返しによる差異がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 1 1.639344 1.639344
30 9 14.754098 16.393443
35 14 22.950820 39.344262
40 17 27.868852 67.213115
45 10 16.393443 83.606557
50 8 13.114754 96.721311
55 2 3.278689 100.000000
>

BGと試料の比較

1. n = 90 BG 繰り返し数 =3
gr1 <- c( 34,39,42,34,27,43,26,36,45,33,44,48,31,36,33,41,34,34,28,36,45,41,32,37,37,40,34,39,41,44 ,
34,47,35,39,34,38,41,24,39,52,47,33,38,34,32,29,40,39,39,24,38,35,51,31,40,24,45,31,46,30 ,
39,29,47,48,33,44,43,34,35,27,27,31,35,36,45,36,39,43,44,28,43,29,39,38,39,44,33,37,42,36)

2. n = 61 栃木県産ホウノキの葉 3.5wg 2.0dg 2012.07.20取得 繰り返し数 =2
gr2 <- c( 31,29,50,56,33,53,36,51,45,40,47,43,38,34,37,46,37,37,37,35,43,41,47,47,43,34,30,31,36,51 ,
40,34,32,47,43,37,53,48,55,43,40,37,52,43,49,47,42,44,54,51,37,41,39,39,39,41,42,43,47,34,42)
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）

データ: gr1 と gr2
t値 = -4.3168, 自由度 = 149, P値 = 2.873e-05
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -6.814527 -2.534835
標本推定値:
平均値x 平均値y
37.17778 41.85246
有意。試料とバックグラウンドに差がある。

> t.test(gr1, gr2)
二標本t検定（Welchの方法）

データ: gr1 と gr2
t値 = -4.2591, 自由度 = 122.827, P値 = 4.043e-05
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -6.847293 -2.502069
標本推定値:
平均値x 平均値y
37.17778 41.85246
有意。試料とバックグラウンドに差がある。

> var.test(gr1, gr2)
二群の等分散性の検定

データ: gr1 と gr2
F = 0.8699, 第1自由度 = 89, 第2自由度 = 60, P値 = 0.5448
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.5396328 1.3726646
標本推定値:
分散比
0.8699348
有意。分散が異なるので（Welchの方法）を使用する。

>

３.補足コメント
41.85246 - 37.17778 = 4.674679 CPM (2.20Bq)
2. n = 61 栃木県産ホウノキの葉 3.5wg 2.0dg 2012.07.20取得 繰り返し数 =2
2.20 * 1000 / 3.5 = 627 Bq/wkg
2.20 * 1000 / 2.0 = 1099 Bq/dkg

ホウノキは、夏場防腐剤として使われることがあります。

１．測定対象
「栃木県産ハリキリの葉 17.5wg 5.6dg 2012.07.12採取」の分析
柏餅用の葉として使用する。防腐剤が含まれているので、夏の保存食に使われる。

２．測定結果
1. n=30 糊台
x <- c( 40, 38, 27, 37, 36, 32, 40, 35, 29, 44, 35, 36, 40, 39, 37, 39, 33, 31, 35, 39, 31, 36, 38, 48, 32, 48, 35, 42, 45, 28)

2. n=30 栃木県産ハリキリ 17.5wg 5.6dg 2012.07.12採取
x <- c( 49, 50, 55, 47, 61, 43, 52, 45, 53, 55, 49, 43, 49, 54, 66, 48, 46, 56, 52, 49, 54, 59, 47, 65, 52, 56, 57, 55, 57, 46)

3. n=30 糊台
x <- c( 38, 35, 36, 46, 41, 29, 33, 35, 34, 29, 31, 47, 39, 46, 37, 51, 38, 34, 44, 45, 54, 45, 37, 51, 38, 25, 26, 34, 29, 45)

4. n=30 栃木県産ハリキリ 17.5wg 5.6dg 2012.07.12採取
x <- c( 47, 60, 55, 46, 55, 64, 48, 54, 64, 77, 48, 47, 44, 59, 52, 63, 55, 61, 64, 54, 38, 45, 73, 45, 58, 56, 60, 56, 50, 54)

5. n=30 糊台
x <- c( 36, 29, 27, 34, 27, 35, 30, 34, 28, 36, 25, 35, 23, 31, 30, 28, 37, 31, 35, 29, 24, 32, 30, 41, 34, 34, 29, 30, 22, 27)

全体の分析
x <- c( 40,38,27,37,36,32,40,35,29,44,35,36,40,39,37,39,33,31,35,39,31,36,38,48,32,48,35,42,45,28 ,
49,50,55,47,61,43,52,45,53,55,49,43,49,54,66,48,46,56,52,49,54,59,47,65,52,56,57,55,57,46 ,
38,35,36,46,41,29,33,35,34,29,31,47,39,46,37,51,38,34,44,45,54,45,37,51,38,25,26,34,29,45 ,
47,60,55,46,55,64,48,54,64,77,48,47,44,59,52,63,55,61,64,54,38,45,73,45,58,56,60,56,50,54 ,
36,29,27,34,27,35,30,34,28,36,25,35,23,31,30,28,37,31,35,29,24,32,30,41,34,34,29,30,22,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 ,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)
> mean(x)
[1] 42.68
> var(x)
[1] 130.9305
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 76.4086, 第1自由度 = 4, 第2自由度 = 145, P値 < 2.2e-16
有意。群間に差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 16.3109, 自由度 = 4, P値 = 0.002629
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 88.0846, 第1自由度 = 4.000, 第2自由度 = 71.577, P値 < 2.2e-16
有意。群間に差異がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 3 2.0000000 2.00000
25 17 11.3333333 13.33333
30 21 14.0000000 27.33333
35 30 20.0000000 47.33333
40 11 7.3333333 54.66667
45 26 17.3333333 72.00000
50 15 10.0000000 82.00000
55 15 10.0000000 92.00000
60 8 5.3333333 97.33333
65 2 1.3333333 98.66667
70 1 0.6666667 99.33333
75 1 0.6666667 100.00000
>

BGと試料の比較
x <- c( 40,38,27,37,36,32,40,35,29,44,35,36,40,39,37,39,33,31,35,39,31,36,38,48,32,48,35,42,45,28 ,
49,50,55,47,61,43,52,45,53,55,49,43,49,54,66,48,46,56,52,49,54,59,47,65,52,56,57,55,57,46 ,
38,35,36,46,41,29,33,35,34,29,31,47,39,46,37,51,38,34,44,45,54,45,37,51,38,25,26,34,29,45 ,
47,60,55,46,55,64,48,54,64,77,48,47,44,59,52,63,55,61,64,54,38,45,73,45,58,56,60,56,50,54 ,
36,29,27,34,27,35,30,34,28,36,25,35,23,31,30,28,37,31,35,29,24,32,30,41,34,34,29,30,22,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 ,
1,1,1,1,1,1,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 ,
1,1,1,1,1,1,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 = 76.4086, 第1自由度 = 4, 第2自由度 = 145, P値 < 2.2e-16
有意。試料とバックグラウンドに差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 16.3109, 自由度 = 4, P値 = 0.002629
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 88.0846, 第1自由度 = 4.000, 第2自由度 = 71.577, P値 < 2.2e-16
有意。試料とバックグラウンドに差異がある。

>

繰り返しによる影響

1. n = 90 BG 繰り返し数 =3
x <- c( 40,38,27,37,36,32,40,35,29,44,35,36,40,39,37,39,33,31,35,39,31,36,38,48,32,48,35,42,45,28 ,
38,35,36,46,41,29,33,35,34,29,31,47,39,46,37,51,38,34,44,45,54,45,37,51,38,25,26,34,29,45 ,
36,29,27,34,27,35,30,34,28,36,25,35,23,31,30,28,37,31,35,29,24,32,30,41,34,34,29,30,22,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 ,
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] 35.33333
> var(x)
[1] 45.75281
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 13.6996, 第1自由度 = 2, 第2自由度 = 87, P値 = 6.725e-06
有意。繰り返しによる差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 8.7621, 自由度 = 2, P値 = 0.01251
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 16.8141, 第1自由度 = 2.000, 第2自由度 = 55.844, P値 = 1.923e-06
有意。繰り返しによる差異がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 3 3.333333 3.333333
25 17 18.888889 22.222222
30 21 23.333333 45.555556
35 29 32.222222 77.777778
40 8 8.888889 86.666667
45 9 10.000000 96.666667
50 3 3.333333 100.000000
>

2. n = 60 栃木県産ハリキリ 17.5wg 5.6dg 2012.07.12採取 繰り返し数 =2
x <- c( 49,50,55,47,61,43,52,45,53,55,49,43,49,54,66,48,46,56,52,49,54,59,47,65,52,56,57,55,57,46 ,
47,60,55,46,55,64,48,54,64,77,48,47,44,59,52,63,55,61,64,54,38,45,73,45,58,56,60,56,50,54)
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)
> mean(x)
[1] 53.7
> var(x)
[1] 55.80678
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 2.0436, 第1自由度 = 1, 第2自由度 = 58, P値 = 0.1582
有意。繰り返しによる差異がある。

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 4.2465, 自由度 = 1, P値 = 0.03933
有意。分散が異なるので先の分散分析は無効。

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 2.0436, 第1自由度 = 1.000, 第2自由度 = 50.949, P値 = 0.1589
有意。繰り返しによる差異がある。

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
35 1 1.666667 1.666667
40 3 5.000000 6.666667
45 17 28.333333 35.000000
50 12 20.000000 55.000000
55 15 25.000000 80.000000
60 8 13.333333 93.333333
65 2 3.333333 96.666667
70 1 1.666667 98.333333
75 1 1.666667 100.000000
>

BGと試料の比較

1. n = 90 BG 繰り返し数 =3
gr1 <- c( 40,38,27,37,36,32,40,35,29,44,35,36,40,39,37,39,33,31,35,39,31,36,38,48,32,48,35,42,45,28 ,
38,35,36,46,41,29,33,35,34,29,31,47,39,46,37,51,38,34,44,45,54,45,37,51,38,25,26,34,29,45 ,
36,29,27,34,27,35,30,34,28,36,25,35,23,31,30,28,37,31,35,29,24,32,30,41,34,34,29,30,22,27)

2. n = 60 栃木県産ハリキリ 17.5wg 5.6dg 2012.07.12採取 繰り返し数 =2
gr2 <- c( 49,50,55,47,61,43,52,45,53,55,49,43,49,54,66,48,46,56,52,49,54,59,47,65,52,56,57,55,57,46 ,
47,60,55,46,55,64,48,54,64,77,48,47,44,59,52,63,55,61,64,54,38,45,73,45,58,56,60,56,50,54)
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）

データ: gr1 と gr2
t値 = -15.622, 自由度 = 148, P値 < 2.2e-16
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -20.68997 -16.04336
標本推定値:
平均値x 平均値y
35.33333 53.70000
有意。試料とバックグラウンドに差異がある。

> t.test(gr1, gr2)
二標本t検定（Welchの方法）

データ: gr1 と gr2
t値 = -15.3137, 自由度 = 117.792, P値 < 2.2e-16
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -20.74178 -15.99155
標本推定値:
平均値x 平均値y
35.33333 53.70000
有意。試料とバックグラウンドに差異がある。

> var.test(gr1, gr2)
二群の等分散性の検定

データ: gr1 と gr2
F = 0.8198, 第1自由度 = 89, 第2自由度 = 59, P値 = 0.3927
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.5070552 1.2960152
標本推定値:
分散比
0.8198432
有意。分散が異なるので（Welchの方法）を使用する。

>

３.補足コメント
栃木県産ハリキリ 17.5wg 5.6dg 2012.07.12採取
柏餅用の葉として使用する。防腐剤が含まれているので、夏の保存食に使われる。
53.70000 - 35.33333 = 18.36667CPM (8.63Bq)
.63 * 1000 / 17.5 = 493 Bq/kg
8.63 * 1000 / 5.6 = 1540 Bq/kg

近所の農家から、タバコの廃棄葉を入手しようとしたら、
契約によりできません
と断られてしまった。
タバコは、重金属を集める性質がある(Pb, Cdの灰中濃度が高い)ので
ウラン系列の濃度を知る指標に使用かと思ったが、うまくいかなかった。
Pb-Bi等ウラン系列の核種を蓄積する性質のある植物を知っていたら知らせて欲しい。

現在は、農産加工品、つまり、塩漬けによる線量低下を測定しようとしているが
はっきり言って下がってくれない。
乾重量あたり1000-6000Bq/kgが300-1000位(1/3-1/10)に落ちるだけで、
湿重量40Bq/kg (地震前の輸入汚染食品を国産品で薄めた後の国内平均濃度)までは下がらない。
線量測定まで手が回らない状態。
結構高線量なので、手が荒れやすい。食品を扱うときには手袋が必須。
素手で船頭やって、乳酸を絞った後、手を洗って拭いて、インスに手を触れたらば
カウントが倍に上がった(30-40→60-70)。
漬物をいじっている間に測定ができれば良いが、
漬物をいじっているときには漬物による汚染対策でインスを使ってしまい測定できない状態。

漏れ一人が測定するのではなく、他の人による調査をキボンヌ。

タバコの「エコー」「わかば」「しんせい」「ゴールデンバット」は200円です。
海外製タバコ葉を前に測ったところ、汚染の疑いがありました。タバコにはもともとポロニウムが含まれています。

ニンジン　ゆで(食塩2%水) 2012.08.13収穫 116Bq/wkg 971Bq/dkg → 26.6Bq/wkg 256Bq/dkg

キャベツ 塩漬け 2012.07.20収穫 193Bq/wkg 1980Bq/dkg → (不明) 868Bq/dkg
　塩漬け 2012.08.03　→ (不明) 593Bq/dkg
　塩漬け 2012.08.13　→ (不明) 51Bq/dkg
これは、キャベツを千切りにして、流水で１時間さらし、塩漬け(10-20%食塩)した。
その後、約１週間放置して、圧搾し漬け汁を絞り廃棄。流水で２４時間さらして、塩漬け(10-20%食塩)、これを繰り返し、水でさらした試料を乾燥し測定したもの。
単に水でさらしてもあまり減らない。圧搾して、高濃度の漬け汁をどれだけ除去するか、が、線量低下に関係する。
単に、水でさらして塩抜きしても線量が下がらない。下がり方のばらつきは、圧搾時に残った水量が原因と考えている。

BGが暴れて、測定にならない日々が続いています。
８月の0.2-0.35uSv/h、９月の0.15-0.20uSv/h程度の新規供給に加えて
粉塵の蓄積があちらこちらに見られます。
気がついた点を少々。
測定場所、栃木県南部のどこか。一説によると、栃木・小山・真岡地区。
測定日時、2012.10.06-08
測定器、インスペクター＋。30秒読み捨てて、その後の３０秒間の最大値と最小値。

特記無しは、おおむね、測定点上空30-50cm。
洗濯機横。部屋の隅。0.17-0.19uSv/h。
土間入り口。0.15-0.19
土間のトタン屋根。0.13-0.15
同雨どい。 0.15-0.19(上空15cm前後)
土間入り口。0.08-0.10
土間上がり口(床直置き、上方向)。0.18-0.20
土間上がり口(横方向)。0.08-0.10
土間室内(水平方向)。0.08-0.10
屋根裏排気口。0.13-0.15

どうも、天井裏に溜まった粉塵が、夜間の気温低下で室内に流入し、
測定時の途中から線量が上昇し、BGが安定しない原因となっているようです。
泥埃が原因かと土間の測定を行ってみたのですが、全体的には低く、掃除が困難な場所が高線量になっている模様です。

BGが安定していないので、ちょっと信頼性が低いのですが
関西産食塩の線量が上昇している可能性があります。
８月以前は、BGに比べて有意差がある状態で、食塩の線量が低くなりました。
しかし、９月に購入した食塩は、瀬戸内産・長崎産を問わず、BGとの有意差が取れませんでした。
７月に関西地区で降下線量の上昇がみられました。降下した粉塵が海に落下し、関西地区の海洋汚染が進んでいる可能性があります。

つまり、瀬戸内産食塩を原料としているやさしおを標準物質として校正している場合に、
やさしおに含まれている瀬戸内産食塩の線量が上昇していると、負の誤差になります。
校正用線源の取り扱いに注意してください。

>>162
土間入り口(屋外)。0.15-0.19
土間入り口(屋内)。0.08-0.10

雨に濡れた衣類の測定結果。
1. n=30 糊台
x <- c( 38,29,38,50,58,42,48,44,29,44,32,45,38,44,21,36,40,44,38,42,32,34,28,52,26,34,41,33,40,43)

2. n=30 2012.10.23 雨の中付近を徘徊した上着(未乾燥)
x <- c( 44,52,37,39,39,51,47,45,36,34,33,48,45,31,34,38,41,45,50,36,45,38,39,37,34,31,30,36,44,39)

3. n=30 糊台
x <- c( 34,48,43,41,40,51,40,38,24,44,45,50,55,31,46,35,46,30,34,32,41,39,43,49,45,26,36,34,38,34)

4. n=30 2012.10.23 雨の中付近を徘徊した上着(未乾燥)
x <- c( 49,41,44,39,44,61,44,54,45,42,51,46,43,43,34,35,33,42,36,40,40,47,34,55,44,33,33,49,56,38)

5. n=30 糊台
x <- c( 39,33,34,51,49,27,35,34,44,27,35,41,54,39,37,43,53,30,36,29,38,46,40,47,37,37,36,44,41,37)

6. n=30 2012.10.23 雨の中付近を徘徊した下着(未乾燥)
x <- c( 35,36,44,43,37,43,42,51,41,44,37,36,48,40,40,41,43,37,46,33,38,42,32,40,49,46,41,31,44,40)

7. n=30 糊台
x <- c( 47,42,42,39,40,45,29,48,31,42,42,37,43,31,33,38,36,34,44,29,37,39,40,31,27,40,39,44,36,37)

8. n=45 2012.10.23 雨の中付近を徘徊した下着(未乾燥)
x <- c( 42,49,40,41,48,47,45,44,44,48,40,51,43,57,52,39,43,41,46,41,43,30,46,47,34,50,52,39,56,40,63,47,42,58,55,35,49,38,50,55,51,37,34,49,33)

9. n=30 糊台
x <- c( 33,42,46,33,37,36,43,35,36,49,38,32,40,44,38,39,37,52,32,43,36,28,37,34,56,32,38,30,31,30)

BG
> mean(x) [1] 38.71333
> var(x) [1] 48.80989
雨の中付近を徘徊した上着(未乾燥)
> mean(x) [1] 41.55
> var(x) [1] 47.98051
雨の中付近を徘徊した下着(未乾燥)
> mean(x) [1] 43.38667
> var(x) [1] 45.53766

単位 CPM インスペクターですので、334で除すると uSv/h になります。
栃木県南部、一説によると、佐野・小山・真岡地区。
2012.10.23 日に小雨の中を徘徊して帰った服(ほぼ乾いていた)をポリ袋に入れて保存。測定。
上着は水が抜ける性質のある布地(濡れたまま吊るしておくと、いくら手で絞っても水滴が下に落下する)
下着は水を保つ性質のある布地。どちらも綿100%。中国製。

中曽根が自分の『天地有情』の中にも『回顧録』の中にも書いてます、俺はＣＩＡのテストを受けた、
英語もあった、論文も書いた、パスした、自分から進んでＣＩＡのテストを受けたちゅーことですね。
それで彼はアメリカに派遣されます、で中曽根は自分で自慢げに書いてますけど、
色んな原子力発電所の法律は俺がみんな作ったってね、野党改進党です、自民党、自由党時代で、
その野党の若造が作れるはずがない、そこの背後にＣＩＡがみんな絡んでます。
そうして日本に原子力発電所ができるんです、東電とか関西電力に作れています。
だけどまったくやってないことなのにやれというわけですから、無理があったわけですね。
無理がありました、だから今日、福島が事故があるのはそこなんですよ。

でブラックボックスといって触ってはいけない、中を見せないような状態なのに原発が日本に来ます。
だから今の福島のある原発はＧＥが開発して間もなくて危険極まりないのに、
彼らは無条件に入れて信じきったわけです、だから無条件に入れてそしてその数年後に、
欠陥商品であるということがわかるわけです、設計者が言うわけですよ、あれは間違ってたって。
だけどアメリカは、その「マーク１」ていうんですけど最初のやつが、まだいっぱいあります。
その設計をし直したり、色んな部品をつけ直したりして、ずうっとやってきてるわけですよ。
日本はまったくやらないまま四十年間、欠陥商品をそのまま使い続けてきたわけです。
日本は四十年間経ってるのに、「危ない」とＧＥも原発を作ったジェネラルエレクトリックも、
報告書を東電に差し上げたのに日本は改善一つしない。
http://www.youtube.com/watch?v=TuVjmXdufS4

でみんなが安全だと思わせるために凄い金を使ったわけですよ。
その金が全部電気代になってる、ということをみなさんは知らないといけないんです。
日本の電気料金は世界の電気料金の三倍なんです、三倍ですよ、でまだこれから上がります。
彼らはみんな太ったわけです、だから自民党の政治家も民主党の政治家も、
ほとんど反対しないじゃないですか、今でさえ反対の声上げないじゃないですか。
ここまで福島の人が苦しんでいるのに何をやってるんだと。
http://www.youtube.com/watch?v=3glGABd52fk

「2012.10.05 使用靴下(足が腫れた)」の分析

1. n=60 糊台
x <- c( 40,39,28,38,43,40,31,45,36,31,31,26,36,38,29,32,44,31,37,35,33,37,47,34,47,20,36,42,43,48,39,47,32,38,37,44,47,43,30,29,42,27,39,38,47,39,46,33,37,39,44,51,41,27,36,42,29,43,46,45)

2. n=30 2012.10.05 使用靴下(足が腫れた)
x <- c( 45,38,38,33,43,36,39,43,35,47,33,40,43,43,40,43,42,46,42,35,39,46,38,44,36,56,48,39,38,36)

3. n=30 糊台
x <- c( 24,30,31,31,37,31,29,27,47,40,36,32,29,40,39,33,35,46,35,46,33,37,52,36,40,40,38,35,42,35)

4. n=30 2012.10.05 使用靴下(足が腫れた)
x <- c( 45,39,45,36,42,37,57,47,21,38,38,38,37,34,39,39,64,47,34,44,39,37,44,39,50,43,25,47,39,36)

5. n=30 糊台
x <- c( 34,39,38,33,40,46,42,27,33,35,39,35,43,43,45,34,29,44,42,41,35,33,27,31,37,36,33,37,41,46)

床に蓄積しているようなので、床の処理を優先してやっています。

BGと試料の比較

1. n = 120 BG 繰り返し数 =3
gr1 <- c( 40,39,28,38,43,40,31,45,36,31,31,26,36,38,29,32,44,31,37,35,33,37,47,34,47,20,36,42,43,48,
39,47,32,38,37,44,47,43,30,29,42,27,39,38,47,39,46,33,37,39,44,51,41,27,36,42,29,43,46,45 ,
24,30,31,31,37,31,29,27,47,40,36,32,29,40,39,33,35,46,35,46,33,37,52,36,40,40,38,35,42,35 ,
34,39,38,33,40,46,42,27,33,35,39,35,43,43,45,34,29,44,42,41,35,33,27,31,37,36,33,37,41,46)

2. n = 60 2012.10.05 使用靴下(足が腫れた) 繰り返し数 =2
gr2 <- c( 45,38,38,33,43,36,39,43,35,47,33,40,43,43,40,43,42,46,42,35,39,46,38,44,36,56,48,39,38,36 ,
45,39,45,36,42,37,57,47,21,38,38,38,37,34,39,39,64,47,34,44,39,37,44,39,50,43,25,47,39,36)
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）

データ: gr1 と gr2
t値 = -3.3565, 自由度 = 178, P値 = 0.0009652
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -5.425442 -1.407891
標本推定値:
平均値x 平均値y
37.31667 40.73333

> t.test(gr1, gr2)
二標本t検定（Welchの方法）

データ: gr1 と gr2
t値 = -3.2851, 自由度 = 111.598, P値 = 0.001363
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -5.477495 -1.355838
標本推定値:
平均値x 平均値y
37.31667 40.73333

> var.test(gr1, gr2)
二群の等分散性の検定

データ: gr1 と gr2
F = 0.8795, 第1自由度 = 119, 第2自由度 = 59, P値 = 0.5499
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.5542763 1.3492261
標本推定値:
分散比
0.8794716

>

床に１−２時間寝ると、喉が痛くなる。ところが、椅子で寝るとおこらない。
気道程度、濡れ壁式集塵機で、線量が落ちるのではないか、ということではじめたのが以下の測定。
生データは余りにも大量すぎるので、掲載不可。大体半日ぐらいの連続測定(記録ソフトが15時間程度で異常停止するため)。

作成した集塵機は、
発泡スチロール製箱(25*15*10程度)で、15cm側を上側約5cmを切断。残る３面にタオルをぶら下げ、蓋をした。
蓋に、直径10cm程度の穴をあけ、卓上扇風機を下から上に風が動くように置く。
これを床面におき、通常の試料測定用机の線量変化を測定した。

扇風機動作前
> mean(x)
[1] 39.34407
> var(x)
[1] 44.45018
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 6 1.0169492 1.016949
25 23 3.8983051 4.915254
30 103 17.4576271 22.372881
35 202 34.2372881 56.610169
40 130 22.0338983 78.644068
45 83 14.0677966 92.711864
50 31 5.2542373 97.966102
55 9 1.5254237 99.491525
60 2 0.3389831 99.830508
65 0 0.0000000 99.830508
70 1 0.1694915 100.000000

扇風機動作後。
> mean(x)
[1] 36.46329
> var(x)
[1] 36.78771
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 27 2.27848101 2.278481
25 114 9.62025316 11.898734
30 317 26.75105485 38.649789
35 372 31.39240506 70.042194
40 248 20.92827004 90.970464
45 83 7.00421941 97.974684
50 18 1.51898734 99.493671
55 5 0.42194093 99.915612
60 1 0.08438819 100.000000

> mean(x)
[1] 36.352
> var(x)
[1] 35.81995
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
15 1 0.08 0.08
20 24 1.92 2.00
25 137 10.96 12.96
30 315 25.20 38.16
35 391 31.28 69.44
40 264 21.12 90.56
45 98 7.84 98.40
50 18 1.44 99.84
55 2 0.16 100.00

> mean(x)
[1] 37.09242
> var(x)
[1] 40.52478
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 11 2.0332717 2.033272
25 47 8.6876155 10.720887
30 140 25.8780037 36.598891
35 147 27.1719039 63.770795
40 132 24.3992606 88.170055
45 50 9.2421442 97.412200
50 10 1.8484288 99.260628
55 4 0.7393715 100.000000

> mean(x)
[1] 36.63696
> var(x)
[1] 35.487
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 13 1.430143 1.430143
25 87 9.570957 11.001100
30 223 24.532453 35.533553
35 303 33.333333 68.866887
40 201 22.112211 90.979098
45 64 7.040704 98.019802
50 17 1.870187 99.889989
55 0 0.000000 99.889989
60 1 0.110011 100.000000
> mean(x)
[1] 37.05952
> var(x)
[1] 39.63607
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 6 1.4285714 1.428571
25 40 9.5238095 10.952381
30 107 25.4761905 36.428571
35 122 29.0476190 65.476190
40 89 21.1904762 86.666667
45 43 10.2380952 96.904762
50 12 2.8571429 99.761905
55 1 0.2380952 100.000000

床に寝たときに、多少喉が楽になったが、放射能の影響か、湿度が上がったためか、は不明。
湿度計の目盛りは、5%(１目盛り)しか動かないので、測定誤差の範囲内。
大体１日で、300-500ccの水が蒸発する。
しばらく使用して、中に入れたタオルの線量を測定すればはっきりするかもしれない。
発泡スチロールの上からの測定では、BGと共に同じ uSv/h値の範囲(３０秒読み捨て、その後１分間の最大値と最小値を比較)。

分量が多いので、一部省略します
1. n=88 糊台
x <- c( 31,45,45,35,48,29,36,39,47,35,36,32,39,45,45,27,34,32,35,37,33,30,43,37,39,36,35,41,35,38,40,50,37,34,40,30,40,40,46,43,36,35,32,
43,32,32,45,35,37,33,59,38,33,45,36,42,34,38,26,41,18,42,34,35,38,37,32,33,36,27,25,38,45,41,35,40,38,29,38,41,28,45,38,24,36,40,48,60)

2. n=30 2012.11.11-2012.12.12の間使用した集塵機のフィルターに使ったタオル
x <- c( 39,33,44,34,38,41,37,38,34,29,42,37,49,46,37,36,32,48,45,40,51,29,36,34,33,44,43,34,42,41)

3. n=45 糊台
x <- c( 30,37,34,32,53,37,38,34,35,32,55,34,34,34,39,39,34,27,47,30,31,41,35,37,42,46,35,40,34,40,46,27,40,39,28,40,41,36,49,37,35,41,40,28,34)

4. n=30 2012.11.11-2012.12.12の間使用した集塵機のフィルターに使ったタオル
x <- c( 25,39,43,51,41,45,34,39,35,35,41,33,54,48,39,46,29,51,33,37,44,39,30,42,43,30,36,36,30,48)

5. n=30 糊台
x <- c( 39,40,31,38,40,37,42,37,46,46,38,33,29,41,33,31,37,41,32,35,41,52,34,36,40,32,30,32,46,36)

6. n=30 2012.11.11-2012.12.12の間使用した集塵機のフィルターに使ったタオル
x <- c( 29,39,30,34,43,42,35,32,47,36,49,40,35,34,43,45,41,36,35,36,39,30,42,35,39,37,34,38,33,35)

7. n=30 糊台
x <- c( 45,31,39,39,22,47,40,28,43,25,32,33,40,29,49,35,47,28,43,32,33,38,34,29,38,31,37,25,46,31)

8. n=30 2012.11.11-2012.12.12の間使用した集塵機のフィルターに使ったタオル
x <- c( 37,41,41,43,38,40,41,35,48,50,39,44,48,38,46,34,38,44,48,43,40,41,36,30,35,45,34,38,36,32)

9. n=30 糊台
x <- c( 47,46,38,37,37,38,43,29,40,34,29,25,40,40,40,31,37,44,33,37,49,34,44,40,41,41,48,33,25,32)

全体の分析
> mean(x)
[1] 37.793
> var(x)
[1] 40.14709
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 1.3354, 第1自由度 = 8, 第2自由度 = 334, P値 = 0.2248

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 10.1848, 自由度 = 8, P値 = 0.2523

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 1.4059, 第1自由度 = 8.000, 第2自由度 = 119.984, P値 = 0.2008

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
15 1 0.2915452 0.2915452
20 2 0.5830904 0.8746356
25 27 7.8717201 8.7463557
30 77 22.4489796 31.1953353
35 106 30.9037901 62.0991254
40 75 21.8658892 83.9650146
45 44 12.8279883 96.7930029
50 8 2.3323615 99.1253644
55 2 0.5830904 99.7084548
60 1 0.2915452 100.0000000
>
BGと試料の比較
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 5.7115, 第1自由度 = 1, 第2自由度 = 341, P値 = 0.01740

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 1.8656, 自由度 = 1, P値 = 0.1720

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 6.1058, 第1自由度 = 1.000, 第2自由度 = 267.874, P値 = 0.01410

>

繰り返しによる影響

1. n = 223 BG 繰り返し数 =5
> mean(x)
[1] 37.19731
> var(x)
[1] 42.5645
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 0.5177, 第1自由度 = 4, 第2自由度 = 218, P値 = 0.7228

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 2.7326, 自由度 = 4, P値 = 0.6035

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 0.4303, 第1自由度 = 4.000, 第2自由度 = 84.391, P値 = 0.7864

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
15 1 0.4484305 0.4484305
20 2 0.8968610 1.3452915
25 22 9.8654709 11.2107623
30 53 23.7668161 34.9775785
35 66 29.5964126 64.5739910
40 45 20.1793722 84.7533632
45 28 12.5560538 97.3094170
50 3 1.3452915 98.6547085
55 2 0.8968610 99.5515695
60 1 0.4484305 100.0000000
>

2. n = 120 2012.11.11-2012.12.12の間使用した集塵機のフィルターに使ったタオル 繰り返し数 =4
> mean(x)
[1] 38.9
> var(x)
[1] 34.07395
> oneway.test(x ~ g, var = T)
一元配置分散分析

データ: x と g
F = 1.0829, 第1自由度 = 3, 第2自由度 = 116, P値 = 0.3592

> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）

データ: x と g
バートレットのK二乗値 = 5.4507, 自由度 = 3, P値 = 0.1416

> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）

データ: x と g
F = 1.4109, 第1自由度 = 3.00, 第2自由度 = 63.95, P値 = 0.2477

> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
25 5 4.166667 4.166667
30 24 20.000000 24.166667
35 40 33.333333 57.500000
40 30 25.000000 82.500000
45 16 13.333333 95.833333
50 5 4.166667 100.000000
>

BGと試料の比較

1. n = 223 BG 繰り返し数 =5
2. n = 120 2012.11.11-2012.12.12の間使用した集塵機のフィルターに使ったタオル 繰り返し数 =4
> t.test(gr1, gr2, v=T)
二標本t検定（分散が等しいと仮定できるとき）

データ: gr1 と gr2
t値 = -2.3899, 自由度 = 341, P値 = 0.01740
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -3.1040585 -0.3013226
標本推定値:
平均値x 平均値y
37.19731 38.90000

> t.test(gr1, gr2)
二標本t検定（Welchの方法）

データ: gr1 と gr2
t値 = -2.471, 自由度 = 267.874, P値 = 0.01410
対立仮説: 母平均の差は，0ではない
95 パーセント信頼区間: -3.0593787 -0.3460025
標本推定値:
平均値x 平均値y
37.19731 38.90000

> var.test(gr1, gr2)
二群の等分散性の検定

データ: gr1 と gr2
F = 1.2492, 第1自由度 = 222, 第2自由度 = 119, P値 = 0.1771
対立仮説: 分散比は，1ではない
95 パーセント信頼区間: 0.9036585 1.7006307
標本推定値:
分散比
1.24918

>

これみればわかるとおり、少し線量が増えたのがタオル。
多少粉塵が減るみたい、という結果。
BGの度数分布の形を見ればわかるとおり、
高原状態から脱して、わずかに違う複数の線量差がある大気が、入れ替わりながら移動するという現象は、抑えられた気分。

インスペクターたすの相場はいくらくらい？

糊台の最近の数値。
暖房入れて換気が悪くなったせいか、bgが暴れて安定していない。
それで、測定ができないでいる。
bgの乱れを参考程度に。
37,34,39,41,38,32,39,39,34,41,46,47,30,38,38,44,46,38,37,30,37,49,48,38,44,48,47,32,41,32,29,39,41,30,39,36,
26,32,33,43,31,34,40,45,31,38,36,44,47,43,47,36,39,39,31,40,35,38,41,28,36,39,32,47,35,36,38,37,25,31,48,34,
39,35,31,40,36,36,37,30,43,34,41,38,32,37,44,40,46,28,31,36,38,34,48,46,39,40,32,43,37,34,40,39,41,31,45,37,
34,44,38,37,34,51,37,33,47,44,38,30,34,38,31,30,40,42,39,26,36,41,41,39,45,37,38,42,30,38,32,57,35,26,35,26,
38,32,37,36,44,41,40,37,33,34,46,44,37,38,31,41,40,42,49,29,30,31,25,37,47,40,40,22,33,38,32,32,43,47,62,26,
42,35,38,40,33,35,37,49,47,43,39,41,36,32,40,40,38,32,36,40,34,44,34,25,34,48,51,29,33,32,39,28,38,46,43,36,
41,38,33,30,37,38,34,35,39,49,35,32,40,42,35,31,40,44,34,32,48,45,35,36,47,33,47,44,35,43,30,44,41,44,33,37,
33,29,30,30,39,31,47,41,40,42,36,36,42,25,31,38,39,35,41,38,34,25,41,43,36,40,30,46,33,42,24,33,33,37,39,37,
39,36,28,28,36,30,45,50,24,23,46,40,23,33,38,36,29,26,42,49,44,23,35,32,31,38,40,30,44,41,42,51,32,33,43,31,
29,25,50,35,33,39,41,39,34,38,40,29,43,35,29,40,29,38,48,51,35,24,39,34,40,26,41,35,34,33,49,42,38,39,40,39,
35,34,43,33,34,40,38,37,28,46,32,39,32,27,44,42,27,42,34,41,31,39,36,32,44,36,30,34,36,38,44,35,45,41,45,32,
36,41,46,38,31,52,27,35,45,53,34,34,37,37,35,29,52,38,34,40,44,43,34,45,37,53,38,33,49,36,20,29,36,28,37,46,
32,39,38,44,34,42,32,35,40,44,37,35,37,50,28,38,33,31,28,32,46,29,33,35,50,39,37,33,28,45,48,30,34,30,34,31,
32,38,38,32,42,44,39,37,36,32,41,33,31,47,33,25,37,37,38,35,24,36,28,41,43,36,37,33,41,35,52,40,34,35,38,35,
53,40,47,37,42,31,48,42,33,30,37,36,32,43,34,38,34,28,26,32,44,38,44,39,37,48,38,34,45,35,27,37,40,30,36,41,
41,29,26,32,41,37,30,37,31,38,45,34,31,36,39,25,40,40,38,37,32,50,35,30,42,42,44,38,42,25,40,39,43,34,30,36,
41,33,42,46,41,35,34,33,32,27,30,36,49,29,38,22,40,36,44,38,28,37,33,33,32,35,33,35,35,46,30,27,37,47,36,37,
40,39,47,40,41,35,27,31,43,47,30,38,42,32,40,32,32,34,44,46,40,41,43,43,43,30,35,39,43,46,31,31,39,46,42,43,
42,47,42,43,35,29,25,32,45,39,47,39,34,37,38,37,31,44,42,37,39,45,39,37,41,38,37,42,38,33,38,46,28,34,43,33,
40,38,40,35,31,36,29,30,34,46,42,38,36,27,37,38,47,44,43,26,29,39,32,36,38,34,42,35,43,25,41,32,34,29,33,38,
34,39,42,35,46,44,35,30,37,37,34,32,31,29,46,29,36,30,28,27,29,55,27,39,41,39,49,33,46,37,36,34,37,45,30,34,
40,31,32,31,38,22,34,38,50,43,27,45,37,43,39,43,31,39,32,32,40,32,31,38,36,40,41,39,39,44,36,41,45,56,45,31,
39,38,44,27,45,50,33,48,43,41,44,26,40,38,34,43,48,39,39,36,48,36,27,24,32,36,37,34,38,40,41,40,38,41,33,28,
49,33,42,38,50,35,37,27,39,32,36,46,30,41,52,44,49,45,33,30,27,30,39,37,41,45,34,34,36,34,30,38,37,36,40,49,
40,36,36,34,39,39,48,43,37,30,37,54,41,35,52,49,30,35,31,40,35,41,35,35,36,33,45,39,33,47,31,40,38,28,41,36,
33,41,33,37,45,39,43,35,36,41,38,31,45,32,36,29,45,30,47,35,48,34,34,34,38,42,39,28,30,29,41,37,33,44,47,34,
37,35,36,36,39,38,30,35,43,42,32,41,25,46,44,28,38,37,36,43,40,37,47,40,37,40,30,43,42,33,30,40,39,32,37,41,
29,28,31,36,43,28,34,35,40,37,36,40,41,29,53,39,42,37,35,27,33,43,39,46,35,39,37,35,36,33,46,45,42,30,40,46,
38,35,35,30,38,37,36,39,43,39,51,27,40,35,30,44,31,43,31,27,33,37,39,43,36,44,34,43,38,39,31,33,36,30,30,36,
35,38,34,30,48,39,43,39,38,27,37,40,37,40,47,39,35,33,38,43,40,33,36,26,39,42,43,44,47,37,41,37,40,36,37,37,
32,37,45,39,41,39,35,39,37,48,32,29,47,32,30,31,32,40,37,42,40,34,37,32,29,38,40,44,43,40,41,45,42,39,41,27,
40,45,36,41,41,52,41,41,44,45,45,49,31,36,33,40,41,32,44,38,47,43,38,34,33,42,39,36,36,46,51,43,50,35,45,37,
41,38,32,37,38,37,28,36,34,40,27,50,43,37,46,25,25,50,40,34,43,38,41,43,48,34,35,50,35,38,27,43,46,40,31,32,
37,36,41,43,44,36,30,39,33,31,31,42,48,37,45,39,32

>>178
インスペクタースレで聞いてくれ。
http://uni.2ch.net/test/read.cgi/radiation/1331598144/l50
こっちなら、並行輸入店・日本正規代理店と取引のある人がいるから。

>163
栃木さんまだ測ってたんですねｗ
関西や瀬戸内海の食塩が汚染しているとのことですが、海は繋がっているんだから時間がたつにつれて、汚染するのは当然だと思います。
シュミレーションでは原発事故から300日後には日本周辺の太平洋側の海がすべて汚染。500日後には日本海側もすべて汚染。1000日後にはアメリカ西海岸全域の海が汚染。
まだ福島原発からは放射性物質が毎日でていますから今後さらにひどくなるでしょう。

参考
http://blog.goo.ne.jp/flyhigh_2012/e/99f9d328ef7dda22c196ee043404caed

自分の生命健康，未来を守りたいのなら選挙に行こう!


高濃度放射性物質拡散反対!

原発再稼動反対!

基本的人権まで無くそうと憲法改正を強行しようとしている売国奴政権は要らない!

徴兵制度及び国防軍化反対!

福島第一原発事故以来今も尚，毎時1000万ベクレル以上の
高濃度放射性物質を事故以来拡散し続けており，
事故の収束すら全くお手上げ状態にもかかわらず
原発再稼動を強行する諸悪の根源の売国奴政権の存続断固反対!

自民・公明・維新の会・みんなの党・民主は
日本国国民及び日本国自体を壊滅させる!

ほう

とちぎーーーーーーwww数字の羅列がちこみしろーーwwwwww
グンマーと抗争中で忙しいんか？wwww

意味もなく age

(栃木県)が、見放されて淋しいのか久しぶりに降臨ｗｗｗ

がんばれ！(栃木県)！

都議 樺山卓司（自殺扱いで変死）

6/30までブログで放射線量測定結果を発表