Semi-Supervised Learning for Detecting Human Trafficking

TL;DR

This paper introduces a semi-supervised learning approach using Laplacian SVM variants to detect potential human trafficking advertisements online, leveraging limited labeled data and unlabeled web content.

Contribution

We develop S3VM-R, an enhanced semi-supervised learning method that incorporates exogenous information for better detection of trafficking-related ads.

Findings

S3VM-R outperforms existing semi-supervised methods in identifying relevant ads.

The approach effectively utilizes limited labeled data with abundant unlabeled data.

Law enforcement experts verified the high-interest advertisements identified by our method.

Abstract

Human trafficking is one of the most atrocious crimes and among the challenging problems facing law enforcement which demands attention of global magnitude. In this study, we leverage textual data from the website "Backpage"- used for classified advertisement- to discern potential patterns of human trafficking activities which manifest online and identify advertisements of high interest to law enforcement. Due to the lack of ground truth, we rely on a human analyst from law enforcement, for hand-labeling a small portion of the crawled data. We extend the existing Laplacian SVM and present S3VM-R, by adding a regularization term to exploit exogenous information embedded in our feature space in favor of the task at hand. We train the proposed method using labeled and unlabeled data and evaluate it on a fraction of the unlabeled data, herein referred to as unseen data, with our expert's further verification. Results from comparisons between our method and other semi-supervised and supervised approaches on the labeled data demonstrate that our learner is effective in identifying advertisements of high interest to law enforcement

Tables

Table 1. Table 1: Different features and their corresponding groups.

No.	Feature Group	Ref.
1	Advertisement Language Pattern	kennedy2012predictive ; Li:2008:IKC:1478784 ; conf/setn/KanarisKS06
	- Third person language
	- First person plural pronouns
	- Kolmogorov complexity
	- n-grams (1)
	- n-grams (2)
	- n-grams (3)
2	Words and Phrases of Interest	hetter2012 ; Lloyd2012 ; Goodman2011
3	Countries of Interest	HT_report2015
4	Multiple Victims Advertised	kennedy2012predictive
5	Victim Weight	tvpa2000 ; weightchart
6	Reference to Website or Spa Massage Therapy	kennedy2012predictive
	- Reference to a website
	- Reference to a Spa Massage Therapy

Table 2. Table 2: Description of the dataset.

Name	Value
Raw	20,822
Filtered	3,543
Unlabeled	3,343
Labeled	Positive	Negative
	70	130

Table 3. Table 3: AUC and accuracy results with 10-fold cross-validation on the labeled data. The best performance is in bold.

Learner	AUC			Accuracy
	${ℱ}_1$	${ℱ}_2$	$\{{ℱ}_1,{ℱ}_2\}$	${ℱ}_1$	${ℱ}_2$	$\{{ℱ}_1,{ℱ}_2\}$
$S^{3}VM-R$	0.91	0.9	0.96	0.91	0.9	0.97
Laplacian SVM	0.9	0.9	0.9	0.91	0.9	0.92
LabelSpreading (RBF)	0.78	0.87	0.84	0.8	0.85	0.86
LabelSpreading (KNN)	0.68	0.80	0.74	0.71	0.8	0.8
Co-Training (SVM)	0.82	0.94	0.92	0.85	0.94	0.93
SVM	0.82	0.9	0.91	0.85	0.92	0.93
KNN	0.76	0.91	0.81	0.79	0.92	0.84
Gaussian NB	0.78	0.91	0.9	0.82	0.9	0.9
Logistic Regression	0.82	0.89	0.88	0.85	0.92	0.92
AdaBoost	0.82	0.85	0.85	0.85	0.88	0.88
Random Forest	0.81	0.89	0.89	0.83	0.91	0.92

Table 4. Table 4: Precision, recall and F1-score for the positive and negative classes using ℱ 1 subscript ℱ 1 \mathcal{F}_{1} . Experiments were run using 10-fold cross-validation on the labeled data. The best performance is in bold.

Learner	Precision		Recall		F1-score
	$clas{s}_p$	$clas{s}_n$	$clas{s}_p$	$clas{s}_n$	$clas{s}_p$	$clas{s}_n$
$S^{3}VM-R$	0.91	0.92	0.91	0.93	0.91	0.92
Laplacian SVM	0.86	0.89	0.88	0.9	0.87	0.88
LabelSpreading (RBF)	0.76	0.78	0.77	0.73	0.8	0.81
LabelSpreading (KNN)	0.65	0.7	0.71	0.68	0.69	0.73
Co-Training (SVM)	0.81	0.84	0.71	0.92	0.73	0.87
SVM	0.86	0.83	0.68	0.96	0.74	0.88
KNN	0.72	0.8	0.63	0.88	0.65	0.83
Gaussian NB	0.79	0.81	0.72	0.85	0.73	0.81
Logistic Regression	0.81	0.85	0.71	0.93	0.74	0.88
AdaBoost	0.86	0.83	0.68	0.95	0.74	0.88
Random Forest	0.77	0.85	0.73	0.89	0.73	0.86

Table 5. Table 5: Precision, recall and F1-score for the positive and negative classes using ℱ 2 subscript ℱ 2 \mathcal{F}_{2} . Experiments were run using 10-fold cross-validation on the labeled data. The best performance is in bold.

Learner	Precision		Recall		F1-score
	$clas{s}_p$	$clas{s}_n$	$clas{s}_p$	$clas{s}_n$	$clas{s}_p$	$clas{s}_n$
$S^{3}VM-R$	0.91	0.9	0.9	0.9	0.91	0.92
Laplacian SVM	0.91	0.9	0.9	0.91	0.89	0.92
LabelSpreading (RBF)	0.8	0.86	0.82	0.83	0.81	0.85
LabelSpreading (KNN)	0.7	0.75	0.73	0.78	0.79	0.77
Co-Training (SVM)	0.96	0.91	0.91	0.97	0.93	0.93
SVM	0.93	0.91	0.84	0.97	0.87	0.93
KNN	0.87	0.92	0.88	0.94	0.87	0.93
Gaussian NB	0.78	0.96	0.94	0.87	0.84	0.91
Logistic Regression	0.98	0.89	0.81	0.98	0.88	0.93
AdaBoost	0.88	0.88	0.75	0.95	0.78	0.91
Random Forest	0.93	0.89	0.81	0.97	0.85	0.93

