from QSP400 import dataPro as data
import numpy as np
from keras.models import Sequential
from keras.regularizers import l2
from keras.layers import LSTM, Dense, Dropout, Activation, initializers, GRU, SimpleRNN, ConvLSTM2D
from keras import optimizers
from sklearn.model_selection import train_test_split
from sklearn.ensemble import RandomForestRegressor
import matplotlib.pyplot as plt
from keras.layers import Conv1D, GlobalAveragePooling1D, MaxPooling1D, BatchNormalization, Activation, Flatten
from keras.layers import Conv2D, MaxPooling2D, Embedding, Bidirectional, LeakyReLU
from sklearn.metrics import roc_curve, auc
import numpy as np
# ac,label=data.ac()
ac_p, label = data.deal()
aac = data.fe()
ctd = data.CTD()
gaac = data.gaac()
X = np.concatenate((aac, gaac, ac_p), axis=1) # X=ac,ctd,kmer
# print(X)
def calculate_performace(test_num, pred_y, labels):
tp = 0
fp = 0
tn = 0
fn = 0
for index in range(test_num):
if labels[index] == 1:
if labels[index] == pred_y[index]:
tp = tp + 1
import shap
# import test as test
import numpy as np
import pandas as pd
from QSP400 import dataPro as data
import matplotlib.pyplot as plt
import numpy as np
from sklearn.model_selection import train_test_split
from keras.models import load_model
from keras.models import model_from_json
from keras.models import model_from_yaml
# model = open('model_architecture_1.json', 'r')
yaml_string = open('model_arthitecture_2.yaml', 'r')
model = model_from_yaml(yaml_string)
ac_p, label, ac = data.deal()
aac = data.fe()
ctd = data.CTD()
gaac = data.gaac()
# model=test.do_model()
# model = load_model('my_model.h5')
X = np.concatenate((aac, gaac, ac_p), axis=1) # X=ac,ctd,kmer
X = np.array(ac_p)
label = np.array(label)
# x_train, x_test, y_train, y_test = train_test_split(X, label, test_size=0.2, random_state=10)
x_train, x_test, y_train, y_test = train_test_split(X,
label,
test_size=0.2,
random_state=10)
# x_train = np.array([x for i, x in enumerate(X) if i % num_cross_val != fold])
# x_test = np.array([x for i, x in enumerate(X) if i % num_cross_val == fold])
# y_train = np.array([x for i, x in enumerate(label) if i % num_cross_val != fold])