# Patent

## Title
Link spam detection using smooth classification function

## Abstract

A collection of web pages is considered as a directed graph in which the pages themselves are nodes and the hyperlinks between the pages are directed edges in the graph. A trusted entity identifies training examples for spam pages and normal pages. A random walk is conducted through the directed graph that includes the collection of web pages and the stationary probabilities, and transitional probabilities, among the nodes in the directed graph are obtained. A classifier training component estimates a classification function that changes slowly on densely connected subgraphs within the directed graph. The classification function assigns a value to each of the nodes in the directed graph and identifies them as spam or normal pages based upon whether the value meets a given function threshold value.

## Background
The present application is based on and claims the benefit of U.S. provisional patent application Ser. No. 60/927,647, filed May 4, 2007, the content of which is hereby incorporated by reference in its entirety.


BACKGROUND
Web search engines are currently in wide use, and are used to return a ranked list of web sites in response to a search query input by a user. It can be very valuable to have a web page returned high in the ranked list of web pages for a wide variety of different queries. This may increase the likelihood that a user will view a given web page.
Therefore, in order to increase web traffic to a given site, the authors of certain cites have tried to artificially manipulate the ranked list returned by search engines such that the web sites authored by those authors are ranked higher then they would normally be ranked. The particular manipulation techniques used by such authors depends on how a given web search engine ranks the pages for a given query. Any of the different manipulation techniques used by such authors are referred to as “spamming” techniques.
Some search engines use link analysis algorithms in order to generate the ranked list of web pages returned in response to a query. In general, link analysis algorithms identify an importance of a given web page, based upon the number of links that point to that web page. It is assumed that related web pages (those that have related content) have links to one another. Therefore, the more links that point to a web page, the more important the web page may be regarded by the search engine.
In order to manipulate this type of search engine, web spammers (those employing spamming techniques) sometimes attempt to create a large number of links to their web pages by having unrelated web pages (web page with unrelated content) linked to their web pages. This can be done using automated techniques to post links to their web sites onto other web pages, or simply by creating a large number of their own web pages and web sites, and then placing links in those web pages and web sites to all the other web pages and web sites which they created. This increases the number of links to any given web page or web site created by the author, regardless of whether it has related content. Similarly, some web sites reciprocally exchange links. When two unrelated web sites exchange links, at least one, and possibly both, of them are very likely to be spam (web sites that receive the benefit of spamming techniques).
It can be seen that spamming techniques can produce spam that misleads a search engine into returning low quality, or even entirely irrelevant, information to a user in response to a query. Therefore, a number of techniques have been developed in order to identify spam so that it can be removed from the ranked search results returned by a search engine. For instance, human experts can generally identify web spam in a very effective manner. However, it is quite easy for a spammer to create a large number of spam pages and to manipulate their link structure. It is thus impractical to detect web spam using only human judges. Therefore, some automatic approaches have been developed to identifying spam. One category of such approaches is referred to as a supervised approach in which some known examples of spam are provided to the system, and the system learns to recognize spam from those examples.
One such technique builds a ranking measure for web pages modeled on a user randomly following hyperlinks through the web pages. This ranking measure is well known as PageRank used by the Google search engine. At each web page, the modeled user either selects an outlink uniformly at random to follow with a certain probability, or jumps to a new web page selected from the whole web uniformly at random with the remaining probability. The stationary probability of a web page in this “random walk” is regarded as the ranking score of the web page. The basic assumption behind such a technique is that a hyperlink from one page to another is a recommendation of the second page by the author of the first page. If this assumption is recursively applied, then a web page is considered to be important if many important web pages point to it.
By using random jumps to uniformly selected pages, this system accommodates the problem that some high quality pages have no out links, although they are pointed to by many other web pages.
This concept of random jumps has also been adopted, in another way, to address the problem of web spam. Basically, the random user described above is allowed to jump to a set of pages (seed pages) which have been judged as being high quality, normal pages, by human experts. Assuming this choice for the random jumps, the stationary probability of a web page is regarded as its trust score, and a web page with a trust score smaller than a given threshold value is considered to be spam.
This type of system can also be understood as follows: initially, only the selected good seed pages have trust scores equal to one, and the trust scores of other web pages are zero. Each seed page then iteratively propagates its trust score to its neighbors, and its neighbors further propagate their received scores to their neighbors. The underlying assumption in this algorithm is that web pages of high quality seldom point to spam pages.
A counterpart to this algorithm allows the random web user to either select an inlink uniformly at random to follow, in reverse, with a certain probability, or jump to a new web page randomly selected from a web page set which has been judged as spam by human experts with the remaining probability. The stationary probability of a web page is, in this system, referred to as its antitrust rank, or antitrust score. A web page will be classified as spam if its score is larger than a chosen threshold value. In terms of the propagation understanding, the scores in this system are propagated in the reverse direction along the inlinks. The basic underlying assumption of this type of system is that a web page pointing to spam pages is likely to be spam, itself.
Another system is referred to as a functional ranking system. It considers a general ranking function that depends on incoming paths of various lengths weighted by some chosen damping function that decreases with distance. In other words, links from pages that are a greater distance from the subject web page are weighted by weight that is damped less than links from closer web pages. That is, spam pages may gain an artificially high score under a system that simply ranks the pages based on the number of links to it, because a spam page may be formed by using a spamming technique to have many incoming links from its immediate neighbor pages. However, spam pages of this type can be demoted using this system by choosing a damping function that ignores the direct contribution of links from pages directly adjacent the given page, and only valuing links that start at least one link away from the subject page.
Yet another technology to be considered is general machine learning technology. In this technology, features must be selected that are useful in detecting spam, and each web page is then represented as a vector having each element described by one type of spam feature. The features can be the number of inlinks, the number of outlinks, scores under any of the above-mentioned algorithms, etc. Then, a classifier is chosen, such as a neural network, a decision tree, a support vector machine (SVM), etc., and it is trained with a set of examples of normal and spam web pages which have been judged by human experts. The trained classifier is then used to predict a given web page as spam or not spam (i.e., as spam or a content page). One difficulty with this methodology is that the efficiency of a spam feature is generally validated only on the web pages which are not sampled from the entire web uniformly at random, but instead from large websites and highly ranked web pages. Consequently, the trained classifier is biased to those selected pages, and it does not generalize well to the entire web.
The discussion above is merely provided for general background information and is not intended to be used as an aid in determining the scope of the claimed subject matter.
SUMMARY
A collection of web pages is considered as a directed graph in which the pages themselves are nodes and the hyperlinks between the pages are directed edges in the graph. A trusted entity identifies training examples for spam pages and normal pages. A random walk is conducted through the directed graph. A classifier built on the random walk estimates a classification function that changes slowly on densely connected subgraphs within the directed graph. The classification function assigns a value to each of the nodes in the directed graph and identifies them as spam or normal pages based upon whether the value meets a given function threshold value.
This Summary is provided to introduce a selection of concepts in a simplified form that are further described below in the Detailed Description. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used as an aid in determining the scope of the claimed subject matter. The claimed subject matter is not limited to implementations that solve any or all disadvantages noted in the background.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram of one illustrative embodiment of a link spam detection system.
FIG. 2 is a flow diagram illustrating one embodiment of the overall operation of the system shown in FIG. 1.
FIG. 3 is a flow diagram illustrating one embodiment of obtaining stationary and transition probabilities for a directed graph.
FIG. 4 is one embodiment of a strongly connected directed web graph.
FIG. 5 illustrates one embodiment of the web graph shown in FIG. 4 after classification.
FIG. 6 is a block diagram of one embodiment of an illustrative computing environment.



