-
Notifications
You must be signed in to change notification settings - Fork 17
/
Copy pathclassifytraj.py
106 lines (79 loc) · 2.38 KB
/
classifytraj.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
from sklearn.linear_model import SGDClassifier
from sklearn.feature_extraction.text import CountVectorizer
from sklearn.feature_extraction.text import TfidfTransformer
import pandas as pd
from sklearn.decomposition import TruncatedSVD
from sklearn import preprocessing
from sklearn.pipeline import Pipeline
import numpy as np
from sklearn.neighbors import KNeighborsClassifier
from sklearn.linear_model import LogisticRegression
from sklearn.ensemble import RandomForestClassifier
import csv
#k-fold cross validation
from sklearn.cross_validation import KFold
#accuracy
from sklearn.metrics import accuracy_score
# my method
from sklearn.ensemble import VotingClassifier
def compute_and_print(): #prints all results
stats = []
for i, (train_index, test_index) in enumerate(kf):
#10-fold cross validation (9 samples for training, 1 for testing)
X_train1, X_test = X_train[train_index], X_train[test_index]
Y_train1, Y_test = Y_train[train_index], Y_train[test_index]
probas_ = pipeline.fit(X_train1,Y_train1).predict(X_test)
stats.append(accuracy_score(Y_test, probas_))
return stats
stats = []
stats1 = []
stats2 = []
stats3 = []
#Read Data
df=pd.read_csv("grids.csv")
#print df
le = preprocessing.LabelEncoder()
le.fit(df["TripId"])
Y_train=le.transform(df["TripId"])
X_train1=df['Grids']
X_train=np.array(X_train1)
vectorizer=CountVectorizer()
transformer=TfidfTransformer()
svd=TruncatedSVD(n_components=300, random_state=42)
kf = KFold(len(X_train), n_folds=10)
#knn
clf=KNeighborsClassifier(n_neighbors=7,n_jobs=-1) #best exper for KNN
pipeline = Pipeline([
('vect', vectorizer),
('tfidf', transformer),
('svd',svd),
('clf', clf)
])
stats1 = compute_and_print()
#randomforest
clf=RandomForestClassifier(n_estimators=50,n_jobs=-1)
pipeline = Pipeline([
('vect', vectorizer),
('tfidf', transformer),
('svd',svd),
('clf', clf)
])
stats2 = compute_and_print()
#logistic regression
clf=LogisticRegression()
pipeline = Pipeline([
('vect', vectorizer),
('tfidf', transformer),
('svd',svd),
('clf', clf)
])
stats3 = compute_and_print()
csv_out = open('EvaluationsMetricAccuracy', 'wb')
clwriter = csv.writer(csv_out)
for i in range(1,11):
stats.append('Fold'+str(i) )
fieldnames = ['Accuracy','KNN','RandomForests','LogisticRegression']
rows = zip(stats, stats1, stats2,stats3)
clwriter.writerow(fieldnames)
clwriter.writerows(rows)
csv_out.close()