{
  "id": 382376,
  "title": "Theoretical Best Score with only auxiliary = False",
  "url": "/competitions/icecube-neutrinos-in-deep-ice/discussion/382376",
  "author_name": "",
  "post_date": "2023-01-30T19:18:33.583299300Z",
  "votes": 32,
  "comment_count": 5,
  "views": 0,
  "content": "<p>I'd like to initiate a post to estimate the best score we could get from using only the sensors tagged with auxiliary = False. </p>\n<p>So I ran a little experiment on the first 30 000 events from the first batch. If you just do a simple line fitting, and evaluate the error of your forecast vs the true value, you should find something similar to this distribution:<br>\n<img src=\"https://i.imgur.com/Gz8Dvqd.png\" alt=\"\"></p>\n<p>We can see clearly a mixture of two distributions. <br>\nThe way I interpreted it: </p>\n<ul>\n<li>The red distribution is looking like a purely random distribution. It is centered on pi/2 and is symmetric, it looks typically like the result of a random fit.</li>\n<li>The green distribution, on the other hand, looks like the part in which the line fitted is imperfect, but meaningful. It is typically the case when a few outliers might change a bit our forecast. </li>\n</ul>\n<p>If we follow the reasoning until the end, by improving the models using only auxiliary = False we can asymptote to a best-case scenario where all the errors from the green distribution are 0 and where the errors from the red distribution are on average 1.5.</p>\n<p>A quick estimation for my 30k samples tells me there are ~70% of points that belong to the random distribution (I really did it quickly, so it might not be totally accurate).</p>\n<p>So in the best-case scenario where the error of the green distribution is always 0, we should tend to a \"best score\" somewhere around 1.1 (0.7 x pi/2)</p>",
  "messages": [
    {
      "id": "2122231",
      "postDate": "01/30/2023 19:18:33",
      "content": "<p>I'd like to initiate a post to estimate the best score we could get from using only the sensors tagged with auxiliary = False. </p>\n<p>So I ran a little experiment on the first 30 000 events from the first batch. If you just do a simple line fitting, and evaluate the error of your forecast vs the true value, you should find something similar to this distribution:<br>\n<img src=\"https://i.imgur.com/Gz8Dvqd.png\" alt=\"\"></p>\n<p>We can see clearly a mixture of two distributions. <br>\nThe way I interpreted it: </p>\n<ul>\n<li>The red distribution is looking like a purely random distribution. It is centered on pi/2 and is symmetric, it looks typically like the result of a random fit.</li>\n<li>The green distribution, on the other hand, looks like the part in which the line fitted is imperfect, but meaningful. It is typically the case when a few outliers might change a bit our forecast. </li>\n</ul>\n<p>If we follow the reasoning until the end, by improving the models using only auxiliary = False we can asymptote to a best-case scenario where all the errors from the green distribution are 0 and where the errors from the red distribution are on average 1.5.</p>\n<p>A quick estimation for my 30k samples tells me there are ~70% of points that belong to the random distribution (I really did it quickly, so it might not be totally accurate).</p>\n<p>So in the best-case scenario where the error of the green distribution is always 0, we should tend to a \"best score\" somewhere around 1.1 (0.7 x pi/2)</p>",
      "rawMarkdown": "I'd like to initiate a post to estimate the best score we could get from using only the sensors tagged with auxiliary = False. \n\nSo I ran a little experiment on the first 30 000 events from the first batch. If you just do a simple line fitting, and evaluate the error of your forecast vs the true value, you should find something similar to this distribution:\n![](https://i.imgur.com/Gz8Dvqd.png)\n\nWe can see clearly a mixture of two distributions. \nThe way I interpreted it: \n- The red distribution is looking like a purely random distribution. It is centered on pi/2 and is symmetric, it looks typically like the result of a random fit.\n- The green distribution, on the other hand, looks like the part in which the line fitted is imperfect, but meaningful. It is typically the case when a few outliers might change a bit our forecast. \n\nIf we follow the reasoning until the end, by improving the models using only auxiliary = False we can asymptote to a best-case scenario where all the errors from the green distribution are 0 and where the errors from the red distribution are on average 1.5.\n\nA quick estimation for my 30k samples tells me there are ~70% of points that belong to the random distribution (I really did it quickly, so it might not be totally accurate).\n\nSo in the best-case scenario where the error of the green distribution is always 0, we should tend to a \"best score\" somewhere around 1.1 (0.7 x pi/2)",
      "votes": null
    },
    {
      "id": "2122250",
      "postDate": "01/30/2023 19:47:17",
      "content": "<p>The green distribution is likely the \"track\" category of events, which will be predicted pretty accurately using line fitting. I think the red distribution will likely be a combination of the other types of events (e.g. cascade &amp; double-bang)</p>\n<p><img src=\"https://i.imgur.com/7Q88JH3.png\" alt=\"\"></p>\n<p>Interesting how the cascade events have a 10 degree angular resolution. This gives us a lower bound on the cascade accuracy.</p>\n<p><a href=\"https://indico.cern.ch/event/472838/contributions/1150248/attachments/1296103/1932596/CAP_2016_IceCube.pdf\" target=\"_blank\">From page 6 of this slidepack</a></p>",
      "rawMarkdown": "The green distribution is likely the \"track\" category of events, which will be predicted pretty accurately using line fitting. I think the red distribution will likely be a combination of the other types of events (e.g. cascade & double-bang)\n\n![](https://i.imgur.com/7Q88JH3.png)\n\nInteresting how the cascade events have a 10 degree angular resolution. This gives us a lower bound on the cascade accuracy.\n\n[From page 6 of this slidepack](https://indico.cern.ch/event/472838/contributions/1150248/attachments/1296103/1932596/CAP_2016_IceCube.pdf)",
      "votes": null
    },
    {
      "id": "2122263",
      "postDate": "01/30/2023 19:58:49",
      "content": "<p>Thanks for sharing, this is in fact very interesting! I join you on your conclusions regarding green/red distributions, and I guess other methodes should be investigated  for cascades and double bangs other than a simple line fitting…</p>",
      "rawMarkdown": "Thanks for sharing, this is in fact very interesting! I join you on your conclusions regarding green/red distributions, and I guess other methodes should be investigated  for cascades and double bangs other than a simple line fitting...",
      "votes": null
    },
    {
      "id": "2123588",
      "postDate": "01/31/2023 15:39:28",
      "content": "<p>My predictions also showed the same distribution as yours.<br>\nAnd I found so many predictions are pointing up; particles are estimated to be going down, while the ground truth is omni directional.<br>\nI think there are many down particles <strong>as noise events</strong> (, which is named \"Atmospheric neutrino\", maybe?)</p>\n<p><img src=\"https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2F979b32b7cd0f73327a81b00db96f06ad%2Fhist.png?generation=1675179897149169&amp;alt=media\" alt=\"histogram of ground truth\"></p>\n<p><img src=\"https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2Fbff56a7bc894b80d7667b0c238224b40%2Fhist2.jpg?generation=1675179942082220&amp;alt=media\" alt=\"histogram of predictions\"></p>",
      "rawMarkdown": "My predictions also showed the same distribution as yours.\nAnd I found so many predictions are pointing up; particles are estimated to be going down, while the ground truth is omni directional.\nI think there are many down particles **as noise events** (, which is named \"Atmospheric neutrino\", maybe?)\n\n![histogram of ground truth](https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2F979b32b7cd0f73327a81b00db96f06ad%2Fhist.png?generation=1675179897149169&alt=media)\n\n![histogram of predictions](https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2Fbff56a7bc894b80d7667b0c238224b40%2Fhist2.jpg?generation=1675179942082220&alt=media)",
      "votes": null
    },
    {
      "id": "2125198",
      "postDate": "02/01/2023 14:15:03",
      "content": "<p>Following on this topic, and thanks to datasaurus comment, I started to work on a classifier that would help detect the \"easy\" samples for which a line fit will be sufficient. <a href=\"https://www.kaggle.com/code/bowaka/icecube-simple-detector-of-muon-neutrino-type\" target=\"_blank\">Notebook Here</a>.</p>\n<p>I used only basic statistics, for now, but I managed to get already promising results. See the roc curve below:<br>\n<img src=\"https://i.imgur.com/GK8kJ5P.png\" alt=\"\"></p>\n<p>Used with a threshold low enough, such an approach can help maki a pre-segmentation of the samples. The illustration below shows the distribution of samples after a simple boosting algo forecast.<br>\n<img src=\"https://i.imgur.com/h76ZZ65.png\" alt=\"\"></p>",
      "rawMarkdown": "Following on this topic, and thanks to datasaurus comment, I started to work on a classifier that would help detect the \"easy\" samples for which a line fit will be sufficient. [Notebook Here](https://www.kaggle.com/code/bowaka/icecube-simple-detector-of-muon-neutrino-type).\n\nI used only basic statistics, for now, but I managed to get already promising results. See the roc curve below:\n![](https://i.imgur.com/GK8kJ5P.png)\n\nUsed with a threshold low enough, such an approach can help maki a pre-segmentation of the samples. The illustration below shows the distribution of samples after a simple boosting algo forecast.\n![](https://i.imgur.com/h76ZZ65.png)",
      "votes": null
    },
    {
      "id": "2126459",
      "postDate": "02/02/2023 08:38:36",
      "content": "<p>I had a look at your interesting notebook <a href=\"https://www.kaggle.com/bowaka\" target=\"_blank\">@bowaka</a>.<br>\nI think there is a delicate trade-off between accuracy and speed here. <br>\nA richer set of input variables to the muon/else classifier leads to a better AUC, but it is also(typically) slower to generate.<br>\nIn the version on your notebook it takes about 3min to generate output for 5k events -&gt; 2h for 1 batch.<br>\nIf we need to pre-classify all test events this might be an issue.<br>\nIn general finding a fast and accurate muon/else classifier seems to me the next major challenge.</p>",
      "rawMarkdown": "I had a look at your interesting notebook @bowaka.\nI think there is a delicate trade-off between accuracy and speed here. \nA richer set of input variables to the muon/else classifier leads to a better AUC, but it is also(typically) slower to generate.\nIn the version on your notebook it takes about 3min to generate output for 5k events -> 2h for 1 batch.\nIf we need to pre-classify all test events this might be an issue.\nIn general finding a fast and accurate muon/else classifier seems to me the next major challenge.",
      "votes": null
    }
  ],
  "comments": [
    {
      "id": 2122250,
      "author_name": "anjum48",
      "author_url": "",
      "post_date": "01/30/2023 19:47:17",
      "content": "<p>The green distribution is likely the \"track\" category of events, which will be predicted pretty accurately using line fitting. I think the red distribution will likely be a combination of the other types of events (e.g. cascade &amp; double-bang)</p>\n<p><img src=\"https://i.imgur.com/7Q88JH3.png\" alt=\"\"></p>\n<p>Interesting how the cascade events have a 10 degree angular resolution. This gives us a lower bound on the cascade accuracy.</p>\n<p><a href=\"https://indico.cern.ch/event/472838/contributions/1150248/attachments/1296103/1932596/CAP_2016_IceCube.pdf\" target=\"_blank\">From page 6 of this slidepack</a></p>",
      "votes": null,
      "replies": [
        {
          "id": 2122263,
          "author_name": "bowaka",
          "author_url": "",
          "post_date": "01/30/2023 19:58:49",
          "content": "<p>Thanks for sharing, this is in fact very interesting! I join you on your conclusions regarding green/red distributions, and I guess other methodes should be investigated  for cascades and double bangs other than a simple line fitting…</p>",
          "votes": null,
          "replies": []
        }
      ]
    },
    {
      "id": 2123588,
      "author_name": "yosshi999",
      "author_url": "",
      "post_date": "01/31/2023 15:39:28",
      "content": "<p>My predictions also showed the same distribution as yours.<br>\nAnd I found so many predictions are pointing up; particles are estimated to be going down, while the ground truth is omni directional.<br>\nI think there are many down particles <strong>as noise events</strong> (, which is named \"Atmospheric neutrino\", maybe?)</p>\n<p><img src=\"https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2F979b32b7cd0f73327a81b00db96f06ad%2Fhist.png?generation=1675179897149169&amp;alt=media\" alt=\"histogram of ground truth\"></p>\n<p><img src=\"https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2Fbff56a7bc894b80d7667b0c238224b40%2Fhist2.jpg?generation=1675179942082220&amp;alt=media\" alt=\"histogram of predictions\"></p>",
      "votes": null,
      "replies": []
    },
    {
      "id": 2125198,
      "author_name": "bowaka",
      "author_url": "",
      "post_date": "02/01/2023 14:15:03",
      "content": "<p>Following on this topic, and thanks to datasaurus comment, I started to work on a classifier that would help detect the \"easy\" samples for which a line fit will be sufficient. <a href=\"https://www.kaggle.com/code/bowaka/icecube-simple-detector-of-muon-neutrino-type\" target=\"_blank\">Notebook Here</a>.</p>\n<p>I used only basic statistics, for now, but I managed to get already promising results. See the roc curve below:<br>\n<img src=\"https://i.imgur.com/GK8kJ5P.png\" alt=\"\"></p>\n<p>Used with a threshold low enough, such an approach can help maki a pre-segmentation of the samples. The illustration below shows the distribution of samples after a simple boosting algo forecast.<br>\n<img src=\"https://i.imgur.com/h76ZZ65.png\" alt=\"\"></p>",
      "votes": null,
      "replies": []
    },
    {
      "id": 2126459,
      "author_name": "cesarjesusvalls",
      "author_url": "",
      "post_date": "02/02/2023 08:38:36",
      "content": "<p>I had a look at your interesting notebook <a href=\"https://www.kaggle.com/bowaka\" target=\"_blank\">@bowaka</a>.<br>\nI think there is a delicate trade-off between accuracy and speed here. <br>\nA richer set of input variables to the muon/else classifier leads to a better AUC, but it is also(typically) slower to generate.<br>\nIn the version on your notebook it takes about 3min to generate output for 5k events -&gt; 2h for 1 batch.<br>\nIf we need to pre-classify all test events this might be an issue.<br>\nIn general finding a fast and accurate muon/else classifier seems to me the next major challenge.</p>",
      "votes": null,
      "replies": []
    }
  ],
  "raw_markdown_by_id": {
    "2122231": "I'd like to initiate a post to estimate the best score we could get from using only the sensors tagged with auxiliary = False. \n\nSo I ran a little experiment on the first 30 000 events from the first batch. If you just do a simple line fitting, and evaluate the error of your forecast vs the true value, you should find something similar to this distribution:\n![](https://i.imgur.com/Gz8Dvqd.png)\n\nWe can see clearly a mixture of two distributions. \nThe way I interpreted it: \n- The red distribution is looking like a purely random distribution. It is centered on pi/2 and is symmetric, it looks typically like the result of a random fit.\n- The green distribution, on the other hand, looks like the part in which the line fitted is imperfect, but meaningful. It is typically the case when a few outliers might change a bit our forecast. \n\nIf we follow the reasoning until the end, by improving the models using only auxiliary = False we can asymptote to a best-case scenario where all the errors from the green distribution are 0 and where the errors from the red distribution are on average 1.5.\n\nA quick estimation for my 30k samples tells me there are ~70% of points that belong to the random distribution (I really did it quickly, so it might not be totally accurate).\n\nSo in the best-case scenario where the error of the green distribution is always 0, we should tend to a \"best score\" somewhere around 1.1 (0.7 x pi/2)",
    "2122250": "The green distribution is likely the \"track\" category of events, which will be predicted pretty accurately using line fitting. I think the red distribution will likely be a combination of the other types of events (e.g. cascade & double-bang)\n\n![](https://i.imgur.com/7Q88JH3.png)\n\nInteresting how the cascade events have a 10 degree angular resolution. This gives us a lower bound on the cascade accuracy.\n\n[From page 6 of this slidepack](https://indico.cern.ch/event/472838/contributions/1150248/attachments/1296103/1932596/CAP_2016_IceCube.pdf)",
    "2122263": "Thanks for sharing, this is in fact very interesting! I join you on your conclusions regarding green/red distributions, and I guess other methodes should be investigated  for cascades and double bangs other than a simple line fitting...",
    "2123588": "My predictions also showed the same distribution as yours.\nAnd I found so many predictions are pointing up; particles are estimated to be going down, while the ground truth is omni directional.\nI think there are many down particles **as noise events** (, which is named \"Atmospheric neutrino\", maybe?)\n\n![histogram of ground truth](https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2F979b32b7cd0f73327a81b00db96f06ad%2Fhist.png?generation=1675179897149169&alt=media)\n\n![histogram of predictions](https://www.googleapis.com/download/storage/v1/b/kaggle-forum-message-attachments/o/inbox%2F2026558%2Fbff56a7bc894b80d7667b0c238224b40%2Fhist2.jpg?generation=1675179942082220&alt=media)",
    "2125198": "Following on this topic, and thanks to datasaurus comment, I started to work on a classifier that would help detect the \"easy\" samples for which a line fit will be sufficient. [Notebook Here](https://www.kaggle.com/code/bowaka/icecube-simple-detector-of-muon-neutrino-type).\n\nI used only basic statistics, for now, but I managed to get already promising results. See the roc curve below:\n![](https://i.imgur.com/GK8kJ5P.png)\n\nUsed with a threshold low enough, such an approach can help maki a pre-segmentation of the samples. The illustration below shows the distribution of samples after a simple boosting algo forecast.\n![](https://i.imgur.com/h76ZZ65.png)",
    "2126459": "I had a look at your interesting notebook @bowaka.\nI think there is a delicate trade-off between accuracy and speed here. \nA richer set of input variables to the muon/else classifier leads to a better AUC, but it is also(typically) slower to generate.\nIn the version on your notebook it takes about 3min to generate output for 5k events -> 2h for 1 batch.\nIf we need to pre-classify all test events this might be an issue.\nIn general finding a fast and accurate muon/else classifier seems to me the next major challenge."
  },
  "source": "meta"
}