DETAILED DESCRIPTION
Link spam detection in the present specification is discussed in terms of a machine learning problem of classification on directed graphs. FIG. 1 is a block diagram of one illustrative embodiment of a link spam detection system 100. System 100 shows trusted entity 102, directed graph 104, random walk component 106, classifier training component 108, spam classifier 110, and spam detection system 112.
In system 100, a collection of web pages 114 is also shown. The collection of web pages 114 is considered a directed graph, in that the web pages themselves in collection 114 are nodes in the graph while hyperlinks between those web pages are directed edges in the directed graph. Of course, it will be appreciated that the present system can be applied at the domain/host level as well, where the domains/hosts are nodes in the graph and hyperlinks among the web pages in the domains/hosts are the directed edges. For purposes of the present discussion, however, and by way of example only, reference will be made to the web pages in collection 113 as the nodes and hyperlinks between those pages as the directed edges.
FIG. 2 is a flow diagram illustrating one illustrative embodiment of the operation of system 100 shown in FIG. 1. FIGS. 1 and 2 will now be described in conjunction with one another.
It will first be noted that, if web page collection 114 does not form a strongly connected graph, then it is first decomposed into strongly connected components, and the present process proceeds with respect to each of the strongly directed components. The precise definition of what makes a strongly connected graph, or strongly connected component, is set out below. Briefly, however, a graph can be considered strongly connected if each vertex (or node) in the graph is connected to every other vertex (or node) in the graph by some path of directed edges. Decomposing the directed graph into strongly connected components is illustrated by block 150 in FIG. 2 and is shown in phantom indicating that it is only performed, if necessary.
The web page collection 114 is also provided to trusted entity 102, such as a human expert, in identifying link spam. Trusted entity 102 then identifies some examples of spam web pages in web page collection 114 as spam training examples 116. Trusted entity 102 also identifies good web pages (or normal web pages) in web page collection 114 as normal training examples 118. Obtaining these examples is indicated by block 152 in FIG. 2.
Random walk component 106 then receives a definition of a random walk on directed graph 104 (or each strongly connected component in directed graph 104), the random walk being defined by translation probabilities (set out below in Eqs. 20-22). Receiving this definition is indicated by block 153 in FIG. 1. Based on the defined random walk, component 106 obtains stationary probabilities associated with each node in directed graph 104. The stationary probabilities are indicated by block 120 in FIG. 1. Obtaining these probabilities for directed graph 104 is indicated by block 154 in FIG. 2, and the stationary probabilities are obtained by conducting the defined random walk through the directed graph 104. This is discussed in greater detail below with respect to FIG. 3.
In any case, once examples 116 and 118 and probabilities 120 and 122 are obtained, classifier training component 108 trains a classifier that can be used in link spam detection. The classifier is shown in FIG. 1 as spam classifier 110. In one embodiment, training a classifier is performed by generating a smooth classification function based on the probabilities 120 and 122 over the detected graph 104. In generating the classification function, the values of the classification function are forced to be close to known values for examples 116 and 118. In other words, the values of the classification function are forced to be close to the values that indicate spam and normal pages at the nodes in the graph that are actually known to be spam and normal pages as identified by the trusted entity 102. For example, assume that the value of −1 indicates that the node is spam, while the value of 1 indicates that the node is a normal content page. Then, the classification function is forced to be at least close to the values of 1 or −1 at the pages known to be normal pages and spam pages, respectively. This is indicated by block 156 in FIG. 2.
The closeness between the classification function and the known values can be measured in a variety of different ways. For example, the closeness can be measured using least square loss, hinge loss, precision/recall measure, the F1-score, the ROC score, or the AUC score, as examples.
In accordance with one embodiment, the classification function is not only close to known values at known nodes, but it is relatively smooth in that it changes relatively slowly on densely connected subgraphs. In other words, the nodes that reside close to one another on the subgraph may likely have values which are relatively close to one another. However, if they are known to be one spam node and one normal node, respectively, then the classification function changes by a large amount between those nodes, but this lack of smoothness is penalized in the chosen cost function that is optimized.
This can provide significant advantages over prior systems. In prior systems, for instance, those pages closely related to spam pages were deemed as spam while all other pages were deemed as normal pages. In another prior system, those pages close to normal pages were deemed normal pages, while all other pages were deemed spam. The present system includes information related to both normal pages and spam pages in classifying a given page under consideration as content or spam. Also, because pages that are relatively close to one another on the directed graph are assumed to be the same type (pages close to a known spam page are likely to be spam pages, while pages close to a known normal page are likely to be normal pages) by making the function smooth and relatively slow changing pages in the directed subgraph that are close to a known normal content page will have classification function values that more likely indicate it to be a normal content page. Similarly, those pages in the directed subgraph that are close to a spam page will have classification function values that are likely to indicate that it will be a spam page. The classification function value can change abruptly, if necessary. Again, however, this is penalized.
In any case, the spam classifier is then used to assign values to all unlabeled nodes in the directed graph 104. A threshold can be set, and those pages that meet the classification threshold may be deemed to be normal pages, while those nodes having a value that does not meet the threshold value may be deemed to be spam pages. In one embodiment, simply the sign of the value calculated using the classification function is used to determine whether the associated node is spam or content. This effectively sets the classification function threshold value at 0. It may, however, be desirable to set the value at a level other than 0. For instance, if it is more desirable in a given application to error on the side of classifying spam pages as normal pages, then the threshold may be set below 0. On the other hand, if a given application deems it more desirable to error on the side of classifying normal pages, as spam pages, then the threshold value can be set above 0, etc. Using the classification function embodied in spam classifier 110 to perform spam detection is indicated by block 158 in FIG. 2.
FIG. 3 is a flow diagram showing one illustrative embodiment of random walk component 106 in obtaining the stationary probabilities 120. In one embodiment, random walk component 106 simply selects, at random, a starting node in directed graph 104. This is indicated by block 180 in FIG. 3. Component 106 then randomly follows links in graph 104 starting at the selected starting node. This is indicated by block 182. It will be noted that, at each step, component 106 can follow inlinks or outlinks from the given web page uniformly at random. If the outlinks are followed, component 106 simply follows links from the current page to another page to which it is linked through an outlink. However, if inlinks are used, then component 106 travels backward along the links that link to the current page, to the page at which the inlink originates. For the present description, following the outlink will be used, although either inlinks or outlinks could be used, as desired.
Component 106 continues to follow the links, uniformly at random, for a sufficient amount of time. This is indicated by block 184. The amount of time will depend upon the size of the collection of web pages 114. As component 106 is performing this random walk, it calculates the stationary probability distribution for the various nodes in graph 104. The “transition probabilities” are the probabilities of transitioning from any given node on graph 104 to another node. The “stationary probability distribution” assumes that component 106 starts from a randomly chosen node in graph 104, and jumps to an adjacent node by choosing an outlink. Assume for the sake of example that this is repeated infinitely many times for the various nodes in graph 104. Then, if graph 104 is connected (that is, using such a random walk, any point can be reached from any other point), then the fraction of time component 106 spends at a given node converges to a fixed number (where the corresponding numbers for all nodes sum to 1), and that fixed number is actually independent of the choice of starting nodes. In other words, the stationary probability distributions are the probabilities of being in any given node on directed graph 104.
Component 106 can use any given metric to determine whether it has performed the random walk sufficiently long enough. For example, where the stationary probabilities do not change by a given amount (that is they are changing very little or very slowly with each given iteration of the jump) then, component 106 may deem that it has performed the random walk long enough. In any case, once the random walk has been performed for sufficiently long time, component 106 calculates the final stationary probabilities 120 that are output to classifier training component 108. This is indicated by block 186 in FIG. 3.
It may seem at first that performing the classification over a large directed graph may take an inordinately large amount of time. It has been found that it can be done quite quickly, using a relatively small number of training examples. For instance, in a directed graph having 20 million web pages connected by directed edges (links) with 10,000 examples of spam web pages and 20,000 examples of content web pages, the classification can be performed in several minutes.
Having thus described transductive detection of spam pages in an intuitive sense, it will now be described in a more formal way. First, a discussion of some specific items of notation will be made.
Let G=(V,E) denote a directed graph, where V is the set of vertices, and E the set of edges. For a given edge eεE; denote the initial vertex of e by e−, and the terminal vertex of e by e+. Also denote by (u, v) an edge from the vertex u to the vertex v. It is clear that an undirected graph can be regarded as a directed graph with each edge being double oriented. A graph G is weighted if it is associated with a function w: E→R+ which assigns a positive number w(e) to each edge e of G: Let G=(V,E,w) denote a weighted directed graph. The function w is called the weight function of G: The in-degree d− and the out-degree d+ of a vertex v ε V are respectively defined as:




  
    
      
        
          
            
              d
              -
            
            ⁡
            
              (
              v
              )
            
          
          =
          
            
              ∑
              
                {
                
                  
                    e
                    ❘
                    
                      e
                      +
                    
                  
                  =
                  v
                
                }
              
            
            ⁢
            
              w
              ⁡
              
                (
                e
                )
              
            
          
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              
                d
                +
              
              ⁡
              
                (
                v
                )
              
            
          
          =
          
            
              ∑
              
                {
                
                  
                    e
                    ❘
                    
                      e
                      -
                    
                  
                  =
                  v
                
                }
              
            
            ⁢
            
              w
              ⁡
              
                (
                e
                )
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        1
      
    
  




