• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            posts - 195,  comments - 30,  trackbacks - 0
            At first, let's define the sample function:

            def sample(dist, num_samples=1):
                """
                Uses the inverse CDF method to return samples drawn from an
                (unnormalized) discrete distribution.

                Arguments:

                dist -- (unnormalized) distribution

                Keyword arguments:

                num_samples -- number of samples to draw
                
            """

                cdf = cumsum(dist)
                r = uniform(size=num_samples) * cdf[-1]

                return cdf.searchsorted(r)
            As we can see, the sample function input two parameters, one is dist, which can be an un-normalized distribution, another is the sample we want to draw.

            Let's see how to generate corpus for Dirichlet--multinomial unigram language model
            def generate_corpus(beta, mean, N):
                """
                Returns a corpus of tokens drawn from a Dirichlet--multinomial
                unigram language model. Each token is an instance of one of V
                unique word types, represented by indices 0, , V - 1.

                Arguments:

                beta -- concentration parameter for the Dirichlet prior
                mean -- V-dimensional mean of the Dirichlet prior
                N -- number of tokens to generate
                
            """

                pass # YOUR CODE GOES HERE
                #print mean
                #print beta 
                #print dot(mean,beta)
                #print dirichlet(mean*beta,size=1)
                temp=sample(dirichlet(beta*array(mean),size=1),N)
                #print temp
                return temp
            please keep in mind the dirichlet function is  “from numpy.random.mtrand import dirichlet"
            and the parameters it receives are corresponding to beta*array(mean). beta is the concentration factor, and mean is the vector which sum to 1.



            another way is to generate corpus is using the property:
            P(D'|D,H)= Nv+beta_nv/N+beta
            def generate_corpus_collapsed(beta, mean, N):
                """
                Returns a corpus of tokens drawn from a Dirichlet--multinomial
                unigram language model using the 'collapsed' generative process
                (i.e., phi is not explicitly represented). Each token is an
                instance of one of V unique word types.

                Arguments:

                beta -- concentration parameter for the Dirichlet prior
                mean -- V-dimensional mean of the Dirichlet prior
                N -- number of tokens to generate
                
            """

                V = len(mean) # vocabulary size

                corpus = zeros(N, dtype=int) # corpus

                Nv = zeros(V, dtype=int) # counts for each word type

                pass # YOUR CODE GOES HERE
                for n in xrange(N):
                    corpus[n]=sample((Nv+beta*array(mean))/(n+beta),1)
                    Nv[corpus[n]]+=1;    
                return corpus

            Let's see how to generate corpus for Mixture of Dirichlet-multinomial unigram language model 

            def generate_corpus(alpha, m, beta, n, D, Nd):
                """
                Returns a grouped corpus drawn from a mixture of
                Dirichlet--multinomial unigram language models.

                Arguments:

                alpha -- concentration parameter for the Dirichlet prior over theta
                m -- T-dimensional mean of the Dirichlet prior over theta
                beta -- concentration parameter for the Dirichlet prior over phis
                n -- V-dimensional mean of the Dirichlet prior over phis
                D -- number of documents to generate
                Nd -- number of tokens to generate per document
                
            """
                corpus = GroupedCorpus()

                pass # YOUR CODE GOES HERE
                #determine the topic the distribution for topic dirichlet(dot(m,alpha),size=1)
                #given the topic, the distribtuion for word dirichlet(dot(n,beta),size=1)
                theta=dirichlet(alpha*array(m),1)
                phis=dirichlet(beta*array(n),len(m))
                for d in range(0,D):
                    [t]=sample(theta,1)
                    #print groupVcab
                    corpus.add(str(d),str(t),[str(x) for x in sample(phis[t,:],Nd)])           
                return corpus
            注意是T個topic (group),  phis=dirichlet(beta*array(n),len(m))  產生了T個 dirichlet distribution,相同的topic t應該取同一個 dirichlet distribution phis[t,:]
            posted on 2012-10-28 10:13 luis 閱讀(630) 評論(0)  編輯 收藏 引用 所屬分類: Python
            <2012年10月>
            30123456
            78910111213
            14151617181920
            21222324252627
            28293031123
            45678910

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            文章分類

            文章檔案

            友情鏈接

            搜索

            •  

            最新評論

            閱讀排行榜

            評論排行榜

            久久久一本精品99久久精品88| 99久久99久久| 热99RE久久精品这里都是精品免费 | 久久精品国产亚洲av影院| 久久国产精品无码HDAV| 国产精品嫩草影院久久| 模特私拍国产精品久久| 99久久99久久久精品齐齐| 国内精品久久久久久中文字幕| 要久久爱在线免费观看| 亚洲午夜久久影院| 精品综合久久久久久97| 精品久久久久久无码中文野结衣| 久久精品国产欧美日韩99热| 日本精品久久久中文字幕| 亚洲伊人久久综合中文成人网| 精品国产乱码久久久久久人妻 | 人妻无码αv中文字幕久久琪琪布| 国内精品久久久久影院优| 久久精品二区| 国内精品伊人久久久久AV影院| 18禁黄久久久AAA片| 爱做久久久久久| 国产精品久久99| 久久夜色精品国产亚洲| 久久亚洲中文字幕精品一区四| 亚洲女久久久噜噜噜熟女| 久久精品综合一区二区三区| 久久99精品国产99久久6男男| 久久免费的精品国产V∧| 综合久久久久久中文字幕亚洲国产国产综合一区首 | 日本免费久久久久久久网站| 久久精品水蜜桃av综合天堂| 亚洲精品无码久久久影院相关影片| 久久天天躁狠狠躁夜夜2020| 国产亚洲色婷婷久久99精品91| 一本色道久久88加勒比—综合| 国产精品岛国久久久久| 久久免费高清视频| 久久婷婷综合中文字幕| 国产精品va久久久久久久|