{"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":"code","source":"import pandas as pd, numpy as np\nfrom catboost import CatBoostClassifier\nimport pickle\nimport sys","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:14.639501Z","iopub.execute_input":"2023-06-27T07:06:14.639897Z","iopub.status.idle":"2023-06-27T07:06:15.65104Z","shell.execute_reply.started":"2023-06-27T07:06:14.63982Z","shell.execute_reply":"2023-06-27T07:06:15.650041Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Load Train Data and Labels","metadata":{"papermill":{"duration":0.004542,"end_time":"2023-02-07T00:59:59.189777","exception":false,"start_time":"2023-02-07T00:59:59.185235","status":"completed"},"tags":[]}},{"cell_type":"code","source":"dtypes = {\"session_id\": 'int64',\n          \"index\": np.int16,\n          \"elapsed_time\": np.int32,\n          \"event_name\": 'category',\n          \"name\": 'category',\n          \"level\": np.int8,\n          \"page\": np.float16,\n          \"room_coor_x\": np.float16,\n          \"room_coor_y\": np.float16,\n          \"screen_coor_x\": np.float16,\n          \"screen_coor_y\": np.float16,\n          \"hover_duration\": np.float32,\n          \"text\": 'category',\n          \"fqid\": 'category',\n          \"room_fqid\": 'category',\n          \"text_fqid\": 'category',\n          \"fullscreen\": np.int8,\n          \"hq\": np.int8,\n          \"music\": np.int8,\n          \"level_group\": 'category'\n          }\nuse_col = ['session_id', 'index', 'elapsed_time', 'event_name', 'name', 'level', 'page',\n           'room_coor_x', 'room_coor_y', 'hover_duration', 'text', 'fqid', 'room_fqid', 'text_fqid', 'level_group']","metadata":{"papermill":{"duration":59.284316,"end_time":"2023-02-07T01:00:58.478743","exception":false,"start_time":"2023-02-07T00:59:59.194427","status":"completed"},"tags":[],"execution":{"iopub.status.busy":"2023-06-27T07:06:15.653741Z","iopub.execute_input":"2023-06-27T07:06:15.654124Z","iopub.status.idle":"2023-06-27T07:06:15.663535Z","shell.execute_reply.started":"2023-06-27T07:06:15.65409Z","shell.execute_reply":"2023-06-27T07:06:15.662557Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"targets = pd.read_csv('/kaggle/input/predict-student-performance-from-game-play/train_labels.csv')\ntargets['session'] = targets.session_id.apply(lambda x: int(x.split('_')[0]) )\ntargets['q'] = targets.session_id.apply(lambda x: int(x.split('_')[-1][1:]) )\nprint( targets.shape )\ntargets.head()","metadata":{"papermill":{"duration":0.598155,"end_time":"2023-02-07T01:00:59.082015","exception":false,"start_time":"2023-02-07T01:00:58.48386","status":"completed"},"tags":[],"execution":{"iopub.status.busy":"2023-06-27T07:06:15.665449Z","iopub.execute_input":"2023-06-27T07:06:15.666149Z","iopub.status.idle":"2023-06-27T07:06:17.308841Z","shell.execute_reply.started":"2023-06-27T07:06:15.666111Z","shell.execute_reply":"2023-06-27T07:06:17.307843Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"feature_df = pd.read_csv('/kaggle/input/featur/feature_sort.csv')","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.311892Z","iopub.execute_input":"2023-06-27T07:06:17.312626Z","iopub.status.idle":"2023-06-27T07:06:17.472979Z","shell.execute_reply.started":"2023-06-27T07:06:17.31257Z","shell.execute_reply":"2023-06-27T07:06:17.471849Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Feature Engineer","metadata":{"papermill":{"duration":0.005196,"end_time":"2023-02-07T01:00:59.092865","exception":false,"start_time":"2023-02-07T01:00:59.087669","status":"completed"},"tags":[]}},{"cell_type":"code","source":"def delt_time_def(df):\n    df.sort_values(by=['session_id', 'elapsed_time'], inplace=True)\n    df['d_time'] = df['elapsed_time'].diff(1)\n    df['d_time'].fillna(0, inplace=True)\n    df['delt_time'] = df['d_time'].clip(0, 103000)\n    df['delt_time_next'] = df['delt_time'].shift(-1)\n    return df","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.474581Z","iopub.execute_input":"2023-06-27T07:06:17.475275Z","iopub.status.idle":"2023-06-27T07:06:17.481951Z","shell.execute_reply.started":"2023-06-27T07:06:17.475238Z","shell.execute_reply":"2023-06-27T07:06:17.48101Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def feature_engineer(train, kol_f):\n    global kol_col, kol_col_max\n    kol_col = 20\n    kol_col_max = 11+kol_f*2\n    col = [i for i in range(0,kol_col_max)]\n    new_train = pd.DataFrame(index=train['session_id'].unique(), columns=col, dtype=np.float16)  \n    # Index\n    new_train[21] = new_train.index # \"session_id\"    \n    # FE.d_time\n    new_train[0] = train.groupby(['session_id'])['d_time'].quantile(q=0.15)\n    new_train[1] = train.groupby(['session_id'])['d_time'].quantile(q=0.25)\n    new_train[2] = train.groupby(['session_id'])['d_time'].quantile(q=0.35)\n    new_train[3] = train.groupby(['session_id'])['d_time'].quantile(q=0.45)\n    new_train[4] = train.groupby(['session_id'])['d_time'].quantile(q=0.65)\n    new_train[5] = train.groupby(['session_id'])['d_time'].quantile(q=0.75)\n    new_train[6] = train.groupby(['session_id'])['d_time'].quantile(q=0.85)\n    new_train[7] = train.groupby(['session_id'])['d_time'].quantile(q=0.95)\n    # FE.hover_duration\n    new_train[8] = train.groupby(['session_id'])['hover_duration'].agg('mean')\n    new_train[9] = train.groupby(['session_id'])['hover_duration'].agg('median')\n    new_train[10] = train.groupby(['session_id'])['hover_duration'].agg('std')\n    new_train[11] = train.groupby(['session_id'])['hover_duration'].agg('min')\n    new_train[12] = train.groupby(['session_id'])['hover_duration'].agg('max')\n    q15 = lambda x: x.quantile(0.15); q15.__name__ = \"q0.15\"\n    q25 = lambda x: x.quantile(0.25); q25.__name__ = \"q0.50\"\n    q75 = lambda x: x.quantile(0.75); q75.__name__ = \"q0.75\"\n    q85 = lambda x: x.quantile(0.85); q85.__name__ = \"q0.85\"\n    new_train[13] = train.groupby(['session_id'])['hover_duration'].agg(q15)\n    new_train[14] = train.groupby(['session_id'])['hover_duration'].agg(q25)\n    new_train[15] = train.groupby(['session_id'])['hover_duration'].agg(q75)\n    new_train[16] = train.groupby(['session_id'])['hover_duration'].agg(q85)\n    # FE.dt\n    new_train[17] = new_train[21].apply(lambda x: int(str(x)[:2])).astype(np.uint8) # \"year\"\n    new_train[18] = new_train[21].apply(lambda x: int(str(x)[2:4])+1).astype(np.uint8) # \"month\"\n    new_train[19] = new_train[21].apply(lambda x: int(str(x)[4:6])).astype(np.uint8) # \"day\"\n    new_train[20] = new_train[21].apply(lambda x: int(str(x)[6:8])).astype(np.uint8) + new_train[21].apply(lambda x: int(str(x)[8:10])).astype(np.uint8)/60\n    new_train[21] = 0\n    new_train = new_train.fillna(-1)\n    \n    return new_train","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.483597Z","iopub.execute_input":"2023-06-27T07:06:17.484385Z","iopub.status.idle":"2023-06-27T07:06:17.523069Z","shell.execute_reply.started":"2023-06-27T07:06:17.484333Z","shell.execute_reply":"2023-06-27T07:06:17.52203Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def feature_next_t(row_f, new_train, train, gran_1, gran_2, i):\n    global kol_col\n    kol_col +=1\n    col1 = row_f['col1']\n    val1 = row_f['val1']\n    maska = (train[col1] == val1)\n    if row_f['kol_col'] == 1:       \n        new_train[kol_col] = train[maska].groupby(['session_id'])['delt_time_next'].sum()\n        if gran_1:\n            kol_col +=1\n            new_train[kol_col] = train[maska].groupby(['session_id'])['delt_time'].mean()\n        if gran_2:\n            kol_col +=1\n            new_train[kol_col] = train[maska].groupby(['session_id'])['index'].count()          \n    elif row_f['kol_col'] == 2: \n        col2 = row_f['col2']\n        val2 = row_f['val2']\n        maska = maska & (train[col2] == val2)        \n        new_train[kol_col] = train[maska].groupby(['session_id'])['delt_time_next'].sum()\n        if gran_1:\n            kol_col +=1\n            new_train[kol_col] = train[maska].groupby(['session_id'])['delt_time'].mean()\n        if gran_2:\n            kol_col +=1\n            new_train[kol_col] = train[maska].groupby(['session_id'])['index'].count()\n    return new_train","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.524865Z","iopub.execute_input":"2023-06-27T07:06:17.525667Z","iopub.status.idle":"2023-06-27T07:06:17.545172Z","shell.execute_reply.started":"2023-06-27T07:06:17.525602Z","shell.execute_reply":"2023-06-27T07:06:17.543725Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def feature_next_t_otvet(row_f, new_train, train, gran_1, gran_2, i):\n    global kol_col\n    kol_col +=1\n    col1 = row_f['col1']\n    val1 = row_f['val1']\n    maska = (train[col1] == val1)\n    if row_f['kol_col'] == 1:      \n        new_train[kol_col] = train[maska]['delt_time_next'].sum()\n        if gran_1:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['delt_time'].mean()\n        if gran_2:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['index'].count()          \n    elif row_f['kol_col'] == 2: \n        col2 = row_f['col2']\n        val2 = row_f['val2']\n        maska = maska & (train[col2] == val2)        \n        new_train[kol_col] = train[maska]['delt_time_next'].sum()\n        if gran_1:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['delt_time'].mean()\n        if gran_2:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['index'].count()\n    return new_train","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.546661Z","iopub.execute_input":"2023-06-27T07:06:17.547287Z","iopub.status.idle":"2023-06-27T07:06:17.568566Z","shell.execute_reply.started":"2023-06-27T07:06:17.547244Z","shell.execute_reply":"2023-06-27T07:06:17.567512Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def experiment_feature_next_t_otvet(row_f, new_train, train, gran_1, gran_2, i):\n    global kol_col\n    kol_col +=1\n    if row_f['kol_col'] == 1: \n        maska = train[row_f['col1']] == row_f['val1']\n        new_train[kol_col] = train[maska]['delt_time_next'].sum()\n        if gran_1:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['delt_time'].mean()\n        if gran_2:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['index'].count()          \n    elif row_f['kol_col'] == 2: \n        col2 = row_f['col2']\n        val2 = row_f['val2']\n        maska = (train[col1] == val1) & (train[col2] == val2)        \n        new_train[kol_col] = train[maska]['delt_time_next'].sum()\n        if gran_1:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['delt_time'].mean()\n        if gran_2:\n            kol_col +=1\n            new_train[kol_col] = train[maska]['index'].count()\n    return new_train","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.570046Z","iopub.execute_input":"2023-06-27T07:06:17.570604Z","iopub.status.idle":"2023-06-27T07:06:17.591138Z","shell.execute_reply.started":"2023-06-27T07:06:17.57057Z","shell.execute_reply":"2023-06-27T07:06:17.589578Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def feature_quest_otvet(new_train, train, quest, kol_f):\n    global kol_col\n    kol_col = 20\n    g1 = 0.7 \n    g2 = 0.3 \n\n    feature_q = feature_df[feature_df['quest'] == quest].copy()\n    feature_q.reset_index(drop=True, inplace=True)\n    \n    gran1 = round(kol_f * g1)\n    gran2 = round(kol_f * g2)    \n    for i in range(0, kol_f):         \n        row_f = feature_q.loc[i]\n        new_train = feature_next_t_otvet(row_f, new_train, train, i < gran1, i <  gran2, i) \n    col = [i for i in range(0,kol_col+1)]\n    return new_train[col]","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.592701Z","iopub.execute_input":"2023-06-27T07:06:17.593283Z","iopub.status.idle":"2023-06-27T07:06:17.60719Z","shell.execute_reply.started":"2023-06-27T07:06:17.593244Z","shell.execute_reply":"2023-06-27T07:06:17.605742Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def feature_engineer_new(new_train, train, feature_q, kol_f):\n    g1 = 0.7 \n    g2 = 0.3 \n    gran1 = round(kol_f * g1)\n    gran2 = round(kol_f * g2)    \n    for i in range(0, kol_f): \n        row_f = feature_q.loc[i]       \n        new_train = feature_next_t(row_f, new_train, train, i < gran1, i <  gran2, i)         \n    return new_train","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.608627Z","iopub.execute_input":"2023-06-27T07:06:17.608992Z","iopub.status.idle":"2023-06-27T07:06:17.627187Z","shell.execute_reply.started":"2023-06-27T07:06:17.60896Z","shell.execute_reply":"2023-06-27T07:06:17.625726Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"def feature_quest(new_train, train, quest, kol_f):\n    global kol_col\n    kol_col = 20\n    feature_q = feature_df[feature_df['quest'] == quest].copy()\n    feature_q.reset_index(drop=True, inplace=True)\n    new_train = feature_engineer_new(new_train, train, feature_q, kol_f)\n    col = [i for i in range(0,kol_col+1)]\n    return new_train[col]","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.62909Z","iopub.execute_input":"2023-06-27T07:06:17.629595Z","iopub.status.idle":"2023-06-27T07:06:17.643151Z","shell.execute_reply.started":"2023-06-27T07:06:17.629557Z","shell.execute_reply":"2023-06-27T07:06:17.641683Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"from xgboost import XGBClassifier\ndef create_model(old_train, quests, models, list_kol_f):\n    \n    kol_quest = len(quests)\n    # ITERATE THRU QUESTIONS\n    for q in quests:\n        print('### quest ', q, end='')\n        new_train = feature_engineer(old_train, list_kol_f[q])\n        train_x = feature_quest(new_train, old_train, q, list_kol_f[q])\n        print (' ---- ', 'train_q.shape = ', train_x.shape)\n           \n        # TRAIN DATA\n        train_users = train_x.index.values\n        train_y = targets.loc[targets.q==q].set_index('session').loc[train_users]\n\n        # ---------------------------------------------------------------- #\n        # Model.CatBoostClassifier\n        # ---------------------------------------------------------------- #\n        print(\"[INFO]Proc CatBoostClassifier\")\n        cat_model = CatBoostClassifier(\n            n_estimators = 300,\n            learning_rate= 0.045,\n            depth = 6\n        )\n        cat_model.fit(train_x.astype('float32'), train_y['correct'], verbose=False)\n        models[f'{q}_catboost'] = cat_model\n        # ---------------------------------------------------------------- #\n        # Model.XGBClassifier\n        # ---------------------------------------------------------------- #\n        print(\"[INFO]Proc XGBClassifier\")\n        xgb_model = XGBClassifier(\n            learning_rate=0.02,\n            max_depth= 4,\n            n_estimators= 600,\n            gamma = 0.01, \n            alpha = 0.05,  \n        )\n        xgb_model.fit(train_x.astype('float32'), train_y['correct'])\n        models[f'{q}_xgboost'] = xgb_model\n\n    \n    return models","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.648422Z","iopub.execute_input":"2023-06-27T07:06:17.64884Z","iopub.status.idle":"2023-06-27T07:06:17.663798Z","shell.execute_reply.started":"2023-06-27T07:06:17.648804Z","shell.execute_reply":"2023-06-27T07:06:17.662215Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"models = {}\nbest_threshold = 0.62","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.668834Z","iopub.execute_input":"2023-06-27T07:06:17.669273Z","iopub.status.idle":"2023-06-27T07:06:17.683763Z","shell.execute_reply.started":"2023-06-27T07:06:17.669231Z","shell.execute_reply":"2023-06-27T07:06:17.682516Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"list_kol_f = {\n    1:140,3:110,\n    4:120, 5:220, 6:130, 7:110, 8:110, 9:100, 10:140, 11:120,\n    14: 160, 15:160, 16:130, 17:140             \n             }","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.685463Z","iopub.execute_input":"2023-06-27T07:06:17.685885Z","iopub.status.idle":"2023-06-27T07:06:17.699323Z","shell.execute_reply.started":"2023-06-27T07:06:17.685845Z","shell.execute_reply":"2023-06-27T07:06:17.698005Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Train","metadata":{}},{"cell_type":"code","source":"quests_0_4 = [1, 3] ","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.70123Z","iopub.execute_input":"2023-06-27T07:06:17.701673Z","iopub.status.idle":"2023-06-27T07:06:17.711323Z","shell.execute_reply.started":"2023-06-27T07:06:17.701622Z","shell.execute_reply":"2023-06-27T07:06:17.710135Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# df0_4 = pd.read_csv('/kaggle/input/featur/train_0_4t.csv', dtype=dtypes) \n# kol_lvl = (df0_4 .groupby(['session_id'])['level'].agg('nunique') < 5)\n# list_session = kol_lvl[kol_lvl].index\n# df0_4  = df0_4 [~df0_4 ['session_id'].isin(list_session)]\n# df0_4 = delt_time_def(df0_4)\n\n# quests_0_4 = [1, 3] \n# # list_kol_f = {1:140,3:110}\n\n# models = create_model(df0_4, quests_0_4, models, list_kol_f)\n# del df0_4","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.712821Z","iopub.execute_input":"2023-06-27T07:06:17.713194Z","iopub.status.idle":"2023-06-27T07:06:17.723737Z","shell.execute_reply.started":"2023-06-27T07:06:17.713161Z","shell.execute_reply":"2023-06-27T07:06:17.722733Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"quests_5_12 = [4, 5, 6, 7, 8, 9, 10, 11] ","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.724981Z","iopub.execute_input":"2023-06-27T07:06:17.725319Z","iopub.status.idle":"2023-06-27T07:06:17.741779Z","shell.execute_reply.started":"2023-06-27T07:06:17.725289Z","shell.execute_reply":"2023-06-27T07:06:17.74058Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# df5_12 = pd.read_csv('/kaggle/input/featur/train_5_12t.csv', dtype=dtypes)\n# kol_lvl = (df5_12.groupby(['session_id'])['level'].agg('nunique') < 8)\n# list_session = kol_lvl[kol_lvl].index\n# df5_12 = df5_12[~df5_12['session_id'].isin(list_session)]\n# df5_12 = delt_time_def(df5_12)\n# quests_5_12 = [4, 5, 6, 7, 8, 9, 10, 11] \n\n# # list_kol_f = {4:110, 5:220, 6:120, 7:110, 8:110, 9:100, 10:140, 11:120}\n\n# models = create_model(df5_12, quests_5_12, models, list_kol_f)\n# del df5_12","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.743022Z","iopub.execute_input":"2023-06-27T07:06:17.743353Z","iopub.status.idle":"2023-06-27T07:06:17.756596Z","shell.execute_reply.started":"2023-06-27T07:06:17.743324Z","shell.execute_reply":"2023-06-27T07:06:17.755692Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"quests_13_22 = [14, 15, 16, 17] ","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.757822Z","iopub.execute_input":"2023-06-27T07:06:17.758163Z","iopub.status.idle":"2023-06-27T07:06:17.771502Z","shell.execute_reply.started":"2023-06-27T07:06:17.758133Z","shell.execute_reply":"2023-06-27T07:06:17.770171Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# df13_22 = pd.read_csv('/kaggle/input/featur/train_13_22t.csv', dtype=dtypes) \n# kol_lvl = (df13_22 .groupby(['session_id'])['level'].agg('nunique') < 10)\n# list_session = kol_lvl[kol_lvl].index\n# df13_22  = df13_22 [~df13_22 ['session_id'].isin(list_session)]\n# df13_22 = delt_time_def(df13_22)\n\n# quests_13_22 = [14, 15, 16, 17] \n# # list_kol_f = {14: 160, 15:160, 16:105, 17:140}\n\n# models = create_model(df13_22, quests_13_22, models, list_kol_f)","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.773194Z","iopub.execute_input":"2023-06-27T07:06:17.773589Z","iopub.status.idle":"2023-06-27T07:06:17.783721Z","shell.execute_reply.started":"2023-06-27T07:06:17.773555Z","shell.execute_reply":"2023-06-27T07:06:17.782722Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"# # Saving a Model\n# for q in quests_0_4 + quests_5_12 + quests_13_22:\n#     models[f\"{q}_catboost\"].save_model(f'cat_model_{q}.bin')\n#     models[f\"{q}_xgboost\"].save_model(f'xgb_model_{q}.bin')","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.785263Z","iopub.execute_input":"2023-06-27T07:06:17.785658Z","iopub.status.idle":"2023-06-27T07:06:17.797951Z","shell.execute_reply.started":"2023-06-27T07:06:17.785594Z","shell.execute_reply":"2023-06-27T07:06:17.796665Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### Load Model","metadata":{}},{"cell_type":"code","source":"import xgboost as xgb\nimport joblib\n\nmodels = {}\ndir = '/kaggle/input/psp-train-cat-xgb-mlp-v0627v1/'\nfor q in quests_0_4 + quests_5_12 + quests_13_22:\n    for i in [\"_catboost\",\"_xgboost\",\"_mlp\"] :\n        # -------------------------------------------------------- #\n        # Load CatBoost\n        # -------------------------------------------------------- #\n        if i == \"_catboost\":\n            models[f\"{q}_catboost\"] = CatBoostClassifier().load_model(dir+f'cat_model_{q}.bin')\n        # -------------------------------------------------------- #\n        # Load XGBooxt\n        # -------------------------------------------------------- #\n        if i == \"_xgboost\":\n            xgb_clf = XGBClassifier()\n            xgb_booster = xgb.Booster()\n            xgb_booster.load_model(dir+f'xgb_model_{q}.bin')\n            xgb_clf._Booster = xgb_booster\n            models[f'{q}_xgboost'] = xgb_clf\n        # -------------------------------------------------------- #\n        # Load MLP\n        # -------------------------------------------------------- #\n        if i == \"_mlp\":\n            models[f'{q}_mlp'] = joblib.load(dir+f'mlp_model_{q}.joblib')","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:17.799814Z","iopub.execute_input":"2023-06-27T07:06:17.800175Z","iopub.status.idle":"2023-06-27T07:06:18.759722Z","shell.execute_reply.started":"2023-06-27T07:06:17.800129Z","shell.execute_reply":"2023-06-27T07:06:18.75871Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"models","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:18.760917Z","iopub.execute_input":"2023-06-27T07:06:18.761468Z","iopub.status.idle":"2023-06-27T07:06:18.850647Z","shell.execute_reply.started":"2023-06-27T07:06:18.761436Z","shell.execute_reply":"2023-06-27T07:06:18.849654Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"**Infer Test Data**","metadata":{}},{"cell_type":"code","source":"import jo_wilder\n\ntry:\n    jo_wilder.make_env.__called__ = False\n    env.__called__ = False\n    type(env)._state = type(type(env)._state).__dict__['INIT']\nexcept:\n    pass\n\nenv = jo_wilder.make_env()\niter_test = env.iter_test()    ","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:18.851848Z","iopub.execute_input":"2023-06-27T07:06:18.852921Z","iopub.status.idle":"2023-06-27T07:06:18.880981Z","shell.execute_reply.started":"2023-06-27T07:06:18.852882Z","shell.execute_reply":"2023-06-27T07:06:18.879704Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"import time","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:18.886048Z","iopub.execute_input":"2023-06-27T07:06:18.887962Z","iopub.status.idle":"2023-06-27T07:06:18.892424Z","shell.execute_reply.started":"2023-06-27T07:06:18.887912Z","shell.execute_reply":"2023-06-27T07:06:18.891352Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"g_end4 = 0\ng_end5 = 0\n\nlist_q = {'0-4':quests_0_4, '5-12':quests_5_12, '13-22':quests_13_22}\nfor (test, sam_sub) in iter_test:\n    sam_sub['question'] = [int(label.split('_')[1][1:]) for label in sam_sub['session_id']]    \n    grp = test.level_group.values[0]   \n    sam_sub['correct'] = 1\n    sam_sub.loc[sam_sub.question.isin([5, 8, 10, 13, 15]), 'correct'] = 0  \n    old_train = delt_time_def(test[test.level_group == grp])\n       \n    for q in list_q[grp]:\n        \n        start4 = time.time()\n        new_train = feature_engineer(old_train, list_kol_f[q])\n        new_train = feature_quest_otvet(new_train, old_train, q, list_kol_f[q])\n#         new_train = feature_quest(new_train, old_train, q, kol_f)\n        \n        end4 = time.time() - start4\n        g_end4 += end4\n        \n        start5 = time.time()        \n\n        # ----------------------------------------------------------------------------- #\n        # Blending\n        # ----------------------------------------------------------------------------- #\n        for i in [\"_catboost\",\"_xgboost\",\"_mlp\"] :\n            clf = models[f'{str(q)+i}']\n            if i == \"_catboost\" :\n                p1 = clf.predict_proba(new_train.astype('float32'))[:,1]\n            elif i == \"_xgboost\":\n                p2 = clf.predict_proba(new_train.astype('float32'))[:,1]\n            elif i == \"_mlp\" :\n                p4 = clf.predict_proba(new_train.astype('float32').fillna(0))[:,1]\n        p = 0.99*p1 + 0.009*p2 + 0.001*p4\n        \n        end5 = time.time() - start5\n        g_end5 += end5\n             \n        \n        mask = sam_sub.question == q \n        x = int(p[0]>best_threshold)\n        sam_sub.loc[mask,'correct'] = x      \n        \n        \n    sam_sub = sam_sub[['session_id', 'correct']]      \n    env.predict(sam_sub)","metadata":{"execution":{"iopub.status.busy":"2023-06-27T07:06:18.894743Z","iopub.execute_input":"2023-06-27T07:06:18.895462Z","iopub.status.idle":"2023-06-27T07:06:34.43037Z","shell.execute_reply.started":"2023-06-27T07:06:18.895423Z","shell.execute_reply":"2023-06-27T07:06:34.429376Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# EDA submission.csv","metadata":{"papermill":{"duration":0.011427,"end_time":"2023-02-07T01:02:45.502331","exception":false,"start_time":"2023-02-07T01:02:45.490904","status":"completed"},"tags":[]}},{"cell_type":"code","source":"# df = pd.read_csv('submission.csv')\n# print( df.shape )\n# df.head(60)","metadata":{"papermill":{"duration":0.027432,"end_time":"2023-02-07T01:02:45.541022","exception":false,"start_time":"2023-02-07T01:02:45.51359","status":"completed"},"tags":[],"execution":{"iopub.status.busy":"2023-06-27T07:06:34.431723Z","iopub.execute_input":"2023-06-27T07:06:34.432301Z","iopub.status.idle":"2023-06-27T07:06:34.435979Z","shell.execute_reply.started":"2023-06-27T07:06:34.432265Z","shell.execute_reply":"2023-06-27T07:06:34.434935Z"},"trusted":true},"execution_count":null,"outputs":[]}]}