A path is a tuple of vertices (v1, v2, . . . , vp) with the property that (vi, vi+1) ε E for 1≦i≦p−1. A directed graph is strongly connected when for every pair of vertices u and v there is a path in which v1=u and vp=v. For a strongly connected graph, there is an integer k≧1 and a unique partition V=V0∪V1∪ . . . ∪Vk−1 such that for all 0≦r≦k−1 each edge (u, v) ε E with u ε Vr has v ε Vr+1, where Vk=V0; and k is maximal, that is, there is no other such partition V=V0′∪ . . . ∪Vk′−1′ with k′>k.
When k=1, the graph is aperiodic; otherwise the graph is periodic.
For a given weighted directed graph, there is a natural random walk on the graph with the transition probability function p:V×V→ defined by:




  
    
      
        
          p
          ⁡
          
            (
            
              u
              ,
              v
            
            )
          
        
        =
        
          
            w
            ⁡
            
              (
              
                u
                ,
                v
              
              )
            
          
          
            
              d
              +
            
            ⁡
            
              (
              u
              )
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        2
      
    
  




for all (u, v) ε E, and 0 otherwise. If the graph is strongly connected, there is a unique function π: V→ which satisfies:




  
    
      
        
          
            
              ∑
              
                u
                ∈
                V
              
            
            ⁢
            
              
                π
                ⁡
                
                  (
                  u
                  )
                
              
              ⁢
              
                p
                ⁡
                
                  (
                  
                    u
                    ,
                    v
                  
                  )
                
              
            
          
          =
          
            π
            ⁡
            
              (
              v
              )
            
          
        
        ,
        
          
            and
            ⁢
            
              
            
            ⁢
            
              
                ∑
                v
              
              ⁢
              
                π
                ⁡
                
                  (
                  v
                  )
                
              
            
          
          =
          1
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        3
      
    
  