Table 6. Table 6: Precision, recall and F1-score for the positive and negative classes using { ℱ 1 , ℱ 2 } subscript ℱ 1 subscript ℱ 2 \{\mathcal{F}_{1},\mathcal{F}_{2}\} . Experiments were run using 10-fold cross-validation on the labeled data. The best performance is in bold.

Learner	Precision		Recall		F1-score
	$clas{s}_p$	$clas{s}_n$	$clas{s}_p$	$clas{s}_n$	$clas{s}_p$	$clas{s}_n$
$S^{3}VM-R$	0.97	0.97	0.95	0.98	0.94	0.95
Laplacian SVM	0.96	0.94	0.91	0.96	0.91	0.93
LabelSpreading (RBF)	0.83	0.86	0.82	0.84	0.81	0.86
LabelSpreading (KNN)	0.71	0.74	0.75	0.78	0.8	0.78
Co-Training (SVM)	0.92	0.9	0.9	0.94	0.91	0.92
SVM	0.96	0.92	0.84	0.97	0.89	0.94
KNN	0.84	0.83	0.67	0.95	0.73	0.88
Gaussian NB	0.77	0.96	0.94	0.87	0.84	0.91
Logistic Regression	0.95	0.9	0.79	0.97	0.85	0.93
AdaBoost	0.88	0.88	0.75	0.95	0.78	0.91
Random Forest	0.93	0.9	0.82	0.97	0.86	0.93

Table 7. Table 7: Blind evaluation of the baselines on the two control groups. The best performance is in bold.

Learner	AUC
	${ℱ}_1$	${ℱ}_2$	$\{{ℱ}_1,{ℱ}_2\}$
Laplacian SVM	0.9	0.92	0.93
LabelSpreading (RBF)	0.75	0.85	0.87
LabelSpreading (KNN)	0.7	0.82	0.79
Co-Training (SVM)	0.8	0.9	0.91
SVM	0.8	0.65	0.69
KNN	0.74	0.62	0.77
Gaussian NB	0.77	0.51	0.52
Logistic Regression	0.76	0.62	0.75
AdaBoost	0.77	0.74	0.74
Random Forest	0.8	0.8	0.8

Table 8. Table 8: Significance of the features in ℱ 1 subscript ℱ 1 \mathcal{F}_{1} . The check-marked features show the top features.

No.	Feature Group	${\chi }^2$	Selected
1	Advertisement Language Pattern
	- Third person language	8.4	✓
	- First person plural pronouns	9.5	✓
	- Kolmogorov complexity	0.7	✓
	- n-grams (1)	0.4
	- n-grams (2)	0.0
	- n-grams (3)	0.4
2	Words and Phrases of Interest	0.0
3	Countries of Interest	59.3	✓
4	Multiple Victims Advertised	14.1	✓
5	Victim Weight	0.2
6	Reference to Website or Spa Massage Therapy
	- Reference to website	0.1
	- Reference to Spa Massage Therapy	33.5	✓

Table 9. Table 9: Classification results (AUC) using 10-fold cross-validation and different subsets of the features on the labeled data.

Name	Value
	$\overline{{ℱ}_1^{\ast }}$	${ℱ}_1^{\ast }$	${ℱ}_1$
$S^{3}VM-R$	0.82	0.87	0.91

Equations

$$K(X)\approx -\sum _{i=1}^{n}P(x_i){\log }_{2}P(x_i)$$

$$\underset{f_{\theta }\in {\mathcal{H}}_k}{\min }\gamma ||f_{\theta }|{|}_k^2+C_l\sum _{i=1}^l{H}_1(y_i{f}_{\theta }(x_i))$$

$$f_{\theta }^{\ast }(x)=\sum _{i=1}^l{\alpha }_i^{\ast }\mathbf{K}(x,x_i)$$

$$\underset{w,b,ϵ}{\min }\frac{1}{2}||w|{|}_2^2+C_l\sum _{i=1}^l{ϵ}_i$$

$$s.t.y_i(w.\mathbf{\Phi }(x_i)+b)\ge 1-{ϵ}_i,i=1,...,l$$

$${ϵ}_i\ge 0,i=1,...,l$$

$${\mathbf{F}}_{ij}=\frac{1}{n_f}(\mathbf{\Phi }(x_i)\cdot \mathbf{\Phi }(x_j))$$

$$\underset{f_{\theta }\in {\mathcal{H}}_k}{\min }\frac{1}{2}\sum _{i=1}^l{\text{F}}_{ij}||f_{\theta }(x_i)-f_{\theta }(x_j)|{|}_2^2={\mathbf{f}}_{\theta }^T{\mathcal{L}}^T{\mathbf{f}}_{\theta }$$

$$\underset{f_{\theta }\in {\mathcal{H}}_k}{\min }\gamma ||f_{\theta }|{|}_k^2+C_l\sum _{i=1}^l{H}_1(y_i{f}_{\theta }(x_i))+C_r{\text{f}}_{\theta }^T\mathcal{L}{\text{f}}_{\theta }+C_s{\text{f}}_{\theta }^T{\mathcal{L}}^{\prime }{\text{f}}_{\theta }$$

$$f_{\theta }^{\ast }(x)=\sum _{i=1}^{l+u}{\alpha }_i^{\ast }\mathbf{K}(x,x_i)$$

$$\underset{\alpha ,b,ϵ}{\min }\frac{1}{2}{\alpha }^T\mathbf{K}\alpha +C_l\sum _{i=1}^l{ϵ}_i+\frac{C_r}{2}{\alpha }^T\mathbf{K}\mathcal{L}\mathbf{K}\alpha +\frac{{\gamma }_I}{2(l+u{)}^2}{\alpha }^T\mathbf{K}{\mathcal{L}}^{\prime }\mathbf{K}\alpha$$

$$s.t.y_i(\sum _{j=1}^{l+u}{\alpha }_j\mathbf{K}(x_i,x_j)+b)\ge 1-{ϵ}_i,i=1,...,l$$

$${ϵ}_i\ge 0,i=1,...,l$$

$$L(\alpha ,ϵ,b,\beta ,\gamma )=\frac{1}{2}{\alpha }^T\mathbf{K}(I+C_r\mathcal{L}+\frac{{\gamma }_I}{(l+u{)}^2}{\mathcal{L}}^{\prime })\alpha +C_l\sum _{i=1}^l{ϵ}_i$$

$$-\sum _{i=1}^l{\beta }_i(y_i(\sum _{j=1}^{l+u}{\alpha }_j\mathbf{K}(x_i,x_j)+b)-1+{ϵ}_i)-\sum _{i=1}^l{\gamma }_i{ϵ}_i$$

$$\frac{\partial L}{\partial b}=0\to \sum _{i=1}^l{\beta }_i{y}_i=0$$

$$\frac{\partial L}{\partial {ϵ}_i}=0\to C_l-{\beta }_i-{\gamma }_i=0\to 0\le {\beta }_i\le C_l$$

$$L^R(\alpha ,\beta )=\frac{1}{2}{\alpha }^T\mathbf{K}(I+C_r\mathcal{L}+\frac{{\gamma }_I}{(l+u{)}^2}{\mathcal{L}}^{\prime })\alpha$$

$$-\sum _{i=1}^l{\beta }_i(y_i(\sum _{j=1}^{l+u}{\alpha }_j\mathbf{K}(x_i,x_j)+b)-1+{ϵ}_i)$$

$$L^R(\alpha ,\beta )=\frac{1}{2}{\alpha }^T\mathbf{K}(I+C_r\mathcal{L}+\frac{{\gamma }_I}{(l+u{)}^2}{\mathcal{L}}^{\prime })\alpha$$

$$-{\alpha }^T\mathbf{K}{\mathbf{J}}^T\mathbf{Y}\beta +\sum _{i=1}^l{\beta }_i$$

$$\mathbf{K}(I+C_r\mathcal{L}+\frac{{\gamma }_I}{(l+u{)}^2}{\mathcal{L}}^{\prime })\alpha -\mathbf{K}{\mathbf{J}}^T\mathbf{Y}\beta =0$$

$${\alpha }^{\ast }=(I+C_r\mathcal{L}+\frac{{\gamma }_I}{(l+u{)}^2}{\mathcal{L}}^{\prime }{)}^{-1}{\mathbf{J}}^T\mathbf{Y}{\beta }^{\ast }$$

$${\beta }^{\ast }={\mathrm{argmax}}_{\beta \in {\mathbb{R}}^l}-\frac{1}{2}{\beta }^T\mathbf{Q}\beta +\sum _{i=1}^l{\beta }_i$$

$$s.t.\sum _{i=1}^l{\beta }_i{y}_i=0$$

$$0\le {\beta }_i\le C_l$$

$$\mathbf{Q}=\mathbf{YJK}(I+(C_r\mathcal{L}+\frac{{\gamma }_I}{(l+u{)}^2}{\mathcal{L}}^{\prime })\mathbf{K}{)}^{-1}{\mathbf{J}}^T\mathbf{Y}$$

Full text

1 Introduction

According to the United Nation unodc , human trafficking is defined as the modern slavery or the trade of humans mostly for the purpose of sexual exploitation and forced labor, via different improper ways including force, fraud and deception. The United States’ Trafficking Victim Protection Act of 2000 (TVPA 2000) tvpa2000 was the first U.S. legislation passed against human trafficking. Human trafficking has ever since received increased national and societal concern HT_report2015 but still demands persistent fight against from all over the globe. No country is immune and the problem is rapidly growing with little to no law enforcement addressing the issue. This problem is amongst the challenging ones facing law enforcement as it is difficult to identify victims and counter traffickers.

Before the advent of the Internet, human traffickers were under risks of being arrested by law enforcement while advertising their victims on streets desplaces92 . However, move to the Internet has made it easier and less dangerous for sex sellers nicholas2012 as they no longer needed to advertise on the streets. There are now a plethora of websites that host and provide sexual services under categories of escort, adult entertainment, massage services, etc., which help sex sellers and buyers maintain their anonymity. Although some services such as the Craiglist’s adult section and myredbook.com were shut down recently, there are still many websites such as the Backpage.com that provide such services and many new are frequently created. Traffickers even use dating and social networking websites such as the Twitter, Facebook, Instagram and Tinder to reach out to sex buyers and their followers. Although the Internet has presented new trafficking related challenges for law enforcement, it has also provided readily and publicly available rich source of information which could be gleaned from online sex advertisements for fighting this crime kennedy2012predictive . However, the problem is we lack the ground truth and obtaining the labels through hand-labeling is indeed tedious and expensive even for a small subset of data– this is the point where the semi-supervised setting comes in handy.

Despite considerable attention which has been devoted to studying supervised, unsupervised and semi-supervised learning settings via different applications mitchell2006learning ; beigi2016signed ; backstrom2011supervised ; alvari2016identifying ; beigi2016exploiting ; mitchell1997machine ; beigi2014leveraging , semi-supervised learning, i.e., learning from labeled and unlabeled examples, is still one of the most interesting yet challenging problems in the machine learning community belkin2006manifold . The idea is simple though– we shall have an approach that makes a better use of unlabeled data to boost performance. This is pretty close to the most natural learning that occurs in the world. For the most part, we as humans are exposed only to a small number of labeled instances; yet we successfully generalize well by effective utilization of a large amount of unlabeled data. This motivates us to use unlabeled samples to improve recognition performance while developing classifiers.

In this article, expanding on our previous work alvari2016non , we use the data crawled from the adult entertainment section of the website Backpage.com and extend the existing Laplacian SVM framework belkin2006manifold to detect escort advertisements of high interest to law enforcement. Here, we merely focus on the online advertisements although the Internet has triggered many other activities including attracting the victims, communicating with customers and rating the escort services. We thus highlight several contributions of the current research as follows.

Based on the literature, we created different groups of features that capture the characteristics of potential human trafficking activities. The less likely human trafficking related posts were then filtered out using these features. We also conducted a feature importance analysis to demonstrate how these features contribute to the proposed learner. 2. 2.

We extended the Laplacian SVM belkin2006manifold and proposed the semi-supervised support vector machine learning algorithm, $S^{3}VM-R$. In particular, we incorporated additional information of our feature space as a regularization term into the standard optimization formulation with regard to the Laplacian SVM. We also used geometry of the underlying data as an intrinsic regularization term in Laplacian SVM. 3. 3.

We trained our model on both of the labeled and unlabeled data and sent back the identified human trafficking related advertisements to an expert from law enforcement for further verification. We then validated our approach on a small subset of the unlabeled data (i.e. unseen data) with further verification of the expert. 4. 4.

We performed comparisons between our approach and several semi-supervised and supervised baselines on both of the labeled and unseen data (so-called blind evaluation). 5. 5.

