{"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":"!pip install lightgbm\nimport numpy as np # linear algebra\nimport pandas as pd # data processing, CSV file I/O (e.g. pd.read_csv)\nimport matplotlib.pyplot as plt\nimport os\nimport lightgbm as lgbm\nfrom sklearn.decomposition import PCA\nfrom collections import defaultdict","metadata":{"_uuid":"8f2839f25d086af736a60e9eeb907d3b93b6e0e5","_cell_guid":"b1076dfc-b9ad-4769-8c92-a6c4dae69d19","execution":{"iopub.status.busy":"2023-02-15T13:53:44.879417Z","iopub.execute_input":"2023-02-15T13:53:44.880381Z","iopub.status.idle":"2023-02-15T13:53:56.143703Z","shell.execute_reply.started":"2023-02-15T13:53:44.880307Z","shell.execute_reply":"2023-02-15T13:53:56.141798Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Load data","metadata":{}},{"cell_type":"code","source":"df = pd.read_csv('/kaggle/input/research-project-01-around-multimodal-singlecell/HVG_cite_train_test_donor_3000.csv')\ntrain_cite_input_df = pd.read_hdf(\"/kaggle/input/open-problems-multimodal/train_cite_inputs.h5\")\ntrain_cite_target_df = pd.read_hdf(\"/kaggle/input/open-problems-multimodal/train_cite_targets.h5\")\nmetadata_df = pd.read_csv(\"/kaggle/input/open-problems-multimodal/metadata.csv\")\n\nday_4_cell_ids = set(metadata_df[metadata_df[\"day\"] == 4][\"cell_id\"])\n\nday_4_test_index = train_cite_input_df.index[train_cite_input_df.index.isin(day_4_cell_ids)]\ntrain_index = train_cite_input_df.index[~train_cite_input_df.index.isin(day_4_cell_ids)]\n\n#new_features = df['seurat_v3']\nnew_feat = train_cite_input_df[train_cite_input_df.columns[train_cite_input_df.columns.isin(df['seurat_v3'])]]","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:55:05.403207Z","iopub.execute_input":"2023-02-15T13:55:05.403601Z","iopub.status.idle":"2023-02-15T13:56:09.746793Z","shell.execute_reply.started":"2023-02-15T13:55:05.403564Z","shell.execute_reply":"2023-02-15T13:56:09.745517Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"new_feat","metadata":{"execution":{"iopub.status.busy":"2023-02-15T11:45:52.500656Z","iopub.execute_input":"2023-02-15T11:45:52.501123Z","iopub.status.idle":"2023-02-15T11:45:52.577662Z","shell.execute_reply.started":"2023-02-15T11:45:52.501087Z","shell.execute_reply":"2023-02-15T11:45:52.576361Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"cell_types = [\n#    \"HSC\",\n    \"EryP\",\n#    \"BP\",\n#    \"MasP\",\n#    \"MkP\",\n#    \"Mop\",\n#    \"NeuP\",\n#    \"hidden\"\n]","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:09.749796Z","iopub.execute_input":"2023-02-15T13:56:09.750214Z","iopub.status.idle":"2023-02-15T13:56:09.755503Z","shell.execute_reply.started":"2023-02-15T13:56:09.750178Z","shell.execute_reply":"2023-02-15T13:56:09.754037Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"target_cds = [\n    \"CD36\",\n#    \"CD41\",\n#    \"CD32\",\n#    \"CD45\",\n#     \"CD88\",\n#     \"CD48\",\n#     \"CD62L\",\n#     \"CD49b\",\n#     \"CD45RA\",\n#     \"CD82\",\n#     \"CD38\",\n#     \"CD71\",\n#     \"CD115\",\n#     \"CD11a\",\n#     \"CD244\",\n]","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:09.771509Z","iopub.execute_input":"2023-02-15T13:56:09.771892Z","iopub.status.idle":"2023-02-15T13:56:09.783763Z","shell.execute_reply.started":"2023-02-15T13:56:09.771861Z","shell.execute_reply":"2023-02-15T13:56:09.782466Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"for cd_it in target_cds:\n    for cell_type_it in cell_types:\n        set_cell_type_it = set(metadata_df[metadata_df[\"cell_type\"] == cell_type_it][\"cell_id\"].tolist())\n        cell_type_test_index = day_4_test_index[day_4_test_index.isin(set_cell_type_it)]\n        cell_type_train_index = train_index[train_index.isin(set_cell_type_it)]    \n\n        del set_cell_type_it","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:09.788222Z","iopub.execute_input":"2023-02-15T13:56:09.788613Z","iopub.status.idle":"2023-02-15T13:56:09.870137Z","shell.execute_reply.started":"2023-02-15T13:56:09.788582Z","shell.execute_reply":"2023-02-15T13:56:09.868900Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"X = new_feat.loc[cell_type_train_index]\ny = train_cite_target_df.loc[cell_type_train_index][cd_it]","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:09.871620Z","iopub.execute_input":"2023-02-15T13:56:09.871962Z","iopub.status.idle":"2023-02-15T13:56:10.024559Z","shell.execute_reply.started":"2023-02-15T13:56:09.871922Z","shell.execute_reply":"2023-02-15T13:56:10.023471Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Define model & train","metadata":{}},{"cell_type":"code","source":"data = X  \nlabel = y\ntrain_data = lgbm.Dataset(data, label=label)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:10.030632Z","iopub.execute_input":"2023-02-15T13:56:10.031375Z","iopub.status.idle":"2023-02-15T13:56:10.036607Z","shell.execute_reply.started":"2023-02-15T13:56:10.031331Z","shell.execute_reply":"2023-02-15T13:56:10.035300Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"param = {\n    'boosting_type':'gbdt',\n    'num_leaves':31,\n    'max_depth': -1, \n    'learning_rate':0.1,\n    'n_estimators': 100, \n    'subsample_for_bin':200000, \n    'objective':None, \n    'class_weight':None, \n    'min_split_gain':0.0, \n    'min_child_weight':0.001, \n    'min_child_samples':20, \n    'subsample':1.0, \n    'subsample_freq':0,\n    'colsample_bytree':1.0, \n    'reg_alpha':0.0, \n    'reg_lambda':0.0, \n    'random_state':42, \n    'n_jobs':None, \n    'importance_type':'split',\n}","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:10.038115Z","iopub.execute_input":"2023-02-15T13:56:10.038594Z","iopub.status.idle":"2023-02-15T13:56:10.053311Z","shell.execute_reply.started":"2023-02-15T13:56:10.038551Z","shell.execute_reply":"2023-02-15T13:56:10.051903Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"num_round = 10\nbst = lgbm.train(param, train_data, num_round)\nbst.save_model('model_default.txt', num_iteration=bst.best_iteration)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:10.055342Z","iopub.execute_input":"2023-02-15T13:56:10.056339Z","iopub.status.idle":"2023-02-15T13:56:31.089406Z","shell.execute_reply.started":"2023-02-15T13:56:10.056271Z","shell.execute_reply":"2023-02-15T13:56:31.088304Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"model_default = lgbm.Booster(model_file='model_default.txt')","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:31.112154Z","iopub.execute_input":"2023-02-15T13:56:31.112911Z","iopub.status.idle":"2023-02-15T13:56:31.128502Z","shell.execute_reply.started":"2023-02-15T13:56:31.112872Z","shell.execute_reply":"2023-02-15T13:56:31.127375Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# LGBM FI","metadata":{}},{"cell_type":"markdown","source":"### FI lgbm trained seurat","metadata":{}},{"cell_type":"code","source":"lgbm.plot_importance(model_default, importance_type = 'split',max_num_features=80, figsize=(12,22) )","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:31.130798Z","iopub.execute_input":"2023-02-15T13:56:31.131668Z","iopub.status.idle":"2023-02-15T13:56:33.310796Z","shell.execute_reply.started":"2023-02-15T13:56:31.131621Z","shell.execute_reply":"2023-02-15T13:56:33.309554Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"lgbm_fi_seurat = model_default.feature_importance(importance_type='split', iteration=10)\nprint(lgbm_fi_seurat)\nlgbm_fi_seurat = pd.DataFrame(lgbm_fi_seurat, index = [x for x in range(len(lgbm_fi_seurat))])\nlgbm_fi_seurat['gene_id'] = X.columns.tolist()\nlgbm_fi_seurat = lgbm_fi_seurat.rename(columns={0: 'fi_value'})\nlgbm_fi_seurat = lgbm_fi_seurat.query('fi_value > 0')\n\nlgbm_fi_seurat = lgbm_fi_seurat.sort_values(by='fi_value', ascending=False)\nlgbm_fi_seurat = lgbm_fi_seurat[:100]\nlgbm_fi_seurat = lgbm_fi_seurat.reset_index(drop=True)\nlgbm_fi_seurat","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:33.312904Z","iopub.execute_input":"2023-02-15T13:56:33.314227Z","iopub.status.idle":"2023-02-15T13:56:33.351175Z","shell.execute_reply.started":"2023-02-15T13:56:33.314175Z","shell.execute_reply":"2023-02-15T13:56:33.349915Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# PCA 'seurat'","metadata":{}},{"cell_type":"code","source":"X","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:33.352677Z","iopub.execute_input":"2023-02-15T13:56:33.353079Z","iopub.status.idle":"2023-02-15T13:56:33.397308Z","shell.execute_reply.started":"2023-02-15T13:56:33.353040Z","shell.execute_reply":"2023-02-15T13:56:33.396080Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"pca = PCA()\ns = pca.fit_transform(X)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:56:53.168350Z","iopub.execute_input":"2023-02-15T13:56:53.168813Z","iopub.status.idle":"2023-02-15T13:57:04.259184Z","shell.execute_reply.started":"2023-02-15T13:56:53.168779Z","shell.execute_reply":"2023-02-15T13:57:04.257816Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"s","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:57:04.273310Z","iopub.execute_input":"2023-02-15T13:57:04.273709Z","iopub.status.idle":"2023-02-15T13:57:04.287968Z","shell.execute_reply.started":"2023-02-15T13:57:04.273674Z","shell.execute_reply":"2023-02-15T13:57:04.285579Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"s = pd.DataFrame(s, columns=X.columns)\ns","metadata":{"execution":{"iopub.status.busy":"2023-02-15T13:57:18.769012Z","iopub.execute_input":"2023-02-15T13:57:18.769512Z","iopub.status.idle":"2023-02-15T13:57:18.808926Z","shell.execute_reply.started":"2023-02-15T13:57:18.769455Z","shell.execute_reply":"2023-02-15T13:57:18.807730Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# LGBM FI seurat","metadata":{}},{"cell_type":"markdown","source":"### Make a dataset file to load to model","metadata":{}},{"cell_type":"code","source":"data = s\nlabel = y\ntrain_data = lgbm.Dataset(data, label=label)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.948307Z","iopub.status.idle":"2023-02-15T09:44:45.948793Z","shell.execute_reply.started":"2023-02-15T09:44:45.948601Z","shell.execute_reply":"2023-02-15T09:44:45.948619Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"param = {\n    'boosting_type':'gbdt',\n    'num_leaves':31,\n    'max_depth': -1, \n    'learning_rate':0.1,\n    'n_estimators': 100, \n    'subsample_for_bin':200000, \n    'objective':None, \n    'class_weight':None, \n    'min_split_gain':0.0, \n    'min_child_weight':0.001, \n    'min_child_samples':20, \n    'subsample':1.0, \n    'subsample_freq':0,\n    'colsample_bytree':1.0, \n    'reg_alpha':0.0, \n    'reg_lambda':0.0, \n    'random_state':42, \n    'n_jobs':None, \n    'importance_type':'split',\n}","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.950973Z","iopub.status.idle":"2023-02-15T09:44:45.951334Z","shell.execute_reply.started":"2023-02-15T09:44:45.951157Z","shell.execute_reply":"2023-02-15T09:44:45.951175Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"num_round = 10\nbst = lgbm.train(param, train_data, num_round)\nbst.save_model('model_pca_seurat.txt', num_iteration=bst.best_iteration)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.952774Z","iopub.status.idle":"2023-02-15T09:44:45.953159Z","shell.execute_reply.started":"2023-02-15T09:44:45.952977Z","shell.execute_reply":"2023-02-15T09:44:45.952995Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"model_pca_seurat = lgbm.Booster(model_file='/kaggle/working/model_pca_seurat.txt')","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.954347Z","iopub.status.idle":"2023-02-15T09:44:45.954762Z","shell.execute_reply.started":"2023-02-15T09:44:45.954567Z","shell.execute_reply":"2023-02-15T09:44:45.954586Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### feature importance pca_seurat","metadata":{}},{"cell_type":"code","source":"lgbm_PCA_fi_seurat = model_default.feature_importance(importance_type='split', iteration=10)\nlgbm_PCA_fi_seurat = pd.DataFrame(lgbm_PCA_fi_seurat, index = [x for x in range(len(lgbm_PCA_fi_seurat))])\nlgbm_PCA_fi_seurat['gene_id'] = X.columns.tolist()\nlgbm_PCA_fi_seurat = lgbm_PCA_fi_seurat.rename(columns={0: 'fi_value'})\nlgbm_PCA_fi_seurat = lgbm_PCA_fi_seurat.query('fi_value > 0')\nlgbm_PCA_fi_seurat = lgbm_PCA_fi_seurat.sort_values(by='fi_value', ascending=False)\nlgbm_PCA_fi_seurat = lgbm_PCA_fi_seurat[:100]\nlgbm_PCA_fi_seurat = lgbm_PCA_fi_seurat.reset_index(drop=True)\nlgbm_PCA_fi_seurat","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.955984Z","iopub.status.idle":"2023-02-15T09:44:45.956653Z","shell.execute_reply.started":"2023-02-15T09:44:45.956394Z","shell.execute_reply":"2023-02-15T09:44:45.956415Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# LGBM FI PCA.invers_transform()","metadata":{}},{"cell_type":"code","source":"ss = pca.inverse_transform(s)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.957930Z","iopub.status.idle":"2023-02-15T09:44:45.958306Z","shell.execute_reply.started":"2023-02-15T09:44:45.958121Z","shell.execute_reply":"2023-02-15T09:44:45.958140Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"ss","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.959157Z","iopub.status.idle":"2023-02-15T09:44:45.959518Z","shell.execute_reply.started":"2023-02-15T09:44:45.959327Z","shell.execute_reply":"2023-02-15T09:44:45.959344Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"data = ss\nlabel = y\ntrain_data = lgbm.Dataset(data, label=label)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.960955Z","iopub.status.idle":"2023-02-15T09:44:45.961345Z","shell.execute_reply.started":"2023-02-15T09:44:45.961148Z","shell.execute_reply":"2023-02-15T09:44:45.961167Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"param = {\n    'boosting_type':'gbdt',\n    'num_leaves':31,\n    'max_depth': -1, \n    'learning_rate':0.1,\n    'n_estimators': 100, \n    'subsample_for_bin':200000, \n    'objective':None, \n    'class_weight':None, \n    'min_split_gain':0.0, \n    'min_child_weight':0.001, \n    'min_child_samples':20, \n    'subsample':1.0, \n    'subsample_freq':0,\n    'colsample_bytree':1.0, \n    'reg_alpha':0.0, \n    'reg_lambda':0.0, \n    'random_state':42, \n    'n_jobs':None, \n    'importance_type':'split',\n}","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.962762Z","iopub.status.idle":"2023-02-15T09:44:45.963686Z","shell.execute_reply.started":"2023-02-15T09:44:45.963400Z","shell.execute_reply":"2023-02-15T09:44:45.963454Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"num_round = 10\nbst = lgbm.train(param, train_data, num_round)\nbst.save_model('model_pca_inv_t_seurat.txt', num_iteration=bst.best_iteration)","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.965080Z","iopub.status.idle":"2023-02-15T09:44:45.965497Z","shell.execute_reply.started":"2023-02-15T09:44:45.965279Z","shell.execute_reply":"2023-02-15T09:44:45.965298Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"model_pca_seurat = lgbm.Booster(model_file='/kaggle/working/model_pca_inv_t_seurat.txt')","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.966849Z","iopub.status.idle":"2023-02-15T09:44:45.967238Z","shell.execute_reply.started":"2023-02-15T09:44:45.967041Z","shell.execute_reply":"2023-02-15T09:44:45.967060Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"lgbm_PCA_inv_t_fi_seurat = model_default.feature_importance(importance_type='split', iteration=10)\nlgbm_PCA_inv_t_fi_seurat = pd.DataFrame(lgbm_PCA_inv_t_fi_seurat, index = [x for x in range(len(lgbm_PCA_inv_t_fi_seurat))])\nlgbm_PCA_inv_t_fi_seurat['gene_id'] = X.columns.tolist()\nlgbm_PCA_inv_t_fi_seurat = lgbm_PCA_inv_t_fi_seurat.rename(columns={0: 'fi_value'})\nlgbm_PCA_inv_t_fi_seurat = lgbm_PCA_inv_t_fi_seurat.query('fi_value > 0')\nlgbm_PCA_inv_t_fi_seurat = lgbm_PCA_inv_t_fi_seurat.sort_values(by='fi_value', ascending=False)\nlgbm_PCA_inv_t_fi_seurat = lgbm_PCA_inv_t_fi_seurat[:100]\nlgbm_PCA_inv_t_fi_seurat = lgbm_PCA_inv_t_fi_seurat.reset_index(drop=True)\nlgbm_PCA_inv_t_fi_seurat","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.968534Z","iopub.status.idle":"2023-02-15T09:44:45.968899Z","shell.execute_reply.started":"2023-02-15T09:44:45.968720Z","shell.execute_reply":"2023-02-15T09:44:45.968737Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"# Intersect FI's LGBM vs PCA_LGBM vs PCA_inv_t_LGBM","metadata":{}},{"cell_type":"markdown","source":"### PCA_LGBM vs PCA_inv_t_LGBM","metadata":{}},{"cell_type":"code","source":"len(set(lgbm_PCA_inv_t_fi_seurat['gene_id'])&set(lgbm_PCA_fi_seurat['gene_id']))","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.970804Z","iopub.status.idle":"2023-02-15T09:44:45.971174Z","shell.execute_reply.started":"2023-02-15T09:44:45.970992Z","shell.execute_reply":"2023-02-15T09:44:45.971011Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### LGBM vs PCA_inv_t_LGBM","metadata":{}},{"cell_type":"code","source":"len(set(lgbm_fi_seurat['gene_id'])&set(lgbm_PCA_inv_t_fi_seurat['gene_id']))","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.972357Z","iopub.status.idle":"2023-02-15T09:44:45.972779Z","shell.execute_reply.started":"2023-02-15T09:44:45.972582Z","shell.execute_reply":"2023-02-15T09:44:45.972601Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### LGBM vs PCA_LGBM","metadata":{}},{"cell_type":"code","source":"len(set(lgbm_fi_seurat['gene_id'])&set(lgbm_PCA_fi_seurat['gene_id']))","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.974042Z","iopub.status.idle":"2023-02-15T09:44:45.974408Z","shell.execute_reply.started":"2023-02-15T09:44:45.974228Z","shell.execute_reply":"2023-02-15T09:44:45.974247Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"### LGBM vs PCA_LGBM vs PCA_inv_t_LGBM","metadata":{}},{"cell_type":"code","source":"len(set(lgbm_fi_seurat['gene_id'])&set(lgbm_PCA_inv_t_fi_seurat['gene_id'])&set(lgbm_PCA_fi_seurat['gene_id']))","metadata":{"execution":{"iopub.status.busy":"2023-02-15T09:44:45.975621Z","iopub.status.idle":"2023-02-15T09:44:45.976683Z","shell.execute_reply.started":"2023-02-15T09:44:45.976395Z","shell.execute_reply":"2023-02-15T09:44:45.976447Z"},"trusted":true},"execution_count":null,"outputs":[]}]}