いろんなデータの分布をR(R言語)で可視化 ヒストグラムと散布図 世界の主な火山

Rコード

# 分布 ヒストグラムと散布図
# 使用するデータは、『データブック・オブ・ザ・ワールド 世界各国要覧と最新統計 2016 Vol.28』
# p.5

# 世界の主な火山
# 火山の名称は、この目的上、必要ないので省略

# 文字化けする場合
par(family = "Osaka")

世界の主な火山 <- data.frame(高度 = c(3676, 3428, 3142, 2891, 2821,
                                    2462, 1731, 1486, 813, 5895,
                                    4100, 3718, 2632, 4750, 3323,
                                    1491, 1281, 6908, 5896, 5452,
                                    5399, 5230, 4392, 4220, 3432,
                                    3187, 2550, 2047, 1397, 4170,
                                    2797, 2518, 1222, 3794))

世界の主な火山


x <- hist(世界の主な火山$高度, col = "steelblue", right = FALSE, main = "世界の主な火山の高度 単位 : m")
x


# 階級幅(区間、区間の幅)を500
y <- hist(世界の主な火山$高度, breaks = seq(0, 7000, 500), col = "steelblue", right = FALSE, main = "世界の主な火山の高度 単位 : m")
y


summary(世界の主な火山$高度)


plot(世界の主な火山$高度, ylab = "世界の主な火山$高度 単位 : m")

R Console

> # 分布 ヒストグラムと散布図
> # 使用するデータは、『データブック・オブ・ザ・ワールド 世界各国要覧と最新統計 2016 Vol.28』
> # p.5
> 
> # 世界の主な火山
> # 火山の名称は、この目的上、必要ないので省略
> 
> # 文字化けする場合
> par(family = "Osaka")
> 
> 世界の主な火山 <- data.frame(高度 = c(3676, 3428, 3142, 2891, 2821,
+                                     2462, 1731, 1486, 813, 5895,
+                                     4100, 3718, 2632, 4750, 3323,
+                                     1491, 1281, 6908, 5896, 5452,
+                                     5399, 5230, 4392, 4220, 3432,
+                                     3187, 2550, 2047, 1397, 4170,
+                                     2797, 2518, 1222, 3794))
> 
> 世界の主な火山
   高度
1  3676
2  3428
3  3142
4  2891
5  2821
6  2462
7  1731
8  1486
9   813
10 5895
11 4100
12 3718
13 2632
14 4750
15 3323
16 1491
17 1281
18 6908
19 5896
20 5452
21 5399
22 5230
23 4392
24 4220
25 3432
26 3187
27 2550
28 2047
29 1397
30 4170
31 2797
32 2518
33 1222
34 3794
> 
> 
> x <- hist(世界の主な火山$高度, col = "steelblue", right = FALSE, main = "世界の主な火山の高度 単位 : m")
> x
$breaks
[1]    0 1000 2000 3000 4000 5000 6000 7000

$counts
[1] 1 6 8 8 5 5 1

$density
[1] 2.941176e-05 1.764706e-04 2.352941e-04 2.352941e-04 1.470588e-04 1.470588e-04 2.941176e-05

$mids
[1]  500 1500 2500 3500 4500 5500 6500

$xname
[1] "世界の主な火山$高度"

$equidist
[1] TRUE

attr(,"class")
[1] "histogram"
> 
> 
> # 階級幅(区間、区間の幅)を500
> y <- hist(世界の主な火山$高度, breaks = seq(0, 7000, 500), col = "steelblue", right = FALSE, main = "世界の主な火山の高度 単位 : m")
> y
$breaks
 [1]    0  500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000

$counts
 [1] 0 1 5 1 2 6 5 3 4 1 3 2 0 1

$density
 [1] 0.000000e+00 5.882353e-05 2.941176e-04 5.882353e-05 1.176471e-04 3.529412e-04 2.941176e-04
 [8] 1.764706e-04 2.352941e-04 5.882353e-05 1.764706e-04 1.176471e-04 0.000000e+00 5.882353e-05

$mids
 [1]  250  750 1250 1750 2250 2750 3250 3750 4250 4750 5250 5750 6250 6750

$xname
[1] "世界の主な火山$高度"

$equidist
[1] TRUE

attr(,"class")
[1] "histogram"
> 
> 
> summary(世界の主な火山$高度)
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
    813    2476    3255    3360    4208    6908 
> 
> 
> plot(世界の主な火山$高度, ylab = "世界の主な火山$高度 単位 : m")

ヒストグラムスクリーンショット

階級幅(区間区間の幅)は、Rデフォルト(スタージェス)

f:id:my_notes:20170605191116p:plain

階級幅(区間区間の幅)は、500

f:id:my_notes:20170605191251p:plain

散布図のスクリーンショット

f:id:my_notes:20170605191329p:plain

参考文献

データブック オブ・ザ・ワールド―世界各国要覧と最新統計〈2016 Vol.28〉

データブック オブ・ザ・ワールド―世界各国要覧と最新統計〈2016 Vol.28〉