We demonstrated the effect of varying different hyperparameters used in our learner on its performance.

The rest of the paper is organized as follows. In Section 2, we review the prior studies on human trafficking. Section 3 covers our data preparation, feature engineering, unsupervised filtering and expert assisted labeling. We detail our semi-supervised learning approach in Section 4 by deriving the required equations. Section 5 provides in-depth explanation of our experiments. Section 6 concludes the paper by providing future research directions.

2 Related Work

Recently, several studies have examined the role of the Internet and related technology in facilitating human trafficking hughes2005demand ; hughes2002use ; latonero2011human . For example, the work of hughes2005demand studied how closely sex trafficking is intertwined with new technologies. According to hughes2002use , sexual exploitation of women and children is a global human right crisis that is being escalated by the use of new technologies. Researchers have studied relationships between new technologies and human trafficking and advantages of the Internet for sex traffickers. For instance, findings from a group of experts from the Council of Europe demonstrated that the Internet and sex industry are closely interlinked and volume and content of the material on the Internet promoting human trafficking are unprecedented latonero2011human .

One of the earliest works which leveraged data mining techniques for online human trafficking was latonero2011human , wherein the authors conducted data analysis on the adult section of the website Backpage.com. Their findings confirmed that female escort post frequency would increase in Dallas, Texas, leading up to the Super Bowl 2011 event. In a similar attempt, other studies dominique2015 ; 2016arXiv160205048M have investigated impact of large public events such as the Super Bowl on sex trafficking by exploring advertisement volume, trends and movement of advertisements along with the scope and volume of demand associated with such events. The work of dominique2015 for instance, concluded that large events such as the Super Bowl which attract significant amount of concentration of people in a relatively short period of time and in a confined urban area, could be a desirable location for sex traffickers to bring their victims for commercial sexual exploitation. Similarly, the data-driven approach of 2016arXiv160205048M showed that in some but not all events, one can see a correlation between occurrence of the event and statistically significant evidence of an influx of sex trafficking activity. Also, certain studies conf/semweb/SzekelyKSPSYKNM15 have tried to build large distributed systems to store and process available online human trafficking data in order to perform entity resolution and create ontological relations between entities.

Beyond these works, the work of 2015arXiv150906659N studied the problem of isolating sources of human trafficking from online advertisements with a pairwise entity resolution approach. Specifically, they used phone number as a strong feature and trained a classifier to predict if two ads are from the same source. This classifier was then used to perform entity resolution using a heuristically learned value for the score of classifier. Another work of kennedy2012predictive used Backpage.com data and extracted most likely human trafficking spatio-temporal patterns with the help of law enforcement. Note that unlike our method, this work did not employ any machine learning methodologies for automatically identifying human trafficking related advertisements. The work of doi:10.1080/23322705.2015.1015342 also deployed machine learning for the advertisement classification problem, by training a supervised learning classifier on labeled data (based on phone numbers of known traffickers) provided by a victim advocacy group. We note that while phone numbers can provide a very precise set of positive labeled data, there are clearly many posts with previously unseen phone numbers.

In contrast, we do not solely rely on phone numbers for labeling our data. Instead, our expert analyze each post’s content to identify whether it is human trafficking related or not. To do so, we first filter out most likely advertisements using several feature groups and pass a small sample to the expert for hand-labeling. Then, we train our semi-supervised learner on both of the labeled and unlabeled data which in turn lets us evaluate our approach on new coming (unseen) data later. We note that our semi-supervised approach can also be used as a complementary method to procedures such as those described in doi:10.1080/23322705.2015.1015342 as we can significantly expand the training set for use with supervised learning.

Finally, note that our current research is different from our previous work alvari2016non and we list the key nuances here:

• •

In this study we experiment with a much larger dataset. To obtain such dataset, we use the same raw data from alvari2016non , but this time with slight modifications of the thresholds that were used for filtering out less likely human trafficking related advertisements.

• •

As opposed to our previous research which deployed only one feature space, in this work, two feature spaces that have complementary roles to each other are used.

• •

In this paper we present a new framework based on the existing Laplacian SVM belkin2006manifold , by adding a regularization term to the standard optimization problem and solving the new optimization equation derived from there. In contrast, alvari2016non utilized the off-the-shelf graph based semi-supervised learner, LabelSpreading method Zhou04learningwith , without any further manipulation of the original approach.

• •

Unlike alvari2016non in which we did not compare our method with other approaches, this work compares our proposed framework against other semi-supervised and supervised learners. Also unlike our previous work in which only one group of human trafficking related advertisements were passed to two experts for validation, here in order to reduce the inconsistency, two control groups of advertisements–those of interest to law enforcement and those of not– are sent to only one expert for verification.

3 Data Preparation

We collected about 20K publicly available listings from the U.S. posted on Backpage.com in March, 2016. Each post includes a title, description, time stamp, poster’s age, poster’s ID, location, image, and sometimes video and audio. The description usually lists the attributes of the individual(s) and contact phone numbers. In this work, we only focus on the textual component of the data. This free-text data required significant cleaning due to a variety of issues common to textual analytics (i.e. misspellings, format of phone numbers, etc.). We also acknowledge that the information in data could be intentionally inaccurate, such as poster’s name, age and even physical appearance (e.g. bra cup size, weight). Figure 1 shows an actual post from Backpage.com. To illustrate geographic diversity of the listings, we use the Tableau111https://www.tableau.com/ software to visualize choropleth map of phone frequency with respect to the different states in Figure 2, wherein darker colors mean higher frequencies.

Next, we will explain most important characteristics of potential human trafficking advertisements which are captured by our feature groups.

3.1 Feature Engineering

Though many advertisements on Backpage.com are posted by posters selling their own services without coercion and intervention of traffickers, some do exhibit many common trafficking triggers. For example, in contrast to Figure 1, Figure 3 shows an advertisement that could be an evidence of human trafficking. This advertisement indicates several potential properties of human trafficking, including advertising for multiple escorts with the first individual coming from Asia and very young. In what follows, such common properties of human trafficking related advertisements are discussed in more detail.

Inspired by the literature, we define and extract 6 groups of features from advertisements (see Table 1). These features could be amongst the strong indicators of human trafficking. Let us now briefly describe each group of features used in our work. Note each feature listed here is ultimately treated as a binary variable.

3.1.1 Advertisement Language Pattern

