Skip to content

Instantly share code, notes, and snippets.

@brendano
Created December 29, 2013 19:29
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save brendano/8173861 to your computer and use it in GitHub Desktop.
Save brendano/8173861 to your computer and use it in GitHub Desktop.
File Function Line
48.0 % /Users/brendano/Desktop/hier_lda/code/treetm.jl cgsIterPath 389 qnewLL,pnewLL = proposePath!(newpath, mm, V,di,word, first ? nothing : oldpath, :simulate)
41.6 % /Users/brendano/Desktop/hier_lda/code/treetm.jl cgsIterPath 390 qoldLL,poldLL = proposePath!(oldpath, mm, V,di,word, first ? nothing : oldpath, :evaluate)
1.1 % /Users/brendano/Desktop/hier_lda/code/treetm.jl cgsIterPath 391 logA = pnewLL-poldLL + qoldLL-qnewLL # (pnew-qnew) - (pold-qnew)
3.0 % /Users/brendano/Desktop/hier_lda/code/treetm.jl cgsIterPath 402 incrementFullpath!(mm.cTopicWord, newpath, word, +1)
3.1 % /Users/brendano/Desktop/hier_lda/code/treetm.jl cgsIterPath 408 incrementFullpath!(mm.cTopicWord, oldpath, word, -1)
2.7 % /Users/brendano/Desktop/hier_lda/code/treetm.jl getindex 82 getindex(c::ClassCountTable, k) = c.counts[k]
2.2 % /Users/brendano/Desktop/hier_lda/code/treetm.jl incrementFullpath! 125 x = b ? x.right : x.left
1.2 % /Users/brendano/Desktop/hier_lda/code/treetm.jl incrementFullpath! 161 x = b ? x.right : x.left
4.0 % /Users/brendano/Desktop/hier_lda/code/treetm.jl incrementFullpath! 162 increment!(x.counts, k, delta)
20.2 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 342 w0 = (nn.counts[wordID] + betaHere/V - on0) / (nn.counts.total + betaHere - int(on_cur))
18.4 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 344 w1 = (nn.counts[wordID] + betaHere/V - on1) / (nn.counts.total + betaHere - int(on_cur))
6.3 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 345 p0 = (cur_docnode.left.count + mm.gammaConc/2 - on0)
6.6 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 346 p1 = (cur_docnode.right.count+ mm.gammaConc/2 - on1)
7.3 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 347 q = p1*w1 / (p1*w1 + p0*w0)
1.5 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 350 dec = bool(rand() < q)
8.5 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 356 proposalLL += dec ? log(q) : log(1-q)
10.5 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 357 unnormLL += dec ? log(p1/(p0+p1)) : log(p0/(p0+p1)) # only the prior branch prob
4.2 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 373 wp = (nn.counts[wordID] + mm.betaConc/V +rem) / (nn.counts.total + mm.betaConc +rem)
1.1 % /Users/brendano/Desktop/hier_lda/code/treetm.jl proposePath! 374 unnormLL += log(wp)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment