{"id":256,"date":"2012-07-09T02:10:08","date_gmt":"2012-07-09T02:10:08","guid":{"rendered":"http:\/\/blog.espol.edu.ec\/nemo\/?p=256"},"modified":"2012-07-09T18:15:09","modified_gmt":"2012-07-09T18:15:09","slug":"la-estadistica-en-el-hallazgo-del-boson-de-higgs-y-el-significado-del-5-sigma","status":"publish","type":"post","link":"https:\/\/blog.espol.edu.ec\/nemo\/2012\/07\/09\/la-estadistica-en-el-hallazgo-del-boson-de-higgs-y-el-significado-del-5-sigma\/","title":{"rendered":"La estad\u00edstica en el hallazgo del Boson de Higgs y el significado del 5-sigma"},"content":{"rendered":"<p>La noticia del hallazgo del Boson de Higgs fue dada en t\u00e9rminos de 5-sigma, y he le\u00eddo en varios post que hay mucha confusi\u00f3n sobre que esto significa (se habla de que 5-sigma significa 95% de confianza, 99,5% y cosas as\u00ed), es por ello que en el presente post espero explicar un poco de que se trata el experimento y que significa dicha forma de presentar el resultado.<\/p>\n<p>Para entender la terminolog\u00eda primero se debe entender el experimento realizado, no soy f\u00edsico as\u00ed que si alguien encuentra un error en lo que voy a describir, que me lo haga saber y lo corregir\u00e9 al instante. Lo primero que debo resaltar es que el Boson de Higgs no se ha detectado directamente en el CERN, el experimento consiste en confirmar su existencia de manera indirecta, seg\u00fan la teor\u00eda una de las formas en que se produce el bos\u00f3n de Higgs es a trav\u00e9s de la fusi\u00f3n de los Gluones, estos producen un Boson de Higgs e\u00a0inmediatamente este\u00a0se desintegra\u00a0en Quarks y lo que se mide es la energ\u00eda producida en este proceso, con estos datos se puede armar un gr\u00e1fico de N\u00famero de Fusiones o Colisiones versus Energ\u00eda producida; en paralelo se tiene la curva para las colisiones que se sabe que no involucran al Boson de Higgs.<\/p>\n<div style=\"width: 416px\" class=\"wp-caption aligncenter\"><a href=\"http:\/\/www.phdcomics.com\/higgs\/index.php?page=10\"><img loading=\"lazy\" decoding=\"async\" class=\"  \" src=\"http:\/\/www.phdcomics.com\/higgs\/darkmatters2_10s.gif\" alt=\"Proceso de detecci\u00f3n del Boson de Higgs (explicado por PHDcomics)\" width=\"416\" height=\"546\" \/><\/a><p class=\"wp-caption-text\">Proceso de detecci\u00f3n del Boson de Higgs (explicado por PHDcomics)<\/p><\/div>\n<p>Ok, entonces entendido (mas o menos) el experimento vamos a la parte estad\u00edstica que se resume en determinar si el n\u00famero de collisiones al nivel de energ\u00eda <em>\"x\"<\/em>es diferente a lo que deber\u00eda ser si no hubiese producci\u00f3n de\u00a0Bosones de Higgs, y aqu\u00ed entra el famoso 5-sigma.<\/p>\n<p>Suponiendo que no hay producci\u00f3n de Bosones de Higgs en todas las fusiones realizadas, se debe cumplir que al nivel\u00a0<em>\"x\"<\/em>\u00a0de energ\u00eda el n\u00famero de colisiones detectadas\u00a0siguen una distribuci\u00f3n normal, as\u00ed\u00a0este enfoque le permite a uno calcular la probabilidad de un falso positivo. Lo que nos dicen es:<\/p>\n<blockquote><p>\u201cObservamos en nuestros datos se\u00f1ales claras de una nueva part\u00edcula, al nivel de 5 sigma, en la regi\u00f3n de masas de alrededor de 126 gigaelectronvoltios [unas 134 veces la masa de un prot\u00f3n]\u201d<\/p><\/blockquote>\n<div>Entonces, \u00bfqu\u00e9 quieren decir los f\u00edsicos presentando resultados \"al nivel 5-sigma\"? Bueno, digamos que es s\u00f3lo otra convenci\u00f3n de informar los valores alfa, sabemos que sigma se refiere a la desviaci\u00f3n est\u00e1ndar de la poblaci\u00f3n, y 5-sigma significa que aceptan eventos como estad\u00edsticamente\u00a0significativos\u00a0si estos caen a una distancia mayor a 5 veces sigma con respecto a la media, es decir que <strong>el 5-sigma no se refiere al nivel de confianza de que lo descubierto sea un Boson de Higgs, sino que es el nivel de confianza de que han encontrado algo diferente de lo que se esperaba en el experimento si no hubiese producci\u00f3n de Boson de Higgs.<\/strong><\/div>\n<div><\/div>\n<div>Se puede obtener el valores de p para el 5-sigma en R haciendo:<\/div>\n<div style=\"overflow: auto\">\n<div class=\"geshifilter\">\n<pre style=\"font-family: monospace\"><span>&gt;<\/span> <span style=\"color: #cc66cc\">1<\/span> <span>-<\/span> <a href=\"http:\/\/inside-r.org\/r-doc\/stats\/pnorm\"><span style=\"color: #003399;font-weight: bold\">pnorm<\/span><\/a><span style=\"color: #009900\">(<\/span><span style=\"color: #cc66cc\">5<\/span><span style=\"color: #009900\">)<\/span>\r\n<span style=\"color: #009900\">[<\/span><span style=\"color: #cc66cc\">1<\/span><span style=\"color: #009900\">]<\/span> 2.866516e<span>-<\/span>07<\/pre>\n<\/div>\n<\/div>\n<p><strong>esto es p = 0.0000003<\/strong>.<\/p>\n<p>&nbsp;<\/p>\n<p>Finalmente, les dejo un gr\u00e1fico para entender un poco el 5-sigma en una Normal Est\u00e1ndar (el c\u00f3digo R para generarlo se muestra debajo del mismo)<\/p>\n<p>&nbsp;<\/p>\n<div id=\"attachment_284\" style=\"width: 457px\" class=\"wp-caption aligncenter\"><a href=\"http:\/\/blog.espol.edu.ec\/nemo\/files\/2012\/07\/5-sigma.png\"><img loading=\"lazy\" decoding=\"async\" aria-describedby=\"caption-attachment-284\" class=\" wp-image-284 \" src=\"http:\/\/blog.espol.edu.ec\/nemo\/files\/2012\/07\/5-sigma.png\" alt=\"5-sigma en una Normal Est\u00e1ndar\" width=\"457\" height=\"361\" srcset=\"https:\/\/blog.espol.edu.ec\/nemo\/files\/2012\/07\/5-sigma.png 762w, https:\/\/blog.espol.edu.ec\/nemo\/files\/2012\/07\/5-sigma-300x237.png 300w\" sizes=\"auto, (max-width: 457px) 100vw, 457px\" \/><\/a><p id=\"caption-attachment-284\" class=\"wp-caption-text\">5-sigma en una Normal Est\u00e1ndar<\/p><\/div>\n<div style=\"overflow: auto\">\n<div class=\"geshifilter\">\n<pre style=\"font-family: monospace\"><a href=\"http:\/\/inside-r.org\/r-doc\/base\/library\"><span style=\"color: #003399;font-weight: bold\">library<\/span><\/a><span style=\"color: #009900\">(<\/span><a href=\"http:\/\/inside-r.org\/packages\/cran\/ggplot2\"><span>ggplot2<\/span><\/a><span style=\"color: #009900\">)<\/span>\r\nx <span>&lt;-<\/span> <a href=\"http:\/\/inside-r.org\/r-doc\/base\/seq\"><span style=\"color: #003399;font-weight: bold\">seq<\/span><\/a><span style=\"color: #009900\">(<\/span><span>-<\/span><span style=\"color: #cc66cc\">6<\/span><span style=\"color: #339933\">,<\/span><span style=\"color: #cc66cc\">6<\/span><span style=\"color: #339933\">,<\/span><a href=\"http:\/\/inside-r.org\/r-doc\/base\/length\"><span style=\"color: #003399;font-weight: bold\">length<\/span><\/a>=<span style=\"color: #cc66cc\">1000<\/span><span style=\"color: #009900\">)<\/span>  <span style=\"color: #666666;font-style: italic\"># valores de \"x\"<\/span>\r\ny <span>&lt;-<\/span> <a href=\"http:\/\/inside-r.org\/r-doc\/stats\/dnorm\"><span style=\"color: #003399;font-weight: bold\">dnorm<\/span><\/a><span style=\"color: #009900\">(<\/span>x<span style=\"color: #009900\">)<\/span>               <span style=\"color: #666666;font-style: italic\"># distribucion Normal<\/span>\r\n<a href=\"http:\/\/inside-r.org\/r-doc\/stats\/df\"><span style=\"color: #003399;font-weight: bold\">df<\/span><\/a> <span>&lt;-<\/span> <a href=\"http:\/\/inside-r.org\/r-doc\/base\/data.frame\"><span style=\"color: #003399;font-weight: bold\">data.frame<\/span><\/a><span style=\"color: #009900\">(<\/span><span style=\"color: #0000ff\">\"x\"<\/span> = x<span style=\"color: #339933\">,<\/span><span style=\"color: #0000ff\">\"y\"<\/span> = y<span style=\"color: #009900\">)<\/span> \r\n\r\n<span style=\"color: #666666;font-style: italic\"># plot<\/span>\r\ntext_block <span>&lt;-<\/span> <span style=\"color: #0000ff\">\"Al nivel 5-sigma significa que <span style=\"color: #000099;font-weight: bold\">\\n<\/span>encontraron algo diferente que <span style=\"color: #000099;font-weight: bold\">\\n<\/span>ca\u00eda a 5-sigma de distancia <span style=\"color: #000099;font-weight: bold\">\\n<\/span>de la media\"<\/span>\r\n<a href=\"http:\/\/inside-r.org\/packages\/cran\/ggplot\"><span>ggplot<\/span><\/a><span style=\"color: #009900\">(<\/span><a href=\"http:\/\/inside-r.org\/r-doc\/stats\/df\"><span style=\"color: #003399;font-weight: bold\">df<\/span><\/a><span style=\"color: #339933\">,<\/span> aes<span style=\"color: #009900\">(<\/span>x<span style=\"color: #339933\">,<\/span>y<span style=\"color: #009900\">)<\/span><span style=\"color: #009900\">)<\/span> <span>+<\/span>\r\n  geom_line<span style=\"color: #009900\">(<\/span>size=<span style=\"color: #cc66cc\">1<\/span><span style=\"color: #009900\">)<\/span> <span>+<\/span>\r\n  annotate<span style=\"color: #009900\">(<\/span><span style=\"color: #0000ff\">\"text\"<\/span><span style=\"color: #339933\">,<\/span> x=<span style=\"color: #cc66cc\">1.7<\/span><span style=\"color: #339933\">,<\/span> y=<span style=\"color: #cc66cc\">0.15<\/span><span style=\"color: #339933\">,<\/span> label=text_block<span style=\"color: #339933\">,<\/span> size=<span style=\"color: #cc66cc\">4<\/span><span style=\"color: #339933\">,<\/span> hjust=<span style=\"color: #cc66cc\">0<\/span><span style=\"color: #009900\">)<\/span> <span>+<\/span>\r\n  annotate<span style=\"color: #009900\">(<\/span><span style=\"color: #0000ff\">\"segment\"<\/span><span style=\"color: #339933\">,<\/span> x = <span style=\"color: #cc66cc\">5<\/span><span style=\"color: #339933\">,<\/span> xend=<span style=\"color: #cc66cc\">5<\/span><span style=\"color: #339933\">,<\/span> y = <span style=\"color: #cc66cc\">0<\/span><span style=\"color: #339933\">,<\/span> yend = <span style=\"color: #cc66cc\">0.05<\/span><span style=\"color: #339933\">,<\/span> linetype=<span style=\"color: #0000ff\">\"dashed\"<\/span><span style=\"color: #009900\">)<\/span> <span>+<\/span>\r\n  annotate<span style=\"color: #009900\">(<\/span><span style=\"color: #0000ff\">\"text\"<\/span><span style=\"color: #339933\">,<\/span> x=<span style=\"color: #cc66cc\">5<\/span><span style=\"color: #339933\">,<\/span> y=<span style=\"color: #cc66cc\">0.05<\/span><span style=\"color: #339933\">,<\/span> label=<span style=\"color: #0000ff\">\"5-sigma\"<\/span><span style=\"color: #339933\">,<\/span> vjust=<span>-<\/span><span style=\"color: #cc66cc\">0.5<\/span><span style=\"color: #009900\">)<\/span><\/pre>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>La noticia del hallazgo del Boson de Higgs fue dada en t\u00e9rminos de 5-sigma, y he le\u00eddo en varios post que hay mucha confusi\u00f3n sobre que esto significa (se habla de que 5-sigma significa 95% de confianza, 99,5% y cosas as\u00ed), es por ello que en el presente post espero explicar un poco de que se trata el experimento y que significa dicha forma de presentar el resultado.<\/p><p><a class=\"more-link btn\" href=\"https:\/\/blog.espol.edu.ec\/nemo\/2012\/07\/09\/la-estadistica-en-el-hallazgo-del-boson-de-higgs-y-el-significado-del-5-sigma\/\">Seguir leyendo<\/a><\/p>\n","protected":false},"author":3803,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1688,771,1534],"tags":[356594,356592,356591,1171628,356596,356593,356595,1171627],"class_list":["post-256","post","type-post","status-publish","format-standard","hentry","category-estadistica","category-graficos","category-r","tag-5-sigma-significado-5-sigma","tag-boson-de-higgs","tag-el-significado-del-5-sigma-en-el-hallazgo-del-boson-de-higgs","tag-estadistica","tag-grafico-de-distribucion-normal","tag-higgs-boson","tag-la-estadistica-en-el-hallazgo-del-boson-de-higgs","tag-r","item-wrap"],"_links":{"self":[{"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/posts\/256","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/users\/3803"}],"replies":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/comments?post=256"}],"version-history":[{"count":25,"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/posts\/256\/revisions"}],"predecessor-version":[{"id":291,"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/posts\/256\/revisions\/291"}],"wp:attachment":[{"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/media?parent=256"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/categories?post=256"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/nemo\/wp-json\/wp\/v2\/tags?post=256"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}