The first group consists of different language related features. For the first and second features, we identify posts which have third person language (more likely to be written by someone other than the escort) and posts which contain first person plural pronouns such as ‘we’ and ‘our’ (more likely to be an organization) kennedy2012predictive .

To ensure their anonymity, traffickers would deploy techniques to generate diverse information and hence make their posts look more complicated. They usually do this to avoid being identified by either human analysts or automated programs. Thus, to obtain the third feature we take an approach from complexity theory, namely Kolmogorov complexity, which is defined as length of shortest program to reproduce a string of characters on a universal machine such as the Turing Machine Li:2008:IKC:1478784 . Since the Kolmogorov complexity is not computable, we approximate the complexity of an advertisement content by first removing stop words and then computing entropy of the content Li:2008:IKC:1478784 . To illustrate this, let $X$ denote the content and $x_{i}$ be a given word in the content. We use the following equation shannon2001mathematical to calculate the entropy of the content and thus approximate the Kolmogorov complexity of $X$:

[TABLE]

We expect higher values of the entropy correspond to human trafficking. Finally, we discretize the result by using the threshold of 4 which was found empirically in our experiments.

For the next features, we use word-level n-grams to find common language patterns of advertisements. This particular choice is because of the fact that character-level n-grams have already shown to be useful in detecting unwanted content for spam detection conf/setn/KanarisKS06 . We set $n=4$ and use the range of (4,4) to compute normalized n-grams (using TF-IDF) of each advertisement content. We ultimately create a matrix whose rows and columns correspond to the advertisements contents and their associated 4-grams, respectively. We rank all elements of this matrix in a descending order and pick the top 3 ones. Finally for each advertisement content, 3 elements with the column numbers associated with the top elements are chosen. This way, 3 more features will be added to our feature set. Overall, we have 6 features related to the language of the advertisement.

3.1.2 Words and Phrases of Interest

Despite the fact that advertisements on Backpage.com do not directly mention sex with children, customers who prefer children know to look for words and phrases such as “sweet, candy, fresh, new in town, new to the game” hetter2012 ; Lloyd2012 ; Goodman2011 . We thus investigate within the posts to see if they contain such words as they could be highly related with human trafficking in general.

3.1.3 Countries of Interest

We identify if the individual being escorted is coming from other countries such as those in Southeast Asia (especially from China, Vietnam, Korea and Thailand, as we observed in our data) HT_report2015 .

3.1.4 Multiple Victims Advertised

Some advertisements advertise for multiple women at the same time. We consider the presence of more than one victim as a potential evidence of organized human trafficking kennedy2012predictive .

3.1.5 Victim Weight

We take into account the weight of the individual being escorted as a feature (if it is available). This information is particularly useful assuming that for the most part, lower body weights (under 115 lbs) correlate with smaller and underage girls tvpa2000 ; weightchart and thereby human trafficking.

3.1.6 Reference to Website or Spa Massage Therapy

The presence of a link in the advertisement either referencing to an outside website (especially infamous ones) or spa massage therapy could be an indicator of more elaborate organization kennedy2012predictive . In particular, in case of spa therapy, we observed many advertisements interrelated with advertising for young Asian girls and their erotic massage abilities. Therefore, the last group of features has two binary features for presence of any website and spa.

Finally, in order to extract all of the above features, we first clean the original data and conduct preprocessing. By applying these features, we draw a random sample of 3,543 instances out of our dataset for further analysis to see if they are evidences of human trafficking– this is described in the next section.

3.2 Unsupervised Filtering

Having detailed our feature set, we now construct a feature vector for each instance by creating a vector of 12 binary features that correspond to the important characteristics of human trafficking. Hereafter, we refer to this feature space, as our first feature space and denote it with $\mathcal{F}_{1}$. As mentioned earlier, we draw 3,543 instances from our raw data by filtering out those that do not posses any of the binary features. We will refer to this as our filtered dataset. For the sake of visualization, a 2-D projection (using the t-SNE transformation maaten2008visualizing ) of the filtered dataset is depicted in Figure 4. The purpose of this figure is to demonstrate how hard it is for basic clustering techniques such as the K-means, to correctly assign labels to unlabeled instances using only few existing labeled ones.

Now, we shall define our second feature space, namely $\mathcal{F}_{2}$, which will be used to compute geometry of the underlying data. Note that our proposed framework will utilize both of the feature spaces in the form of regularization terms, to detect advertisements of high interest to law enforcement. After conducting standard preprocessing techniques on the filtered dataset, we build $\mathcal{F}_{2}$ by transforming the filtered data into a 3,543$\times$3,543 matrix of TF-IDF similarity features. Each entry in this matrix simply shows the similarity between a pair of advertisements in our filtered dataset.

Note that since we lack the ground truth, we would rely on a human analyst (expert) for labeling the listings as either ‘of interest’ or ‘of not interest’ to law enforcement. In the next section, we select a smaller yet finer grain subset of this data to be sent to the expert. This alleviates the burden of the tedious work of hand-labeling.

3.3 Expert Assisted Labeling

We first obtain a sample of 200 listings from the filtered dataset. This set of listings was labeled by our expert from law enforcement who is specialized in this type of crime. From this subset, the law enforcement professional identified 70 instances to be of interest to law enforcement and the rest to be not human trafficking related. However, we are still left with a large amount of the unlabeled examples (3,343 instances) in our dataset. The ratio of the labeled to unlabeled instances in our dataset is very small (about 0.06). The statistics of our dataset is summarized in Table 2.

4 Semi-Supervised Learning Framework

Here, we first introduce some preliminary notations necessary for the rest of the discussion and then outline our proposed semi-supervised approach, $ُS^3VM-R$,fordetectingonlinehumantrafficking.Noteassaidearlier,ourframeworkisanextensiontotheexistingLaplacianSVM~{}\cite[cite]{@@bibref{Authors Phrase1YearPhrase2}{belkin2006manifold}{@@citephrase{(}}{@@citephrase{)}}}.Inparticular,weincorporatedanotherregularizationtermintothestandardLaplacianSVMtoleveragetheadditionalinformationofourfirstfeaturespaceandthensolvedtheassociatedoptimizationproblem.Consequently,similarnotationisadoptedthroughoutthefollowingsection.Furthermore,weshallonceagainnotethatourcurrentresearchdoesnotutilizeanyoff-the-shelfgraphbasedsemi-supervisedleanerincontrasttoourpreviousresearch~{}\cite[cite]{@@bibref{Authors Phrase1YearPhrase2}{alvari2016non}{@@citephrase{(}}{@@citephrase{)}}}.\par$

