-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathAggragate.py
More file actions
133 lines (103 loc) · 4.34 KB
/
Aggragate.py
File metadata and controls
133 lines (103 loc) · 4.34 KB
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
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
import pickle
from sklearn import tree
import warnings
from sklearn import preprocessing
from sklearn.ensemble import RandomForestClassifier
from sklearn import datasets
from skimage.io import imread
from sklearn.neighbors import KNeighborsClassifier
from sklearn import svm
import numpy as np
import pandas as pd
import math
from os import listdir
from os.path import isfile,join
data = pd.read_csv('Phones_gyroscope.csv')
target = data['gt']
lab_enc = preprocessing.LabelEncoder()
target = lab_enc.fit_transform(target)
testtarget = target[600000:3205431]#size of test set. Previously set to 5000. I tried to see how long 2.5 million records would take, it took too long.
del data['Index']
del data['User']
del data['Model']
del data['Device']
del data['gt']
del data['Arrival_Time']
del data['Creation_Time']
print("Data imported")
result=""
def train_test(size):
with warnings.catch_warnings():
warnings.simplefilter("ignore")
#set up lists to hold prediction
testset=data[600000:3205431]
new_model=0
dt_pre_B=[]
rtree_pre_A=[]
dt_pre_A =[]
rtree_pre_B=[]
dt_pre_all=[]
rtree_pre_all=[]
trainset=data[:size]
#find our mid point
mid = int(math.floor(size/2))
#create subset A
Adata=data[:mid]
Atarget=target[:mid]
#create subset B
Bdata = data[mid:size]
Btarget = target[mid:size]
nn = KNeighborsClassifier(n_neighbors=3)
dt = tree.DecisionTreeClassifier(max_depth=5)
rtree = RandomForestClassifier(max_depth=5)
print("Begin first round training")
#first round training
dt.fit(Adata,Atarget)
rtree.fit(Bdata,Btarget)
#now we have them predict the other half
for i in range(0,mid):
dt_pre_B.append(dt.predict(Bdata.iloc[i])[0])
rtree_pre_A.append(rtree.predict(Adata.iloc[i])[0])
#now we swap the subsets and classifiers
dt = tree.DecisionTreeClassifier(max_depth=5)
rtree = RandomForestClassifier(max_depth=5)
print("Begin second round training")
#second round of training
dt.fit(Bdata,Btarget)
rtree.fit(Adata,Atarget)
#again, we predict the other half
for i in range(0,mid):
dt_pre_A.append(dt.predict(Adata.iloc[i])[0])
rtree_pre_B.append(rtree.predict(Bdata.iloc[i])[0])
dt_pre_all=dt_pre_A
rtree_pre_all=rtree_pre_A
for i in range(0,(mid)):
dt_pre_all.append(dt_pre_B[i])
rtree_pre_all.append(rtree_pre_B[i])
rtree_pre_all = np.asarray(rtree_pre_all)
dt_pre_all = np.asarray(dt_pre_all)
trainset['rtree'] = pd.Series(rtree_pre_all, index=trainset.index)
trainset['dt'] = pd.Series(dt_pre_all, index=trainset.index)
print("Agraget Training set created")
nn.fit(trainset,target[:size])
print("Agraget Model Trained")
dt_pre_test=[]
rtree_pre_test = []
#begin providing test data to measure accuracy
for i in range(0,2605431):#this needs to match the size of the test set
current = i
dt_pre_test.append(dt.predict(testset.iloc[current])[0])#just like in training we predict with dt and rtree
rtree_pre_test.append(rtree.predict(testset.iloc[current])[0])
dt_pre_test = np.asarray(dt_pre_test)
rtree_pre_test = np.asarray(rtree_pre_test)
testset['rtree'] = pd.Series(rtree_pre_test, index=testset.index)
testset['dt'] = pd.Series(dt_pre_test, index=testset.index)#add the predictions as features
print("Agraget test set generated")
for i in range(0,2605431):#now we feed it into the final algorithm
if testtarget[current] == nn.predict(testset.iloc[current]):
new_model = new_model +1
print("Training set Size: "+str(size))
print("Correct: " + str(new_model))
return(str(size) + "," + str(new_model))
result= ""+ str(train_test(500000)) #+ str(train_test(750000)) + str(train_test(1000000)) + str(train_test(1250000)) + str(train_test(1500000)) + str(train_test(1750000)) + str(train_test(2000000)) + str(train_test(2250000)) + str(train_test(2500000)) + str(train_test(2750000)) +str(train_test(3000000))
print(result)