［R］サンプルのデータセット「面積が1万平方マイルを超える陸塊」

datasetsパッケージのislandsは、面積が1万平方マイルを超える陸塊をまとめたサンプルのデータセット。引用元は以下のとおり。

The World Almanac and Book of Facts, 1975, page 406

> length(islands)
[1] 48
> head(islands)
      Africa   Antarctica         Asia    Australia Axel Heiberg 
       11506         5500        16988         2968           16 
      Baffin 
         184 
> tail(islands)
          Taiwan         Tasmania Tierra del Fuego            Timor 
              14               26               19               13 
       Vancouver         Victoria 
              12               82 
> names(islands)
 [1] "Africa"           "Antarctica"       "Asia"            
 [4] "Australia"        "Axel Heiberg"     "Baffin"          
 [7] "Banks"            "Borneo"           "Britain"         
[10] "Celebes"          "Celon"            "Cuba"            
[13] "Devon"            "Ellesmere"        "Europe"          
[16] "Greenland"        "Hainan"           "Hispaniola"      
[19] "Hokkaido"         "Honshu"           "Iceland"         
[22] "Ireland"          "Java"             "Kyushu"          
[25] "Luzon"            "Madagascar"       "Melville"        
[28] "Mindanao"         "Moluccas"         "New Britain"     
[31] "New Guinea"       "New Zealand (N)"  "New Zealand (S)" 
[34] "Newfoundland"     "North America"    "Novaya Zemlya"   
[37] "Prince of Wales"  "Sakhalin"         "South America"   
[40] "Southampton"      "Spitsbergen"      "Sumatra"         
[43] "Taiwan"           "Tasmania"         "Tierra del Fuego"
[46] "Timor"            "Vancouver"        "Victoria"        

全部で48個納められている。単位は平方マイルであることに注意。これを平方メートルに変換するには、1マイルは約1609mであるので、それぞれ1609^2=2588881、平方キロメールに変換するには1.609^2=2.588881で割ればよい。北海道、本州、九州で試してみる。

> s <- c("Hokkaido", "Honshu", "Kyushu")
> islands[s]
Hokkaido   Honshu   Kyushu 
      30       89       14 
> islands[s] / 1609 ^ 2
    Hokkaido       Honshu       Kyushu 
1.158802e-05 3.437779e-05 5.407742e-06 
> islands[s] / 1.609 ^ 2
 Hokkaido    Honshu    Kyushu 
11.588018 34.377787  5.407742 

ニュージーランドの北島と南島でも試してみる。

> s <- grep("New Zealand", names(islands))
> islands[s]
New Zealand (N) New Zealand (S) 
             44              58 
> islands[s] / 1609 ^ 2
New Zealand (N) New Zealand (S) 
   1.699576e-05    2.240350e-05 
> islands[s] / 1.609 ^ 2
New Zealand (N) New Zealand (S) 
       16.99576        22.40350