The first equation in Equation 3 is called the balance equation, and π is called the Perron vector. For a general directed graph, there is no closed form solution for π. If the graph is both strongly connected and aperiodic, the random walk defined by Eq. 2 converges to the Perron vector π. Unless stated otherwise, the directed graphs considered are always assumed to be strongly connected. One embodiment of a strongly connected graph is shown in FIG. 4. The nodes (or vertices) are labeled 1-9 while the edges are shown as arrows.
Now, a number of discrete operators on directed graphs are defined. The operators are discrete analogs of the corresponding differential operators on Riemannian manifolds. As discussed below, the discrete operators are then used to develop a discrete analog of classical regularization theory. Consequently, as in other regularization based machine learning algorithms in vectorial spaces (for instance, support vector machines (SVMs)) the present classification algorithm for directed graphs is derived from the discrete regularization.
In any case, let F(V) denote the set of all real-valued functions on V; and F(E) the set of all real-valued functions on E. The function set F(V μ)can be regarded as a Hilbert space H(V μ) with the inner product defined by:




  
    
      
        
          
            〈
            
              φ
              ,
              φ
            
            〉
          
          
            H
            ⁡
            
              (
              V
              )
            
          
        
        =
        
          
            ∑
            
              v
              ∈
              V
            
          
          ⁢
          
            
              φ
              ⁡
              
                (
                v
                )
              
            
            ⁢
            
              φ
              ⁡
              
                (
                v
                )
              
            
            ⁢
            
              π
              ⁡
              
                (
                v
                )
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        4
      
    
  




where φ,φεF(V). Let c(e)=π(e−)p(e). The number c(e) is called the ergodic flow on e. It is easy to check that the ergodic flow is a circulation, that is:




  
    
      
        
          
            
              ∑
              
                {
                
                  
                    e
                    ❘
                    
                      e
                      -
                    
                  
                  =
                  v
                
                }
              
            
            ⁢
            
              c
              ⁡
              
                (
                e
                )
              
            
          
          =
          
            
              ∑
              
                {
                
                  
                    e
                    ❘
                    
                      e
                      +
                    
                  
                  =
                  v
                
                }
              
            
            ⁢
            
              c
              ⁡
              
                (
                e
                )
              
            
          
        
        ,
        
          ∀
          
            v
            ∈
            V
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        5
      
    
  




A Hilbert space H(E) over F(E) can be constructed with the inner product defined by:




  
    
      
        
          
            〈
            
              θ
              ,
              ψ
            
            〉
          
          
            H
            ⁡
            
              (
              E
              )
            
          
        
        =
        
          
            ∑
            
              e
              ∈
              E
            
          
          ⁢
          
            
              θ
              ⁡
              
                (
                e
                )
              
            
            ⁢
            
              ψ
              ⁡
              
                (
                e
                )
              
            
            ⁢
            
              c
              ⁡
              
                (
                e
                )
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        6
      
    
  




Where ψεF(E).
The discrete gradient ∇:H(V)→H(E) is defined as an operator:

(∇φ)(e):=φ(e+)−φ(e−), ∀φεH(V)  Eq. 7

For simplicity, (∇φ)(e) is also denoted as ∇eφ. For gaining an intuition of this definition, one may imagine a set of buckets, and some of them are connected by tubes. Assume a tube e which connects buckets e− and e+, and the quantities of fluid in buckets e− and e+ to be φ(e−) and φ(e+). Then the flow through the tube should be proportional to the pressure difference and hence to φ(e+)−φ(e−). When the fluid distributes itself uniformly among buckets, that is, φ is constant, the pressure differences will disappear and consequently there will be no flow in tubes any more, that is, ∇φ vanishes everywhere.
As in the continuous case, the discrete divergence div:H(E) →H(V) can be defined as the dual of −∇ that is:

∇φ,ψH(E)=φ,−divψH(V)  Eq. 8

