{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"name":"python","version":"3.10.14","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"},"kaggle":{"accelerator":"none","dataSources":[{"sourceId":81933,"databundleVersionId":9643020,"sourceType":"competition"}],"dockerImageVersionId":30786,"isInternetEnabled":false,"language":"python","sourceType":"notebook","isGpuEnabled":false}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"code","source":"# This Python 3 environment comes with many helpful analytics libraries installed\n# It is defined by the kaggle/python Docker image: https://github.com/kaggle/docker-python\n# For example, here's several helpful packages to load\n\nimport numpy as np # linear algebra\nimport pandas as pd # data processing, CSV file I/O (e.g. pd.read_csv)\n\n# Input data files are available in the read-only \"../input/\" directory\n# For example, running this (by clicking run or pressing Shift+Enter) will list all files under the input directory\n\nimport os\nfor dirname, _, filenames in os.walk('/kaggle/input'):\n    for filename in filenames:\n        print(os.path.join(dirname, filename))\n\n# You can write up to 20GB to the current directory (/kaggle/working/) that gets preserved as output when you create a version using \"Save & Run All\" \n# You can also write temporary files to /kaggle/temp/, but they won't be saved outside of the current session","metadata":{"_uuid":"8f2839f25d086af736a60e9eeb907d3b93b6e0e5","_cell_guid":"b1076dfc-b9ad-4769-8c92-a6c4dae69d19","trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:53.609711Z","iopub.execute_input":"2024-12-05T12:49:53.610109Z","iopub.status.idle":"2024-12-05T12:49:58.658625Z","shell.execute_reply.started":"2024-12-05T12:49:53.610071Z","shell.execute_reply":"2024-12-05T12:49:58.657438Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"import pandas as pd\nfrom sklearn.decomposition import PCA\nfrom sklearn.preprocessing import StandardScaler\nfrom sklearn.model_selection import train_test_split\nfrom sklearn.neural_network import MLPRegressor\nfrom sklearn.metrics import mean_squared_error, r2_score","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:58.660260Z","iopub.execute_input":"2024-12-05T12:49:58.660698Z","iopub.status.idle":"2024-12-05T12:49:59.519328Z","shell.execute_reply.started":"2024-12-05T12:49:58.660664Z","shell.execute_reply":"2024-12-05T12:49:59.518438Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"def encode_string_columns(df):\n    for column in df.columns:\n        if df[column].dtype == object:\n            codes, _ = pd.factorize(df[column])\n            df[column] = codes\n    return df\n\ndf = pd.read_csv(\"/kaggle/input/child-mind-institute-problematic-internet-use/train.csv\")\ndf_test = pd.read_csv(\"/kaggle/input/child-mind-institute-problematic-internet-use/test.csv\")\nid = df_test['id']\ndf = df.drop(columns=['id', 'sii'])\ndf_test = df_test.drop(columns=['id'])\ndf = encode_string_columns(df)\ndf.fillna(0, inplace=True)\ndf_test = encode_string_columns(df_test)\ndf_test.fillna(0, inplace=True)\n\ncommon_columns = df.columns.intersection(df_test.columns)\nprint(common_columns)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:59.520693Z","iopub.execute_input":"2024-12-05T12:49:59.521502Z","iopub.status.idle":"2024-12-05T12:49:59.640294Z","shell.execute_reply.started":"2024-12-05T12:49:59.521445Z","shell.execute_reply":"2024-12-05T12:49:59.639151Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Set up features and target\nX = df[common_columns]\ny = df.drop(columns=common_columns)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:59.642276Z","iopub.execute_input":"2024-12-05T12:49:59.642600Z","iopub.status.idle":"2024-12-05T12:49:59.650785Z","shell.execute_reply.started":"2024-12-05T12:49:59.642570Z","shell.execute_reply":"2024-12-05T12:49:59.649719Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Split data into training and validation sets\nX_train, X_val, y_train, y_val = train_test_split(X, y, test_size=0.2, random_state=42)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:59.652138Z","iopub.execute_input":"2024-12-05T12:49:59.653156Z","iopub.status.idle":"2024-12-05T12:49:59.667196Z","shell.execute_reply.started":"2024-12-05T12:49:59.653102Z","shell.execute_reply":"2024-12-05T12:49:59.666197Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Standardize features\nscaler = StandardScaler()\nX_train = scaler.fit_transform(X_train)\nX_val = scaler.transform(X_val)\nX_test = scaler.transform(df_test)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:59.668556Z","iopub.execute_input":"2024-12-05T12:49:59.668888Z","iopub.status.idle":"2024-12-05T12:49:59.701339Z","shell.execute_reply.started":"2024-12-05T12:49:59.668845Z","shell.execute_reply":"2024-12-05T12:49:59.700220Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Define the neural network model for regression\nmodel = MLPRegressor(hidden_layer_sizes=(10, 5), max_iter=1000, random_state=42)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:59.702568Z","iopub.execute_input":"2024-12-05T12:49:59.702891Z","iopub.status.idle":"2024-12-05T12:49:59.707675Z","shell.execute_reply.started":"2024-12-05T12:49:59.702860Z","shell.execute_reply":"2024-12-05T12:49:59.706708Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Train the model\nmodel.fit(X_train, y_train)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:49:59.708911Z","iopub.execute_input":"2024-12-05T12:49:59.709259Z","iopub.status.idle":"2024-12-05T12:50:08.356801Z","shell.execute_reply.started":"2024-12-05T12:49:59.709228Z","shell.execute_reply":"2024-12-05T12:50:08.355772Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Predict on validation set\ny_pred = model.predict(X_val)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.357995Z","iopub.execute_input":"2024-12-05T12:50:08.358367Z","iopub.status.idle":"2024-12-05T12:50:08.368346Z","shell.execute_reply.started":"2024-12-05T12:50:08.358335Z","shell.execute_reply":"2024-12-05T12:50:08.365253Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Predict on the test set\ny_test_predict = model.predict(X_test)\ny_test_predict_df = pd.DataFrame(y_test_predict, columns=y.columns)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.371952Z","iopub.execute_input":"2024-12-05T12:50:08.372427Z","iopub.status.idle":"2024-12-05T12:50:08.385847Z","shell.execute_reply.started":"2024-12-05T12:50:08.372382Z","shell.execute_reply":"2024-12-05T12:50:08.384564Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Concatenate the predicted values with the original test DataFrame\ndf_test_with_predictions = pd.concat([df_test.reset_index(drop=True), y_test_predict_df], axis=1)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.387451Z","iopub.execute_input":"2024-12-05T12:50:08.388171Z","iopub.status.idle":"2024-12-05T12:50:08.407322Z","shell.execute_reply.started":"2024-12-05T12:50:08.388118Z","shell.execute_reply":"2024-12-05T12:50:08.406178Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Reorder columns to match train DataFrame\ndf_test_with_predictions = df_test_with_predictions[df.columns]","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.409349Z","iopub.execute_input":"2024-12-05T12:50:08.410409Z","iopub.status.idle":"2024-12-05T12:50:08.428314Z","shell.execute_reply.started":"2024-12-05T12:50:08.410362Z","shell.execute_reply":"2024-12-05T12:50:08.426759Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"print(df_test_with_predictions.head())","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.430828Z","iopub.execute_input":"2024-12-05T12:50:08.432524Z","iopub.status.idle":"2024-12-05T12:50:08.464612Z","shell.execute_reply.started":"2024-12-05T12:50:08.432476Z","shell.execute_reply":"2024-12-05T12:50:08.463413Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"# Load the CSV file into a DataFrame\ndf = pd.read_csv(\"/kaggle/input/child-mind-institute-problematic-internet-use/train.csv\")\ndf_test = df_test_with_predictions","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.465773Z","iopub.execute_input":"2024-12-05T12:50:08.466205Z","iopub.status.idle":"2024-12-05T12:50:08.526650Z","shell.execute_reply.started":"2024-12-05T12:50:08.466149Z","shell.execute_reply":"2024-12-05T12:50:08.525475Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"df = pd.DataFrame(df)\ndf_test = pd.DataFrame(df_test)\n\nprint(df.shape)\ncol = df.columns\n\n\n\n\nprint(id)\ndf = df.drop(columns = [\"id\",'sii'])\n\nstring_columns = df.select_dtypes(include='object').columns\nprint(df)\nprint(df_test)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.528417Z","iopub.execute_input":"2024-12-05T12:50:08.528780Z","iopub.status.idle":"2024-12-05T12:50:08.564296Z","shell.execute_reply.started":"2024-12-05T12:50:08.528747Z","shell.execute_reply":"2024-12-05T12:50:08.563147Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"\n# Encode the string columns in DataFrame\nencoded_df = encode_string_columns(df)\n# print(encoded_df.head())\nencoded_df_test = encode_string_columns(df_test)\nencoded_df_test.fillna(0, inplace=True)\nencoded_df.fillna(0, inplace=True)\n\nprint(encoded_df.shape)\nprint(encoded_df_test.shape)\ndef choose_optimal_components(explained_variance_ratio, threshold=0.95):\n    cumulative_variance = explained_variance_ratio.cumsum()\n    optimal_components = (cumulative_variance >= threshold).argmax() + 1\n    return optimal_components, cumulative_variance\n\n\ndef perform_pca(data, data_2, threshold=0.95):\n    # Excluding the last column\n    features = data.iloc[:, :-1]\n    features_test = data_2.iloc[:,:-1]\n\n    # Standardize the features\n    scaler = StandardScaler()\n    standardized_data = scaler.fit_transform(features)\n    standardized_data_test = scaler.transform(features_test)\n\n    # Perform PCA without specifying the number of components\n    pca = PCA()\n    pca.fit(standardized_data)\n\n    # Determine the optimal number of components\n    optimal_components, cumulative_variance = choose_optimal_components(pca.explained_variance_ratio_, threshold)\n\n    # Perform PCA with the optimal number of components\n    pca = PCA(n_components=optimal_components)\n    principal_components = pca.fit_transform(standardized_data)\n\n    # Create a DataFrame with the principal components\n    columns = [f\"PC{i+1}\" for i in range(optimal_components)]\n    pca_df = pd.DataFrame(data=principal_components, columns=columns)\n    \n    principle_components_test= pca.transform(standardized_data_test)\n    columns_test = [f\"PC{i+1}\" for i in range(optimal_components)]\n    pca_test_df = pd.DataFrame(data = principle_components_test, columns = columns_test)\n\n    return pca_df, optimal_components, cumulative_variance, pca_test_df\n\n\n    \ndata = encoded_df\ndata_2 = encoded_df_test\n\n    # Perform PCA excluding the last column\nthreshold = 0.95  # Desired level of explained variance\npca_result, optimal_components, cumulative_variance, pca_test_df= perform_pca(data, data_2, threshold)\n\n   \n    \nprint(pca_result)\n\npca_result.fillna(0, inplace=True)\nprint(pca_result)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.565790Z","iopub.execute_input":"2024-12-05T12:50:08.566134Z","iopub.status.idle":"2024-12-05T12:50:08.943751Z","shell.execute_reply.started":"2024-12-05T12:50:08.566099Z","shell.execute_reply":"2024-12-05T12:50:08.941773Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"from sklearn.model_selection import train_test_split\nfrom sklearn.neural_network import MLPClassifier\nfrom sklearn.preprocessing import StandardScaler\nfrom sklearn.metrics import accuracy_score\nfrom sklearn.metrics import mean_squared_error, r2_score","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.944874Z","iopub.execute_input":"2024-12-05T12:50:08.945236Z","iopub.status.idle":"2024-12-05T12:50:08.956446Z","shell.execute_reply.started":"2024-12-05T12:50:08.945196Z","shell.execute_reply":"2024-12-05T12:50:08.955163Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"X = pca_result\nX.fillna(method = 'ffill', inplace = True)\ndf_test = pca_test_df\ndf_test.fillna(method = 'ffill', inplace = True)\n\ny = pd.read_csv(\"/kaggle/input/child-mind-institute-problematic-internet-use/train.csv\")\n# Features and target\n\ny = y['sii']\ny.fillna(0, inplace = True)\nprint(X)\nprint(y)\n\n# Split data into training and validation sets\nX_train, X_val, y_train, y_val = train_test_split(X, y, test_size=0.2, random_state=42)\nX_test = df_test\n\n# Standardize features\nscaler = StandardScaler()\nX_train = scaler.fit_transform(X_train)\nX_val = scaler.transform(X_val)\nX_test = scaler.transform(X_test)\n\n# Define the neural network model\nmodel = MLPClassifier(hidden_layer_sizes=(10, 5), max_iter=1000, random_state=42)\n\n# Train the model\nmodel.fit(X_train, y_train)\n\n# Predict on validation set\ny_pred = model.predict(X_val)\n\nmse = mean_squared_error(y_val, y_pred)\nr2 = r2_score(y_val, y_pred)\n\n# Evaluate the model\naccuracy = accuracy_score(y_val, y_pred)\nprint(f'Validation Accuracy: {accuracy:.2f}')\nprint(mse)\nprint(r2)\n\n\ny_test_predict = model.predict(X_test)\nprint(y_test_predict)\n\ny_test_predict = pd.DataFrame(y_test_predict.astype(int), columns = [\"sii\"])","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:08.957515Z","iopub.execute_input":"2024-12-05T12:50:08.957903Z","iopub.status.idle":"2024-12-05T12:50:12.179417Z","shell.execute_reply.started":"2024-12-05T12:50:08.957860Z","shell.execute_reply":"2024-12-05T12:50:12.178082Z"}},"outputs":[],"execution_count":null},{"cell_type":"code","source":"submission = pd.concat([id, y_test_predict], axis=1)\n\nprint(submission.head())\n\nsubmission.to_csv(\"/kaggle/working/submission.csv\", index=False)","metadata":{"trusted":true,"execution":{"iopub.status.busy":"2024-12-05T12:50:12.181153Z","iopub.execute_input":"2024-12-05T12:50:12.182639Z","iopub.status.idle":"2024-12-05T12:50:12.201819Z","shell.execute_reply.started":"2024-12-05T12:50:12.182585Z","shell.execute_reply":"2024-12-05T12:50:12.200369Z"}},"outputs":[],"execution_count":null}]}