Upgrade to Pro — share decks privately, control downloads, hide ads and more …

ベイズ・EMアルゴリズムの紹介

yasubei
September 23, 2019

 ベイズ・EMアルゴリズムの紹介

初心者向けの勉強会のため前半はデータサイエンスに関する紹介をしています。
後半からベイズ統計学とEMアルゴリズムの紹介を行っています。
気になるところだけ見たほうが効率よいと思います。

動画はこちらです。
https://www.youtube.com/watch?v=DM2zyS4aaRg&t=710

yasubei

September 23, 2019
Tweet

More Decks by yasubei

Other Decks in Science

Transcript

  1. ࠓ ि ͷ Θ ͔ Β Μ ઌ ੜ ϕΠζ౷ܭֶ

    ͕ Θ͔ΒΜ 1ZUIPO#FHJOOFSTԭೄΦʔΨφΠβʔ ΍͢΂͐
  2. ࣗݾ঺հ :BTVCFJ ΍͢΂͐ ZBTVCFJUXJ ZBTVIBSVTV[VLJ IUUQTIVHLVODPN ΤϯδχΞͰ͢ɻ ࡳຈͰήʔϜ։ൃ೥ɹˠɹԭೄͰ։ൃ෼ੳͳͲΛ೥ )VHLVOͱ)VHLVO%BUBͱ͍͏ձࣾΛ΍͍ͬͯ·͢ɻ اۀ͞Μͱֶੜ͞ΜͱҰॹʹ࢓ࣄͯ͠8JO8JO8JOͳؔ܎Λங͘ͷ͕޷͖Ͱ͢ɻ

    ࠷ۙ͸ҿΈձͷ૿Ճʹ೰ΜͰ͓ΓɺςΩύΩؾ͕ར͘पΓͷਓͨͪͱൺֱͯ͠ɺ૬มΘΒͣؾ ͕ར͔ͳ͍ࣗ෼ʹࣦ๬͍ͯ͠Δɻ ؤுͬͯτʔΫΛ͢ΔΑ͏ʹ৺͕͚͍ͯΔ͕ɺͨͿΜۭؾ͸ಡΊ͍ͯͳ͍ͱࢥΘΕΔɻ ࠷ۙ޷͖ͳϚϯΨ͸౳෼ͷՖՇɺ੨ͷΦʔέετϥɺༀ԰ͷͻͱΓ͝ͱɻ
  3. ࣾձత੹೚ɾࣄྫ ಺ఆࣙୀ཰໰୊͚ͩͰ͸ͳ͍ʂʮϦΫφϏγϣοΫʯͰਓࡐۀքશମ͕਒্͕͑Δཧ༝ɹIUUQTBTDJJKQFMFN ϦΫϧʔτΩϟϦΞʮϦΫφϏब׆ࣄ݅ʯ͕ϓϥΠόγʔ໰୊ʹ౤͔͚͛ͨ΋ͷIUUQTOFXTZBIPPDPKQCZMJOFZBNBNPUPJDIJSP ࢀߟ63- ࢀߟɿ৬ۀ҆ఆ๏ ୈޒेҰ৚ɹ৬ۀ঺հࣄۀऀɺٻਓऀɺ࿑ಇऀͷืूΛߦ͏ऀɺืूडୗऀɺ࿑ಇऀڙڅࣄۀऀٴͼ࿑ಇऀڙڅΛड͚Α͏ͱ͢ΔऀʢҎԼ͜ͷ৚ʹ͓͍ͯʮ৬ۀ঺հࣄۀऀ౳ʯͱ͍͏ɻʣฒͼʹ͜ΕΒͷ୅ཧਓɺ࢖༻ਓͦͷଞͷैۀऀ͸ɺਖ਼౰ͳཧ ༝ͳ͘ɺͦͷۀ຿্औΓѻͭͨ͜ͱʹ͍ͭͯ஌ΓಘͨਓͷൿີΛ࿙Βͯ͠͸ͳΒͳ͍ɻ৬ۀ঺հࣄۀऀ౳ٴͼ͜ΕΒͷ୅ཧਓɺ࢖༻ਓͦͷଞͷैۀऀͰͳ͘ͳͭͨޙʹ͓͍ͯ΋ɺಉ༷ͱ͢Δɻ ˓̎ɹ৬ۀ঺հࣄۀऀ౳ٴͼ͜ΕΒͷ୅ཧਓɺ࢖༻ਓͦͷଞͷैۀऀ͸ɺલ߲ͷൿີͷ΄͔ɺͦͷۀ຿ʹؔͯ͠஌Γಘͨݸਓ৘ใͦͷଞްੜ࿑ಇলྩͰఆΊΔऀʹؔ͢Δ৘ใΛɺΈͩΓʹଞਓʹ஌Βͤͯ͸ͳΒͳ͍ɻ৬ۀ঺հࣄۀऀ౳ٴͼ͜ΕΒͷ ୅ཧਓɺ࢖༻ਓͦͷଞͷैۀऀͰͳ͘ͳͭͨޙʹ͓͍ͯ΋ɺಉ༷ͱ͢Δɻ ɾब৬αΠτʮϦΫφϏʯͷར༻ཤྺΛ%.1ϑΥϩʔͱ͍͏ผαʔϏεͰར༻ɻ಺ఆࣙୀ཰

    Λఏڙ͓ͯ͠Γɺब׆ੜͷ஌Βͳ͍ͱ͜ΖͰब৬׆ಈ͕ෆརʹͳΔՄೳੑ͕͋ͬͨɻ ʢސ٬֤ࣾ͸࠾༻ɾෆ࠾༻ͷ൑அʹӨڹͨ͠ࣄ࣮Λ൱ఆ͍ͯ͠Δʣ ɾ৬ۀ҆ఆ๏ҧ൓ୈ৚߲ʹҧ൓ Լهࢀߟ ɻް࿑ল͔Βߦ੓ࢦಋɻ ɾຊਓಉҙ͕͋ͬͯ΋Ξ΢τʢ΄΅Չ઎ঢ়ଶͷͨΊʣ ɾϚΠφϏͷྨࣅαʔϏεʮ13BJ0ʯ΋಺ఆࣙୀ཰Λࢉग़͍ͯ͠Δ͕ɺೖྗιʔε͕Τϯτ ϦʔγʔτͷͨΊ໰୊ʹͳͬͯ͸͍ͳ͍ɻʢͱ͸ֶ͍͑ੜ͕ϚΠφϏʹొ࿥ͨ͠σʔλΛ໨ త֎Ͱ࢖༻͍ͯ͠ΔՄೳੑ͕͋ΓɺάϨʔͰ͸͋Δʣ ɾʮࣗ෼ͷݸਓ৘ใ͸͋Δఔ౓औΒΕͯ΋ྑ͍ʯͱߟ͑Δਓ͕͍Δ͕ɺָ؍తʹߟ͑͗͢Δ ͷ͸/(Ͱ͋Δɻ ϦΫφϏʮ಺ఆࣙୀ཰ʯఏڙ໰୊
  4. ϕΠζ౷ܭֶͱ͸ʁ هड़౷ܭֶͱਪܭ౷ܭֶͷҧ͍ɹIUUQTUPLFJOFUCBTJDJOEVDUJWFTUBUJTUJDT ϕΠζ౷ܭֶͱ͸ʁॳ৺ऀ޲͚ͷ΍͍͞͠ղઆɹIUUQTUPLFJOFUCBTJDCBZFTJBOTUBUJTUJDT ϕΠζ౷ܭͱ͸ʁී௨ͷ౷ܭͱԿ͕ҧ͏ʁపఈղઆʂɹIUUQTVEFNZCFOFTTFDPKQBJCBZFTJBOTUBUJTUJDTIUNM ࢀߟ63- छྨ ֓ཁ ྫ هड़౷ܭֶ σʔλͷಛ௃ΛΘ͔Γ΍͘͢දݱ͢Δֶ໰

    தԝ஋ɾඪ४ภࠩ ՄࢹԽ ਪܭ౷ܭֶ ඪຊΛ෼ੳͯ͠฼ूஂʹ͍ͭͯਪଌ͢Δֶ໰ ̓ݕఆɾРೋ৐ݕఆɾ૬ؔ෼ੳ ճؼɾ෼ྨɾΫϥελϦϯά ϕΠζ౷ܭֶ ৘ใ͕ೖΔͨͼʹ֬཰Λมߋ͍͖ͯ͠ σʔλͷશମ૾Λ௫ΜͰֶ͍͘໰ ϕΠζͷఆཧ ࠷໬ਪఆɾ."1ਪఆɾϕΠζਪఆ ਖ਼ن෼෍ͷਪఆɾΧςΰϦΧϧ෼෍ͷਪఆ ࣝผϞσϧɾੜ੒Ϟσϧ (..ɾ̓෼෍ɾҼࢠ෼ੳɾ&.ΞϧΰϦζϜ Ψ΢εաఔɾ࠷ٸ߱Լ๏ɾϒʔεςΟϯά ϕΠδΞϯωοτϫʔΫɾϚϧίϑ࿈࠯ϞϯςΧϧϩ๏
  5. ϕΠζ౷ܭֶͱ͸ʁ 1ZUIPO#FHJOOFSTެࣜTMBDLɹCFHJMBC@CBZFTɹ ࢀߟ ֬཰ɾ౷ܭ ߴߍେֶ ϕΠζ౷ܭ "*σʔλαΠΤϯε ස౓ओٛ ओ؍ओٛ ֬཰͸ܾఆࡁΈ

    શମΛ෦෼͔Βਪଌ ʮநग़ʯ ֬཰͸ະܾఆ શମͰ͸ͳ͘ະདྷ ʮΞοϓσʔτʯ શମ฼ूஂ σʔλ͸੩త ෦෼ඪຊ ؍ଌऀΑΓશମ աڈʹࣄલ֬཰ σʔλ͸ಈతʹ؍ଌ ະདྷʹࣄޙ֬཰ શମΑΓ؍ଌऀ ౷ܭͷߟ͑ํΛࠜఈ͔Β ม͑Δඞཁ͋Γ
  6. ϕΠζͷఆཧ P(B|A) = P(A|B)P(B) P(A) ࣄޙ֬཰ ໬౓ ࣄલ֬཰ ϕΠζ౷ܭ͸ɺσʔλ͕ෆे෼Ͱ΋ɺ z͋Δࣄଶ͕ൃੜ͢Δ֬཰zʢʹࣄલ֬཰ʣΛ࠷ॳʹઃఆͨ͠ޙɺ

    ͞ΒͳΔ৘ใ͕ಘΒΕΔ౓ʹz͋Δࣄଶ͕ൃੜ͢Δ֬཰zʢʹࣄޙ֬཰ʣΛߋ৽͍͖ͯ͠ɺ ຊདྷى͜ΔͰ͋Ζ͏ࣄ৅ͷ֬཰ʢओ؍֬཰ʣΛಋ͖ग़͠·͢ɻ
  7. ౳෼ͷϕΠζ౷ܭֶ ౳෼ͷՖՇΛ ಡΜͩ͜ͱ͕ͳ͍ ຖिϚΨδϯͰ ಡΜͰΔ ίϛοΫεશרॴ࣋ Կ౓΋ϧʔϓ ࣄલ֬཰ ΄΅θϩʹ͍ۙ ҰՖγϣʔτ

    ೋ೫େ͖ΊϦϘϯ ࡾۢϔουϑΥϯ ࢛༿͏͗͞ϦϘϯ ޒ݄੕ΞΫη ҰՖෲࠇˍϚΠϖʔε ೋ೫πϯσϨˍՈ଒ࢥ͍ ࡾۢେਓ͍͠ˍরΕ԰ ࢛༿ݩؾˍԕྀ͕ͪ ޒ݄ਅ໘໨ ໬౓ ΄΅θϩʹ͍ۙ ߴ͍ ඇৗʹߴ͍ ࣄޙ֬཰ θϩʹ͍ۙ͸ͣ ʢ·͙Ε΋͋Δʣ ͓ͦΒ͘ߴ͍ ͓ͦΒ͘ඇৗʹߴ͍ ʢޱบ΍࢓૲Ͱ༧ଌՄೳʣ
  8. ਖ਼ղ ʆŋТŋ ʎ ݱࡏ ະདྷ աڈ    

               ˞͜Ε͸ͪΐͬͱࣗ৴ͳ͍
  9. ౳෼ͷϕΠζ౷ܭֶ ౳෼ͷՖՇΛ ಡΜͩ͜ͱ͕ͳ͍ ຖिϚΨδϯͰ ಡΜͰΔ ίϛοΫεશרॴ࣋ Կ౓΋ϧʔϓ ࣄલ֬཰ ΄΅θϩʹ͍ۙ ҰՖγϣʔτ

    ೋ೫େ͖ΊϦϘϯ ࡾۢϔουϑΥϯ ࢛༿͏͗͞ϦϘϯ ޒ݄੕ΞΫη ҰՖෲࠇˍϚΠϖʔε ೋ೫πϯσϨˍՈ଒ࢥ͍ ࡾۢେਓ͍͠ˍরΕ԰ ࢛༿ݩؾˍԕྀ͕ͪ ޒ݄ਅ໘໨ ໬౓ ΄΅θϩʹ͍ۙ ߴ͍ ඇৗʹߴ͍ ࣄޙ֬཰ θϩʹ͍ۙ͸ͣ ʢ·͙Ε΋͋Δʣ ͓ͦΒ͘ߴ͍ ͓ͦΒ͘ඇৗʹߴ͍ ʢޱบ΍࢓૲Ͱ༧ଌՄೳʣ ࠓ౓͸ೋ೫͚ͩ౰ͯͯΈ·͠ΐ͏ ʆŋТŋ
  10. 8JLJQFEJBΑΓൈਮ  Ψ΢ε෼෍ ਖ਼ن෼෍ ͱ͸ʁ Ψ΢ε෼෍ ਖ਼ن෼෍ ͱ͸ɺฏۉ஋ͷ෇ۙʹूੵ͢ΔΑ͏ͳσʔλͷ෼෍Λදͨ͠࿈ ଓతͳม਺ʹؔ͢Δ֬཰෼෍Ͱ͋Δɻ த৺ۃݶఆཧʹΑΓɺಠཱͳଟ਺ͷҼࢠͷ࿨ͱͯ͠ද͞ΕΔ֬཰ม਺͸ਖ਼ن෼෍ʹै͏ɻ

    ͜ͷ͜ͱʹΑΓਖ਼ن෼෍͸౷ܭֶ΍ࣗવՊֶɺࣾձՊֶͷ༷ʑͳ৔໘Ͱෳࡶͳݱ৅Λ؆୯ʹද͢Ϟσϧͱͯ͠༻͍ΒΕ͍ͯΔɻ ͨͱ͑͹࣮ݧʹ͓͚Δଌఆͷޡࠩ͸ਖ਼ن෼෍ʹैͬͯ෼෍͢ΔͱԾఆ͞Εɺෆ͔֬͞ͷධՁ͕ܭࢉ͞Ε͍ͯΔɻ ֬཰ີ౓ؔ਺ ඪ४ਖ਼ن෼෍͕࣋ͭ֬཰ີ౓ؔ਺ͷάϥϑ
  11. ʲ࠶ܝʳϕΠζ౷ܭֶͱ͸ʁ ֬཰ɾ౷ܭ ߴߍେֶ ϕΠζ౷ܭ "*σʔλαΠΤϯε ස౓ओٛ ओ؍ओٛ ֬཰͸ܾఆࡁΈ શମΛ෦෼͔Βਪଌ ʮநग़ʯ

    ֬཰͸ະܾఆ શମͰ͸ͳ͘ະདྷ ʮΞοϓσʔτʯ શମ฼ूஂ σʔλ͸੩త ෦෼ඪຊ ؍ଌऀΑΓશମ աڈʹࣄલ֬཰ σʔλ͸ಈతʹ؍ଌ ະདྷʹࣄޙ֬཰ શମΑΓ؍ଌऀ ฏۉ஋ͱू߹͕ऩଋ͞ΕΔ·Ͱ ଟ͘ͷσʔλ͕ඞཁ গͳ͍σʔλͰ΋ ਖ਼ن෼෍౳ʹԾఆͯ͠ஔ͖׵͑Δ͜ͱͰ ϞσϧΛߋ৽͍ͯ͘͜͠ͱ͕Մೳ
  12. ख़࿅ͷݚڀऀͷʮצͱܦݧʯΛ୭Ͱ΋؆୯ʹ࠶ݱʙՊֶٕज़ৼڵػߏ IUUQTXXXKTUHPKQQSBOOPVODFJOEFYIUNM ࢀߟ63- ࣄྫᶃ ख़࿅ͷݚڀऀͷʮצͱܦݧʯΛ୭Ͱ΋؆୯ʹ࠶ݱ ໨తʢ΍Γ͍ͨ͜ͱʣ ୯݁থߏ଄ղੳͷߴਫ਼౓Խͱߴޮ཰Խ ՝୊ʢϘτϧωοΫʣ ద੾ͳࢼྉͷબผͱܭଌ৚݅ͷઃఆΛߦ͏ඞཁ͕͋Δ ˠɹݚڀऀͷʮצͱܦݧʯʹཔ͍ͬͯͨ

    ˠɹਖ਼֬ͳࢉग़ʹ͸਺͔͔࣌ؒΔ ະܦݧͷ ݚڀऀ ख़࿅ͷ ݚڀऀ ຊݚڀάϧʔϓͷ։ൃͨ͠γεςϜ ࣄલ֬཰ θϩ צͱܦݧ ਺෼ͷ༧උܭଌͰಘΒΕΔ਺ඦఔ౓ͷσʔλ͔Β༧ଌ ໬౓ θϩ ߴ͍ ख़࿅ͷܦݧऀΑΓ΋͞Βʹߴ͍ ࣄޙ֬཰ θϩ ߴ͍ ख़࿅ͷܦݧऀΑΓ΋͞Βʹߴ͍
  13. 1Z$PO+1ػցֶशʹ͓͚ΔϋΠύʔύϥϝʔλ࠷దԽͷཧ࿦ͱ࣮ફ IUUQTTQFBLFSEFDLDPNONBTBIJSPKJKJFYVFYJOJPLFSVIBJIBIBSBNFUB[VJTIJIVBGBMTFMJMVOUPTIJKJBO ࢀߟ63- ࣄྫᶄ "VUP.-ͷϋΠύʔύϥϝʔλ࠷దԽ ໨తʢ΍Γ͍ͨ͜ͱʣ ϋΠύʔύϥϝʔλ࠷దԽ ՝୊ʢϘτϧωοΫʣ ద੾ͳύϥϝʔλΛઃఆ͢Δਓࡐ͕ݶΒΕΔ ύϥϝʔλΛߟ͑Δͷʹ΋͕͔͔࣌ؒΔ

    ˠɹϕΠζ౷ܭֶ͕༷ʑͳΞϧΰϦζϜͰ࢖༻͞Ε͍ͯΔ ී௨ͷ ΤϯδχΞ ख़࿅ͷ ΤϯδχΞ "VUP.- ࣄલ֬཰ ϥϯμϜ ஌ࣝɾܦݧ (1&*ɾ51&ɾ#0)#ɾ.5#0ͳͲ֤छΞϧΰϦζϜ ໬౓ ௿͍ ߴ͍ ߴ͍ʢख़࿅ͷΤϯδχΞΑΓߴ͍͔Ͳ͏͔͸ෆ໌ʣ ࣄޙ֬཰ ௿͍ ߴ͍ ߴ͍ʢख़࿅ͷΤϯδχΞΑΓߴ͍͔Ͳ͏͔͸ෆ໌ʣ
  14. ϕΠζ౷ܭͷద༻෼໺ ֬཰ɾ౷ܭɾਂ૚ֶश ϕΠζҎ֎ͷख๏ ϕΠζ౷ܭ ਺ཧϞσϧ Ϋϥ΢υͰͷֶश ը૾ೝࣝɾࣗવݴޠॲཧ౳ʣ (16౳ͷϦονͳϦιʔε͕࢖༻Ͱ͖ΔͳΒ ਂ૚ֶशͰߴ͍ਫ਼౓Λग़͢͜ͱ͕Մೳ ਺ཧϞσϧͷͨΊઆ໌͠΍͍͢

    ਫ਼౓Ͱ͸ྼΔ Τοδ"* ΤοδΞφϦςΟΫε ਂ૚ֶशʢਪ࿦ʣ͕࢖͑ΔͳΒߴ͍ਫ਼౓ ˞J1IPOFͳͲిྗʹ༨༟͕͋ΔσόΠεͰ࢖༻Մೳ গͳ͍Ϧιʔεɾগͳ͍σʔλͰ΋ ͋Δఔ౓ੑೳΛग़͢͜ͱ͕Մೳ "VUP.- ʢػցֶशͷࣗಈԽʣ গͳ͍&QPDIͰֶशͰ͖Δ͜ͱΛ ໨ࢦ͢ͷʹ͸গʑෆ޲͖ (1&*ɾ51&ɾ#0)#ɾ.5#0ͳͲ ༷ʑͳܗͰ༻͍ΒΕ͍ͯΔ গͳ͍ࢼߦճ਺Λޮ཰తʹར༻Մೳ ਂ૚ֶशΑΓઆ໌͠΍͍͢
  15. ਖ਼ن෼෍͸୯७ͳϞσϧ NJO  IUUQTXXXVEFNZDPNDPVSTFDPNQVUFSWJTJPOMFBSOMFDUVSF ਖ਼ن෼෍ͷ̏ͭͷ֦ு NJO  IUUQTXXXVEFNZDPNDPVSTFDPNQVUFSWJTJPOMFBSOMFDUVSF જࡏม਺ͱ&.ΞϧΰϦζϜ NJO

     IUUQTXXXVEFNZDPNDPVSTFDPNQVUFSWJTJPOMFBSOMFDUVSF ࠞ߹ਖ਼ن෼෍ NJO  IUUQTXXXVEFNZDPNDPVSTFDPNQVUFSWJTJPOMFBSOMFDUVSF 6EFNZϕΠζਪఆͱάϥϑΟΧϧϞσϧ ਖ਼ن෼෍͸୯७ͳϞσϧʙࠞ߹ਖ਼ن෼෍ 6EFNZϕΠζਪఆͱάϥϑΟΧϧϞσϧ
  16. 1ZUIPO#FHJOOFSTԭೄɾऩࢧ ߲໨ ֹۚ ิ଍ ऩೖ =  Ͱ༗ྉԽ͞Ε͔ͯΒͷੵཱۚʢճ෼  ߹ಉษڧձͰ͸ܦඅΛ1Z%BUB0LJOBXB͔Βෛ୲ͯ͠௖

    ͍ͨͨΊɺ1Z#FHJͰ͸௃ऩ͍ͯ͠·ͤΜɻ ࢧग़ = ·ͬͨ͘࢖͍ͬͯͳ͍ঢ়ଶͰ͢ɻ ࣍ճ1ZUIPO#FHJOOFSTͰొஃऀΛݝ֎͔ΒݺͿࡍ ͷ౉ߤඅʹ౰ͯΔ༧ఆͰ͢ɻ ऩࢧ =  ͪΐͬͱཷ·͖ͬͯ·ͨ͠Ͷʔɻ