where φεH(V),ψεH(E). By a straightforward computation, the following is obtained:




  
    
      
        
          
            (
            
              div
              ⁢
              
                
              
              ⁢
              ψ
            
            )
          
          ⁢
          
            (
            v
            )
          
        
        =
        
          
            1
            
              π
              ⁡
              
                (
                v
                )
              
            
          
          ⁢
          
            (
            
              
                
                  ∑
                  
                    {
                    
                      
                        e
                        ❘
                        
                          e
                          -
                        
                      
                      =
                      v
                    
                    }
                  
                
                ⁢
                
                  
                    c
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                  ⁢
                  
                    ψ
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                
              
              -
              
                
                  ∑
                  
                    {
                    
                      
                        e
                        ❘
                        
                          e
                          +
                        
                      
                      =
                      v
                    
                    }
                  
                
                ⁢
                
                  
                    c
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                  ⁢
                  
                    ψ
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                
              
            
            )
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        9
      
    
  




By following the above fluid model, the divergence measures the net flows at buckets. Now the concept of circulation can be generalized in terms of divergence. A function ψεH(E) is called a circulation if and only if div ψ=0.
The discrete Laplacian Δ: H(V)→H(V)is defined by:




  
    
      
        
          Δ
          :=
          -
        
        ⁣
        
          
            1
            2
          
          ⁢
          
            div
            ·
            ∇
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        10
      
    
  




Compared with its counterpart in the continuous case, the additional factor in Eq. 10 is due to edges being oriented. From Eq. 10:




  
    
      
        
          
            〈
            
              Δφ
              ,
              φ
            
            〉
          
          
            H
            ⁡
            
              (
              V
              )
            
          
        
        =
        
          
            
              1
              2
            
            ⁢
            
              
                〈
                
                  
                    ∇
                    φ
                  
                  ,
                  
                    ∇
                    φ
                  
                
                〉
              
              
                H
                ⁡
                
                  (
                  E
                  )
                
              
            
          
          =
          
            
              〈
              
                φ
                ,
                Δφ
              
              〉
            
            
              H
              ⁡
              
                (
                V
                )
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        11
      
    
  




Note that the first equation in Eq. 11 is a discrete analog of Green's formula. In addition, Eq. 11 implies that Δ is self-adjoint. In particular, when φ=φ, then:




  
    
      
        
          
            〈
            
              Δφ
              ,
              φ
            
            〉
          
          
            H
            ⁡
            
              (
              V
              )
            
          
        
        =
        
          
            
              1
              2
            
            ⁢
            
              
                〈
                
                  
                    ∇
                    φ
                  
                  ,
                  
                    ∇
                    φ
                  
                
                〉
              
              
                H
                ⁡
                
                  (
                  E
                  )
                
              
            
          
          =
          
            
              1
              2
            
            ⁢
            
              
                
                
                  ∇
                  φ
                
                
              
              
                H
                ⁡
                
                  (
                  E
                  )
                
              
              2
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        12
      
    
  




which implies that Δ is positive semi-definite. By substituting Eqs. 7 and 9 into Eq. 10:




  
    
      
        
          
            (
            Δφ
            )
          
          ⁢
          
            (
            v
            )
          
        
        -
        
          
            1
            
              2
              ⁢
              
                π
                ⁡
                
                  (
                  v
                  )
                
              
            
          
          ·
          
            (
            
              
                
                  ∑
                  
                    {
                    
                      
                        e
                        ❘
                        
                          e
                          +
                        
                      
                      =
                      v
                    
                    }
                  
                
                ⁢
                
                  
                    c
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                  ⁢
                  
                    φ
                    ⁡
                    
                      (
                      
                        e
                        -
                      
                      )
                    
                  
                
              
              +
              
                
                  ∑
                  
                    {
                    
                      
                        e
                        ❘
                        
                          e
                          -
                        
                      
                      =
                      v
                    
                    }
                  
                
                ⁢
                
                  
                    c
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                  ⁢
                  
                    φ
                    ⁡
                    
                      (
                      
                        e
                        +
                      
                      )
                    
                  
                
              
            
            )
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        13
      
    
  




when the graph is undirected, that is, each edge being double oriented, Eq. 13 reduces to:




  
    
      
        
          
            (
            Δφ
            )
          
          ⁢
          
            (
            v
            )
          
        
        =
        
          
            φ
            ⁡
            
              (
              v
              )
            
          
          -
          
            
              1
              
                d
                ⁡
                
                  (
                  v
                  )
                
              
            
            ⁢
            
              
                ∑
                
                  u
                  -
                  v
                
              
              ⁢
              
                
                  w
                  ⁡
                  
                    (
                    
                      u
                      ,
                      v
                    
                    )
                  
                
                ⁢
                
                  φ
                  ⁡
                  
                    (
                    v
                    )
                  
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        14
      
    
  