4.1 Technical Preliminaries

We assume a set of $l$ labeled pairs $\{(x_{i},y_{i})\}_{i=1}^{l}$ and an unlabeled set of $u$ instances $\{x_{l+i}\}_{i=1}^{u}$, where $x_{i}\in\mathbb{R}^{n}$ and $y_{i}\in\{+1,-1\}$. Recall for the standard soft-margin support vector machine, the following optimization problem is solved:

[TABLE]

In the above equation, $f_{\theta}(\cdot)$ is a decision function of the form $f_{\theta}(\cdot)=w.\mathbf{\Phi}(\cdot)+b$ where $\theta=(w,b)$ are the parameters of the model, and $\mathbf{\Phi(\cdot)}$ is the feature map which is usually implemented using the kernel trick cortes1995support . Also, the function $H_{1}(\cdot)=\max(0,1-\cdot)$ is the Hinge Loss function.

The classical Representer theorem belkin2005manifold suggests that solution to the optimization problem exists in a Hilbert space $\mathcal{H}_{k}$ and is of the following form:

[TABLE]

where $\mathbf{K}$ is the $l\times l$ Gram matrix over labeled samples. Equivalently, the above problem can be written as:

[TABLE]

[TABLE]

Next, we will use the above optimization equation as our basis to derive the formulations for our proposed semi-supervised learner.

4.2 The proposed Method

The basic assumption behind semi-supervised learning methods is to leverage unlabeled instances in order to restructure hypotheses during the learning process. In this paper, exogenous information extracted from both of our feature spaces is further exploited to make a better use of the unlabeled examples. To do so, we first introduce matrix $\mathbf{F}$ in $\mathcal{F}_{1}$ and over both of the labeled and unlabeled samples with $\mathbf{F}_{ij}$ defined as follows:

[TABLE]

where $n_{f}$ is the number of features in $\mathcal{F}_{1}$ (here, $n_{f}=12$). We force the instances $x_{i}$ and $x_{j}$ in our dataset to have same label if they both possess same features. To account for this, a regularization term is added to the standard equation and the following optimization is solved:

[TABLE]

where $\mathbf{f}=[f(x_{1}),...,f(x_{l+u})]^{T}$ and $\mathcal{L}$ is the Laplacian matrix based on $\mathbf{F}$ given by $\mathcal{L}=\mathbf{D}-\mathbf{F}$, and $\mathbf{D}_{ii}=\sum_{j=1}^{l+u}\mathbf{F}_{ij}$. The intuition here is that any two instances which are composed of same features are more likely to have same labels than others. Next, by solving a similar optimization problem, we are able to capture data geometry in $\mathcal{F}_{2}$ as $\mathbf{f}^{T}_{\theta}\mathcal{L}^{\prime T}\mathbf{f}_{\theta}$ (also referred to as the intrinsic smoothness penalty term belkin2006manifold ). Here, $\mathcal{L}^{\prime}$ is the Laplacian of matrix $\mathbf{A}$ associated with the data adjacency graph $\mathbf{G}$ in $\mathcal{F}_{2}$.

We construct $\mathbf{G}$ with $(l+u)$ nodes in $\mathcal{F}_{2}$, and by adding an edge between each pair of nodes $\langle i,j\rangle$, if the edge weight $W_{ij}$ exceeds a given threshold. For computing the edge weights, we use the heat kernel grigor2006heat as a function of the Euclidean distance between two samples in $\mathcal{F}_{2}$, hence we set $W_{ij}=\exp^{-||x_{i}-x_{j}||^{2}/4t}$.

Following the notations used in belkin2006manifold and by including our regularization term as well as the intrinsic smoothness penalty term, we would extend the standard equation by solving the following optimization:

[TABLE]

Note one typical value for the smoothness penalty coefficient $C_{s}$ is $\frac{\gamma_{I}}{(l+u)^{2}}$, where $\frac{1}{(l+u)^{2}}$ is a natural scale factor for empirical estimate of the Laplace operator and $\gamma_{I}$ is a regularization term belkin2006manifold . Again, solution in $\mathcal{H}_{k}$ would be in the following form:

[TABLE]

Here $\mathbf{K}$ is the $(l+u)\times(l+u)$ Gram matrix over all samples. The equation 8 could be then written as follows:

[TABLE]

[TABLE]

With introduction of the Lagrangian multipliers $\beta$ and $\gamma$, we write the Lagrangian function of the above equation as follows:

[TABLE]

Obtaining the dual representation, requires taking the following steps:

[TABLE]

[TABLE]

With the above equations, we formulate the reduced Lagrangian as a function of only $\alpha$ and $\beta$ as follows:

[TABLE]

This equation is further simplified as follows:

[TABLE]

In the above equation, $\mathbf{J}=[\mathbf{I}~{}\mathbf{0}]$ is a $l\times(l+u)$ matrix, $\mathbf{I}$ is the $l\times l$ identity matrix and $\mathbf{Y}$ is a diagonal matrix consisting of the labels of the labeled examples.

In the followings, we first take the derivative of $L^{R}$ with respect to $\alpha$ and then set $\frac{\partial L^{R}(\alpha,\beta)}{\partial\alpha}=0$:

[TABLE]

Accordingly, we obtain $\alpha^{*}$ by solving the following equation:

[TABLE]

Next, we obtain the dual problem in the form of a quadratic programming problem by substituting $\alpha$ back in the reduced Lagrangian function:

[TABLE]

[TABLE]

where $\beta=[\beta_{1},...,\beta_{l}]^{T}\in\mathbb{R}^{l}$ are the Lagrangian multipliers and $\mathbf{Q}$ is obtained as follows:

[TABLE]

We summarize the proposed semi-supervised framework in Algorithm 1. Our optimization problem is very similar to the standard optimization problem solved for SVMs, hence we use a standard optimizer for SVMs to solve our problem.

5 Experimental Study

In this section, we provide a comprehensive analysis of the proposed framework by designing a series of experiments on the filtered dataset. First, we explain several approaches used in this study. Next, various results are discussed: (1) comparisons on the labeled data were made between our method and other approaches, (2) experiments were performed on a fraction of the unlabeled data (i.e., unseen data), and the results were further verified by our expert to see what fraction is of interest to law enforcement, (3) blind evaluation was conducted to examine other approaches on the unseen data, and finally, (4) experiments were designed to analyze effect of varying different hyperparameters on our method as well as impact of different groups of features in $\mathcal{F}_{1}$ on our approach.

