インスペクター系総合 ３【inspector+,Alert】
332 ：名無しに影響はない(栃木県)：
プラスチックによる遮蔽に必要な厚さの測定。
机の上に、CDケースに入れた、５円玉の穴の中にやさしおを充填し紙製シールで密封したものを置き
この上に、CDケースを何枚か置くことによって、やさしおから出るβ線の遮蔽に必要なプラスチックの厚さを求めた。
分量が多いので、一部分だけ。
繰り返しによる差異
BG
x <- c(41, 35, 40, 41, 42, 46, 32, 37, 46, 31, 36, 36, 41, 35, 48, 37, 40, 37, 49, 42, 35, 26, 42, 37,
32, 41, 45, 44, 32, 35, 48, 36, 45, 45, 37, 36, 39, 35, 42, 26, 43, 45, 48, 37, 37, 37, 47, 42, 34, 43,
37, 45, 35, 42, 35, 48, 52, 39, 39, 28, 45, 43, 43, 32, 31, 33, 41, 28, 43, 28, 34, 22, 31, 42, 33, 37, 35, 31, 33, 27, 37, 51, 37, 37, 38, 46, 36, 41, 32, 44 ,
45, 32, 43, 42, 55, 49, 25, 45, 28, 48, 48, 37, 39, 51, 36, 47, 50, 49, 46, 32, 35, 38, 44, 30, 40, 53,
41, 37, 51, 37, 41, 40, 31, 39, 38, 46, 35, 43, 30, 47, 39, 39, 44, 37, 36, 39, 36, 40, 42, 37, 50 ,
39, 52, 32, 38, 36, 31, 43, 45, 47, 30, 33, 32, 21, 44, 35, 52, 35, 40, 30, 41, 45, 50, 37, 42, 49, 40,
47, 41, 36, 35, 46, 35, 42, 52, 43, 41, 40, 36, 50, 35, 38, 39, 35, 42, 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,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,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 ,
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)
> mean(x)
[1] 39.4086
> var(x)
[1] 42.51322
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.4755, 第1自由度 = 2, 第2自由度 = 183, P値 = 0.08693
有意。繰り返しによる差異があり
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.8018, 自由度 = 2, P値 = 0.6697
有意。分散が一様でなく先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.4383, 第1自由度 = 2.000, 第2自由度 = 96.759, P値 = 0.09266
有意。繰り返しによる差異があり
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 2 1.0752688 1.075269
25 8 4.3010753 5.376344
30 25 13.4408602 18.817204
35 61 32.7956989 51.612903
40 46 24.7311828 76.344086
45 31 16.6666667 93.010753
50 12 6.4516129 99.462366
55 1 0.5376344 100.000000
>
机の上に、CDケースに入れた、５円玉の穴の中にやさしおを充填し紙製シールで密封したものを置き
この上に、CDケースを何枚か置くことによって、やさしおから出るβ線の遮蔽に必要なプラスチックの厚さを求めた。
分量が多いので、一部分だけ。
繰り返しによる差異
BG
x <- c(41, 35, 40, 41, 42, 46, 32, 37, 46, 31, 36, 36, 41, 35, 48, 37, 40, 37, 49, 42, 35, 26, 42, 37,
32, 41, 45, 44, 32, 35, 48, 36, 45, 45, 37, 36, 39, 35, 42, 26, 43, 45, 48, 37, 37, 37, 47, 42, 34, 43,
37, 45, 35, 42, 35, 48, 52, 39, 39, 28, 45, 43, 43, 32, 31, 33, 41, 28, 43, 28, 34, 22, 31, 42, 33, 37, 35, 31, 33, 27, 37, 51, 37, 37, 38, 46, 36, 41, 32, 44 ,
45, 32, 43, 42, 55, 49, 25, 45, 28, 48, 48, 37, 39, 51, 36, 47, 50, 49, 46, 32, 35, 38, 44, 30, 40, 53,
41, 37, 51, 37, 41, 40, 31, 39, 38, 46, 35, 43, 30, 47, 39, 39, 44, 37, 36, 39, 36, 40, 42, 37, 50 ,
39, 52, 32, 38, 36, 31, 43, 45, 47, 30, 33, 32, 21, 44, 35, 52, 35, 40, 30, 41, 45, 50, 37, 42, 49, 40,
47, 41, 36, 35, 46, 35, 42, 52, 43, 41, 40, 36, 50, 35, 38, 39, 35, 42, 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,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,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 ,
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)
> mean(x)
[1] 39.4086
> var(x)
[1] 42.51322
> oneway.test(x ~ g, var = T)
一元配置分散分析
データ: x と g
F = 2.4755, 第1自由度 = 2, 第2自由度 = 183, P値 = 0.08693
有意。繰り返しによる差異があり
> bartlett.test(x, g)
分散の一様性の検定（バートレット検定）
データ: x と g
バートレットのK二乗値 = 0.8018, 自由度 = 2, P値 = 0.6697
有意。分散が一様でなく先の分散分析が無効。
> oneway.test(x ~ g)
一元配置分散分析（等分散を仮定しない場合）
データ: x と g
F = 2.4383, 第1自由度 = 2.000, 第2自由度 = 96.759, P値 = 0.09266
有意。繰り返しによる差異があり
> dosuu.bunpu(x, 5)
freq pcnt cum.pcnt
20 2 1.0752688 1.075269
25 8 4.3010753 5.376344
30 25 13.4408602 18.817204
35 61 32.7956989 51.612903
40 46 24.7311828 76.344086
45 31 16.6666667 93.010753
50 12 6.4516129 99.462366
55 1 0.5376344 100.000000
>
|
|
|
333 ：名無しに影響はない(栃木県)：2012/04/04(水) 19:48:09.78 ID:QmVsyECa
回帰分析
y <- c(38, 33, 39, 39, 37, 30, 44, 43, 37, 40, 47, 35, 37, 43, 40, 40, 37, 43, 44, 42, 39, 47, 46, 40, 41, 43, 40, 43, 36, 33 ,
52, 47, 35, 51, 28, 35, 34, 39, 38, 41, 38, 36, 43, 27, 39, 46, 42, 39, 37, 41, 41, 44, 37, 37, 45, 40, 42, 26, 25, 44 ,
41, 42, 33, 41, 41, 43, 36, 31, 39, 30, 33, 38, 44, 38, 42, 28, 36, 35, 38, 54, 41, 44, 53, 34, 49, 50, 37, 30, 39, 34 ,
52, 39, 31, 36, 42, 46, 42, 39, 64, 38, 46, 42, 36, 30, 32, 42, 29, 36, 50, 40, 35, 35, 37, 30, 29, 40, 43, 39, 51, 45,
47, 33, 50, 30, 48, 44, 34, 28, 35, 31, 35, 31, 35, 39, 46, 45, 47, 30, 31, 44, 29, 31, 27, 45, 33, 33, 33, 41, 41, 35,
31, 34, 32, 40, 36, 40, 27, 39, 34, 33, 39, 43, 29, 46, 34, 32, 35, 36, 43, 30, 43, 35, 35, 41, 42, 37, 47, 38, 33, 44,
43, 31, 52, 40, 33, 42, 44, 41, 50, 47, 36, 29, 48, 53, 31, 45, 52, 25, 35, 32, 44, 35, 43, 43, 40, 48, 44, 36, 36, 37,
38, 40, 42, 45, 42, 45, 41, 42, 35, 27, 36, 43, 46, 35, 44, 32, 42, 31, 37, 37, 47, 40, 47, 50, 31, 34, 40, 31, 32, 35,
44, 28, 52, 32, 28, 44, 33, 43, 32, 42, 31, 23, 51, 39, 38, 42, 41, 43, 35, 44, 31, 38, 34, 44, 34, 42, 40, 54, 42, 39, 38, 43, 40, 41, 44, 34, 43, 48, 37, 43, 42, 30, 36, 34, 48, 39 ,
39, 34, 44, 32, 36, 39, 35, 30, 28, 38, 40, 47, 36, 33, 44, 57, 37, 46, 38, 40, 57, 43, 43, 45, 35, 40, 51, 30, 37, 44, 38, 33, 45, 39 ,
53, 41, 37, 55, 56, 55, 42, 47, 48, 70, 51, 33, 53, 51, 62, 65, 46, 53, 49, 37, 49, 42, 56, 52, 53, 59, 59, 40, 54, 43, 53, 42 ,
46, 42, 49, 52, 39, 31, 36, 32, 43, 37, 51, 39, 46, 60, 37, 44, 47, 38, 44, 44, 49, 57, 46, 44, 49, 38, 37, 30, 40, 53, 42, 28, 43, 47, 48, 39, 40, 43, 30, 47, 51, 46, 45 ,
32, 40, 39, 44, 27, 44, 38, 26, 35, 32, 39, 35, 40, 36, 29, 42, 42, 31, 46, 55, 47, 51, 42, 34, 48, 29, 45, 37, 39, 35, 48, 38, 32, 34, 42, 48, 45, 43, 35, 40, 30, 49, 34, 35, 38, 50, 46, 43, 53, 30 ,
51, 45, 45, 44, 37, 32, 43, 36, 36, 36, 33, 34, 33, 47, 38, 34, 48, 50, 34, 37, 39, 33, 42, 44, 33, 58, 32, 33, 46, 38, 49, 36, 41, 43, 42, 34, 37, 43, 23, 40, 36, 36, 39, 42, 46, 36, 43, 37, 38, 39, 38)
y <- c(38, 33, 39, 39, 37, 30, 44, 43, 37, 40, 47, 35, 37, 43, 40, 40, 37, 43, 44, 42, 39, 47, 46, 40, 41, 43, 40, 43, 36, 33 ,
52, 47, 35, 51, 28, 35, 34, 39, 38, 41, 38, 36, 43, 27, 39, 46, 42, 39, 37, 41, 41, 44, 37, 37, 45, 40, 42, 26, 25, 44 ,
41, 42, 33, 41, 41, 43, 36, 31, 39, 30, 33, 38, 44, 38, 42, 28, 36, 35, 38, 54, 41, 44, 53, 34, 49, 50, 37, 30, 39, 34 ,
52, 39, 31, 36, 42, 46, 42, 39, 64, 38, 46, 42, 36, 30, 32, 42, 29, 36, 50, 40, 35, 35, 37, 30, 29, 40, 43, 39, 51, 45,
47, 33, 50, 30, 48, 44, 34, 28, 35, 31, 35, 31, 35, 39, 46, 45, 47, 30, 31, 44, 29, 31, 27, 45, 33, 33, 33, 41, 41, 35,
31, 34, 32, 40, 36, 40, 27, 39, 34, 33, 39, 43, 29, 46, 34, 32, 35, 36, 43, 30, 43, 35, 35, 41, 42, 37, 47, 38, 33, 44,
43, 31, 52, 40, 33, 42, 44, 41, 50, 47, 36, 29, 48, 53, 31, 45, 52, 25, 35, 32, 44, 35, 43, 43, 40, 48, 44, 36, 36, 37,
38, 40, 42, 45, 42, 45, 41, 42, 35, 27, 36, 43, 46, 35, 44, 32, 42, 31, 37, 37, 47, 40, 47, 50, 31, 34, 40, 31, 32, 35,
44, 28, 52, 32, 28, 44, 33, 43, 32, 42, 31, 23, 51, 39, 38, 42, 41, 43, 35, 44, 31, 38, 34, 44, 34, 42, 40, 54, 42, 39, 38, 43, 40, 41, 44, 34, 43, 48, 37, 43, 42, 30, 36, 34, 48, 39 ,
39, 34, 44, 32, 36, 39, 35, 30, 28, 38, 40, 47, 36, 33, 44, 57, 37, 46, 38, 40, 57, 43, 43, 45, 35, 40, 51, 30, 37, 44, 38, 33, 45, 39 ,
53, 41, 37, 55, 56, 55, 42, 47, 48, 70, 51, 33, 53, 51, 62, 65, 46, 53, 49, 37, 49, 42, 56, 52, 53, 59, 59, 40, 54, 43, 53, 42 ,
46, 42, 49, 52, 39, 31, 36, 32, 43, 37, 51, 39, 46, 60, 37, 44, 47, 38, 44, 44, 49, 57, 46, 44, 49, 38, 37, 30, 40, 53, 42, 28, 43, 47, 48, 39, 40, 43, 30, 47, 51, 46, 45 ,
32, 40, 39, 44, 27, 44, 38, 26, 35, 32, 39, 35, 40, 36, 29, 42, 42, 31, 46, 55, 47, 51, 42, 34, 48, 29, 45, 37, 39, 35, 48, 38, 32, 34, 42, 48, 45, 43, 35, 40, 30, 49, 34, 35, 38, 50, 46, 43, 53, 30 ,
51, 45, 45, 44, 37, 32, 43, 36, 36, 36, 33, 34, 33, 47, 38, 34, 48, 50, 34, 37, 39, 33, 42, 44, 33, 58, 32, 33, 46, 38, 49, 36, 41, 43, 42, 34, 37, 43, 23, 40, 36, 36, 39, 42, 46, 36, 43, 37, 38, 39, 38)
334 ：名無しに影響はない(栃木県)：2012/04/04(水) 19:48:40.99 ID:QmVsyECa
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 , 1,1,1,1,1,1,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,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,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,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3 ,
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,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 ,
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 ,
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,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)
Call:
lm(formula = y ~ x)
Coefficients:
(Intercept) x
41.3532 -0.6165
> plot(x,y,main="Regression line")
> abline(result,col="red")
>
BGが39.4086だから
(41.3532-39.4086)/0.6165 =3.15
よって、CDケース3.15枚でバックグラウンドと同じ線量に落とせるから、
５円玉の穴に詰め込んだやさしおの線量を遮蔽するのに必要なプラスチックの分量は 3.15*2+1=7.3mmとなる。
>>327-329 より、「検知部が受光可能な範囲を中心部(検出部中心線上の、検出器から2.5cm離れた点)から6.6cm離れた円周上とする。」
今回の測定よりβ線の受光可能深さが 7.3mmだから、深さ3.65mm受光可能な範囲を底面積とする円柱と、β線を検出可能な試料の大きさを近似する(本来は下面が球面となるはず)。
机面からの距離(cm) 受光範囲の半径(cm) 円柱の容積
0.0 2.5 7.16
0.5 3.32 12.63
1.0 4.14 19.65
1.5 4.95 28.09
2.0 5.78 38.30
2.5 6.63 50.40
3.0 7.42 63.13
円柱体積の計算は http://www.benricho.org/calculate/Cylinder.html による。
2,2,2,2,2,2,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,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,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 ,
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,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 ,
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 ,
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,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)
Call:
lm(formula = y ~ x)
Coefficients:
(Intercept) x
41.3532 -0.6165
> plot(x,y,main="Regression line")
> abline(result,col="red")
>
BGが39.4086だから
(41.3532-39.4086)/0.6165 =3.15
よって、CDケース3.15枚でバックグラウンドと同じ線量に落とせるから、
５円玉の穴に詰め込んだやさしおの線量を遮蔽するのに必要なプラスチックの分量は 3.15*2+1=7.3mmとなる。
>>327-329 より、「検知部が受光可能な範囲を中心部(検出部中心線上の、検出器から2.5cm離れた点)から6.6cm離れた円周上とする。」
今回の測定よりβ線の受光可能深さが 7.3mmだから、深さ3.65mm受光可能な範囲を底面積とする円柱と、β線を検出可能な試料の大きさを近似する(本来は下面が球面となるはず)。
机面からの距離(cm) 受光範囲の半径(cm) 円柱の容積
0.0 2.5 7.16
0.5 3.32 12.63
1.0 4.14 19.65
1.5 4.95 28.09
2.0 5.78 38.30
2.5 6.63 50.40
3.0 7.42 63.13
円柱体積の計算は http://www.benricho.org/calculate/Cylinder.html による。