Eq. 14 has been widely used to define the Laplacian for an undirected graph. Now, define a family of functions {δv}vεV with δv(u)=Iu=v, which is clearly a basis of H(V). The matrix form of Δ with respect to this basis has the following components:




  
    
      
        
          
            Δ
            am
          
          ⁡
          
            (
            
              u
              ,
              v
            
            )
          
        
        =
        
          {
          
            
              
                
                  -
                  
                    ⁣
                    
                      
                        
                          c
                          ⁡
                          
                            (
                            
                              u
                              ,
                              v
                            
                            )
                          
                        
                        +
                        
                          c
                          ⁡
                          
                            (
                            
                              v
                              ,
                              u
                            
                            )
                          
                        
                      
                      
                        2
                        ⁢
                        
                          π
                          ⁡
                          
                            (
                            u
                            )
                          
                        
                      
                    
                  
                
              
              
                
                  
                    u
                    ≠
                    v
                  
                  ,
                
              
            
            
              
                1
              
              
                
                  u
                  =
                  v
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        15
      
    
  




This matrix is not symmetric. However, if another basis {π−1/2(v)δv}vεV is chosen, then Δ can be represented as a symmetric matrix:




  
    
      
        
          
            Δ
            sm
          
          ⁡
          
            (
            
              u
              ,
              v
            
            )
          
        
        =
        
          {
          
            
              
                
                  -
                  
                    ⁣
                    
                      
                        
                          c
                          ⁡
                          
                            (
                            
                              u
                              ,
                              v
                            
                            )
                          
                        
                        +
                        
                          c
                          ⁡
                          
                            (
                            
                              v
                              ,
                              u
                            
                            )
                          
                        
                      
                      
                        2
                        ⁢
                        
                          
                            
                              π
                              ⁡
                              
                                (
                                u
                                )
                              
                            
                            ⁢
                            
                              π
                              ⁡
                              
                                (
                                v
                                )
                              
                            
                          
                        
                      
                    
                  
                
              
              
                
                  
                    u
                    ≠
                    v
                  
                  ,
                
              
            
            
              
                1
              
              
                
                  u
                  =
                  v
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        16
      
    
  




This matrix has been used to define Laplacian for directed graphs.
Now, learning on directed graphs using the above analysis is discussed. Given a directed graph G=(V,E,w), and a discrete label set L={−1, 1}, the vertices in a subset S⊂V have labels in L. The task is to predict the labels of those unclassified vertices in Sc, the complement of S. The present link spam detection problem can be cast into classification on a directed graph. For instance, FIG. 5 shows the vertices (or nodes) in the graph shown in FIG. 4 classified as spam or normal pages. The solid nodes are classified as normal pages while those shown in phantom are spam pages.
Define a function y with y(v)=1 or −1 if vεS, and 0 if vεSc. For classifying those unclassified vertices in Sc, define a discrete regularization:




  
    
      
        
          argmin
          
            φ
            ∈
            
              H
              ⁡
              
                (
                V
                )
              
            
          
        
        ⁢
        
          {
          
            
              
                
                
                  ∇
                  φ
                
                
              
              
                H
                ⁡
                
                  (
                  E
                  )
                
              
              2
            
            +
            
              C
              ⁢
              
                
                  
                  
                    φ
                    -
                    y
                  
                  
                
                
                  H
                  ⁡
                  
                    (
                    V
                    )
                  
                
                2
              
            
          
          }
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        17
      
    
  




where C>0 is the regularization parameter. In the objective function, the first term forces the classification function to be relatively smooth, and perhaps as smooth as possible and the second term forces the classification function to fit the given labels as well as possible.
When choosing the basis {δv}vεV, Eq. 17 can be written as:




  
    
      
        
          argmin
          
            φ
            ∈
            
              H
              ⁡
              
                (
                V
                )
              
            
          
        
        ⁢
        
          {
          
            
              
                ∑
                
                  e
                  ∈
                  E
                
              
              ⁢
              
                
                  π
                  ⁡
                  
                    (
                    
                      e
                      -
                    
                    )
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    e
                    )
                  
                
                ⁢
                
                  
                    (
                    
                      
                        φ
                        ⁡
                        
                          (
                          
                            e
                            +
                          
                          )
                        
                      
                      -
                      
                        φ
                        ⁡
                        
                          (
                          
                            e
                            -
                          
                          )
                        
                      
                    
                    )
                  
                  2
                
              
            
            +
            
              C
              ⁢
              
                
                  ∑
                  
                    v
                    ∈
                    V
                  
                
                ⁢
                
                  
                    π
                    ⁡
                    
                      (
                      v
                      )
                    
                  
                  ⁢
                  
                    
                      (
                      
                        
                          φ
                          ⁡
                          
                            (
                            v
                            )
                          
                        
                        -
                        
                          y
                          ⁡
                          
                            (
                            v
                            )
                          
                        
                      
                      )
                    
                    2
                  
                
              
            
          
          }
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        18
      
    
  




Again, the first term makes the function relatively smooth over all nodes while the second term forces the function to fit the labeled nodes to a desired closeness. If each function in H(V) is scaled with a factor π−1/2 (in other words, choose another basis {π−1/2(v)δv}vεV,), then Eq. 18 will be transformed into:




  
    
      
        
          argmin
          
            φ
            ∈
            
              H
              ⁡
              
                (
                V
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        
          {
          
            
              
                ∑
                
                  e
                  ∈
                  E
                
              
              ⁢
              
                
                  π
                  ⁡
                  
                    (
                    
                      e
                      -
                    
                    )
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    e
                    )
                  
                
                ⁢
                
                  
                    (
                    
                      
                        
                          φ
                          ⁡
                          
                            (
                            
                              e
                              +
                            
                            )
                          
                        
                        
                          
                            π
                            ⁡
                            
                              (
                              
                                e
                                +
                              
                              )
                            
                          
                        
                      
                      -
                      
                        
                          φ
                          ⁡
                          
                            (
                            
                              e
                              -
                            
                            )
                          
                        
                        
                          
                            π
                            ⁡
                            
                              (
                              
                                e
                                
                                  -
                                  1
                                
                              
                              )
                            
                          
                        
                      
                    
                    )
                  
                  2
                
              
            
            +
            
              C
              ⁢
              
                
                  ∑
                  
                    v
                    ∈
                    V
                  
                
                ⁢
                
                  
                    (
                    
                      
                        φ
                        ⁡
                        
                          (
                          v
                          )
                        
                      
                      -
                      
                        y
                        ⁡
                        
                          (
                          v
                          )
                        
                      
                    
                    )
                  
                  2
                
              
            
          
          }
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        19
      
    
  




However, it can be seen that Eq. 18 is much more natural than Eq. 19.
A random walk over a given directed graph can be defined in many different ways. Three exemplary types of random walk used in spam detection are:
1. Following outlinks uniformly at random.

    
    
        Formally, define a random walk with:
    
    






  
    
      
        
          p
          ⁡
          
            (
            
              u
              ,
              v
            
            )
          
        
        =
        
          
            w
            ⁡
            
              (
              
                u
                ,
                v
              
              )
            
          
          
            
              d
              +
            
            ⁡
            
              (
              u
              )
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        20
      
    
  




This is the one discussed above with respect to FIG. 3.
2. Following links uniformly at random regardless of directionality. Formally, define a random walk with:




  
    
      
        
          p
          ⁡
          
            (
            
              u
              ,
              v
            
            )
          
        
        =
        
          
            
              w
              ⁡
              
                (
                
                  u
                  ,
                  v
                
                )
              
            
            +
            
              w
              ⁡
              
                (
                
                  v
                  ,
                  u
                
                )
              
            
          
          
            
              
                d
                +
              
              ⁡
              
                (
                u
                )
              
            
            +
            
              
                d
                -
              
              ⁡
              
                (
                u
                )
              
            
          
        
      
    
    
      
        Eq
        ,
        
          
        
        ⁢
        21
      
    
  




3. Following inlinks uniformly at random. Formally, define a random walk with:




  
    
      
        
          p
          ⁡
          
            (
            
              u
              ,
              v
            
            )
          
        
        =
        
          
            w
            ⁡
            
              (
              
                v
                ,
                u
              
              )
            
          
          
            
              d
              _
            
            ⁡
            
              (
              u
              )
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        22
      
    
  




Other choices of random walks can be used as well.
Assigning values to the nodes in directed graph 104 basically requires selection of a random walk definition (transition probabilities) and solving Eq. 18 above for each of the nodes. This is set out above with respect to FIG. 2. Solving for Eq. 18 is set out more formally in pseudocode in Table 1 below for the random walk that inversely follows the links. To solve the optimization problem in Eq. 18, differentiate the objective function with respect to φ and then obtain:

Δamφ+C(φ−y)=0  Eq. 23

where the first term on the left hand side is derived from Eq. 11 via the differential rule on inner products. The above equation can be written as:

(CI+Δam)φ=Cy  Eq. 24

where I is the identity matrix. This linear system has the closed-form solution:

φ=C(CI+Δam)−1y  Eq. 25

although it may be more efficient to solve the linear system directly, rather than computing the inverse.
In the algorithm in Table 1 below, a parameter αε]0,1[ is used instead of Cε]0,∞[. The relationship between αand C can be expressed as:




  
    
      
        α
        =
        
          1
          
            1
            +
            C
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        26
      
    
  




In the last step in Table 1, the classification is based on the sign of the function value on each vertex. As mentioned above with respect to FIG. 2, this is equivalent to setting the classification threshold to 0.







TABLE 1





TRANSDUCTIVE LINK SPAM DETECTION















Given a web graph G = (V,E), some web pages S⊂V have been


manually labeled as content or spam. The graph is strongly connected.


Otherwise, it is decomposed into strongly connected components. The


remaining unclassified web pages in V may be classified as follows:








1.
Define a random walk which chooses an inlink uniformly at random



to follow. Formally, this random walk has the transition probabilities:


 







  
    
      p
      ⁡
      
        (
        
          u
          ,
          v
        
        )
      
    
    =
    
      
        w
        ⁡
        
          (
          
            v
            ,
            u
          
          )
        
      
      
        
          d
          -
        
        ⁡
        
          (
          u
          )
        
      
    
  
  ,






 



for any u, v in V. Let π denote the vector which satisfies:


 







  
    
      ∑
      
        u
        ∈
        V
      
    
    ⁢
    
      
    
    ⁢
    
      
        π
        ⁡
        
          (
          u
          )
        
      
      ⁢
      
        p
        ⁡
        
          (
          
            u
            ,
            v
          
          )
        
      
    
  
  =
  
    
      π
      ⁡
      
        (
        v
        )
      
    
    .
  






 


2.
Denote by P the matrix with the elements p(u,v), and Π the diagonal



matrix with the diagonal elements being stationary probabilities π (u)



and zeros everywhere else. Form the matrix:


 







  L
  =
  
    Π
    -
    
      α
      ⁢
      
        
          
            Π
            ⁢
            
              
            
            ⁢
            P
          
          +
          
            
              P
              T
            
            ⁢
            Π
          
        
        2
      
    
  






 



where α is a parameter in ]0,1[.


3.
Define a function y on V with y(v) = 1 or −1 if the web page v is



labeled as content or spam, and 0 if v is unlabeled. Solve the linear



system:



Lφ = Πy,



and classify each unlabeled web page v as sign φ(v).









FIG. 6 illustrates an example of a suitable computing system environment 300 on which embodiments may be implemented. The computing system environment 300 is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality of the claimed subject matter. Neither should the computing environment 300 be interpreted as having any dependency or requirement relating to any one or combination of components illustrated in the exemplary operating environment 300.
Embodiments are operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well-known computing systems, environments, and/or configurations that may be suitable for use with various embodiments include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, telephony systems, distributed computing environments that include any of the above systems or devices, and the like.
Embodiments may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Some embodiments are designed to be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules are located in both local and remote computer storage media including memory storage devices.
With reference to FIG. 6, an exemplary system for implementing some embodiments includes a general-purpose computing device in the form of a computer 310. Components of computer 310 may include, but are not limited to, a processing unit 320, a system memory 330, and a system bus 321 that couples various system components including the system memory to the processing unit 320. The system bus 321 may be any of several types of bus structures including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus also known as Mezzanine bus.
Computer 310 typically includes a variety of computer readable media. Computer readable media can be any available media that can be accessed by computer 310 and includes both volatile and nonvolatile media, removable and non-removable media. By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes both volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic. storage devices, or any other medium which can be used to store the desired information and which can be accessed by computer 310. Communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media. The term “modulated data signal” means a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal. By way of example, and not limitation, communication media includes wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above should also be included within the scope of computer readable media.
The system memory 330 includes computer storage media in the form of volatile and/or nonvolatile memory such as read only memory (ROM) 331 and random access memory (RAM) 332. A basic input/output system 333 (BIOS), containing the basic routines that help to transfer information between elements within computer 310, such as during start-up, is typically stored in ROM 331. RAM 332 typically contains data and/or program modules that are immediately accessible to and/or presently being operated on by processing unit 320. By way of example, and not limitation, FIG. 6 illustrates operating system 334, application programs 335, other program modules 336, and program data 337. Any part of system 100 can be in programs 335, modules 336, or anywhere else, as desired.
The computer 310 may also include other removable/non-removable volatile/nonvolatile computer storage media. By way of example only, FIG. 6 illustrates a hard disk drive 341 that reads from or writes to non-removable, nonvolatile magnetic media, a magnetic disk drive 351 that reads from or writes to a removable, nonvolatile magnetic disk 352, and an optical disk drive 355 that reads from or writes to a removable, nonvolatile optical disk 356 such as a CD ROM or other optical media. Other removable/non-removable, volatile/nonvolatile computer storage media that can be used in the exemplary operating environment include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like. The hard disk drive 341 is typically connected to the system bus 321 through a non-removable memory interface such as interface 340, and magnetic disk drive 351 and optical disk drive 355 are typically connected to the system bus 321 by a removable memory interface, such as interface 350.
The drives and their associated computer storage media discussed above and illustrated in FIG. 6, provide storage of computer readable instructions, data structures, program modules and other data for the computer 310. In FIG. 6, for example, hard disk drive 341 is illustrated as storing operating system 344, application programs 345, other program modules 346, and program data 347. Note that these components can either be the same as or different from operating system 334, application programs 335, other program modules 336, and program data 337. Operating system 344, application programs 345, other program modules 346, and program data 347 are given different numbers here to illustrate that, at a minimum, they are different copies.
A user may enter commands and information into the computer 310 through input devices such as a keyboard 362, a microphone 363, and a pointing device 361, such as a mouse, trackball or touch pad. Other input devices (not shown) may include a joystick, game pad, satellite dish, scanner, or the like. These and other input devices are often connected to the processing unit 320 through a user input interface 360 that is coupled to the system bus, but may be connected by other interface and bus structures, such as a parallel port, game port or a universal serial bus (USB). A monitor 391 or other type of display device is also connected to the system bus 321 via an interface, such as a video interface 390. In addition to the monitor, computers may also include other peripheral output devices such as speakers 397 and printer 396, which may be connected through an output peripheral interface 395.
The computer 310 is operated in a networked environment using logical connections to one or more remote computers, such as a remote computer 380. The remote computer 380 may be a personal computer, a hand-held device, a server, a router, a network PC, a peer device or other common network node, and typically includes many or all of the elements described above relative to the computer 310. The logical connections depicted in FIG. 6 include a local area network (LAN) 371 and a wide area network (WAN) 373, but may also include other networks. Such networking environments are commonplace in offices, enterprise-wide computer networks, intranets and the Internet.
When used in a LAN networking environment, the computer 310 is connected to the LAN 371 through a network interface or adapter 370. When used in a WAN networking environment, the computer 310 typically includes a modem 372 or other means for establishing communications over the WAN 373, such as the Internet. The modem 372, which may be internal or external, may be connected to the system bus 321 via the user input interface 360, or other appropriate mechanism. In a networked environment, program modules depicted relative to the computer 310, or portions thereof, may be stored in the remote memory storage device. By way of example, and not limitation, FIG. 6 illustrates remote application programs 385 as residing on remote computer 380. It will be appreciated that the network connections shown are exemplary and other means of establishing a communications link between the computers may be used.
Although the subject matter has been described in language specific to structural features and/or methodological acts, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims.

## Claims

1. A computer implemented method of grouping a collection of web pages into a plurality of groups, comprising:
obtaining training examples of the web pages in the plurality of groups, each training example having a target function value,
wherein obtaining training examples comprises obtaining a plurality of web pages in their collection that are labeled by a trusted entity as spam training pages, and also further comprises obtaining a plurality of web pages in the collection that are labeled by the trusted entity as normal training pages;
performing a random walk over a directed graph representing the collection of web pages to obtain stationary probabilities associated with nodes in the directed graph, the nodes representing the web pages in the collection;
identifying a classification function that assigns a classification function value to each of the nodes in the directed graph based on the stationary probabilities associated with the nodes, based on classification function values assigned to other nodes in the directed graph, and based on the target function values; and
labeling nodes in the directed graph as belonging to one of the plurality of groups based on classification function values assigned to the nodes by the identified classification function.

2. The method of claim 1 wherein identifying a classification function comprises:
optimizing a cost function that assigns a penalty based on a size of a difference in classification function values between nodes in the directed graph.

3. The method of claim 2 wherein identifying a classification function comprises:
optimizing the cost function that assigns a penalty based on a difference between classification function values assigned to nodes representing the training examples and the target function values for those nodes.

4. The method of claim 3 wherein performing a random walk, comprises:
obtaining a random walk definition defining the random walk over the nodes in the directed graph.

5. The method of claim 4 wherein the random walk definition specifies transitionary probabilities associated with the directed graph.

6. The method of claim 5 wherein labeling comprises:
classifying the nodes in the directed graph as representing spam pages or normal pages.

7. The method of claim 1 and further comprising:
prior to performing the random walk, breaking the directed graph into strongly connected graph components.