Discovering Topical Interactions in Text-based Cascades using Hidden Markov Hawkes Processes
|
|
- Junior Carson
- 5 years ago
- Views:
Transcription
1 Discovering Topical Interactions in Text-based Cascades using Hidden Markov Hawkes Processes Jayesh Choudhari, Anirban Dasgupta IIT Gandhinagar, India {choudhari.jayesh, Indrajit Bhattacharya TCS Research, India Srikanta Bedathur IIT Delhi, India arxiv: v1 [cs.lg] 12 Sep 2018 Abstract Social media conversations unfold based on complex interactions between users, topics and time. While recent models have been proposed to capture network strengths between users, users topical preferences and temporal patterns between posting and response times, interaction patterns between topics has not been studied. We argue that social media conversations naturally involve interacting rather than independent topics. Modeling such topical interaction patterns can additionally help in inference of latent variables in the data such as diffusion parents and topics of events. We propose the Hidden Markov Hawkes Process (HMHP) that incorporates topical Markov Chains within Hawkes processes to jointly model topical interactions along with useruser and user-topic patterns. We propose a Gibbs sampling algorithm for HMHP that jointly infers the network strengths, diffusion paths, the topics of the posts as well as the topictopic interactions. We show using experiments on real and semisynthetic data that HMHP is able to generalize better and recover the network strengths, topics and diffusion paths more accurately that state-of-the-art baselines. More interestingly, HMHP finds insightful interactions between topics in real tweets which no existing model is able to do. This can potentially lead to actionable insights enabling, e.g., user targeting for influence maximization. INTRODUCTION A popular area of recent research has been the study of information diffusion cascades, where information spreads over a social network when a parent event from one infected node influences a child event at neighboring node [5], [11], [19], [6], [10]. The action of propagating information between two neighboring nodes depends on various factors, such as the strength of influence between the nodes, the topical content of the parent event and the extent of interest of the child node towards that topic. Explosion of social media data has made it possible to analyze and evaluate different models that seek to explain such information cascades. However, many relevant variables such as the network influence strengths, the identity of influencing or parent event for any event, and the actual topics are typically unobserved for most social network data. Therefore, these need to be inferred before any analysis of the diffusion patterns can be performed. A recent body of work in this area has proposed increasingly sophisticated models for information cascades. To account for the temporal burstiness of events, Hawkes processes, which are self-exciting point processes, have been proposed to model their time stamps [16]. Also, influences travel more quickly over stronger social ties. This is modeled by the Network Hawkes process [13] by incorporating connection strengths between users into the intensity function of the Hawkes process. When events additionally have associated textual data, parent and child events in a cascade are typically similar in their topical content. The Hawkes Topic Model [11] captures this by combining the Network Hawkes process with a LDAlike generative process for the textual content where the topic mixture for a child event is close to that for a parent event. Given partially observed data, inference algorithms for both the Network Hawkes and the Hawkes Topic models recover latent parent identities for events and the network strengths. In addition, the Hawkes Topic model enables recovery of the latent topics of the events. The main strength of this inference algorithm is its collective nature - the recovery of topics, parents and network strengths reinforce each other. Despite its strengths, the Hawkes Topic Model fails to capture one important aspect of the richness of information cascades. Typically, there are sequential patterns in the textual content of different events within a cascade. In terms of topical representation of the events, topics of parent and child events display repeating patterns. Consider the following pair of parent-child tweets detected by our model: THOUGH SEEM- INGL DELICATE, MONARCH #BUTTERFLIES ARE REMARK- ABL RESILIENT & THEIR DECLINE CAN BE REVERSED #PESTICIDES and EPA SET TO REVEAL TOUGH NEW SULFUR EMISSIONS RULE #CLIMATECHANGE #SUSTAIN- ABILIT #POLLUTION. It is not obvious how to place these two tweets in the same topic. Our model places them in different topics (say pesticides and butterflies and climate change and pollution) and detects that this topic pair appears frequently in other parent-child tweets in our data. Adding to this their posting time information, it has sufficient evidence to detect a parent-child relation between this pair. However, just like topic distributions, topic interaction patterns are also typically latent, and need to be inferred. Interestingly, inferring topic interaction patterns in turn benefits from more accurate parent and topic assignment to events. This calls for joint inference of topic interaction patterns and the other latent variables. We propose a generative model for textual cascades that captures topical interactions in addition to temporal burstiness, network strengths, and user topic preferences. Temporal and network burstiness is captured using a Hawkes Process over the social network [13]. Topical interactions are modeled by combing this process with a Markov Chain over event topics
2 along a path of the diffusion cascade. Such topical Markov Chains have appeared in the literature [8], [9], [1], [2], [4], but in the very different context of modeling word sequences in a document. Our model effectively integrates topical Markov Chains with the Network Hawkes process and we call it the Hidden Markov Hawkes Process (HMHP). Making use of conjugate priors, we derive a simple and efficient collapsed Gibbs Sampling algorithm for inference using our model. This algorithm is significantly more scalable than the variational algorithm for the Hawkes Topic Model, and allows us to analyze large collections of real tweets. We validate using experiments that HMHP fits information cascades better compared to models that do not consider topical interactions. With the aid of semi-synthetic data, we show that modeling topical interactions indeed leads to more accurate identification of event parents and topics. More importantly, HMHP is able to identify interesting topical interactions in real tweets, which is beyond the capability of any existing information diffusion model. One of the underlying aims of modeling information cascades is to develop accurate models that can predict the growth and virality of cascades on specific topics started by specific users such models can then be used in practice to select a set of users to incentivize to start a conversation about the topics desired. Our modeling of the topic-topic interaction can broaden this approach by providing access to users who generate content about related set of topics, from which, ultimately, the conversations tend to drift to the topics of interest. Mechanisms that can estimate such possible related candidate topics to target could thus be practically very useful, in addition to providing unique insights about how conversations evolve on social networks. MODEL We consider a set of nodes V = {1,..., n}, representing content producers, and a set of directed edges E among them, representing the underlying graph using which information propagates. For every edge (u, v), W uv denotes the weight of the edge, capturing the extent of influence between content producers or users u and v. While we assume the underlying graph structure to be observed, the actual weights W uv are not. This captures the intuition that while the graph structure may be known, the actual measure of influence that one user has on another is not. For each event e, representing e.g. a social media post, we observe its posting time t e, the user c e who creates the post and the document d e that is associated with the post. The posting time follows a Hawkes Process [17] incorporating user-user edge weights [13]. Event e could be spontaneous, or a diffusion event, meaning that it is triggered by some recent event, which we call its parent z e, created by one of the followees of user c e. The document d e for e is drawn using a topic-model over a vocabulary of size W. For spontaneous events, the topic choice η e depends on the topic preference of the user c e. We deviate from existing literature in modeling the topic of a diffusion event. Instead of being identical or close to the topic of the triggering event [11], the diffusion event may be on a related topic. We model this transition between related topics using a Markov Chain over topics, involving a topic transition matrix T. This enables us to capture repeating patterns in topical transitions between parent and child events. The generated event e is then consumed by each of the followers of the user c e, thereby triggering them in turn to create multiple events and producing an information cascade. Since the topical sequence is hidden and observed only indirectly through the content of the events, we call our model the Hidden Markov Hawkes Process (HMHP). We next describe the details of our generative model. This can be broken up into two distinct stages - generating the cascade events, and generating the event documents. Generating Cascade Events We define an event e by a tuple (t e, c e, z e, d e ) where t e indicates the time for event e, c e the id of the creator node and z e the unique parent event that triggered the creation of event e, set to 0 if the event e is spontaneous, and d e is the textual document. The generative model for HMHP consists of two phases. We first generate (t e, c e, z e ) for all events using a multivariate Hawkes process (MHP) following existing models [11], [13], and then use a hidden Markov process based topic model to generate all the associated documents d e. A multivariate Hawkes process models the fact that the users can mutually excite each other to produce events. A Hawkes process can also be represented as a superposition of Poisson processes [17]. For each node v, we define λ v (t), a rate at time t, as a superposition of the base intensity µ v for user v, and the impulse responses for each event e n that has happened at time t n at a followee node c n. λ v (t) = µ v (t) + H t n=1 h cn,v(t t n ) (1) where h cn,v(t t n ) is the impulse response of user c n on the user v, and H t is the history of events upto time t. The impulse response can be decomposed as the product of the influence W uv of user u on v, and a time-kernel as follows: h u,v ( t) = W u,v f( t) (2) We note that while there has been a number of recent works on modeling the time-kernel effectively [13], [5], [11], we use a simple exponential kernel f( t) = exp( t), as this is not the main thrust of our work. Following [17], we generate the events using a level-wise generation process. Level 0, denoted as Π 0, contains all the spontaneous events, generated using the base rates of the users. The events Π l at level l, are generated according to the following non-homogenous Poisson process Π l P oisson h cn, (t t n ) (3) (t n,c n,z n) Π l 1 The above process can be simulated by the following algorithm for each event e n = (t n, c n, z n ) Π l 1, and for
3 each neighbor v of cn, we draw timestamps according to the non-homogeneous Poisson process P oisson(hcn,v (t tn )); we generate an event e for each of these timestamps, and set ze = en (parent) and ce = v (producer node). V Dirichlet hyperparameter for user-topic preference sha1_base64="furcnptcwyyrxvbeloam9fsohsc=">aaacbnicbvdlssnafj3uv42vqesrbkvbvule0juu3lis0bc0iuwmk3bozbjmjkijwbnxv9y4umst3+dov3hszqgtf45nhmv99wtpixkzdvfrm1tfwnzq75t7uzu7r9h0d9mwqckx5owckgazkeuu56iipghqkgka4gqtt21ifpbahack7apsl0zjtiokkdkub532zabjcyus1h/ursjncgi5d3cnxa+1bbb9rzgknaq0abvdxzryw0tnmwek8yqlcphtpwxi6eozqqw3uysfoepgporhhzfrhr5/iwcnjutwigr+nef5+zvirzfsrspo0uxclkryf+0uaaiay+npm0u4xixkmovaksm4ehfqqrntmaug1v4gnsccsdhkmdsfzpnkv9c9ajt1y7i8b7zsqjjo4awfghdjgcrtbheiahsdgetydv/bmpbkvxrvxswitgdxmmfhtxucpgewzjq==</latexit> User-topic preference v sha1_base64="1rt/pkn6onhko1wqd1dwfyviqim=">aaab/hicbvdns8mwhe3n15xf1r29bifgabqi6ekgxjxocb+wlpkm6rawjivjb6xmf8wlb0w8+od4878x3xrqzqchj/d+p/lywprrpr3n26ptbg5t79r3g3v7b4dh9vfjx4lmtldggk5djeijhls01qzmkwlquniyccc3px+eakooi/6jwlfolgnmui22kwg56owcryhnzfv46ocfshtgtp+0sanejw5ewqnan7c8vejhlcneiavgrpnqv0bsu8zivofliqqit9gjazlkchklxbh5/dckbgmhtsha7hqf28ukfflpjozid1rq14p/uenmh3f+axlaajx8uh4oxblwdzbiyojfiz3bcejtvziz4giba2ftvmce7ql9dj/7ltom334arvua3qqintcaugauuqqfcgy7oaqxy8axewzv1zl179bhcrrmvttn8afw5w+jrjvi</latexit> Topic-topic vector sha1_base64="4krulav1mc1ov3+lzermtpj9h8=">aaab/hicbvblswmxgmz6rpw12qoxbe8lv0r9cqflx4r2ad0l5lnztvqpjkkyxl/stepcji1r/izx9jtt2dtg6eddpfrytpxq43nfztr6xubwdm2nvru3f3dohh33tmwujl0smvsdcgncqcbdqw0jg1qrxcng+th0tvt7j0rpkswdyvmscjqwnkegsun3eqsrbrnnurcmaiczqbuu2v5c0bv4lfksao0bm5x0esccajmjghre+l5qwqmpqzmishmsapahp0zgmlrwiex0w8/azegavgczs2smmnku/nwredznptnjkjnrzk8x/vgfmkuuwoclndbf48vcsmwgkljuamvueg5zbgrcinivee6qqnravui3bx/7ykuldthyv5d9fnts3vr01cajowtnwwrvogzvqav2aqq6ewst4c56cf+fd+vimrjnvtgp8gfp5a4aqlu8=</latexit> K Dirichlet hyperparameter for topic-topic vectors sha1_base64="jckw2oyx2xdrrp5bgtbddaxqss=">aaab6hicbvbns8naej3ur1q/qh69lbbbu0le0jmuvahewraf0iay2u7atztn2n0ijfqxepggifd/kjf/jds2b219mpb4b4azeueiudau++0u1t3nrek26wd3b39g/lhuuvhqwlzlgivsegggwx2dtccowkcmkucgwh49uz335cpxksh8wkqt+iq8ldzqixuuo+x664vxcoskq8nfqgr71f/uonpzgka0tvouu5ybgz6gynamclnqpxosymr1i11jji9r+nj90ss6smibhrgxjq+bq74mmrlpposb2rtsm9li3e//zuqkjr/2myyq1knliuzgkmiy+5omuejmxmqsyhs3txi2oooy7mp2rc85zdxseui6rlvr3fzqd3kcrthbe7hhdy4ghrcqr2awadhgv7hzxl0xpx352prwndymwp4a+fzb6drjms=</latexit> Tk sha1_base64="9nkp7f/eq54eygm2wt3iozc3v=">aaaca3icbvdlssnafj34rpuvdaebwsk4kokiupkcg5cv+oimhmlk0g6dziszivbcwi2/4safim79cxf+jzm2c229mmzhnhu5554wzvrpx/m2vlbx1jc2a1v17z3dvx374lcnrcx6wlbhbyesbfgoelqqhkzpjkgjgskh05us73/qksignf0ncv+gkacxhqjbajapvzcwsi1tcyxewns4x3imcsrhdafpzaoua7ccdvbvo7c/vejglcfc4augrpoqv0csu0xi0xdyxrjez6gerkayffclj/pbijgmweigatphtdwxv6eyfgispums3spfrws/e8bzjq+9npk00wtjuel4oxblwazciyojfizqqeis2q8qjxgemftqubenzfk5db76lpok33/rlruqniqietcarogquuqavcgtboagwewtn4bw/wk/vivvsf89vq5o5an/k+vwbomwgg==</latexit> sha1_base64="s+f2w9slbjudzkvew2ohcs8o0/0=">aaab+3icbvdlssnafj34rpuv69lnbfcluqexunbjcsk9gfnkjpjttt0mgkze7ge/iobf4q49ufc+tdo2iy09cawh3puzc6ciovmacf5ttbwnza3tms79d29/nd+6jru0kmkxrpwhm5cigczgr0ndmcbqkeegcc+sh0tvt7jyavs8sdnqxgx2qswmqo0ua2q0vshiozrg5ci8atqr3xrazhx4lbgvaaiknzh95ujzwiqmnki1nb1uu3nrgpgorr1l1oqejolxgakkgmys/n2qt8zpqqr4k0r2g8v39v5crwztwzgrm9uctekf7ndtmdxfs5e2mmqddfq1hgsu5wwqqomqsq+cwqqiuzwtgdeemonnxvtqnu8pdxse+i5tot9/6y2b6p6qihe3skzpglrlab3aeo6ikkntazekvvvmg9wo/wx2j0zap2jtefwj8/veeu3g==</latexit> sha1_base64="nv4rdzpte2wl5a6vknquqvmdjc=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0m86euoepekfe0htkfstpn26wtdjdccf0jxjwo4tvf5m1/47bnqvsfddzem2fmxpgkro3nftultfwnza3ydmvnd2//wd08aukkuwyblbgj6oruo+asm4bgz1ui1dge1wfdpz20+one/ko5mkgmr0khnegtvwerhgv+9wvzo3b1klfkgqukdrd796g4rlmurdbnw663upcxkqdgccp5vepjglbeyh2lvu0hh1km9pnzizqwxilchb0pc5+nsip7hwkzi0nte1i73szct/vg5moqsg5zlndeq2wbrlgpiezp4ma66qgtgxhdlf7a2ejaiiznh0kjef/nlvdk6qplezb/3qvw7io4ynmapnimpl1chw2haexgm4rle4c0rzovz7nwswktomxmmf+b8/gc9o1z</latexit> sha1_base64="nkvkz+558ccv1p97d30fpbaiyz4=">aaab8nicbvbnswmxem3wr1q/qh69bivgqws96euoepekfewhtmustbntadzzklmhlv0zxjwo4tvf481/9ruqvsfddzem2fmxprkgqb6+0tr6xuvxeruzs7u0fva+p2lznhvew01kbbkqtl0lxfgiqvjsatpni8k40vpn5nudurndqaspdxi6vciwjiktek9hzkn/iq8xcas1uidz4fxif6sgcjtd6ld/ofmwcavmumt7pkkhykkbwssfvvqz5sllzrkpucvtbgn8vnju3zmlagotxglam/v3xm5taydjjhrtcim7li3e//zehnev0euvjobv2yxkm4kbo1n/+obmjybndhcmrhuvsxg1faglqwkc8fffnmvtc/qpqn796twucvikkmtdirok8uuqpdoizqi0ekav6m0d78v79z4wrswvmdlgf+b9/gdzcjbj</latexit> sha1_base64="nkvkz+558ccv1p97d30fpbaiyz4=">aaab8nicbvbnswmxem3wr1q/qh69bivgqws96euoepekfewhtmustbntadzzklmhlv0zxjwo4tvf481/9ruqvsfddzem2fmxprkgqb6+0tr6xuvxeruzs7u0fva+p2lznhvew01kbbkqtl0lxfgiqvjsatpni8k40vpn5nudurndqaspdxi6vciwjiktek9hzkn/iq8xcas1uidz4fxif6sgcjtd6ld/ofmwcavmumt7pkkhykkbwssfvvqz5sllzrkpucvtbgn8vnju3zmlagotxglam/v3xm5taydjjhrtcim7li3e//zehnev0euvjobv2yxkm4kbo1n/+obmjybndhcmrhuvsxg1faglqwkc8fffnmvtc/qpqn796twucvikkmtdirok8uuqpdoizqi0ekav6m0d78v79z4wrswvmdlgf+b9/gdzcjbj</latexit> sha1_base64="xt2aaugnu8p84wwaal0etl0aoza=">aaab6nicbva9swnbej2lxzf+rs1tfongffzstazweleewpjef2c8msvb1jd08ir36cjuitv4io/+nm+qktxww8hhvhpl5sqfszr+e6w193nrfj2zwd3b/+genjunkmmobz4ihpdczlbkrs2rlaso6lgfocsh8px9cx/fejtrkie7ctfigzdjslbmxxspfb9frvg63qoskr8gtsgqlnf/eonep7fqcyxzjiut1mb5exbwsvok73mmr4ma2x66himzogn586jwdogzao0a6ujxp190tommmceg62zhztmbif953cxgv0euvjpzvhyxkmoksqmz/u0gqio3cuii41q4wwkfmc24deluxaj+8surph1r92ndv6o1xm0rrxlo4btowdlamannkefhibwdk/w5knvxxv3phatja+o/8d5/apeljzu=</latexit> sha1_base64="srrkghm85upp5nldeepu+/bunq=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0n0omecf09s0x5ag8pmo2mxbjzhdyou0j/gxmixv1f3vw3btsctpxbwoo9gwbmhang2njet1naw9/3cpvv3z29/p3mojlk4yxbdjepgotkg1ci6xabgr2ekv0jgu2a7hnzo//rk80q+mkmkquyhkkecuwolb+xf9t2qv/pmikvel0gvcjt67ldvklasrmmofp3fs81qu6v4uzgtnllnkaujekqu5zkgqmo8vmpu3jmlqgjemvlgjjxf0/knnz6eoe2m6zmpje9mfif181mdb3kxkazqckwi6jmejoq2d9kwbuyiyawuka4vzwwevwugztoxbgl7+8sloxnd+r+fdetx5xxfggezifc/dhcupwcw1oaomhpmmrvdncexhen9fa8kpzo7hd5zph/qtjzc=</latexit> sha1_base64="xlqszcdubydm9adntjnpx3fdk8=">aaab6nicbvbns8naej3ur1q/qh69lbbbu0m86lhgxznutb/qhrlzbtqlm03nrrk6e/w4kerr/4ib/4bt2ko2vpg4pheddpzgkqkg6777zq2nre2d8q7lb39g8oj6vfj28spzrzfhnrbkanl0lxfgquvjtotqna8k4wuv34nsnxrstqcwcj9ym6uiiujkkvhqcdb1ctuxu3b1knxkfquka5qh71hzfli66qswpmz3mt9doqutdj55v+anhc2soem9srsnu/cw/du4urdikaxtkss5+nsio5exsyiwnrhfsvn1fuj/xi/f8mbpheps5iotf4wpjbitxd9kkdrnkgewukafvzwwmdwuou2nkpwvl9ej+2ruufwvqe31rgv4ijdgzzdjxhwdq24gya0gmeinuev3hzpvdjvzseytequm6fwb87ndwsaja=</latexit> 1 C Topic assignment 2 sha1_base64="631x844qjobqmznedfhqodsquq=">aaab7xicbvdlsgnbeoz1gemr6thlba8hd0g6ekcxjxgma9ilja7mu3gzm4sm71ccpkhlx4u8er/epnvncr70mschqkqm+6ukjxcou9/e2vrg5tb24wd4u7e/sfh6ei4axvmgg8wlbvpr9rykrrvoedj26nhnikkb0wj25nfeulgcq0ecjzymkedjwlbkdqp2evie9veqexx/dnikglyuoc9v7pq9vxleu4qiaptz3atzgcuiocst4tdjplu8pgdma7jiqacbto5tdoyblt+itwxpvcmld/t0xou04ivxnqnfol72z+j/xytc+didcprlyxral4kws1gt2ouklwxnkssougefujwxidwxoaiq6eilll1djs1oj/epwf1mu3erxfoauzuacaricgtxbhrra4bge4rxepo29eo/ex6j1zctntuapvm8fmoko2q==</latexit> sha1_base64="jtid92jlc7gbxl3t8d0od4arfte=">aaab7xicbvdlsgnbeoynrxhfu9ebopgkeykoccjepewtwgwclszdzmzuzzpqkieqfvhhqxkv/482/czlsqrmlgoqqbrq7olqki77/7rxw1jc2t4rbpz3dvf2d8ufr0+rmmn5gwmrtjqjluijeqigst1pdarjj3opgtzo/9csnfvo94djluihssscuxrss8ur9ojeuejx/tnikglyuoec9v75q9vxleu4qiaptz3atzgcuiocst4tdtplu8pgdma7jiqacbto5tdoyzlt+itwxpvcmld/t0xou04ivxnqnfol72z+j/xytc+didcprlyxral4kws1gt2ouklwxnkssougefujwxidwxoaiq5eilll1dj86ia+nxg/rjsu8njkmijnmi5bhafnbidojsawsm8wyu8edp78d69j0vrwctnjuepvm8fl16o2a==</latexit> sha1_base64="97chz67ywynvdv5k76jmcxxumba=">aaab7xicbva9swnbej2lxzf+rs1tfongfe600djgyurzackr9jbzcvr9nap3t0hhpwhgwtfbp0/dv4bn8kvmvhg4pheddpzolrw33/2yusrw9sbhw3szu7e/sh5cojplgzzthgsijdjqhbwsu2llcc26lgmkqcw9hozua3nlabrusdhacjnqgecwztu5qdths3mwvxpgr/hxklqq5qucoeq/81e0rliuolrpume7gpzacug05ezgtdtodkwujoscoo5imamlj/nopoxnkn8rku5kwznxfexoagdnoitezuds0y95m/m/rzda+didcpplfyral4kwqq8jsddlngpkv0co09zdstiqasqsc6jkqgiwx14lztq4fede79su8vjkmijnmi5bhafnbifojsawsm8wyu8ecp78d69j0vrwctnjuepvm8fnfao5g==</latexit> sha1_base64="zjcil+fnvv32rvhpf2yzkz167q=">aaab7xicbvdlsgnbeoynrxhfu9ebopgkeykomeaf08swtwgwclspdczmzuzzmwkieqfvhhqxkv/482/czlsqrmlgoqqbrq7olrw33/2yusrw9sbhw3szu7e/sh5cojplgzzthgsijdjqhbwsu2llcc26lgmkqcw9hozua3nlabrusdhacjnqgecwztu5qdths3mwvxpgr/hxklqq5qucoeq/81e0rliuolrpume7gpzacug05ezgtdtodkwujoscoo5imamlj/nopoxnkn8rku5kwznxfexoagdnoitezuds0y95m/m/rzda+didcpplfyral4kwqq8jsddlngpkv0co09zdstiqasqsc6jkqgiwx14lztq4fede79su8vjkmijnmi5bhafnbifojsawsm8wyu8ecp78d69j0vrwctnjuepvm8fn3qo5w==</latexit> sha1_base64="nv4rdzpte2wl5a6vknquqvmdjc=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0m86euoepekfe0htkfstpn26wtdjdccf0jxjwo4tvf5m1/47bnqvsfddzem2fmxpgkro3nftultfwnza3ydmvnd2//wd08aukkuwyblbgj6oruo+asm4bgz1ui1dge1wfdpz20+one/ko5mkgmr0khnegtvwerhgv+9wvzo3b1klfkgqukdrd796g4rlmurdbnw663upcxkqdgccp5vepjglbeyh2lvu0hh1km9pnzizqwxilchb0pc5+nsip7hwkzi0nte1i73szct/vg5moqsg5zlndeq2wbrlgpiezp4ma66qgtgxhdlf7a2ejaiiznh0kjef/nlvdk6qplezb/3qvw7io4ynmapnimpl1chw2haexgm4rle4c0rzovz7nwswktomxmmf+b8/gc9o1z</latexit> sha1_base64="nv4rdzpte2wl5a6vknquqvmdjc=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0m86euoepekfe0htkfstpn26wtdjdccf0jxjwo4tvf5m1/47bnqvsfddzem2fmxpgkro3nftultfwnza3ydmvnd2//wd08aukkuwyblbgj6oruo+asm4bgz1ui1dge1wfdpz20+one/ko5mkgmr0khnegtvwerhgv+9wvzo3b1klfkgqukdrd796g4rlmurdbnw663upcxkqdgccp5vepjglbeyh2lvu0hh1km9pnzizqwxilchb0pc5+nsip7hwkzi0nte1i73szct/vg5moqsg5zlndeq2wbrlgpiezp4ma66qgtgxhdlf7a2ejaiiznh0kjef/nlvdk6qplezb/3qvw7io4ynmapnimpl1chw2haexgm4rle4c0rzovz7nwswktomxmmf+b8/gc9o1z</latexit> sha1_base64="viatvdjju2nj8ubf07csrizaj8=">aaab9hicbvdlsgnbeoynrxhfu9ebopgkezgg16egbdpese8ifmw2ulvmmt24cxsic75di8efphqx3jzb5wke9degoaiqpvulj8rxgnb/rka+sbm1vf7dlo7t7+qfnwqkxivdjssljesunthjh2nrcc+wkemnoc2z7o5uz3x6jvdyohvqkqtekg4ghnfftjpfjy9c7mjjrgl7nk1fsqj0hwsvotiqqo+gvv3r9mkuhrpojqltxsrptzlrqzgros71uulzia6wa2heq1runj96ss6m0idble1fmszv3xmzdzwahl7pdkkeqmvvjv7ndvmdxlkzj5ju8qwi4jueb2twqkkzyuylsaguca5uzwwizwuazntytgll+8slq1qmnxnxu7ur/l4yjcczzcothwcxw4hq0gcejpmmrvflj68v6tz4wrqurnzmgp7a+fwb0xze/</latexit> sha1_base64="viatvdjju2nj8ubf07csrizaj8=">aaab9hicbvdlsgnbeoynrxhfu9ebopgkezgg16egbdpese8ifmw2ulvmmt24cxsic75di8efphqx3jzb5wke9degoaiqpvulj8rxgnb/rka+sbm1vf7dlo7t7+qfnwqkxivdjssljesunthjh2nrcc+wkemnoc2z7o5uz3x6jvdyohvqkqtekg4ghnfftjpfjy9c7mjjrgl7nk1fsqj0hwsvotiqqo+gvv3r9mkuhrpojqltxsrptzlrqzgros71uulzia6wa2heq1runj96ss6m0idble1fmszv3xmzdzwahl7pdkkeqmvvjv7ndvmdxlkzj5ju8qwi4jueb2twqkkzyuylsaguca5uzwwizwuazntytgll+8slq1qmnxnxu7ur/l4yjcczzcothwcxw4hq0gcejpmmrvflj68v6tz4wrqurnzmgp7a+fwb0xze/</latexit> sha1_base64="viatvdjju2nj8ubf07csrizaj8=">aaab9hicbvdlsgnbeoynrxhfu9ebopgkezgg16egbdpese8ifmw2ulvmmt24cxsic75di8efphqx3jzb5wke9degoaiqpvulj8rxgnb/rka+sbm1vf7dlo7t7+qfnwqkxivdjssljesunthjh2nrcc+wkemnoc2z7o5uz3x6jvdyohvqkqtekg4ghnfftjpfjy9c7mjjrgl7nk1fsqj0hwsvotiqqo+gvv3r9mkuhrpojqltxsrptzlrqzgros71uulzia6wa2heq1runj96ss6m0idble1fmszv3xmzdzwahl7pdkkeqmvvjv7ndvmdxlkzj5ju8qwi4jueb2twqkkzyuylsaguca5uzwwizwuazntytgll+8slq1qmnxnxu7ur/l4yjcczzcothwcxw4hq0gcejpmmrvflj68v6tz4wrqurnzmgp7a+fwb0xze/</latexit> sha1_base64="fz38oznm5eba3bojzrefco7c6d0=">aaab9hicbva9swnbej3zm8avqkxnhcswl0absajzvemb+qhmfezpis2ds7d/cc8cjvslfqxnf+e/czncokpbh7vztazl0we18z1v52193nre3ctnf3b//gshr03nrxqhg2wcxi1q6prselngw3atujqhqfalvh6gbmt8aoni/lg5kk6ed0ihmfm2qs5d8fgqbvkbkmghhbqexw3dnikvfyuoc9ad01e3fli1qgiao1h3ptyfuwu4ezgtdloncwujoscopzjgqp1sfvsunfulr/qxsiunmau/jziaat2jqtszutpuy95m/m/rpkz/5wdcjqlbyral+qkgjiazbeipk2rgtcyhthf7k2fdqigznqeidcfbfnmvnksvz dpfhubtoiml8oasanaldwgag0d4hld4c8boi/pufcxa15x85gt+wpn8axfokt0=</latexit> w Generating Event Documents w sha1_base64="gme7o1v+ghubfih0ri+oz4+6sq8=">aaab6hicbvbns8naej3ur1q/qh69lbbbu0le0jmuvhhswx5ag8pmo2nxbjzhd6ou0f/gxmixv1j3vw3btsctpxbwoo9gwbmbng2rjut1nw9/3cpul3z29/pyodhlr2nimgtxsjwnbqffxi03ajsjmopfegsb2mb2d++xgv5rg8n5me/gojq85o8zkjad+uejw3tnikvfyuoec9x75qzeiwrqhnexqrbuemxg/o8pwjnba6quae8rgdihdsywnupvz/napobpkgisxsiunmau/jziaat2jatszutpsy95m/m/rpia89jmuk9sgzitfsqiicnsazlgcpkre0sou9zestiiksqmzazkq/cwx14lruq51a9xmwldpphuqtoivz8oakanahdwgca4rneiu358f5cd6dj0vrwclnjuepnm8f44gm9w==</latexit> w w sha1_base64="gme7o1v+ghubfih0ri+oz4+6sq8=">aaab6hicbvbns8naej3ur1q/qh69lbbbu0le0jmuvhhswx5ag8pmo2nxbjzhd6ou0f/gxmixv1j3vw3btsctpxbwoo9gwbmbng2rjut1nw9/3cpul3z29/pyodhlr2nimgtxsjwnbqffxi03ajsjmopfegsb2mb2d++xgv5rg8n5me/gojq85o8zkjad+uejw3tnikvfyuoec9x75qzeiwrqhnexqrbuemxg/o8pwjnba6quae8rgdihdsywnupvz/napobpkgisxsiunmau/jziaat2jatszutpsy95m/m/rpia89jmuk9sgzitfsqiicnsazlgcpkre0sou9zestiiksqmzazkq/cwx14lruq51a9xmwldpphuqtoivz8oakanahdwgca4rneiu358f5cd6dj0vrwclnjuepnm8f44gm9w==</latexit> N1 sha1_base64="gme7o1v+ghubfih0ri+oz4+6sq8=">aaab6hicbvbns8naej3ur1q/qh69lbbbu0le0jmuvhhswx5ag8pmo2nxbjzhd6ou0f/gxmixv1j3vw3btsctpxbwoo9gwbmbng2rjut1nw9/3cpul3z29/pyodhlr2nimgtxsjwnbqffxi03ajsjmopfegsb2mb2d++xgv5rg8n5me/gojq85o8zkjad+uejw3tnikvfyuoec9x75qzeiwrqhnexqrbuemxg/o8pwjnba6quae8rgdihdsywnupvz/napobpkgisxsiunmau/jziaat2jatszutpsy95m/m/rpia89jmuk9sgzitfsqiicnsazlgcpkre0sou9zestiiksqmzazkq/cwx14lruq51a9xmwldpphuqtoivz8oakanahdwgca4rneiu358f5cd6dj0vrwclnjuepnm8f44gm9w==</latexit> sha1_base64="gme7o1v+ghubfih0ri+oz4+6sq8=">aaab6hicbvbns8naej3ur1q/qh69lbbbu0le0jmuvhhswx5ag8pmo2nxbjzhd6ou0f/gxmixv1j3vw3btsctpxbwoo9gwbmbng2rjut1nw9/3cpul3z29/pyodhlr2nimgtxsjwnbqffxi03ajsjmopfegsb2mb2d++xgv5rg8n5me/gojq85o8zkjad+uejw3tnikvfyuoec9x75qzeiwrqhnexqrbuemxg/o8pwjnba6quae8rgdihdsywnupvz/napobpkgisxsiunmau/jziaat2jatszutpsy95m/m/rpia89jmuk9sgzitfsqiicnsazlgcpkre0sou9zestiiksqmzazkq/cwx14lruq51a9xmwldpphuqtoivz8oakanahdwgca4rneiu358f5cd6dj0vrwclnjuepnm8f44gm9w==</latexit> sha1_base64="sqr3px+kxlongsmku4og89wlvm=">aaab6nicbvbns8naej3ur1q/qh69lbbbu0leqmecf0+lov2anptndtiu3wzc7koot/biwdfvpqlvplv3l5aoudgcd7m8zmcxlbtxhdb6ewsbm1vvpcle3thxwel9p2jpofcmwi0wsughvkljelufgddrsknace3m79zhmqzwp5akj+hedsr5yro2vhhqd60g54lbdbcg68xjsgrznqfmrp4xzgqe0tfcte56bgd+jynamcfbqpxotyiz0hd1lj1q+9ni1bm5smqqhlgyjq1zql8nmhppp0c2xlrm9ar3lz8z+uljrzxmy6t1kbky0vhkoijyfxvmuqkmrftsyht3n5k2jgqyoxnp2rd8fzfxiftq6rnvr17t1jv5heu4qzo4ri8qeed7qajlwawgmd4htdhoc/ou/oxbc04+cwp/ihz+qpsj2b</latexit> sha1_base64="vhflzed5lzk9eggf4bxnvcqoz9e=">aaab6nicbvbns8naej3ur1q/qh69lbbbu0leqccpepekfe0htkfstpt26wtdidccf0jxjwo4tvf5m1/47bnqvsfddzem2fmxpbidb1v53c2vrg5lzxu7szu7d/ud48apk41w3wsxj3qmo4vio3ksbkncszwkusn4oxjczv/3etrgxesrjwv2idpuibanope7vtcvv9yqowdzjv5okpcj0s9/9qxsyoukelqtndze/qzqlewyaelxmp4qtmdnnxukujbvxsfuqunfllqmj21ji5urvixgxkyiwhzgfedm2zuj/3ndfmmrpxmqszertlguppjgtgz/k4hqnkgcwekzfvzwwkzuu42nzinwvt+ezw0lqqew/xulyv16zyoipzakzydbzwowy00oakmhvamr/dmsoffexc+fq0fj585hj9wpn8ayngncg==</latexit> sha1_base64="ocsr/ttg5yk/zw3rrtsde4ht=">aaab6nicbva9swnbej2lxzf+rs1tfongfe7smdjgxuimg9ijrc3muuw7o0du3tcopitbcwusfux2flv3crxaokdgcd7m8zmcxlbtxhdb6ewtb2zu1fclx0chh2fle/pojpofcm2i0wseghvkljetufgc9rsknade3i787hmqzwp5agj+hedsx5yro2vhprd2rbccavuemstedmpqi7wspw1gmusjvaajqjwfc9njj9rztgtoc8nuo0jzvm6xr6lkkao/wx56pxcwwvewljzkos1d8tg20nkwb7yomeh1byh+5/vte9b9jmsknsjzalgcmjisvibjlhczstmesout7csnqgkmmptkdkqvpwxn0mnvvxcqnfvvhrnpi4ixmalximhn9cao2hbgxim4rle4c0rzovz7nyswgtopnmof+b8/gdph1/</latexit> sha1_base64="vr4nd5thwpyjneewnuzry7onql8=">aaab6nicbva9swnbej2lxzf+rs1tfongfe60igxaxipenb+qhgfvm5cs2ds7dveecoqn2fgousvsvpfuemu0mqha4/3zpizfysca+o6305h3nre6e4w9rbpzg8kh+fthwckotfotdqoquxcjlconwg6ikeabwe4wuz37nsdumsfy0uwt9cm6kjzkjborptqg14nyxa26c5b14uwkajmag/jxfxiznejpmkba9zw3mx5glefm4kzutzumle3ochuwshqh9rpfqtnyzuhcwnlsxqyuh9pzdtsehoftjoizqxxvbn4n9dltxjjz1wmquhjlovcvbatk/nfzmgvmiomllcmul2vsdfvlbmbtsmg4k2+ve7av1xprxr3bqxeyomowhmcwyv4uim63eetwsbgbm/wcm+ocf6cd+dj2vpw8plt+apn8wfrc2a</latexit> sha1_base64="shvvhrrcq35aw8copgerzfcxo1g=">aaab6nicbva9swnbej2lxzf+rs1tfongfe7smdjgyurzqckr9jbzcvl9vao3b1aopitbcwusfux2flv3crxaokdgcd7m8zmcxlbtxhdb6ewtb2zu1fclx0chh2fle/p2jpofcmwi0wsughvkljelufgddrsknace3c78zhsv5rf8mrme/ioja85o8zkj9nbbvcuufv3cbjjvjxuiedzup7qd2owrigne1trnucmxs+ompwjnjf6qcaesgkdc9ssspufr8du6urdikaxssuow6u+jjezaz6ladkbujpw6txd/83qpcet+xmwsgprstshmbtexwfxnhlwhm2jmcwwk21sjg1nfmbhplgwi3vrlm6rdq3pu1xtwk437pi4ixmalximhn9cao2hccxim4ble4c0rzovz7nyswgtopnmof+b8/gamhi2n</latexit> sha1_base64="k+btnaj44a3cv/9mtdgwzrwdzg=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0l60wpbiyepad+gdwwznbrln5uwuxfk6e/w4kerr/4ib/4bt20o2vpg4pheddpzwlrwbtzv2yltbg5t75r3k3v7b4dh7vfjwyezthiiuhun6qabzfmtwi7kkarwk7istm7nfeuklesifzttfikjyspoqlhsaw7qa7fq1bwfydrxc1kfas2b+9ufjiyluromqn930tnkfnlobm4q/qzjsllezrcnqwsxqidfhhqjfxzuiirnmshizu3xm5jbwexqhtjkkz61vvlv7n9titxqc5l2lmullloigtxcrk/jczcoxmikkllclubyvstbvlxqztssh4qy+vk3a95ns1/96rnu6kompwbudwct5cqqnuoqktdccz3ifn0c4l86787fsltnfzcn8gfp5a/kpjz=</latexit> sha1_base64="qk3tk1ubhnvhfmukmhk38a8ttvo=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0le0gpbiyepad+gdwwznbrln5uwuxfk6e/w4kerr/4ib/4bt20o2vpg4pheddpzwlrwbtzv2ymtrw9sbpw3kzu7e/sh7ufrsyezthkiuhuj6qabzfnnwi7kqkarwkbifjm5nffkklesifzstfikzdyspoqlhsa/v+27vq3lzkfxif6qkbrp996s3sfgwozrmuk27vpeaikfkcczwwullglpkxnsixusljveh+fzuktmzyobeibildzmrvydygms9iupbgvmz0svetpzp62mug5yltpmogslrvemieni7g8y4aqzernlkfpc3kricrkje2nkpwl19eja2lmu/v/huvwr8r4ijdczzcofhwbxw4hq0gceqnuev3hzhvdjvzseitequm8fwb87nd/wxjzg=</latexit> Observed word The main focus of our model is to capture repeating patterns in topical transitions between parent and child events. We posit the existence of a fixed number of topics K. The topics, denoted {ζ k }, are assumed to be probability distributions over words (vocabulary with size W) and are generated from a Dirichlet distribution. Since our data of interest is tweets (short documents), we model a document at event e as having a single hidden topic ηe. We also assume the existence of a topictopic interaction matrix T, where T k, again sampled from a Dirichlet distribution, denotes the distribution over child topics for a parent topic k. In order to generate the document de and event e, we first follow a Markov process for sampling the topic ηe of the document conditioned on the topic of its parent event ηze, followed by sampling the words according to the chosen topic. For spontaneous events occurring at any user node u, the topic for its document is sampled randomly from the preferred distribution over topics φu for node u. An important consideration in the design of our model is the use of conjugate priors. As we will see in the Inference section such priors play a crucial role in the design of efficient and simple sampling-based inference algorithms. Models such as HTM [11], which have to sacrifice conjugacy to model data complexity, have to resort to more complex variational inference strategies. Figure 1 shows the plate diagram of the generative model. We summarize the entire generative process below. 1) Generate (te, ce, ze ) for all events according to the process described in previous sub-section. 2) For each topic k: sample ζ k DirW (α) 3) For each topic k: sample T k DirK (β) 4) For each node v: sample φv DirK (γ) 5) For each event e at node ce = v: a) i) if ze = 0 (level 0 event): draw a topic ηe DiscreteK (φv ) ii) else: draw a topic ηe DiscreteK (T ηze ) b) Sample document length Ne P oisson(λ) c) For w = 1... Ne : draw word xe,w DiscreteW (ζ ηe ) sha1_base64="nv4rdzpte2wl5a6vknquqvmdjc=">aaab6nicbvbns8naej3ur1q/oh69lbbbu0m86euoepekfe0htkfstpn26wtdjdccf0jxjwo4tvf5m1/47bnqvsfddzem2fmxpgkro3nftultfwnza3ydmvnd2//wd08aukkuwyblbgj6oruo+asm4bgz1ui1dge1wfdpz20+one/ko5mkgmr0khnegtvwerhgv+9wvzo3b1klfkgqukdrd796g4rlmurdbnw663upcxkqdgccp5vepjglbeyh2lvu0hh1km9pnzizqwxilchb0pc5+nsip7hwkzi0nte1i73szct/vg5moqsg5zlndeq2wbrlgpiezp4ma66qgtgxhdlf7a2ejaiiznh0kjef/nlvdk6qplezb/3qvw7io4ynmapnimpl1chw2haexgm4rle4c0rzovz7nwswktomxmmf+b8/gc9o1z</latexit> sha1_base64="nkvkz+558ccv1p97d30fpbaiyz4=">aaab8nicbvbnswmxem3wr1q/qh69bivgqws96euoepekfewhtmustbntadzzklmhlv0zxjwo4tvf481/9ruqvsfddzem2fmxprkgqb6+0tr6xuvxeruzs7u0fva+p2lznhvew01kbbkqtl0lxfgiqvjsatpni8k40vpn5nudurndqaspdxi6vciwjiktek9hzkn/iq8xcas1uidz4fxif6sgcjtd6ld/ofmwcavmumt7pkkhykkbwssfvvqz5sllzrkpucvtbgn8vnju3zmlagotxglam/v3xm5taydjjhrtcim7li3e//zehnev0euvjobv2yxkm4kbo1n/+obmjybndhcmrhuvsxg1faglqwkc8fffnmvtc/qpqn796twucvikkmtdirok8uuqpdoizqi0ekav6m0d78v79z4wrswvmdlgf+b9/gdzcjbj</latexit> sha1_base64="mm2je9w3bkjtlfphq9hzg0rqpza=">aaab8nicbvbnswmxem3wr1q/qh69bivgqwrf0itq8ojjkthaajclm862odlksbjcxfozvhhqxku/xpv/xrtdg7+ghi8n8pmvcgv3fhcvr3syura+kz5s7k1vbo7v90/abuvaqtpotsngaefxcy3irojnqoekk4ceaxu/9h0fqhit5b8cpbakdsb5zrq2tuk9hduh5bf9helzrpe5mwmvel0gnfwig1a9ex7esawmzomz0fzlaikfacizguullbllkrnqaxucltcae+ezkct5xsh/hsrusfs/u3xm5twzj5hrtkgdmkvvkv7ndtmbxw5l2lmqbl5ojgt2co8/r/3uqzmxdgryjr3t2i2pjoy61kqubd8xzexsfus7po6f0dqjdsijji6qsfofpnoajxqdwqifmjiowf0it hvlxkftoh6a+8zx/4c5bm</latexit> k Topic-word distribution Dirichlet hyperparameter for topic-word distribution sha1_base64="f4rbolcfgxon0oitki69iqchmw=">aaab/hicbvdlssnafj34rpuv7dlnbfcluqexunbjcsk9gfnkdetstt0mgkzeyge+ituxcji1g9x5984abpq1gpdhm65lzlzgpqzpr3n21pb39jc2q7t1hf39g8o7apjnkoyswixjdyrgwau5uzqrmaa00eqkcqbp/1gelv6/ucqfuveg85t6scwfixiblsrrnbdcxieqjw2v+ebtycwg9lnp+xmgvejw5emqtaz2v9emjaspkitdkonxsfvfgfsm8lpro5liqzapjcmq0mfxft5xtz8dj8zjcrris0rgs/v3xsfxkrmzyzj0bo17jxif94w09g1xzcrzpoksngoyjjwcs6bwcgtlgiegwjempmvkwliinr0vtclumtfxiw9i5brtnz7y2b7pqqjhk7qktphlrpcbxshoqilcmrrm3pfb9at9wk9wx+l0twr2mmgp7a+fwclq5vs</latexit> sha1_base64="wkfk2ydhywuan2umovzp7lr62z0=">aaab/3icbvbps8mwhe3nvzn/vquvxojd8draefqkay8ej7hnwetj03qls5ospmkspfhvvhhqxktfw5vfxntrqtcfhdze+/3iywttrpv2ng+rtrs8srpwx29sbg5t79i7ez0lmoljfwsm5f2ifgguk66mmpg7vbkuhiz0w/fv6ffvivru8fs9smfocgnmcvigymwd7xqsehnenpl3gprkmjhrrhtafltaexivurjqjqcewvlxi4swjxmcglbq6taj9hulpmsnhwmkvshmdosaagcpqq5eft/au8nkoeyhn4rpo1d8boupugdfmjkip1lxxiv95g0zhf35oezppwvhsothjuatlgejkgnwbgiiwpkarbcpkerm8oapgr3/sulphfacp2we3pwbf9wddtbitgcj8af56anrkehdaegj+azvii368l6sd6tj9lozap29sefwj8/qkow4a==</latexit> K sha1_base64="jckw2oyx2xdrrp5bgtbddaxqss=">aaab6hicbvbns8naej3ur1q/qh69lbbbu0le0jmuvahewraf0iay2u7atztn2n0ijfqxepggifd/kjf/jds2b219mpb4b4azeueiudau++0u1t3nrek26wd3b39g/lhuuvhqwlzlgivsegggwx2dtccowkcmkucgwh49uz335cpxksh8wkqt+iq8ldzqixuuo+x664vxcoskq8nfqgr71f/uonpzgka0tvouu5ybgz6gynamclnqpxosymr1i11jji9r+nj90ss6smibhrgxjq+bq74mmrlpposb2rtsm9li3e//zuqkjr/2myyq1knliuzgkmiy+5omuejmxmqsyhs3txi2oooy7mp2rc85zdxseui6rlvr3fzqd3kcrthbe7hhdy4ghrcqr2awadhgv7hzxl0xpx352prwndymwp4a+fzb6drjms=</latexit> Figure 1: Graphical Model for HMHP (shown for 2 nodes and only 4 sample events) The resultant joint likelihood can we written as follows: P (E, Φ, T, ζ, η, z α, β, γ, W, µ) = P (φv γ) v V P (ηe T ηze )δze,e0 P (ηe φv )δze,0 # P (xe,w ηe, ζ ηe ) w=1 Z exp! T µv (τ )dτ 0 v V " Z hce,v ( τ )dτ te # µv (te )δce,v δze,0 e E #! T exp "N e " e E v V P (T k β) k=1 e E K e :te0 <te e E P (ζ k α) k=1 0 K hce,ce0 ( (t )) e0 δc e0,v δz 0,e e e0 E (4) Here, δce,v is an indicator for the event e being triggered at node v, δze0,e is an indicator for event e being parent of event e0, τ is (τ te ), (te0 ) is (te0 te ) and T is the time horizon. The first line correspond to drawing user topic preference vectors, vectors for word distributions over topics, and topic-topic probability vectors from the corresponding Dirichlet. The second line corresponds to the probability of selecting the event topic depending on the topic of the parent event. The third line is for generating the event words given the event topic and topic distribution. The fourth line corresponds to the base intensity of the Hawkes processes and the fifth line captures the impulse response. The last term can be interpreted as the impulse response of each event e on all events e0 triggered at node v.
4 INFERENCE Given model definition, the underlying graph structure and the observable features of the events E, our task is to identify the latent variables associated with each event and also estimate the hidden parameters of the model. The latent event variables are the topic η e and the diffusion parent z e for identifying the cascade diffusion structure. For each event e, we will either decide the event to be spontaneous (z e = 0), or identify the unique parent event e that triggered the creation of e (z e = e ). The process parameters to be estimated include the user-user influence values W uv, the topic transition matrix T, and the user-topic preferences. Since exact inference is intractable, we perform inference using Gibbs Sampling, where the strategy is to iteratively sample the value of each hidden variable from its conditional distribution given the current values of all the other variables, and continue this process until convergence. This process can be made more simple and efficient for our model. By making use of the conjugate priors, we can perform collapsed Gibbs Sampling, where we integrate out some of the continuous valued parameter variables from the likelihood function in Eqn.4. Specifically, we integrate out the topic distributions over words ζ, the topic interaction distributions T and the user-topic preference distributions Φ. The continuous valued latent variables that remain are the connection strengths W u,v. The resulting collapsed Gibbs sampling algorithm iteratively samples individual topic and parent assignments and the connection strengths from their conditional distributions given current assignments to all other variables until convergence. The parameter variables that were integrated out are estimated from the samples of the remaining variables upon convergence. We next describe the conditional distributions for sampling the different variables in the algorithm. The overall algorithm is described in Algo.1. Topic assignment: The conditional probability for topic η e of a diffusion event e (z e 0) being k when the current parent has topic k is the following: P (η e = k {x e }, η ze = k, η \ze, {z e }) β k + N ( (ze,e)) k,k ( l β l) + N ( (ze,e)) k w d e N w e 1 K N (Ce) k,l 1 l =1 i=0 (β l + N ( Ce) k,l + i) (Ce) N k 1 i=0 (( l β Ce l ) + N k + i) i=0 (α w + T e k,w + i) Ne 1 i=0 (( w W α w) + T e k + i) (5) denotes the number of parent-child event pairs with Here N (C) k,k topics k and k, ( (z e, e)) denotes all edges excluding (z e, e), C e is the set of edges from event e to its child events ( C e being its complement), and T e k,w is the number of occurrences of word w under topic k in events other than e. Further, N C k = k N C k,k, and N e w is count of w in d e. The first term is the conditional probability of transitioning from parent topic k to this event s topic k, the second term is that of transitioning from this topic k to each child event s topic l, and the third term is that of observing the words in the event document given topic k. It is important to observe how this conditional distribution pools together three different sources of evidence for an event s topic. Even when the document words do not provide sufficient evidence for the topic, the parent and children topics taken together can significantly reduce the uncertainty. For spontaneous events (z e = 0), the conditional probability looks very similar to Eqn. 5. Only the first term changes to γ k + U ( e) v,k ( k γ k) + U ( e) v, where U ( e) v,k is the number of events by user v with topic k discounting event e. This term captures the probability of node v picking topic k from among its preferred topics. Parent assignment: The conditional probability of event event e being the parent z e of an event e looks as follows: P (z e = e E t, z e, W, µ) (β k + N k,k 1) (( K k=1 β k) + N k 1) h u e,u e (t e t e ) Here the first term is the transition probability from topic k of the proposed parent event e to this event s topic k. The second term is the probability of t e being the time of this event e given the occurrence time t e of the proposed parent event e. As with topic identification, we see that evidence for the parent now comes from two different sources. When there is uncertainty about the parent based on the event time, common patterns of topic transitions between existing parentchild events helps in identifying the right parent. The conditional probability of event e being a spontaneous event with no parent is given as: P (z e = 0 E t, z e, W, µ) (γ k + U ue,k 1) (( K k=1 γ k) + U ue 1) µ u e (t e ) Here again, we have two terms related to the topic and the event time. But the topic term captures the probability of user u e spontaneously picking topic k for an event, and the time term captures the probability of the same user spontaneously generating an event at time t e given its base intensity. In theory, every preceding event is a candidate parent for an event e. We limit the number of possible parent candidates for an event by the time interval between events (1 day) to a maximum of 100 candidates. Updating network strengths: For the network strength W u,v, using a Gamma prior Gamma(α, β), the posterior distribution can be approximated as follows: P (W u,v = x E (u,v) t, z) x α1 exp( xβ 1 ) where α 1 = (N u,v +α 1) and β 1 = (N u + 1 β ) 1, and N u,v is the number of parent-child events pairs between nodes u and v, and N u is the number of events at node u. This is again a Gamma distribution Gamma(α 1, β 1 ). Instead of sampling W uv, we set it to be the mean of the corresponding Gamma
5 distribution. Note that in each iteration, we update W uv only for edges that have at least one parent-child influence. A practical issue with estimating W uv is that in real datasets, most edges have very few (typically just one) influence propagation event. This makes statistical estimation of their strengths infeasible. To get around this problem, we share parameters across edges. While there may be many ways to group similar edges, we group together edges that have the same value for the tuple (out-degree(source), indegree(destination)). The intuition is that the influence of the edge (u, v) is determined uniquely by the the popularity (outdegree) of u and the number of different influencers (indegree) of v. We then pool data from all edges in a group and estimate a single connection strength for a group. Once the Markov chain has converged, the parameters that were integrated out are estimated using the samples: ˆζ k,r = T k,r + α r U v,k + γ k ˆφv,k W r =1 T k,r + α = K r k =1 U v,k + γ k ˆT k,k = N k,k + β k K t=1 N k,t + β t Here, T k,r is the number of occurrences of word r in events with topic k, U v,k is the number of events posted by user v with topic k, and N k,k is the number of parent-child events with topics k and k respectively. Algorithm 1 Gibbs Sampler Initialize η e for all events Initialize z e for all events Initialize W u,v for all u, v in the followers map for iter = 0; iter! = maxiter; iter + + do for e allevents do Sampling topic η e P (η e = k η e, z, {X}, W, α, β, γ, µ) for e allevents do Sampling Parent z e P (z e = e E t, z e, η, α, β, γ, W, µ) for (u, v) Edges do Estimate user-user influence W u,v = Mean(Gamma(N uv + α, N u + β )) EXPERIMENTS In this section, we empirically validate the strengths of HMHP against competitive baselines over a large collection of real tweets as well as semi-synthetic data. We first discuss the baseline algorithms for comparison, the datasets on which we evaluate these algorithms, the tasks and the evaluation measures, and finally the experimental results. Evaluated Models Recall that our model captures network structure, textual content and timestamp of the posts / tweets, and identifies the topics and parents of the posts, in addition to reconstructing the network connection strengths. Considering this, we evaluate and compare performance for the following models: HMHP: This is the full version of our model with Gibbs sampling-based inference. Recall that this performs all reconstruction tasks mentioned above jointly, while assuming a single topic for a post and topical interaction patterns. HTM: This is the state-of-the-art model [11] closest to ours that addresses the same tasks. The key modeling differences are two-fold: HTM assumes a topical admixture for each document instead of a mixture, and secondly it models parent and child events to be close in terms of their topical admixture. As a result, it cannot capture any sequential pattern in the topical content of cascades. Additionally, it uses a reasonably complex variational inference strategy due to the absence of conjugacy in the prior distributions. We used an implementation of HTM kindly provided to us by the authors [11]. HWK+Diag: This is a simplification of our model where the topic-topic interaction is restricted to be diagonal. In other words, each topic interacts only with itself. However, the reconstruction tasks are still performed jointly. This model helps in our evaluation in two different ways. First, it helps in assessing the importance of topical interactions. Secondly, this serves as a crude approximation of HTM with only one topic per document. Since the algorithm for the original HTM did not scale for our larger datasets, we use this model as a surrogate for evaluation. HWK LDA: This is motivated by the Network Hawkes model [13], which jointly models the network strengths and the event timestamps, and therefore performs parent assignment and user-user influence estimation jointly. However, it does not model the textual content of the events. Therefore, we augment it with an independent LDA mixture model [15] (LDAMM) to model the textual content. While one view of this model is as an augmentation of the Network-Hawkes model, the other view is that of a simplification of our model, where the topic assignment is decoupled from parent assignment and network reconstruction. Thus, comparison against this model helps in analyzing the importance of doing topic, parent and network strength estimation jointly. The network reconstruction component of HWK LDA is similar to the Network Hawkes model [13], with the only difference being that the hyper-parameters of the model are not estimated from the data. Datasets We perform experiments on two datasets, that we name Twitter and SemiSynth. The Twitter dataset was created by crawling 7M seed nodes from Twitter for the months of March-April We restrict ourselves to 500K tweets corresponding to top 5K hashtags from the most prolific 1M users generated in a contiguous part of March For each tweet, we have the time stamp, creator-id and tweet-text. Note that gold-standard for the parameters or the event labels that we look to estimate is not available in this Twitter dataset. We do not know the true network connection strengths, event topics or cascade structure of the tweets. While retweet
6 information is available, it is important to point out that retweets form a very small fraction of the parent-child relations that we hope to infer. Since we need gold-standard labels to evaluate performance of the models, we additionally create a SemiSynth dataset using the generative process of our model while preserving statistics of the real data to the extent possible. From a sample of the Twitter data we retain the underlying set of nodes and the follower graph. Then we estimate all the parameters of our model (the base rate per user, user-user influence matrix, the topic distributions, resulting topical interactions and the user-topic preferences) from Twitter data. The document lengths are randomly drawn from P oisson(7), since 7 was the average length of the tweets in the Twitter dataset. We finally generate 5 different samples of 1M events using our generative model. All empirical evaluations are averaged over the 5 generated samples (together termed as SemiSynth). The details of the parameter estimations are as follows. We assign as the parent for a tweet e from user u the closest in time preceding tweet from the followees of user u for the last 1 day. If this set is empty, then e is marked as a spontaneous event. Topics to every tweet are assigned by fitting a latent Dirichlet analysis based mixture model (LDAMM [15]) with 100 topics. From the parent and topic assignments, we get the network strengths, user-topic preferences, topic-topic interactions and topic-word distributions. For W uv estimation (as well as subsequent generation) we applied the edge grouping described at the end of Inference section and then estimate the W uv for an edge as the smoothed estimate for the group. Tasks and experimental results We address three tasks. (A) Reconstruction accuracy, (B) Generalization performance and (C) Discovery and analysis of topic interactions. In task (A), we use the SemiSynth dataset to compare the ability of different algorithms to recover the topic and parent of each event and also the network connection strengths. In task (B), we use the Twitter dataset to evaluate how well the algorithms fit held-out data after their parameters have been estimated using training data. In task (C), we investigate the topical interaction matrix generated by HMHP to find interesting and actionable insights about information cascades in the Twitter dataset. Finally, we also briefly demonstrate the scalability of the inference algorithm for HMHP. We next describe the experimental setup for each of the tasks and then the results. Reconstruction Accuracy: We address three reconstruction tasks: network reconstruction, parent identification and topic identification. Since the gold-standard values are not known for the Twitter data, we address these tasks only on the SemiSynth data. For network reconstruction, we measure distance between estimated and true W uv values. We report the error in terms of the median Average Percentage Error (APE), defined as W uv Ŵuv u,v W uv, where, Ŵ uv is the estimated W uv value. For comparison with HTM, we use Total Absolute Error (TAE) (the sum of the absolute differences between the true and estimated W uv values) which was used in the HMHP HWK+Diag HWK LDA Accuracy Recall@ Recall@ Recall@ Recall@ Table I: Cascade Reconstruction: accuracy and the recall with top-k (1,3,5,7) candidate parents. HMHP HWK+Diag HWK LDA Mean APE Median APE Mean APE (N u 100) Median APE (N u 100) Table II: Network reconstruction error (as fraction). original paper. For cascade reconstruction, given a ranked list of predicted parents for each event, we measure accuracy and recall of the parent prediction. The accuracy metric is calculated as the ratio of number events for which the parent is identified correctly to the total number of events. For topic identification, we take a (roughly balanced) sample of event pairs, and then measure (using precision/recall/f1) whether the model accurately predicts the event pairs to be from the same or different topics. The model that is closest to ours is HTM [11]. However, the HTM assigns each event to a topic distribution instead of a single topic. Also, the available HTM code cannot scale to the size of our SemiSynth dataset. Therefore, we compared our model with HTM separately. We first present comparisons with the other baselines on the SemiSynth dataset. All the presented numbers are averaged over five different semisynthetic datasets generated independently. Table I presents the accuracy and the recall@k values for the task of predicting parents, and thereby the cascades. For all the algorithms, the recall improves significantly by the time the top-3 candidates are considered. For recall@3, HMHP performs about 32% better than both the other baselines, whereas recall@1, as well as the accuracy of HMHP are at least 57% better than corresponding number of the baselines. For the network reconstruction task, Table II describes the various APE values for models HMHP, HWK+Diag, and HWK LDA. Given that in our simulation, event generation was truncated by limiting the number of events, the last level of events contributes to negative bias affecting the APE, since the algorithms are still trying to assign children to these events. We report both the mean and the median errors for this task. The mean error for HMHP is about 18% lower than other baselines, and the median error is about 10% lower. We also present the results for the set of nodes where the number of generated children (N u ) is at least 100 (this is an easier task). Here also the mean APE for HMHP is at least 20% lower. Table III presents the result of the topic modeling task, measured by calculating the precision-recall over event pairs as defined earlier. The results show that HMHP performs much better than HWK+Diag and at least 5 6% better than
STA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
Study Notes on the Latent Dirichlet Allocation
Study Notes on the Latent Dirichlet Allocation Xugang Ye 1. Model Framework A word is an element of dictionary {1,,}. A document is represented by a sequence of words: =(,, ), {1,,}. A corpus is a collection
Text Mining for Economics and Finance Latent Dirichlet Allocation
Text Mining for Economics and Finance Latent Dirichlet Allocation Stephen Hansen Text Mining Lecture 5 1 / 45 Introduction Recall we are interested in mixed-membership modeling, but that the plsi model
13: Variational inference II
10-708: Probabilistic Graphical Models, Spring 2015 13: Variational inference II Lecturer: Eric P. Xing Scribes: Ronghuo Zheng, Zhiting Hu, Yuntian Deng 1 Introduction We started to talk about variational
Time-Sensitive Dirichlet Process Mixture Models
Time-Sensitive Dirichlet Process Mixture Models Xiaojin Zhu Zoubin Ghahramani John Lafferty May 25 CMU-CALD-5-4 School of Computer Science Carnegie Mellon University Pittsburgh, PA 523 Abstract We introduce
Learning with Temporal Point Processes
Learning with Temporal Point Processes t Manuel Gomez Rodriguez MPI for Software Systems Isabel Valera MPI for Intelligent Systems Slides/references: http://learning.mpi-sws.org/tpp-icml18 ICML TUTORIAL,
STA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 7 Approximate
A graph contains a set of nodes (vertices) connected by links (edges or arcs)
BOLTZMANN MACHINES Generative Models Graphical Models A graph contains a set of nodes (vertices) connected by links (edges or arcs) In a probabilistic graphical model, each node represents a random variable,
CSCI-567: Machine Learning (Spring 2019)
CSCI-567: Machine Learning (Spring 2019) Prof. Victor Adamchik U of Southern California Mar. 19, 2019 March 19, 2019 1 / 43 Administration March 19, 2019 2 / 43 Administration TA3 is due this week March
Click Prediction and Preference Ranking of RSS Feeds
Click Prediction and Preference Ranking of RSS Feeds 1 Introduction December 11, 2009 Steven Wu RSS (Really Simple Syndication) is a family of data formats used to publish frequently updated works. RSS
Latent Dirichlet Allocation
Latent Dirichlet Allocation 1 Directed Graphical Models William W. Cohen Machine Learning 10-601 2 DGMs: The Burglar Alarm example Node ~ random variable Burglar Earthquake Arcs define form of probability
Collapsed Variational Bayesian Inference for Hidden Markov Models
Collapsed Variational Bayesian Inference for Hidden Markov Models Pengyu Wang, Phil Blunsom Department of Computer Science, University of Oxford International Conference on Artificial Intelligence and
Chris Bishop s PRML Ch. 8: Graphical Models
Chris Bishop s PRML Ch. 8: Graphical Models January 24, 2008 Introduction Visualize the structure of a probabilistic model Design and motivate new models Insights into the model s properties, in particular
Linear Dynamical Systems
Linear Dynamical Systems Sargur N. srihari@cedar.buffalo.edu Machine Learning Course: http://www.cedar.buffalo.edu/~srihari/cse574/index.html Two Models Described by Same Graph Latent variables Observations
Data Mining Techniques
Data Mining Techniques CS 622 - Section 2 - Spring 27 Pre-final Review Jan-Willem van de Meent Feedback Feedback https://goo.gl/er7eo8 (also posted on Piazza) Also, please fill out your TRACE evaluations!
ICML Scalable Bayesian Inference on Point processes. with Gaussian Processes. Yves-Laurent Kom Samo & Stephen Roberts
ICML 2015 Scalable Nonparametric Bayesian Inference on Point Processes with Gaussian Processes Machine Learning Research Group and Oxford-Man Institute University of Oxford July 8, 2015 Point Processes
Deep Poisson Factorization Machines: a factor analysis model for mapping behaviors in journalist ecosystem
000 001 002 003 004 005 006 007 008 009 010 011 012 013 014 015 016 017 018 019 020 021 022 023 024 025 026 027 028 029 030 031 032 033 034 035 036 037 038 039 040 041 042 043 044 045 046 047 048 049 050
Stance classification and Diffusion Modelling
Dr. Srijith P. K. CSE, IIT Hyderabad Outline Stance Classification 1 Stance Classification 2 Stance Classification in Twitter Rumour Stance Classification Classify tweets as supporting, denying, questioning,
Computational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
A Bivariate Point Process Model with Application to Social Media User Content Generation
1 / 33 A Bivariate Point Process Model with Application to Social Media User Content Generation Emma Jingfei Zhang ezhang@bus.miami.edu Yongtao Guan yguan@bus.miami.edu Department of Management Science
Evaluation Methods for Topic Models
University of Massachusetts Amherst wallach@cs.umass.edu April 13, 2009 Joint work with Iain Murray, Ruslan Salakhutdinov and David Mimno Statistical Topic Models Useful for analyzing large, unstructured
CS839: Probabilistic Graphical Models. Lecture 7: Learning Fully Observed BNs. Theo Rekatsinas
CS839: Probabilistic Graphical Models Lecture 7: Learning Fully Observed BNs Theo Rekatsinas 1 Exponential family: a basic building block For a numeric random variable X p(x ) =h(x)exp T T (x) A( ) = 1
27 : Distributed Monte Carlo Markov Chain. 1 Recap of MCMC and Naive Parallel Gibbs Sampling
10-708: Probabilistic Graphical Models 10-708, Spring 2014 27 : Distributed Monte Carlo Markov Chain Lecturer: Eric P. Xing Scribes: Pengtao Xie, Khoa Luu In this scribe, we are going to review the Parallel
Bayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
Gentle Introduction to Infinite Gaussian Mixture Modeling
Gentle Introduction to Infinite Gaussian Mixture Modeling with an application in neuroscience By Frank Wood Rasmussen, NIPS 1999 Neuroscience Application: Spike Sorting Important in neuroscience and for
Graphical Models and Kernel Methods
Graphical Models and Kernel Methods Jerry Zhu Department of Computer Sciences University of Wisconsin Madison, USA MLSS June 17, 2014 1 / 123 Outline Graphical Models Probabilistic Inference Directed vs.
STA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 11 Project
Latent Dirichlet Allocation Introduction/Overview
Latent Dirichlet Allocation Introduction/Overview David Meyer 03.10.2016 David Meyer http://www.1-4-5.net/~dmm/ml/lda_intro.pdf 03.10.2016 Agenda What is Topic Modeling? Parametric vs. Non-Parametric Models
STA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
Kernel Density Topic Models: Visual Topics Without Visual Words
Kernel Density Topic Models: Visual Topics Without Visual Words Konstantinos Rematas K.U. Leuven ESAT-iMinds krematas@esat.kuleuven.be Mario Fritz Max Planck Institute for Informatics mfrtiz@mpi-inf.mpg.de
Scaling Neighbourhood Methods
Quick Recap Scaling Neighbourhood Methods Collaborative Filtering m = #items n = #users Complexity : m * m * n Comparative Scale of Signals ~50 M users ~25 M items Explicit Ratings ~ O(1M) (1 per billion)
Topic Models. Brandon Malone. February 20, Latent Dirichlet Allocation Success Stories Wrap-up
Much of this material is adapted from Blei 2003. Many of the images were taken from the Internet February 20, 2014 Suppose we have a large number of books. Each is about several unknown topics. How can
STA 414/2104: Machine Learning
STA 414/2104: Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistics! rsalakhu@cs.toronto.edu! http://www.cs.toronto.edu/~rsalakhu/ Lecture 9 Sequential Data So far
Variational Inference (11/04/13)
STA561: Probabilistic machine learning Variational Inference (11/04/13) Lecturer: Barbara Engelhardt Scribes: Matt Dickenson, Alireza Samany, Tracy Schifeling 1 Introduction In this lecture we will further
Lecture 13 : Variational Inference: Mean Field Approximation
10-708: Probabilistic Graphical Models 10-708, Spring 2017 Lecture 13 : Variational Inference: Mean Field Approximation Lecturer: Willie Neiswanger Scribes: Xupeng Tong, Minxing Liu 1 Problem Setup 1.1
13 : Variational Inference: Loopy Belief Propagation and Mean Field
10-708: Probabilistic Graphical Models 10-708, Spring 2012 13 : Variational Inference: Loopy Belief Propagation and Mean Field Lecturer: Eric P. Xing Scribes: Peter Schulam and William Wang 1 Introduction
Gaussian Mixture Model
Case Study : Document Retrieval MAP EM, Latent Dirichlet Allocation, Gibbs Sampling Machine Learning/Statistics for Big Data CSE599C/STAT59, University of Washington Emily Fox 0 Emily Fox February 5 th,
1 Random walks and data
Inference, Models and Simulation for Complex Systems CSCI 7-1 Lecture 7 15 September 11 Prof. Aaron Clauset 1 Random walks and data Supposeyou have some time-series data x 1,x,x 3,...,x T and you want
Conquering the Complexity of Time: Machine Learning for Big Time Series Data
Conquering the Complexity of Time: Machine Learning for Big Time Series Data Yan Liu Computer Science Department University of Southern California Mini-Workshop on Theoretical Foundations of Cyber-Physical
Variational Learning : From exponential families to multilinear systems
Variational Learning : From exponential families to multilinear systems Ananth Ranganathan th February 005 Abstract This note aims to give a general overview of variational inference on graphical models.
16 : Approximate Inference: Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models 10-708, Spring 2017 16 : Approximate Inference: Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Yuan Yang, Chao-Ming Yen 1 Introduction As the target distribution
Unified Modeling of User Activities on Social Networking Sites
Unified Modeling of User Activities on Social Networking Sites Himabindu Lakkaraju IBM Research - India Manyata Embassy Business Park Bangalore, Karnataka - 5645 klakkara@in.ibm.com Angshu Rai IBM Research
Construction of Dependent Dirichlet Processes based on Poisson Processes
1 / 31 Construction of Dependent Dirichlet Processes based on Poisson Processes Dahua Lin Eric Grimson John Fisher CSAIL MIT NIPS 2010 Outstanding Student Paper Award Presented by Shouyuan Chen Outline
p L yi z n m x N n xi
y i z n x n N x i Overview Directed and undirected graphs Conditional independence Exact inference Latent variables and EM Variational inference Books statistical perspective Graphical Models, S. Lauritzen
Probabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 295-P, Spring 213 Prof. Erik Sudderth Lecture 11: Inference & Learning Overview, Gaussian Graphical Models Some figures courtesy Michael Jordan s draft
Using first-order logic, formalize the following knowledge:
Probabilistic Artificial Intelligence Final Exam Feb 2, 2016 Time limit: 120 minutes Number of pages: 19 Total points: 100 You can use the back of the pages if you run out of space. Collaboration on the
Advanced Introduction to Machine Learning
10-715 Advanced Introduction to Machine Learning Homework 3 Due Nov 12, 10.30 am Rules 1. Homework is due on the due date at 10.30 am. Please hand over your homework at the beginning of class. Please see
Mixtures of Gaussians. Sargur Srihari
Mixtures of Gaussians Sargur srihari@cedar.buffalo.edu 1 9. Mixture Models and EM 0. Mixture Models Overview 1. K-Means Clustering 2. Mixtures of Gaussians 3. An Alternative View of EM 4. The EM Algorithm
Hidden Markov Models Part 1: Introduction
Hidden Markov Models Part 1: Introduction CSE 6363 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 Modeling Sequential Data Suppose that
DRAFT: A2.1 Activity rate Loppersum
DRAFT: A.1 Activity rate Loppersum Rakesh Paleja, Matthew Jones, David Randell, Stijn Bierman April 3, 15 1 Summary On 1st Jan 14, the rate of production in the area of Loppersum was reduced. Here we seek
19 : Bayesian Nonparametrics: The Indian Buffet Process. 1 Latent Variable Models and the Indian Buffet Process
10-708: Probabilistic Graphical Models, Spring 2015 19 : Bayesian Nonparametrics: The Indian Buffet Process Lecturer: Avinava Dubey Scribes: Rishav Das, Adam Brodie, and Hemank Lamba 1 Latent Variable
Variational Principal Components
Variational Principal Components Christopher M. Bishop Microsoft Research 7 J. J. Thomson Avenue, Cambridge, CB3 0FB, U.K. cmbishop@microsoft.com http://research.microsoft.com/ cmbishop In Proceedings
Introduction to Bayesian inference
Introduction to Bayesian inference Thomas Alexander Brouwer University of Cambridge tab43@cam.ac.uk 17 November 2015 Probabilistic models Describe how data was generated using probability distributions
Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets
Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets Abhirup Datta 1 Sudipto Banerjee 1 Andrew O. Finley 2 Alan E. Gelfand 3 1 University of Minnesota, Minneapolis,
Sparse Stochastic Inference for Latent Dirichlet Allocation
Sparse Stochastic Inference for Latent Dirichlet Allocation David Mimno 1, Matthew D. Hoffman 2, David M. Blei 1 1 Dept. of Computer Science, Princeton U. 2 Dept. of Statistics, Columbia U. Presentation
Modeling User Rating Profiles For Collaborative Filtering
Modeling User Rating Profiles For Collaborative Filtering Benjamin Marlin Department of Computer Science University of Toronto Toronto, ON, M5S 3H5, CANADA marlin@cs.toronto.edu Abstract In this paper
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation. Volker Tresp Summer 2016
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation Volker Tresp Summer 2016 1 Introduction So far we were mostly concerned with supervised learning: we predicted one or several
Petr Volf. Model for Difference of Two Series of Poisson-like Count Data
Petr Volf Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodárenskou věží 4, 182 8 Praha 8 e-mail: volf@utia.cas.cz Model for Difference of Two Series of Poisson-like
Content-based Recommendation
Content-based Recommendation Suthee Chaidaroon June 13, 2016 Contents 1 Introduction 1 1.1 Matrix Factorization......................... 2 2 slda 2 2.1 Model................................. 3 3 flda 3
Lecture 21: Spectral Learning for Graphical Models
10-708: Probabilistic Graphical Models 10-708, Spring 2016 Lecture 21: Spectral Learning for Graphical Models Lecturer: Eric P. Xing Scribes: Maruan Al-Shedivat, Wei-Cheng Chang, Frederick Liu 1 Motivation
Lecture 6: Graphical Models: Learning
Lecture 6: Graphical Models: Learning 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering, University of Cambridge February 3rd, 2010 Ghahramani & Rasmussen (CUED)
Sequential Monte Carlo and Particle Filtering. Frank Wood Gatsby, November 2007
Sequential Monte Carlo and Particle Filtering Frank Wood Gatsby, November 2007 Importance Sampling Recall: Let s say that we want to compute some expectation (integral) E p [f] = p(x)f(x)dx and we remember
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation. Volker Tresp Summer 2014
Bayesian Networks: Construction, Inference, Learning and Causal Interpretation Volker Tresp Summer 2014 1 Introduction So far we were mostly concerned with supervised learning: we predicted one or several
STA 4273H: Sta-s-cal Machine Learning
STA 4273H: Sta-s-cal Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 2 In our
28 : Approximate Inference - Distributed MCMC
10-708: Probabilistic Graphical Models, Spring 2015 28 : Approximate Inference - Distributed MCMC Lecturer: Avinava Dubey Scribes: Hakim Sidahmed, Aman Gupta 1 Introduction For many interesting problems,
RECOVERING NORMAL NETWORKS FROM SHORTEST INTER-TAXA DISTANCE INFORMATION
RECOVERING NORMAL NETWORKS FROM SHORTEST INTER-TAXA DISTANCE INFORMATION MAGNUS BORDEWICH, KATHARINA T. HUBER, VINCENT MOULTON, AND CHARLES SEMPLE Abstract. Phylogenetic networks are a type of leaf-labelled,
Machine Learning for Data Science (CS4786) Lecture 24
Machine Learning for Data Science (CS4786) Lecture 24 Graphical Models: Approximate Inference Course Webpage : http://www.cs.cornell.edu/courses/cs4786/2016sp/ BELIEF PROPAGATION OR MESSAGE PASSING Each
Basic Sampling Methods
Basic Sampling Methods Sargur Srihari srihari@cedar.buffalo.edu 1 1. Motivation Topics Intractability in ML How sampling can help 2. Ancestral Sampling Using BNs 3. Transforming a Uniform Distribution
Sequential Monte Carlo Algorithms
ayesian Phylogenetic Inference using Sequential Monte arlo lgorithms lexandre ouchard-ôté *, Sriram Sankararaman *, and Michael I. Jordan *, * omputer Science ivision, University of alifornia erkeley epartment
6.867 Machine learning, lecture 23 (Jaakkola)
Lecture topics: Markov Random Fields Probabilistic inference Markov Random Fields We will briefly go over undirected graphical models or Markov Random Fields (MRFs) as they will be needed in the context
11 : Gaussian Graphic Models and Ising Models
10-708: Probabilistic Graphical Models 10-708, Spring 2017 11 : Gaussian Graphic Models and Ising Models Lecturer: Bryon Aragam Scribes: Chao-Ming Yen 1 Introduction Different from previous maximum likelihood
Learning Sequence Motif Models Using Expectation Maximization (EM) and Gibbs Sampling
Learning Sequence Motif Models Using Expectation Maximization (EM) and Gibbs Sampling BMI/CS 776 www.biostat.wisc.edu/bmi776/ Spring 009 Mark Craven craven@biostat.wisc.edu Sequence Motifs what is a sequence
Latent Dirichlet Allocation (LDA)
Latent Dirichlet Allocation (LDA) D. Blei, A. Ng, and M. Jordan. Journal of Machine Learning Research, 3:993-1022, January 2003. Following slides borrowed ant then heavily modified from: Jonathan Huang
Introduction to Probabilistic Graphical Models
Introduction to Probabilistic Graphical Models Sargur Srihari srihari@cedar.buffalo.edu 1 Topics 1. What are probabilistic graphical models (PGMs) 2. Use of PGMs Engineering and AI 3. Directionality in
17 : Markov Chain Monte Carlo
10-708: Probabilistic Graphical Models, Spring 2015 17 : Markov Chain Monte Carlo Lecturer: Eric P. Xing Scribes: Heran Lin, Bin Deng, Yun Huang 1 Review of Monte Carlo Methods 1.1 Overview Monte Carlo
Machine Learning Summer School
Machine Learning Summer School Lecture 3: Learning parameters and structure Zoubin Ghahramani zoubin@eng.cam.ac.uk http://learning.eng.cam.ac.uk/zoubin/ Department of Engineering University of Cambridge,
Deep Learning Srihari. Deep Belief Nets. Sargur N. Srihari
Deep Belief Nets Sargur N. Srihari srihari@cedar.buffalo.edu Topics 1. Boltzmann machines 2. Restricted Boltzmann machines 3. Deep Belief Networks 4. Deep Boltzmann machines 5. Boltzmann machines for continuous
Gibbs Sampling in Endogenous Variables Models
Gibbs Sampling in Endogenous Variables Models Econ 690 Purdue University Outline 1 Motivation 2 Identification Issues 3 Posterior Simulation #1 4 Posterior Simulation #2 Motivation In this lecture we take
Recent Advances in Bayesian Inference Techniques
Recent Advances in Bayesian Inference Techniques Christopher M. Bishop Microsoft Research, Cambridge, U.K. research.microsoft.com/~cmbishop SIAM Conference on Data Mining, April 2004 Abstract Bayesian
Lecture 16 Deep Neural Generative Models
Lecture 16 Deep Neural Generative Models CMSC 35246: Deep Learning Shubhendu Trivedi & Risi Kondor University of Chicago May 22, 2017 Approach so far: We have considered simple models and then constructed
On the Problem of Error Propagation in Classifier Chains for Multi-Label Classification
On the Problem of Error Propagation in Classifier Chains for Multi-Label Classification Robin Senge, Juan José del Coz and Eyke Hüllermeier Draft version of a paper to appear in: L. Schmidt-Thieme and
ECE521 Tutorial 11. Topic Review. ECE521 Winter Credits to Alireza Makhzani, Alex Schwing, Rich Zemel and TAs for slides. ECE521 Tutorial 11 / 4
ECE52 Tutorial Topic Review ECE52 Winter 206 Credits to Alireza Makhzani, Alex Schwing, Rich Zemel and TAs for slides ECE52 Tutorial ECE52 Winter 206 Credits to Alireza / 4 Outline K-means, PCA 2 Bayesian
Benchmarking and Improving Recovery of Number of Topics in Latent Dirichlet Allocation Models
Benchmarking and Improving Recovery of Number of Topics in Latent Dirichlet Allocation Models Jason Hou-Liu January 4, 2018 Abstract Latent Dirichlet Allocation (LDA) is a generative model describing the
Bayesian Nonparametrics for Speech and Signal Processing
Bayesian Nonparametrics for Speech and Signal Processing Michael I. Jordan University of California, Berkeley June 28, 2011 Acknowledgments: Emily Fox, Erik Sudderth, Yee Whye Teh, and Romain Thibaux Computer
CSCE 478/878 Lecture 6: Bayesian Learning
Bayesian Methods Not all hypotheses are created equal (even if they are all consistent with the training data) Outline CSCE 478/878 Lecture 6: Bayesian Learning Stephen D. Scott (Adapted from Tom Mitchell
Recurrent Latent Variable Networks for Session-Based Recommendation
Recurrent Latent Variable Networks for Session-Based Recommendation Panayiotis Christodoulou Cyprus University of Technology paa.christodoulou@edu.cut.ac.cy 27/8/2017 Panayiotis Christodoulou (C.U.T.)
Generative Clustering, Topic Modeling, & Bayesian Inference
Generative Clustering, Topic Modeling, & Bayesian Inference INFO-4604, Applied Machine Learning University of Colorado Boulder December 12-14, 2017 Prof. Michael Paul Unsupervised Naïve Bayes Last week
Modeling of Growing Networks with Directional Attachment and Communities
Modeling of Growing Networks with Directional Attachment and Communities Masahiro KIMURA, Kazumi SAITO, Naonori UEDA NTT Communication Science Laboratories 2-4 Hikaridai, Seika-cho, Kyoto 619-0237, Japan
Machine Learning Overview
Machine Learning Overview Sargur N. Srihari University at Buffalo, State University of New York USA 1 Outline 1. What is Machine Learning (ML)? 2. Types of Information Processing Problems Solved 1. Regression
Introduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation. EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016
Introduction to Bayesian Statistics and Markov Chain Monte Carlo Estimation EPSY 905: Multivariate Analysis Spring 2016 Lecture #10: April 6, 2016 EPSY 905: Intro to Bayesian and MCMC Today s Class An
Some Mathematical Models to Turn Social Media into Knowledge
Some Mathematical Models to Turn Social Media into Knowledge Xiaojin Zhu University of Wisconsin Madison, USA NLP&CC 2013 Collaborators Amy Bellmore Aniruddha Bhargava Kwang-Sung Jun Robert Nowak Jun-Ming
Introduction to Probabilistic Machine Learning
Introduction to Probabilistic Machine Learning Piyush Rai Dept. of CSE, IIT Kanpur (Mini-course 1) Nov 03, 2015 Piyush Rai (IIT Kanpur) Introduction to Probabilistic Machine Learning 1 Machine Learning
Forecasting Wind Ramps
Forecasting Wind Ramps Erin Summers and Anand Subramanian Jan 5, 20 Introduction The recent increase in the number of wind power producers has necessitated changes in the methods power system operators
PILCO: A Model-Based and Data-Efficient Approach to Policy Search
PILCO: A Model-Based and Data-Efficient Approach to Policy Search (M.P. Deisenroth and C.E. Rasmussen) CSC2541 November 4, 2016 PILCO Graphical Model PILCO Probabilistic Inference for Learning COntrol
MEI: Mutual Enhanced Infinite Community-Topic Model for Analyzing Text-augmented Social Networks
MEI: Mutual Enhanced Infinite Community-Topic Model for Analyzing Text-augmented Social Networks Dongsheng Duan 1, Yuhua Li 2,, Ruixuan Li 2, Zhengding Lu 2, Aiming Wen 1 Intelligent and Distributed Computing
Probabilistic Time Series Classification
Probabilistic Time Series Classification Y. Cem Sübakan Boğaziçi University 25.06.2013 Y. Cem Sübakan (Boğaziçi University) M.Sc. Thesis Defense 25.06.2013 1 / 54 Problem Statement The goal is to assign
Note for plsa and LDA-Version 1.1
Note for plsa and LDA-Version 1.1 Wayne Xin Zhao March 2, 2011 1 Disclaimer In this part of PLSA, I refer to [4, 5, 1]. In LDA part, I refer to [3, 2]. Due to the limit of my English ability, in some place,
Massachusetts Institute of Technology
Massachusetts Institute of Technology 6.867 Machine Learning, Fall 2006 Problem Set 5 Due Date: Thursday, Nov 30, 12:00 noon You may submit your solutions in class or in the box. 1. Wilhelm and Klaus are
Applying LDA topic model to a corpus of Italian Supreme Court decisions
Applying LDA topic model to a corpus of Italian Supreme Court decisions Paolo Fantini Statistical Service of the Ministry of Justice - Italy CESS Conference - Rome - November 25, 2014 Our goal finding
Topic Models and Applications to Short Documents
Topic Models and Applications to Short Documents Dieu-Thu Le Email: dieuthu.le@unitn.it Trento University April 6, 2011 1 / 43 Outline Introduction Latent Dirichlet Allocation Gibbs Sampling Short Text