5.1 Approaches

We present results for the following methods:

• •

Semi-Supervised: $S^{3}VM-R$, Laplacian support vector machines belkin2006manifold , graph inference based label spreading approach Zhou04learningwith with radial basis function (RBF) and K-nearest neighbors (KNN) kernels, and co-training learner blum1998combining with two support vector machines classifiers (SVM).

• •

Supervised: SVM, KNN, Gaussian naïve Bayes, logistic regression, adaboost and random forest.

For the sake of fair comparison, all algorithms were implemented and run in Python. More specifically, the Python package CVXOPT222http://cvxopt.org/ was used to implement $S^{3}VM-R$ and Laplacian support vector machines, and all other approaches were implemented with the help of the Scikit-learn333http://scikit-learn.org/stable/ package in Python. Note for those methods that require special tuning of parameters, we performed grid search to choose the best set of parameters. Before going any further, we first define main parameters used in each method and then demonstrate their best values picked by our grid search. The discussion on the effect of varying the hyperparameters on our learner is provided in the section 5.3.

• •

$S^{3}VM-R$: we set the penalty parameter as $C_{l}=0.6$ and the regularization parameters $C_{r}=0.2$ and $C_{s}=0.2$. Linear kernel was used in our approach.

• •

Laplacian SVM: we used linear kernel and set the parameters $C_{l}=0.6$ and $C_{s}=0.6$.

• •

LabelSpreading (RBF): RBF Kernel was used and $\gamma$ was set to the default value of 20.

• •

LabelSpreading (KNN): KNN kernel was used and the number of neighbors was set to 5.

• •

Co-training (SVM): we followed the algorithm introduced in blum1998combining and used two SVM as our classifiers. For both SVMs we set the tolerance for stopping criteria to 0.001 and the penalty parameter $C=1$.

• •

SVM: tolerance for stopping criteria was set to the default value of 0.001. Penalty parameter $C$ was set to 1 and linear kernel was used.

• •

KNN: number of neighbors was set to 5.

• •

Gaussian NB: there were no specific parameter to tune.

• •

Logistic regression: we used the ‘l2’ penalty. We also set the parameter $C=1$ (the inverse of regularization strength) and tolerance for stopping criteria to 0.01.

• •

Adaboost: number of estimators was set to 200 and we also set the learning rate to 0.01.

• •

Random forest: we used 200 estimators and the ‘entropy’ criterion was used.

5.2 Classification Results

Here, we first evaluate the entire set of approaches on a small portion of the data for which we already know the labels, i.e., the labeled examples. We note that expert-generated judgmental labeling might be error-prone, though it is served as a surrogate to the ground truth problem.

We used 10-fold cross-validation on the labeled data in the following way. We first divided the set of the labeled samples into 10 different sets of approximately equal size. Each time we held one set out for validation (by removing their labels and adding them to the unlabeled samples) and used the remaining along with the unlabeled samples for the training–this was performed for all approaches for the sake of fair comparison. Finally, we reported the average of 10 different runs, using different combinations of the feature spaces and various evaluation metrics, including the area under curve (AUC), accuracy, precision, recall and F1-score. In table 3, we reported the average AUC and accuracy for each method and each feature space. On the other hand, for precision, recall and F1-score, we reported separate results for each feature space, in tables 4-6, respectively. Note, each of these tables includes separate scores for the positive and negative classes. In general, we observed the followings:

• •

Overall, our approach achieved highest performance on $\mathcal{F}_{1}$ (tables 3 and 4) and $\{\mathcal{F}_{1},\mathcal{F}_{2}\}$ (table 6), in terms of all metrics. However it did not perform well using solely $\mathcal{F}_{2}$ (table 5), i.e. when $C_{r}=0$. This clearly demonstrates the importance of using $C_{r}$ over $C_{s}$.

• •

When the feature space used is $\mathcal{F}_{2}$, Co-training (SVM) is the best method. Next best methods are supervised learners KNN and Gaussian NB. Three remarks can be made here. First, our approach could not always defeat supervised learners as it is seen from tables 3 and 5. This is not surprising and in fact lies at the inherent difference between semi-supervised and supervised methods– unlabeled examples could make the trained model susceptible to error propagation and thus wrong estimation. Second, as it is seen in tables 4-6, achieving very high recall on the negative examples and low score on the positive ones shall not be treated as a potent property, otherwise a trivial classifier which always assigns negative labels to all samples would be the best learner. Third, using $C_{r}$ always improves the performance over $C_{s}$. One point that needs to be clarified is, our ultimate goal is not to achieve high performance on the labeled data, but rather to detect the suspicious (unlabeled) advertisements which could be human trafficking related– this will be explained in more details in 5.2.1.

• •

Compared to the other semi-supervised approaches, our approach either achieved higher or comparable AUC scores. The reason we performed exactly the same as the Laplacian SVM, is because by setting $C_{r}=0$, the two approaches are inherently the same.

• •

For the Laplacian SVM to be able to run on $\mathcal{F}_{1}$, the Laplacian $\mathcal{L}^{\prime}$ has to be constructed using $\mathcal{F}_{1}$ while inherently is supposed to be made using $\mathcal{F}_{2}$. This is because $C_{r}$ is essentially associated with $\mathcal{F}_{1}$, and $C_{s}$ corresponds to $\mathcal{L}^{\prime}$ and correspondingly $\mathcal{F}_{2}$. The same holds for $\{\mathcal{F}_{1},\mathcal{F}_{2}\}$, where we need to construct a new feature space by concatenating $\mathcal{F}_{1}$ and $\mathcal{F}_{2}$ as the Laplacian SVM does not inherently use $\mathcal{F}_{1}$ at all. The new feature space is then used to construct the Laplacian $\mathcal{L}^{\prime}$.

• •

Since our approach inherently incorporates both of the Laplacian matrices corresponding to the two feature spaces $\mathcal{F}_{1}$ and $\mathcal{F}_{2}$, all other baselines were also run using the concatenation of these two feature spaces for the sake of fair comparison. Unlike our approach which used the wise combination of $\mathcal{F}_{1}$ and $\mathcal{F}_{2}$, other methods do not gain high AUC by simply combining the feature spaces.

