注目メッシュの隣接メッシュを検索&縦持ち化
# 注目したい場所の緯度経度
coord_center <- tibble(
id = 1:4,
X = c(140.6041855, 130.4140618, 134.0464118, 139.7016358),
Y = c(35.9323516, 33.1568919, 34.31453200000001, 35.6580339)
)
# 注目場所の中心・周囲のメッシュコードを検索
mesh_center_and_neighbor <- coord_center %>%
mutate(
mesh_center = map2_chr(
.x = X,
.y = Y,
.f = ~ {
jpgrid::coords_to_grid(
X = .x,
Y = .y,
grid_size = '1km'
) %>%
as.character()
}
)
) %>%
mutate(
mesh_neighbor = map(
.x = mesh_center,
.f = ~ {
jpgrid::grid_neighbor(
jpgrid::parse_grid(.x),
n = 1L
)
}
)
) %>%
unnest_longer(
mesh_neighbor
) %>%
unnest_longer(
mesh_neighbor
) %>%
mutate(
mesh_neighbor = as.character(mesh_neighbor)
)
# メッシュの中心緯度経度を検索
target_meshes_center <- mesh_center_and_neighbor %>%
distinct(
id,
mesh_code_3rd = mesh_center
) %>%
mutate(
is_center = 1
)
target_meshes_neighbor <- mesh_center_and_neighbor %>%
distinct(
id,
mesh_code_3rd = mesh_neighbor
) %>%
mutate(
is_center = 0
)
target_meshes <- bind_rows(
target_meshes_center,
target_meshes_neighbor
)
回遊データをラインポリゴン化する方法
# 回遊データをラインポリゴン化
kaiyuu_line_polygons <- postings_with_kaiyuu_id %>%
sf::st_as_sf(
coords = c('longitude', 'latitude')
) %>%
group_by( # group_by + summariseをすると個々の投稿座標がMULTIPOINTでまとまる
kaiyuu_id
) %>%
summarise(
dummy = n(),
do_union = F # 元のデータの順番でMULTIPOINTが作られる
) %>%
sf::st_cast(
'LINESTRING'
) %>%
sf::st_set_crs(4326) %>% # マイマップはWGS84なので
select(
!dummy # 不要な変数の削除
)
共通の属性をもつポイントデータ同士を1つのラインポリゴンに統合する方法
postings_with_kaiyuu_id %>%
sf::st_as_sf(
coords = c('longitude', 'latitude')
) %>%
group_by( # group_by + summariseをすると個々の投稿座標がMULTIPOINTでまとまる
kaiyuu_id
) %>%
summarise(
dummy = n(),
do_union = F # 元のデータの順番でMULTIPOINTが作られる
) %>%
sf::st_cast(
'LINESTRING'
) %>%
sf::st_set_crs(4326) %>% # マイマップはWGS84なので
select(
!dummy # 不要な変数の削除
)
教育は遺伝に勝てるか?
面白かった豆知識
- 一卵性双生児と二卵性双生児のそれぞれにおいて、双子間の形質の相関を比較することで、形質の発現に対する遺伝と環境の影響の重みを分けられるらしい
- 一卵性双生児は同じ環境で育ち同じ遺伝子をもつが、二卵性双生児は同じ環境で育ち部分的に同じ遺伝子をもつ。そのため、一卵性双生児間の相関と二卵性双生児間の相関の差分は、双子間の遺伝の共有度合いの違いによって説明できる
- さらに一卵性双生児で双子間の相関が1に達しないかった部分は、双子間で共有されていない要因、すなわち非共有環境(例:家庭以外の環境)が発現に関わっていると推測できる
- そして残りの部分が共有環境(例:家庭環境など)が発現に関わっていることが推測できる
- これらを組み合わせると、形質発現に対する遺伝・共有環境・非共有環境の影響度合いを分離できる
形質発現に対する遺伝・共有環境・非共有環境の影響度合い
ggplotでNSE処理
for (var in vars) {
# var_1 <- rlang::enquo(var)
p <- ggplot(data = target_gridcells_with_landcover) + geom_sf(aes(fill = !!sym(var)))
ggsave( p, filename = str_c( 'DIR', var, '.jpg' ), width = 8, height = 8, dpi = 300 ) }