{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"pygments_lexer":"ipython3","nbconvert_exporter":"python","version":"3.6.4","file_extension":".py","codemirror_mode":{"name":"ipython","version":3},"name":"python","mimetype":"text/x-python"}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"markdown","source":"## Cài đặt các thư viện cần thiết","metadata":{}},{"cell_type":"code","source":"!pip install ../input/python-datatable/datatable-0.11.0-cp37-cp37m-manylinux2010_x86_64.whl","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:01:02.440874Z","iopub.execute_input":"2022-05-27T07:01:02.441546Z","iopub.status.idle":"2022-05-27T07:01:32.538371Z","shell.execute_reply.started":"2022-05-27T07:01:02.441508Z","shell.execute_reply":"2022-05-27T07:01:32.537109Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import pandas as pd\nimport datatable as dt\nimport numpy as np\nfrom matplotlib import pyplot as plt\nimport seaborn as sns\nimport os\n\nimport riiideducation\n\nfrom sklearn.metrics import roc_auc_score,mean_squared_error,roc_curve\nfrom sklearn.preprocessing import LabelEncoder\nfrom sklearn.preprocessing import StandardScaler\n\nimport tensorflow as tf\nfrom tensorflow.keras import Sequential\nfrom tensorflow.keras.layers import BatchNormalization,Dropout,Dense,Flatten,Conv1D, MaxPooling1D,LSTM,GRU\nfrom tensorflow.keras.optimizers import Adam\nfrom keras.metrics import BinaryAccuracy,MeanSquaredError\nfrom keras import backend as K","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:06:47.849415Z","iopub.execute_input":"2022-05-27T07:06:47.849818Z","iopub.status.idle":"2022-05-27T07:06:47.858919Z","shell.execute_reply.started":"2022-05-27T07:06:47.849779Z","shell.execute_reply":"2022-05-27T07:06:47.857606Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import warnings\nwarnings.filterwarnings(\"ignore\")","metadata":{"_uuid":"8f2839f25d086af736a60e9eeb907d3b93b6e0e5","_cell_guid":"b1076dfc-b9ad-4769-8c92-a6c4dae69d19","execution":{"iopub.status.busy":"2022-05-27T07:02:24.945147Z","iopub.execute_input":"2022-05-27T07:02:24.945504Z","iopub.status.idle":"2022-05-27T07:02:24.949863Z","shell.execute_reply.started":"2022-05-27T07:02:24.945474Z","shell.execute_reply":"2022-05-27T07:02:24.949086Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Kiểm tra dữ liệu đang có","metadata":{}},{"cell_type":"code","source":"os.listdir('../input/riiid-test-answer-prediction')","metadata":{"_uuid":"d629ff2d2480ee46fbb7e2d37f6b5fab8052498a","_cell_guid":"79c7e3d0-c299-4dcb-8224-4455121ee9b0","execution":{"iopub.status.busy":"2022-05-27T07:02:27.405738Z","iopub.execute_input":"2022-05-27T07:02:27.406540Z","iopub.status.idle":"2022-05-27T07:02:27.414060Z","shell.execute_reply.started":"2022-05-27T07:02:27.406487Z","shell.execute_reply":"2022-05-27T07:02:27.413242Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"lectures_csv = dt.fread(\"../input/riiid-test-answer-prediction/lectures.csv\").to_pandas()\nexample_test_csv = dt.fread(\"../input/riiid-test-answer-prediction/example_test.csv\").to_pandas()\ntrain_csv = pd.read_csv(\"../input/riiid-test-answer-prediction/train.csv\", nrows=10**6, low_memory=False)\nquestions_csv = dt.fread(\"../input/riiid-test-answer-prediction/questions.csv\").to_pandas()","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:02:30.015971Z","iopub.execute_input":"2022-05-27T07:02:30.016691Z","iopub.status.idle":"2022-05-27T07:02:31.623109Z","shell.execute_reply.started":"2022-05-27T07:02:30.016638Z","shell.execute_reply":"2022-05-27T07:02:31.622045Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## Pre-Processing","metadata":{"trusted":true}},{"cell_type":"code","source":"# 0 nếu sự kiện là câu hỏi cho người dùng, 1 nếu sự kiện là người dùng đang xem bài giảng. Chỉ cần lấy câu hỏi\ntrain_csv = train_csv[train_csv.content_type_id == 0]\n# giá trị -1 là bài giảng\ntrain_csv = train_csv[train_csv.answered_correctly != -1]","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:02:44.426933Z","iopub.execute_input":"2022-05-27T07:02:44.427315Z","iopub.status.idle":"2022-05-27T07:02:44.568422Z","shell.execute_reply.started":"2022-05-27T07:02:44.427283Z","shell.execute_reply":"2022-05-27T07:02:44.567477Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"train_csv = train_csv.sort_values(['timestamp'], ascending=True).reset_index(drop = True)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:02:46.778417Z","iopub.execute_input":"2022-05-27T07:02:46.778841Z","iopub.status.idle":"2022-05-27T07:02:46.991064Z","shell.execute_reply.started":"2022-05-27T07:02:46.778803Z","shell.execute_reply":"2022-05-27T07:02:46.990128Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"train_csv.head(5)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:02:49.278810Z","iopub.execute_input":"2022-05-27T07:02:49.279434Z","iopub.status.idle":"2022-05-27T07:02:49.298964Z","shell.execute_reply.started":"2022-05-27T07:02:49.279381Z","shell.execute_reply":"2022-05-27T07:02:49.297899Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"content_mean_final = train_csv[['content_id','answered_correctly']].groupby(['content_id']).agg(['mean'])\ncontent_mean_final.columns = [\"answered_correctly_content_mean\"]","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:21.307318Z","iopub.execute_input":"2022-05-27T07:03:21.308169Z","iopub.status.idle":"2022-05-27T07:03:21.382525Z","shell.execute_reply.started":"2022-05-27T07:03:21.308124Z","shell.execute_reply":"2022-05-27T07:03:21.381407Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"user_mean_final = train_csv[['user_id','answered_correctly']].groupby(['user_id']).agg(['mean', 'sum', 'count'])\nuser_mean_final.columns = [\"answered_correctly_user_mean\", 'sum_correct', 'count']","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:23.818808Z","iopub.execute_input":"2022-05-27T07:03:23.819238Z","iopub.status.idle":"2022-05-27T07:03:23.873101Z","shell.execute_reply.started":"2022-05-27T07:03:23.819199Z","shell.execute_reply":"2022-05-27T07:03:23.872203Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"#dữ liệu cho fillna\nelapsed_time_mean_final = train_csv.prior_question_elapsed_time.mean()","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:27.056228Z","iopub.execute_input":"2022-05-27T07:03:27.056588Z","iopub.status.idle":"2022-05-27T07:03:27.067610Z","shell.execute_reply.started":"2022-05-27T07:03:27.056557Z","shell.execute_reply":"2022-05-27T07:03:27.066392Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"train_csv.drop(['timestamp', 'content_type_id'], axis=1, inplace=True)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:31.274166Z","iopub.execute_input":"2022-05-27T07:03:31.274557Z","iopub.status.idle":"2022-05-27T07:03:31.320491Z","shell.execute_reply.started":"2022-05-27T07:03:31.274522Z","shell.execute_reply":"2022-05-27T07:03:31.319482Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## Validation/Train datasets","metadata":{}},{"cell_type":"code","source":"validation = pd.DataFrame()\nfor i in range(4):\n    last_records = train_csv.drop_duplicates('user_id', keep = 'last')\n    train_csv = train_csv[~train_csv.index.isin(last_records.index)]\n    validation = validation.append(last_records)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:33.411871Z","iopub.execute_input":"2022-05-27T07:03:33.412387Z","iopub.status.idle":"2022-05-27T07:03:33.759978Z","shell.execute_reply.started":"2022-05-27T07:03:33.412354Z","shell.execute_reply":"2022-05-27T07:03:33.759151Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"X = pd.DataFrame()\nfor i in range(15):\n    last_records = train_csv.drop_duplicates('user_id', keep = 'last')\n    train_csv = train_csv[~train_csv.index.isin(last_records.index)]\n    X = X.append(last_records)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:35.430591Z","iopub.execute_input":"2022-05-27T07:03:35.431497Z","iopub.status.idle":"2022-05-27T07:03:36.718581Z","shell.execute_reply.started":"2022-05-27T07:03:35.431446Z","shell.execute_reply":"2022-05-27T07:03:36.717727Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"results_c = train_csv[['content_id','answered_correctly']].groupby(['content_id']).agg(['mean'])\nresults_c.columns = [\"answered_correctly_content_mean\"]\n\nresults_u = train_csv[['user_id','answered_correctly']].groupby(['user_id']).agg(['mean', 'sum', 'count'])\nresults_u.columns = [\"answered_correctly_user_mean\", 'sum_correct', 'count']","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:37.759013Z","iopub.execute_input":"2022-05-27T07:03:37.759696Z","iopub.status.idle":"2022-05-27T07:03:37.968500Z","shell.execute_reply.started":"2022-05-27T07:03:37.759643Z","shell.execute_reply":"2022-05-27T07:03:37.967557Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"result_time_mean = train_csv.prior_question_elapsed_time.mean()","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:39.984824Z","iopub.execute_input":"2022-05-27T07:03:39.985210Z","iopub.status.idle":"2022-05-27T07:03:39.991019Z","shell.execute_reply.started":"2022-05-27T07:03:39.985179Z","shell.execute_reply":"2022-05-27T07:03:39.990075Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"X = pd.merge(X, results_u, on=['user_id'], how=\"left\")\nX = pd.merge(X, results_c, on=['content_id'], how=\"left\")","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:42.042652Z","iopub.execute_input":"2022-05-27T07:03:42.043223Z","iopub.status.idle":"2022-05-27T07:03:42.075077Z","shell.execute_reply.started":"2022-05-27T07:03:42.043187Z","shell.execute_reply":"2022-05-27T07:03:42.073971Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"validation = pd.merge(validation, results_u, on=['user_id'], how=\"left\")\nvalidation = pd.merge(validation, results_c, on=['content_id'], how=\"left\")","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:45.395828Z","iopub.execute_input":"2022-05-27T07:03:45.396534Z","iopub.status.idle":"2022-05-27T07:03:45.414562Z","shell.execute_reply.started":"2022-05-27T07:03:45.396481Z","shell.execute_reply":"2022-05-27T07:03:45.413681Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"y = X['answered_correctly']\nX = X.drop(['answered_correctly'], axis=1)\n\ny_val = validation['answered_correctly']\nX_val = validation.drop(['answered_correctly'], axis=1)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:48.511824Z","iopub.execute_input":"2022-05-27T07:03:48.512237Z","iopub.status.idle":"2022-05-27T07:03:48.524955Z","shell.execute_reply.started":"2022-05-27T07:03:48.512200Z","shell.execute_reply":"2022-05-27T07:03:48.523842Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"X.columns","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:50.860967Z","iopub.execute_input":"2022-05-27T07:03:50.861370Z","iopub.status.idle":"2022-05-27T07:03:50.869284Z","shell.execute_reply.started":"2022-05-27T07:03:50.861338Z","shell.execute_reply":"2022-05-27T07:03:50.868051Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"lencoder = LabelEncoder()\n\nX['prior_question_had_explanation'].fillna(False, inplace = True)\nX['prior_question_had_explanation_enc'] = lencoder.fit_transform(X['prior_question_had_explanation'])\nX['answered_correctly_user_mean'].fillna(0.5,  inplace=True)\nX['answered_correctly_content_mean'].fillna(0.5,  inplace=True)\nX['sum_correct'].fillna(0, inplace = True)\nX['count'].fillna(0, inplace = True)\nX['prior_question_elapsed_time'].fillna(result_time_mean, inplace = True)\n\nX_val['prior_question_had_explanation'].fillna(False, inplace = True)\nX_val['prior_question_had_explanation_enc'] = lencoder.fit_transform(X_val['prior_question_had_explanation'])\nX_val['answered_correctly_user_mean'].fillna(0.5,  inplace=True)\nX_val['answered_correctly_content_mean'].fillna(0.5,  inplace=True)\nX_val['sum_correct'].fillna(0, inplace = True)\nX_val['count'].fillna(0, inplace = True)\nX_val['prior_question_elapsed_time'].fillna(result_time_mean, inplace = True)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:53.656337Z","iopub.execute_input":"2022-05-27T07:03:53.656694Z","iopub.status.idle":"2022-05-27T07:03:53.690681Z","shell.execute_reply.started":"2022-05-27T07:03:53.656664Z","shell.execute_reply":"2022-05-27T07:03:53.689747Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"X = X[['answered_correctly_user_mean', 'answered_correctly_content_mean', 'sum_correct', 'count',\n       'prior_question_elapsed_time','prior_question_had_explanation_enc']]\nX_val = X_val[['answered_correctly_user_mean', 'answered_correctly_content_mean', 'sum_correct', 'count',\n       'prior_question_elapsed_time','prior_question_had_explanation_enc']]","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:56.354340Z","iopub.execute_input":"2022-05-27T07:03:56.354720Z","iopub.status.idle":"2022-05-27T07:03:56.365714Z","shell.execute_reply.started":"2022-05-27T07:03:56.354689Z","shell.execute_reply":"2022-05-27T07:03:56.364462Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"scaler = StandardScaler()\nX = scaler.fit_transform(X)\nX_val = scaler.transform(X_val)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:03:59.155616Z","iopub.execute_input":"2022-05-27T07:03:59.155977Z","iopub.status.idle":"2022-05-27T07:03:59.176579Z","shell.execute_reply.started":"2022-05-27T07:03:59.155946Z","shell.execute_reply":"2022-05-27T07:03:59.175635Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### cnn","metadata":{}},{"cell_type":"code","source":"from tensorflow.keras.layers import Conv1D\nK.clear_session()\nX_train = X.reshape(X.shape[0], X.shape[1], 1)\nX_test = X_val.reshape(X_val.shape[0], X_val.shape[1], 1)\n    \nmodel=Sequential()\n\nmodel.add(Conv1D(32, 2, activation='sigmoid', input_shape=X_train[0].shape))\n\nmodel.add(Conv1D(24, 2, activation='relu', padding='causal'))\n\nmodel.add(Conv1D(128, 2, activation='relu', padding='causal'))\nmodel.add(Conv1D(256, 2, activation='relu', padding='causal'))\n\nmodel.add(Dense(64, activation='relu'))\nmodel.add(Dense(32, activation='relu'))\nmodel.add(Dense(28, activation='relu'))\nmodel.add(Flatten())\nmodel.add(Dense(1, activation='sigmoid'))\n\nmodel.compile(optimizer=Adam(learning_rate=0.01), loss='binary_crossentropy', metrics=[tf.keras.metrics.BinaryAccuracy()])","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:04:01.288745Z","iopub.execute_input":"2022-05-27T07:04:01.289466Z","iopub.status.idle":"2022-05-27T07:04:01.562352Z","shell.execute_reply.started":"2022-05-27T07:04:01.289413Z","shell.execute_reply":"2022-05-27T07:04:01.561227Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"history = model.fit(X_train, y, epochs=30, verbose=2, batch_size=50000)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:04:03.891567Z","iopub.execute_input":"2022-05-27T07:04:03.892142Z","iopub.status.idle":"2022-05-27T07:05:21.156428Z","shell.execute_reply.started":"2022-05-27T07:04:03.892093Z","shell.execute_reply":"2022-05-27T07:05:21.155490Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"y_pred = model.predict(X_test)\ny_true = np.array(y_val)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:05:35.564571Z","iopub.execute_input":"2022-05-27T07:05:35.564966Z","iopub.status.idle":"2022-05-27T07:05:36.553404Z","shell.execute_reply.started":"2022-05-27T07:05:35.564933Z","shell.execute_reply":"2022-05-27T07:05:36.552435Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"roc_auc_score(y_true, y_pred)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:06:00.804145Z","iopub.execute_input":"2022-05-27T07:06:00.804701Z","iopub.status.idle":"2022-05-27T07:06:00.820342Z","shell.execute_reply.started":"2022-05-27T07:06:00.804667Z","shell.execute_reply":"2022-05-27T07:06:00.819416Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"fpr, tpr, _ = roc_curve(y_true, y_pred)\n\nplt.plot(fpr,tpr)\nplt.ylabel('True Positive Rate')\nplt.xlabel('False Positive Rate')\nplt.show()","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:07:04.375458Z","iopub.execute_input":"2022-05-27T07:07:04.375857Z","iopub.status.idle":"2022-05-27T07:07:04.609531Z","shell.execute_reply.started":"2022-05-27T07:07:04.375825Z","shell.execute_reply":"2022-05-27T07:07:04.608376Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"mean_squared_error(y_true, y_pred)","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:07:07.215293Z","iopub.execute_input":"2022-05-27T07:07:07.215697Z","iopub.status.idle":"2022-05-27T07:07:07.223148Z","shell.execute_reply.started":"2022-05-27T07:07:07.215668Z","shell.execute_reply":"2022-05-27T07:07:07.221947Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"## Prediction","metadata":{}},{"cell_type":"code","source":"env = riiideducation.make_env()","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:07:09.155282Z","iopub.execute_input":"2022-05-27T07:07:09.155887Z","iopub.status.idle":"2022-05-27T07:07:09.161541Z","shell.execute_reply.started":"2022-05-27T07:07:09.155839Z","shell.execute_reply":"2022-05-27T07:07:09.160693Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"iter_test = env.iter_test()","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:07:10.605911Z","iopub.execute_input":"2022-05-27T07:07:10.606348Z","iopub.status.idle":"2022-05-27T07:07:10.611732Z","shell.execute_reply.started":"2022-05-27T07:07:10.606311Z","shell.execute_reply":"2022-05-27T07:07:10.610905Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"for (test_df, sample_prediction_df) in iter_test:\n    test_df = pd.merge(test_df, user_mean_final, on=['user_id'],  how=\"left\")\n    test_df = pd.merge(test_df, content_mean_final, on=['content_id'],  how=\"left\")\n    \n    test_df['answered_correctly_user_mean'].fillna(0.5,  inplace=True)\n    test_df['answered_correctly_content_mean'].fillna(0.5,  inplace=True)\n    test_df['sum_correct'].fillna(0, inplace=True)\n    test_df['count'].fillna(0, inplace=True)\n    test_df['prior_question_elapsed_time'].fillna(elapsed_time_mean_final, inplace = True)\n    test_df['prior_question_had_explanation'].fillna(False, inplace=True)\n    test_df[\"prior_question_had_explanation_enc\"] = lencoder.transform(test_df[\"prior_question_had_explanation\"])\n\n    # fit transform cnn\n    X = scaler.transform(test_df[['answered_correctly_user_mean', 'answered_correctly_content_mean', 'sum_correct', 'count',\n                                  'prior_question_elapsed_time', 'prior_question_had_explanation_enc']])\n    test_df['answered_correctly'] = model.predict(X.reshape(X.shape[0], X.shape[1], 1))\n    \n    env.predict(test_df.loc[test_df['content_type_id'] == 0, ['row_id', 'answered_correctly']])","metadata":{"execution":{"iopub.status.busy":"2022-05-27T07:07:19.800012Z","iopub.execute_input":"2022-05-27T07:07:19.800672Z","iopub.status.idle":"2022-05-27T07:07:20.463730Z","shell.execute_reply.started":"2022-05-27T07:07:19.800634Z","shell.execute_reply":"2022-05-27T07:07:20.462830Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"","metadata":{"trusted":true},"execution_count":null,"outputs":[]}]}