5.2.1 Blind Evaluation

For the next set of experiments, we first run our method on the entire filtered dataset and without cross-validation. Recall from the previous sections that this is to make a better use of the unlabeled examples. Then the following control experiment was conducted. Our learner was tested on the whole set of the unlabeled examples. Out of 3,343 instances, our approach identified two sets of positive and negative instances. The positive set contained 394 advertisements which were likely to be of interest to law enforcement, whereas the negative set included the remaining 2,962 unlabeled advertisements of probably less interest to law enforcement. Next, to precisely determine the correctly identified fractions of these two sets, we randomly picked two subsets (control groups) of 100 examples from each set for further validation by our expert.

We passed these two control groups to our expert for further verification. The expert-validated results demonstrated that all of the examples in the positive group were of interest to law enforcement, while only two examples from the negative group were not correctly classified as of not being of any interest to law enforcement. Thus, both results support the effectiveness of our framework in identifying highly human trafficking advertisements. Using the same two control groups and AUC metric, we now perform so-called blind evaluation (see table 7) of other baselines. Note, we call this blind since actual labels are not provided and the expert-generated labels might convey uninformative information. In general, supervised methods failed to achieve good results in the blind evaluation compared to most of the semi-supervised methods.

5.3 Hyperparameter Sensitivity

Here, we discuss how altering the hyperparameters $C_{l}$, $C_{r}$ and $C_{s}$ may affect the performance of $S^{3}VM-R$. We start off by fixing the value of $C_{l}$ to 0.6, which was empirically found to work well in our experiments. Also, recall from the previous sections that one typical choice for $C_{s}$ is $\frac{\gamma_{I}}{(l+u)^{2}}$ belkin2006manifold . Here, we set $C_{s}=0.2$ and varied the values of $C_{r}$ as $\{0,0.0002,0.0006,0.2,1.0\}$ and plotted the results in Figure 5. We used the same 10-fold cross-validation setting from the previous section.

We made the following observation. With the slight increase of $C_{r}$, the performance of our approach increased, peaked and then stabilized, i.e., further increase of $C_{r}$ did not change the performance. This suggests significance of deploying the additional information from our first feature space $\mathcal{F}_{1}$, over $\mathcal{F}_{2}$ and its corresponding smoothness penalty parameter $C_{s}$ which is used by $S^{3}VM-R$ and the standard Laplacian SVM.

Next, to see the impact of $C_{l}$ on the performance, we set $C_{r}=0.2$ and varied $C_{l}$ as $\{0.2,0.4,0.6,0.8,1.0\}$. The results are depicted in Figure 5. We note that setting $C_{l}=0$ is meaningless and thus we do not have any performance corresponding to that– otherwise each $\beta_{i}$ in Eq. 19 would be zero. In general, the performance was not particular sensitive to this parameter– varying by 0.2 for values of 0.4 and greater.

Finally, having fixed $C_{l}=0.6$ and $C_{r}=0.2$, we also tried other values for $C_{s}$ including $\sum_{i,j=1}^{l+u}W_{ij}$ suggested by belkin2006manifold and depicted the results in Figure 5. The results suggest that our approach is less sensitive to this parameter compared to $C_{r}$ and $C_{l}$.

5.4 Significance of Features

To examine how much discriminative our feature groups in $\mathcal{F}_{1}$ are, we further conducted an analysis using the labeled examples and the standard feature selection measure $\chi^{2}$ to find the top features– only half of the features with scores greater than a given threshold (0.5) were selected (see table 8 for the complete set of features and their corresponding $\chi^{2}$ scores).

From this list, we noticed that ‘countries of interest’ and ‘reference to spa massage therapy’ were the most discriminative feature groups, while ‘advertisement language pattern’ group (with 3 important features) appeared to be the most dominant feature group.

Figure 6 compares the top features against the less important subset of the features (denoted by $\overline{\mathcal{F}^{*}_{1}}$) in the filtered dataset, in terms of frequency values. Note for clarity, we have removed from this figure, the features with frequency less than 20. According to this figure, our most discriminative features are not necessarily those that appear more often.

To further investigate the importance of each of the top features, we performed classification using the labeled examples and the previous setting, on basis of these two subsets of the features and their combination, i.e., $\mathcal{F}^{*}_{1}$, $\overline{\mathcal{F}^{*}_{1}}$ and $\mathcal{F}_{1}$. The classification results are shown in Table 9. We made the following observations:

• •

Considering only the feature space $\mathcal{F}_{1}$, our approach achieved higher performance compared to all other baselines by either using the whole feature space or the most discriminative features $\mathcal{F}^{*}_{1}$.

• •

Deploying only the features from $\mathcal{F}^{*}_{1}$, we were able to achieve comparable results as if we used the whole feature space $\mathcal{F}_{1}$.

6 Conclusion

Readily available online data from escort advertisements could be leveraged in favor of fight against human trafficking. In this study, having focused on textual information from the available data crawled from Backpage.com, we identified if an escort advertisement can be reflective of human trafficking activities. In particular, we first proposed an unsupervised filtering approach to filter out the data which are more likely involved in human trafficking. We then proposed a semi-supervised learner, namely $S^{3}VM-R$, and trained it on a small portion of the data which was hand-labeled by a human trafficking expert. We used the trained model to identify labels of unseen data. Results suggested our approach is effective at identifying potential human trafficking related advertisements.

Our future plans include replicating the study by integrating more interesting features especially those supported by the criminology literature. Also, since hand-labeling unlabeled examples is expensive, an interesting research direction would be to deploy active learning to enable iterative supervised learning to actively query the user for labels. We also note that real-world data is often more imbalanced compared to our data, and the reason is that number of negative samples usually outweigh positive ones. We would thus like to apply the proposed framework on a more realistic dataset which contains much less suspicious posts than normal posts.

7 Competing Interests

The authors declare that they have no competing interests.

8 Authors’ Contributions

HA developed and implemented the human trafficking detection approach and drafted the manuscript. PS provided guidance through the whole project and revised the manuscript. JS was in contact with an expert from law enforcement who was responsible for hand-labeling portions of the data. All authors read and approved the final manuscript.

Acknowledgment

This work was funded by the Find Me Group, a 501(c)3 dedicated to bring resolution and closure to families of missing persons. The authors would like to thank anonymous reviewers for their valuable suggestions to improve the quality of the paper.