$30 off During Our Annual Pro Sale. View Details »

階層クラスタリングにおける仮説検定

 階層クラスタリングにおける仮説検定

saltcooky

May 23, 2020
Tweet

More Decks by saltcooky

Other Decks in Science

Transcript

  1. ϒʔτετϥοϓ๏ʹΑΔݕఆ ͭͷର৅ ͷҨ఻ࢠ ʹΑΓܗ੒͞ΕΔथܗਤ ܥ౷थ ͸௨Γ 5   

      ը૾Ҿ༻(Ұ෦վม) http://stat.sys.i.kyoto-u.ac.jp/titech/multiboot-j.html
  2. 3Ͱ΍ͬͯΈΔ 12 3Ͱ"6ݕఆΛߦ͏ͨΊʹ͸QWDMVTUύοέʔδΛར༻͢Δ ɹlibrary(pvclust) ɹlibrary(parallel) ɹcl <- makeCluster(detectCores()) #ฒྻԽͷ͓·͡ͳ͍ ɹ

    ɹ# આ໌ม਺͕ߦɺର৅͕ྻͷঢ়ଶʹ͢ΔͨΊʹసஔ ɹmnist_df_t <- ɹɹɹɹmnist_df %>% ɹɹɹɹdplyr::select(-label) %>% ɹɹɹɹt() ɹcolnames(mnist_df_t) <- mnist_df$label ɹsa <- 9^seq(-1,1,length=13)ɹ# Ϛϧνεέʔϧϒʔτετϥοϓͷεέʔϧ഑ྻΛੜ੒ ɹmnist_boot <- pvclust(data = mnist_df_t, ɹ r = 1/sa, ɹ nboot = 2000, ɹ method.hclust = "ward.D2", ɹɹɹɹɹɹɹɹɹɹɹɹɹparallel = cl)
  3. 3Ͱ΍ͬͯΈΔ 13 ݁ՌΛ֬ೝ ɹ> mnist_boot ɹCluster method: ward.D2 ɹDistance :

    euclidean ɹEstimates on edges: ɹ si au bp se.si se.au se.bp v c pchi ɹ1 0.999 1.000 0.999 0.000 0.000 0.000 -3.175 0.156 0.000 ɹ2 0.979 0.990 0.971 0.001 0.001 0.001 -2.117 0.221 0.000 ɹ3 0.579 0.801 0.759 0.004 0.002 0.001 -0.775 0.073 0.000 ɹ4 0.392 0.718 0.655 0.003 0.002 0.001 -0.487 0.089 0.000 ɹ5 0.771 0.900 0.852 0.004 0.002 0.002 -1.151 0.105 0.000 ɹ6 0.112 0.634 0.455 0.003 0.002 0.001 -0.119 0.231 0.000 ɹ7 0.213 0.795 0.323 0.004 0.002 0.002 -0.182 0.641 0.000 ɹ8 0.532 0.770 0.755 0.004 0.002 0.001 -0.716 0.025 0.000 ɹ9 0.266 0.779 0.394 0.004 0.002 0.002 -0.250 0.520 0.000 ɹ10 0.915 0.969 0.878 0.003 0.001 0.002 -1.518 0.354 0.183 ɹ11 0.223 0.757 0.397 0.004 0.002 0.001 -0.211 0.471 0.000 ɹ12 0.229 0.704 0.475 0.003 0.002 0.001 -0.237 0.299 0.000 ɹ13 0.246 0.752 0.409 0.004 0.002 0.001 -0.235 0.464 0.017 ɹ14 0.229 0.691 0.492 0.003 0.002 0.001 -0.239 0.258 0.000