Visualización y Análisis en

class: center, middle, inverse, title-slide

# Visualización y Análisis en <img height="90" width="90" class="fa" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgcm9sZT0iaW1nIgogICB2aWV3Qm94PSIwIDAgMjQgMjQiCiAgIHZlcnNpb249IjEuMSIKICAgaWQ9InN2Zzg0OSIKICAgc29kaXBvZGk6ZG9jbmFtZT0icnN0dWRpby5zdmciCiAgIGlua3NjYXBlOnZlcnNpb249IjEuMCAoNDAzNWE0ZiwgMjAyMC0wNS0wMSkiPgogIDxtZXRhZGF0YQogICAgIGlkPSJtZXRhZGF0YTg1NSI+CiAgICA8cmRmOlJERj4KICAgICAgPGNjOldvcmsKICAgICAgICAgcmRmOmFib3V0PSIiPgogICAgICAgIDxkYzpmb3JtYXQ+aW1hZ2Uvc3ZnK3htbDwvZGM6Zm9ybWF0PgogICAgICAgIDxkYzp0eXBlCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vcHVybC5vcmcvZGMvZGNtaXR5cGUvU3RpbGxJbWFnZSIgLz4KICAgICAgPC9jYzpXb3JrPgogICAgPC9yZGY6UkRGPgogIDwvbWV0YWRhdGE+CiAgPGRlZnMKICAgICBpZD0iZGVmczg1MyIgLz4KICA8c29kaXBvZGk6bmFtZWR2aWV3CiAgICAgaW5rc2NhcGU6ZG9jdW1lbnQtcm90YXRpb249IjAiCiAgICAgcGFnZWNvbG9yPSIjZmZmZmZmIgogICAgIGJvcmRlcmNvbG9yPSIjNjY2NjY2IgogICAgIGJvcmRlcm9wYWNpdHk9IjEiCiAgICAgb2JqZWN0dG9sZXJhbmNlPSIxMCIKICAgICBncmlkdG9sZXJhbmNlPSIxMCIKICAgICBndWlkZXRvbGVyYW5jZT0iMTAiCiAgICAgaW5rc2NhcGU6cGFnZW9wYWNpdHk9IjAiCiAgICAgaW5rc2NhcGU6cGFnZXNoYWRvdz0iMiIKICAgICBpbmtzY2FwZTp3aW5kb3ctd2lkdGg9IjEyMzAiCiAgICAgaW5rc2NhcGU6d2luZG93LWhlaWdodD0iOTMxIgogICAgIGlkPSJuYW1lZHZpZXc4NTEiCiAgICAgc2hvd2dyaWQ9ImZhbHNlIgogICAgIGlua3NjYXBlOnpvb209IjMwLjYyNSIKICAgICBpbmtzY2FwZTpjeD0iMTIiCiAgICAgaW5rc2NhcGU6Y3k9IjEyIgogICAgIGlua3NjYXBlOndpbmRvdy14PSIwIgogICAgIGlua3NjYXBlOndpbmRvdy15PSIyMyIKICAgICBpbmtzY2FwZTp3aW5kb3ctbWF4aW1pemVkPSIwIgogICAgIGlua3NjYXBlOmN1cnJlbnQtbGF5ZXI9InN2Zzg0OSIgLz4KICA8dGl0bGUKICAgICBpZD0idGl0bGU4NDUiPlJTdHVkaW8gaWNvbjwvdGl0bGU+CiAgPHBhdGgKICAgICBkPSJNMTIuMTc4LjAwMmExMi4wMDIgMTIuMDAyIDAgMCAwLTguNjYyIDMuNTE1IDEyLjAwMiAxMi4wMDIgMCAwIDAgMCAxNi45NyAxMi4wMDIgMTIuMDAyIDAgMCAwIDE2Ljk3IDAgMTIuMDAyIDEyLjAwMiAwIDAgMCAwLTE2Ljk3QTEyLjAwMiAxMi4wMDIgMCAwIDAgMTIuMTc5LjAwMnpNNy43NyA1Ljk5NWMuNTYyLjEyOCAxLjA1LjIxNyAxLjY2My4yMTcuOTIxIDAgMS44NjMtLjIxNyAyLjc4Ni0uMjE3IDEuNzkgMCAzLjQ1LjgxNCAzLjQ1IDIuOCAwIDEuNTQtLjkyMSAyLjUxNy0yLjM1IDIuOTNsMi43ODggNC4xMDdoMS4zMDF2MS4wMWgtMS45ODZsLTMuMjkzLTQuOTM0aC0xLjc1N3YzLjkyNGgxLjcxOHYxLjAxSDcuNzd2LTEuMDFoMS40ODNWNy4xMzRMNy43NyA2Ljk1MXYtLjk1N3ptNC40NjYgMS4wMTJjLS41OTYgMC0xLjIxMy4wNTMtMS44NjQuMTI3djMuNzk4bC45NDEuMDJjMi4yOTguMDM0IDMuMTgzLS44NSAzLjE4My0yLjAyNiAwLTEuMzc2LS45OTctMS45MTktMi4yNi0xLjkxOXoiCiAgICAgaWQ9InBhdGg4NDciCiAgICAgc3R5bGU9ImZpbGw6I2Q5NjMxZTtmaWxsLW9wYWNpdHk6MSIgLz4KPC9zdmc+Cg==" />
## Sesión II
### Javier Tamayo Leiva <a href="mailto:j.tamayo.leiva@gmail.com"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgaW5rc2NhcGU6dmVyc2lvbj0iMS4wICg0MDM1YTRmLCAyMDIwLTA1LTAxKSIKICAgc29kaXBvZGk6ZG9jbmFtZT0iZW52ZWxvcGUtc29saWQuc3ZnIgogICBpZD0ic3ZnODM1IgogICB2ZXJzaW9uPSIxLjEiCiAgIHZpZXdCb3g9IjAgMCA1MTIgNTEyIgogICByb2xlPSJpbWciCiAgIGNsYXNzPSJzdmctaW5saW5lLS1mYSBmYS1lbnZlbG9wZSBmYS13LTE2IgogICBkYXRhLWljb249ImVudmVsb3BlIgogICBkYXRhLXByZWZpeD0iZmFzIgogICBmb2N1c2FibGU9ImZhbHNlIgogICBhcmlhLWhpZGRlbj0idHJ1ZSI+CiAgPG1ldGFkYXRhCiAgICAgaWQ9Im1ldGFkYXRhODQxIj4KICAgIDxyZGY6UkRGPgogICAgICA8Y2M6V29yawogICAgICAgICByZGY6YWJvdXQ9IiI+CiAgICAgICAgPGRjOmZvcm1hdD5pbWFnZS9zdmcreG1sPC9kYzpmb3JtYXQ+CiAgICAgICAgPGRjOnR5cGUKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9wdXJsLm9yZy9kYy9kY21pdHlwZS9TdGlsbEltYWdlIiAvPgogICAgICA8L2NjOldvcms+CiAgICA8L3JkZjpSREY+CiAgPC9tZXRhZGF0YT4KICA8ZGVmcwogICAgIGlkPSJkZWZzODM5IiAvPgogIDxzb2RpcG9kaTpuYW1lZHZpZXcKICAgICBpbmtzY2FwZTpjdXJyZW50LWxheWVyPSJzdmc4MzUiCiAgICAgaW5rc2NhcGU6d2luZG93LW1heGltaXplZD0iMCIKICAgICBpbmtzY2FwZTp3aW5kb3cteT0iMjMiCiAgICAgaW5rc2NhcGU6d2luZG93LXg9IjAiCiAgICAgaW5rc2NhcGU6Y3k9IjI1NiIKICAgICBpbmtzY2FwZTpjeD0iMjU2IgogICAgIGlua3NjYXBlOnpvb209IjEuMzkwNjI1IgogICAgIHNob3dncmlkPSJmYWxzZSIKICAgICBpZD0ibmFtZWR2aWV3ODM3IgogICAgIGlua3NjYXBlOndpbmRvdy1oZWlnaHQ9Ijk1NyIKICAgICBpbmtzY2FwZTp3aW5kb3ctd2lkdGg9Ijc3NiIKICAgICBpbmtzY2FwZTpwYWdlc2hhZG93PSIyIgogICAgIGlua3NjYXBlOnBhZ2VvcGFjaXR5PSIwIgogICAgIGd1aWRldG9sZXJhbmNlPSIxMCIKICAgICBncmlkdG9sZXJhbmNlPSIxMCIKICAgICBvYmplY3R0b2xlcmFuY2U9IjEwIgogICAgIGJvcmRlcm9wYWNpdHk9IjEiCiAgICAgYm9yZGVyY29sb3I9IiM2NjY2NjYiCiAgICAgcGFnZWNvbG9yPSIjZmZmZmZmIiAvPgogIDxwYXRoCiAgICAgc3R5bGU9ImZpbGw6IzhhYjBiZjtmaWxsLW9wYWNpdHk6MSIKICAgICBpZD0icGF0aDgzMyIKICAgICBkPSJNNTAyLjMgMTkwLjhjMy45LTMuMSA5LjctLjIgOS43IDQuN1Y0MDBjMCAyNi41LTIxLjUgNDgtNDggNDhINDhjLTI2LjUgMC00OC0yMS41LTQ4LTQ4VjE5NS42YzAtNSA1LjctNy44IDkuNy00LjcgMjIuNCAxNy40IDUyLjEgMzkuNSAxNTQuMSAxMTMuNiAyMS4xIDE1LjQgNTYuNyA0Ny44IDkyLjIgNDcuNiAzNS43LjMgNzItMzIuOCA5Mi4zLTQ3LjYgMTAyLTc0LjEgMTMxLjYtOTYuMyAxNTQtMTEzLjd6TTI1NiAzMjBjMjMuMi40IDU2LjYtMjkuMiA3My40LTQxLjQgMTMyLjctOTYuMyAxNDIuOC0xMDQuNyAxNzMuNC0xMjguNyA1LjgtNC41IDkuMi0xMS41IDkuMi0xOC45di0xOWMwLTI2LjUtMjEuNS00OC00OC00OEg0OEMyMS41IDY0IDAgODUuNSAwIDExMnYxOWMwIDcuNCAzLjQgMTQuMyA5LjIgMTguOSAzMC42IDIzLjkgNDAuNyAzMi40IDE3My40IDEyOC43IDE2LjggMTIuMiA1MC4yIDQxLjggNzMuNCA0MS40eiIKICAgICBmaWxsPSJjdXJyZW50Q29sb3IiIC8+Cjwvc3ZnPgo=" /></a> <a href="https://orcid.org/0000-0003-2610-6957"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgaW5rc2NhcGU6dmVyc2lvbj0iMS4wICg0MDM1YTRmLCAyMDIwLTA1LTAxKSIKICAgc29kaXBvZGk6ZG9jbmFtZT0ib3JjaWQuc3ZnIgogICBpZD0ic3ZnODM3IgogICB2ZXJzaW9uPSIxLjEiCiAgIHZpZXdCb3g9IjAgMCAyNCAyNCIKICAgcm9sZT0iaW1nIj4KICA8bWV0YWRhdGEKICAgICBpZD0ibWV0YWRhdGE4NDMiPgogICAgPHJkZjpSREY+CiAgICAgIDxjYzpXb3JrCiAgICAgICAgIHJkZjphYm91dD0iIj4KICAgICAgICA8ZGM6Zm9ybWF0PmltYWdlL3N2Zyt4bWw8L2RjOmZvcm1hdD4KICAgICAgICA8ZGM6dHlwZQogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3B1cmwub3JnL2RjL2RjbWl0eXBlL1N0aWxsSW1hZ2UiIC8+CiAgICAgIDwvY2M6V29yaz4KICAgIDwvcmRmOlJERj4KICA8L21ldGFkYXRhPgogIDxkZWZzCiAgICAgaWQ9ImRlZnM4NDEiIC8+CiAgPHNvZGlwb2RpOm5hbWVkdmlldwogICAgIGlua3NjYXBlOmN1cnJlbnQtbGF5ZXI9InN2ZzgzNyIKICAgICBpbmtzY2FwZTp3aW5kb3ctbWF4aW1pemVkPSIwIgogICAgIGlua3NjYXBlOndpbmRvdy15PSIyMyIKICAgICBpbmtzY2FwZTp3aW5kb3cteD0iMCIKICAgICBpbmtzY2FwZTpjeT0iMTIiCiAgICAgaW5rc2NhcGU6Y3g9IjEyIgogICAgIGlua3NjYXBlOnpvb209IjI5LjY2NjY2NyIKICAgICBzaG93Z3JpZD0iZmFsc2UiCiAgICAgaWQ9Im5hbWVkdmlldzgzOSIKICAgICBpbmtzY2FwZTp3aW5kb3ctaGVpZ2h0PSI5MjgiCiAgICAgaW5rc2NhcGU6d2luZG93LXdpZHRoPSI3NzYiCiAgICAgaW5rc2NhcGU6cGFnZXNoYWRvdz0iMiIKICAgICBpbmtzY2FwZTpwYWdlb3BhY2l0eT0iMCIKICAgICBndWlkZXRvbGVyYW5jZT0iMTAiCiAgICAgZ3JpZHRvbGVyYW5jZT0iMTAiCiAgICAgb2JqZWN0dG9sZXJhbmNlPSIxMCIKICAgICBib3JkZXJvcGFjaXR5PSIxIgogICAgIGJvcmRlcmNvbG9yPSIjNjY2NjY2IgogICAgIHBhZ2Vjb2xvcj0iI2ZmZmZmZiIKICAgICBpbmtzY2FwZTpkb2N1bWVudC1yb3RhdGlvbj0iMCIgLz4KICA8dGl0bGUKICAgICBpZD0idGl0bGU4MzMiPk9SQ0lEIGljb248L3RpdGxlPgogIDxwYXRoCiAgICAgc3R5bGU9ImZpbGw6IzhhYjBiZjtmaWxsLW9wYWNpdHk6MSIKICAgICBpZD0icGF0aDgzNSIKICAgICBkPSJNMTIgMEM1LjM3MiAwIDAgNS4zNzIgMCAxMnM1LjM3MiAxMiAxMiAxMiAxMi01LjM3MiAxMi0xMlMxOC42MjggMCAxMiAwek03LjM2OSA0LjM3OGMuNTI1IDAgLjk0Ny40MzEuOTQ3Ljk0N3MtLjQyMi45NDctLjk0Ny45NDdhLjk1Ljk1IDAgMCAxLS45NDctLjk0N2MwLS41MjUuNDIyLS45NDcuOTQ3LS45NDd6bS0uNzIyIDMuMDM4aDEuNDQ0djEwLjA0MUg2LjY0N1Y3LjQxNnptMy41NjIgMGgzLjljMy43MTIgMCA1LjM0NCAyLjY1MyA1LjM0NCA1LjAyNSAwIDIuNTc4LTIuMDE2IDUuMDI1LTUuMzI1IDUuMDI1aC0zLjkxOVY3LjQxNnptMS40NDQgMS4zMDN2Ny40NDRoMi4yOTdjMy4yNzIgMCA0LjAyMi0yLjQ4NCA0LjAyMi0zLjcyMiAwLTIuMDE2LTEuMjg0LTMuNzIyLTQuMDk3LTMuNzIyaC0yLjIyMnoiIC8+Cjwvc3ZnPgo=" /></a> <a href="https://scholar.google.com/citations?user=Rr-4cmQwXX4C&hl=es"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgd2lkdGg9IjM4NCIKICAgaGVpZ2h0PSI1MTIiCiAgIHZpZXdCb3g9IjAgMCAzODQgNTEyIgogICB2ZXJzaW9uPSIxLjEiCiAgIGlkPSJzdmc4IgogICBpbmtzY2FwZTp2ZXJzaW9uPSIxLjAgKDQwMzVhNGYsIDIwMjAtMDUtMDEpIgogICBzb2RpcG9kaTpkb2NuYW1lPSJnb29nbGUtc2Nob2xhci5zdmciPgogIDx0aXRsZQogICAgIGlkPSJ0aXRsZTgzMyI+R29vZ2xlIFNjaG9sYXI8L3RpdGxlPgogIDxkZWZzCiAgICAgaWQ9ImRlZnMyIiAvPgogIDxzb2RpcG9kaTpuYW1lZHZpZXcKICAgICBpZD0iYmFzZSIKICAgICBwYWdlY29sb3I9IiNmZmZmZmYiCiAgICAgYm9yZGVyY29sb3I9IiM2NjY2NjYiCiAgICAgYm9yZGVyb3BhY2l0eT0iMS4wIgogICAgIGlua3NjYXBlOnBhZ2VvcGFjaXR5PSIwLjAiCiAgICAgaW5rc2NhcGU6cGFnZXNoYWRvdz0iMiIKICAgICBpbmtzY2FwZTp6b29tPSIwLjc1IgogICAgIGlua3NjYXBlOmN4PSIxOTIiCiAgICAgaW5rc2NhcGU6Y3k9IjI1NiIKICAgICBpbmtzY2FwZTpkb2N1bWVudC11bml0cz0icHgiCiAgICAgaW5rc2NhcGU6Y3VycmVudC1sYXllcj0ibGF5ZXI2IgogICAgIGlua3NjYXBlOmRvY3VtZW50LXJvdGF0aW9uPSIwIgogICAgIHNob3dncmlkPSJ0cnVlIgogICAgIHVuaXRzPSJweCIKICAgICBzaG93Z3VpZGVzPSJ0cnVlIgogICAgIGlua3NjYXBlOmd1aWRlLWJib3g9InRydWUiCiAgICAgaW5rc2NhcGU6d2luZG93LXdpZHRoPSIxNjgwIgogICAgIGlua3NjYXBlOndpbmRvdy1oZWlnaHQ9Ijk1OCIKICAgICBpbmtzY2FwZTp3aW5kb3cteD0iMCIKICAgICBpbmtzY2FwZTp3aW5kb3cteT0iMjMiCiAgICAgaW5rc2NhcGU6d2luZG93LW1heGltaXplZD0iMSI+CiAgICA8aW5rc2NhcGU6Z3JpZAogICAgICAgdHlwZT0ieHlncmlkIgogICAgICAgaWQ9ImdyaWQxMCIKICAgICAgIGVtcHNwYWNpbmc9IjQiIC8+CiAgICA8c29kaXBvZGk6Z3VpZGUKICAgICAgIHBvc2l0aW9uPSIzODQuMDAwMDEsNjQiCiAgICAgICBvcmllbnRhdGlvbj0iMCwxIgogICAgICAgaWQ9Imd1aWRlMjMiCiAgICAgICBpbmtzY2FwZTpsYWJlbD0iYmFzZWxpbmUiCiAgICAgICBpbmtzY2FwZTpsb2NrZWQ9InRydWUiCiAgICAgICBpbmtzY2FwZTpjb2xvcj0icmdiKDAsMCwyNTUpIiAvPgogIDwvc29kaXBvZGk6bmFtZWR2aWV3PgogIDxtZXRhZGF0YQogICAgIGlkPSJtZXRhZGF0YTUiPgogICAgPHJkZjpSREY+CiAgICAgIDxjYzpXb3JrCiAgICAgICAgIHJkZjphYm91dD0iIj4KICAgICAgICA8ZGM6Zm9ybWF0PmltYWdlL3N2Zyt4bWw8L2RjOmZvcm1hdD4KICAgICAgICA8ZGM6dHlwZQogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3B1cmwub3JnL2RjL2RjbWl0eXBlL1N0aWxsSW1hZ2UiIC8+CiAgICAgICAgPGRjOnRpdGxlPkdvb2dsZSBTY2hvbGFyPC9kYzp0aXRsZT4KICAgICAgICA8Y2M6bGljZW5zZQogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9PRkwiIC8+CiAgICAgICAgPGRjOnN1YmplY3Q+CiAgICAgICAgICA8cmRmOkJhZz4KICAgICAgICAgICAgPHJkZjpsaT5Hb29nbGUgU2Nob2xhcjwvcmRmOmxpPgogICAgICAgICAgPC9yZGY6QmFnPgogICAgICAgIDwvZGM6c3ViamVjdD4KICAgICAgPC9jYzpXb3JrPgogICAgICA8Y2M6TGljZW5zZQogICAgICAgICByZGY6YWJvdXQ9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvT0ZMIj4KICAgICAgICA8Y2M6cGVybWl0cwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML1JlcHJvZHVjdGlvbiIgLz4KICAgICAgICA8Y2M6cGVybWl0cwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML0Rpc3RyaWJ1dGlvbiIgLz4KICAgICAgICA8Y2M6cGVybWl0cwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML0VtYmVkZGluZyIgLz4KICAgICAgICA8Y2M6cGVybWl0cwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML0Rlcml2YXRpdmVXb3JrcyIgLz4KICAgICAgICA8Y2M6cmVxdWlyZXMKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvcHViL09GTC9Ob3RpY2UiIC8+CiAgICAgICAgPGNjOnJlcXVpcmVzCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vc2NyaXB0cy5zaWwub3JnL3B1Yi9PRkwvQXR0cmlidXRpb24iIC8+CiAgICAgICAgPGNjOnJlcXVpcmVzCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vc2NyaXB0cy5zaWwub3JnL3B1Yi9PRkwvU2hhcmVBbGlrZSIgLz4KICAgICAgICA8Y2M6cmVxdWlyZXMKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvcHViL09GTC9EZXJpdmF0aXZlUmVuYW1pbmciIC8+CiAgICAgICAgPGNjOnJlcXVpcmVzCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vc2NyaXB0cy5zaWwub3JnL3B1Yi9PRkwvQnVuZGxpbmdXaGVuU2VsbGluZyIgLz4KICAgICAgPC9jYzpMaWNlbnNlPgogICAgPC9yZGY6UkRGPgogIDwvbWV0YWRhdGE+CiAgPGcKICAgICBpbmtzY2FwZTpncm91cG1vZGU9ImxheWVyIgogICAgIGlkPSJsYXllcjYiCiAgICAgaW5rc2NhcGU6bGFiZWw9Imljb24iPgogICAgPHBhdGgKICAgICAgIHN0eWxlPSJmaWxsOiM4YWIwYmY7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlLXdpZHRoOjAuMDY0Nzg3MSIKICAgICAgIGQ9Ik0gMzQzLjc1ODY4LDEwNi42NjI0MyBWIDc5LjQzMDIwNSBMIDM2My41MjM2NSw2My45OTk5OTcgSCAxNDkuNjMzNTQgTCAyMC40NzYzNDUsMTc2LjI3MzYgaCA4NS42NTYwNzUgYyAtMC4xNTUzNCwyLjEyNDk0IC0wLjIxOTE0LDQuMDQ2NDQgLTAuMjE5MTQsNi4yMjU2MyAwLDIwLjg0NDcyIDcuMjE5MiwzOC4wODY2MiAyMS42NzIwMyw1MS44NjA4OSAxNC40NTI4NCwxMy43OTcwMiAzMi4yNTEyNCwyMC42NDc4NCA1My4zMjY1MSwyMC42NDc4NCA0LjkyMzE5LDAgOS43NTA1OSwtMC4zNjc5NCAxNC40Mzg0MiwtMS4wMjQxOSAtMi45MDcyMiw2LjUwMDgyIC00LjM3NDU3LDEyLjUyMzAyIC00LjM3NDU3LDE4LjE0MjI4IDAsOS44NzUyNiA0LjQ5OTI0LDIwLjQzMDQgMTMuNDY3MTUsMzEuNjQxOCAtMzkuMjMzNzcsMi42NzA1IC02OC4wNjExMiw5LjczMjY0IC04Ni40MzcwMiwyMS4xNjMyMiAtMTAuNTMxMDgsNi40OTkwNyAtMTkuMDAwMjA3LDE0LjcwMzk2IC0yNS4zOTAzNDksMjQuNTMxMSAtNi4zOTA1NjksOS44OTkzMyAtOS41Nzc3NTQsMjAuNTE1MjUgLTkuNTc3NzU0LDMxLjk2MTYgMCw5LjY0ODIyIDIuMDYyMzc1LDE4LjMzNjExIDYuMjE5MDcsMjYuMDYyMzMgNC4xNTY2OTQsNy43MjYzIDkuNTc3NzU3LDE0LjA3MDQ3IDE2LjMxMjIyMywxOC45ODQwOCA2LjcxODI1LDQuOTY3ODEgMTQuNDY4OTksOS4xMDA4OCAyMy4yMTksMTIuNDY4NzQgOC43MzQyOSwzLjM0Mzc4IDE3LjQwNjQzLDUuNzE4NTggMjYuMDYxMDYsNy4wNjI1OCA4LjYyNzA3LDEuMzQyMjIgMTcuMjA0NzEsMS45OTg1IDI1LjcwNTc5LDEuOTk4NSAxMy40Njg4NywwIDI2Ljk1MzUzLC0xLjczNDI4IDQwLjU0NzExLC01LjE4NzA3IDEzLjU2MTY1LC0zLjQ4NDYxIDI2LjI4MDIyLC04LjY0MTQzIDM4LjE3MTA1LC0xNS40OTI3IDExLjg1OTM1LC02LjgwNDg4IDIxLjUxNTQ1LC0xNi4wODY1IDI4LjkyMTksLTI3LjcxODMgNy4zOTAyNCwtMTEuNjc5OTggMTEuMDk0NTcsLTI0LjgwNDk5IDExLjA5NDU3LC0zOS4zMzYxMyAwLC0xMS4wMTU4NCAtMi4yNDk2NCwtMjEuMDM4NTIgLTYuNzUwMiwtMzAuMTQwNzMgLTQuNDY4NjQsLTkuMDcyMDIgLTkuOTM3ODUsLTE2LjU0MTAyIC0xNi40NTI3MSwtMjIuMzQ0MDMgLTYuNTAwOCwtNS44MTI2MyAtMTIuOTk5ODcsLTExLjE1NTM5IC0xOS41MTUxMiwtMTUuOTY3OSAtNi41MDA4MywtNC44NDQ4OCAtMTIuMDAwMjEsLTkuNzUwNTggLTE2LjQ2ODg0LC0xNC44MTI5IC00LjQ4NDgsLTUuMDQ2NTcgLTYuNzM0NDQsLTEwLjA1NDE5IC02LjczNDQ0LC0xNC45ODM5NSAwLC00LjkyMTQ1IDEuNzM0MjIsLTkuNjcxODMgNS4yMTU4OCwtMTQuMjY1NTkgMy40NTQ1MSwtNC42MDk1IDcuNjczNzYsLTkuMDQ3OTUgMTIuNjA5NjcsLTEzLjMwNTcxIDQuOTM3NTYsLTQuMjQ5NDQgOS44NzUyMywtOC45Njc4OCAxNC43OTY2NSwtMTQuMTMzMDIgNC45MjE0NywtNS4xNDcxOSA5LjE0MDcyLC0xMS44MjczOSAxMi42MDk3MSwtMjAuMDA4MjIgMy40ODQ2NywtOC4xNzkwNyA1LjIwMzE4LC0xNy40NDQ4OSA1LjIwMzE4LC0yNy43NTY3OSAwLC0xMy40NTI3IC0yLjU0NzE0LC0yNC40NjA2NSAtNy41NDczNSwtMzMuMzEzNDggLTAuNTkzNjksLTEuMDIyNDMgLTEuMjE3NTcsLTEuODAzMzggLTEuODc1MTEsLTMuMDIyMjUgbCA1Ni45MDc0NSwtNDYuNjcyMTM2IHYgMTcuMTE4NTI2IGMgLTcuMzkzNzMsMC45Mjk2OSAtNi42MjQyMiw1LjM0NTgyIC02LjYyNDIyLDEwLjYzNTIgdiAxMjguNjY3MTkgYyAwLDUuOTU4MzIgNC44NzUxLDEwLjgzMzgyIDEwLjgzMzg2LDEwLjgzMzgyIGggMy45ODg2OSBjIDUuOTU4MzUsMCAxMC44MzM4NiwtNC44NzUwNiAxMC44MzM4NiwtMTAuODMzODIgViAxMTcuMjkyODIgYyAwLC01LjI3NjY5IDAuNzc3NDEsLTkuNjg4MDEgLTYuNTYxNjcsLTEwLjYzMDM5IHogTSAyMzYuMzk4NjUsMzI5LjE0MTE0IGMgMS4xNDA5OSwwLjc1MDMgMy43MDM5LDIuNzgwNzUgNy43MTg0LDYuMDM4MzggNC4wNDk1LDMuMjQzMTkgNi43OTcsNS42OTU4MiA4LjI2NTY3LDcuNDE0MzIgMS40Mzg1MSwxLjY2MzgxIDMuNTc5Miw0LjE2NTAxIDYuMzc2MTcsNy41NDczNCAyLjgxMjY4LDMuMzc0NCA0LjcxODQsNi4zMDM5NCA1LjcxODUzLDguNzM0MjUgMS4wMDAxNiwyLjQ3NjcgMi4wMTYwMyw1LjQ2MDg5IDMuMDQ2MzYsOC45NDU1NiAwLjk4NTY3LDMuNDQ0ODggMS40ODQ4Niw2Ljk3NTk1IDEuNDg0ODYsMTAuNTYxNjkgMCwxNy4wNDgxMyAtNi41NjMzOCwyOS42ODAwNyAtMTkuNjU2MDQsMzcuODU5MTUgLTEzLjEyNSw4LjE4MDgzIC0yOC43NjY1MSwxMi4yNzM2OCAtNDYuOTM3NjcsMTIuMjczNjggLTkuMTg3MDksMCAtMTguMjAzMSwtMS4wOTI4OSAtMjcuMDYyNDcsLTMuMTk1MSAtOC44NDMyMiwtMi4xMTY2NSAtMTcuMzExOTIsLTUuMzM2MiAtMjUuMzkwMzUsLTkuNjAxODUgLTguMDc4NDYsLTQuMjU3NzEgLTE0LjU3NzU0LC0xMC4yMDMzNyAtMTkuNTAwNzIsLTE3Ljc5NjU5IC00LjkzNzY0LC03LjY0MDEyIC03LjQwNjQ1LC0xNi40MTQ2NCAtNy40MDY0NSwtMjYuMjQ5NjIgMCwtMTAuMzIwMjIgMi43OTY5MiwtMTkuMjg5ODcgOC40MjIzMywtMjYuOTA1ODggNS41OTM0MywtNy42MjU2NCAxMi45Mzc3NCwtMTMuMzkxOSAyMi4wMzIwOCwtMTcuMzE1NCA5LjA2MjQsLTMuOTQ1ODIgMTguMjQ5NDYsLTYuNzQyMzIgMjcuNTYxNjYsLTguMzk4MjcgOS4zMTIyMSwtMS43MDIzIDE4Ljc5Njc5LC0yLjU1NSAyOC40Mzg0MiwtMi41NTUgNC40Njg2MiwwIDcuOTM1ODIsMC4yNTExNSAxMC40MDQ2NSwwLjY5NjA3IDAuNDU0NTYsMC4yMTkxOCAzLjAzMTg4LDIuMDcwMjUgNy43MzQ1Niw1LjU2MzI2IDQuNzA0MDEsMy40NjIzNyA3LjYyNTY1LDUuNTk1MTkgOC43NTA0Nyw2LjM4NDAxIHogbSAtMy4zNTgyMywtMTAwLjU3NzkgYyAtNy40MDY0OCw4Ljg1OTM4IC0xNy43MzQ1NCwxMy4yODgyIC0zMC45NTM2MywxMy4yODgyIC0xMS44NTkzMywwIC0yMi4yOTc2NiwtNC43NjQ4MiAtMzEuMjY1NTQsLTE0LjMxMTk1IC04Ljk5OTg0LC05LjUyMzA5IC0xNS40MjIzNSwtMjAuMzI4MDMgLTE5LjM0NDA4LC0zMi40MzA2MSAtMy45Mzc1MiwtMTIuMTA4NzEgLTUuOTA1ODUsLTIzLjk4NDIzIC01LjkwNTg1LC0zNS42NDggMCwtMTMuNjk0MiAzLjU5NTQyLC0yNS4zNTE4NCAxMC43ODA5LC0zNC45NzU5OCA3LjE4NzI3LC05LjY0OTUyIDE3LjQ5OTE1LC0xNC40ODQ3NyAzMC45Mzc4NiwtMTQuNDg0NzcgMTEuODc1MDcsMCAyMi4zNzQyMyw1LjAzODI1IDMxLjQzNzA0LDE1LjE1Njc3IDkuMDk0MzQsMTAuMDg0ODIgMTUuNjA5NjEsMjEuNDEzMDMgMTkuNTE2OSwzMy45Njc5OSAzLjkyMTc2LDEyLjUzOTIgNS44NzM0NSwyNC41Mjk3OSA1Ljg3MzQ1LDM1Ljk4Mzk5IDAsMTMuNDQ2NTggLTMuNzAyNTYsMjQuNjA5ODQgLTExLjA3NjYzLDMzLjQ1NDM2IHoiCiAgICAgICBpZD0icGF0aDIiIC8+CiAgPC9nPgo8L3N2Zz4K" /></a> <a href="https://www.researchgate.net/profile/Javier-Tamayo"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgd2lkdGg9IjM4NCIKICAgaGVpZ2h0PSI1MTIiCiAgIHZpZXdCb3g9IjAgMCAzODQgNTEyIgogICB2ZXJzaW9uPSIxLjEiCiAgIGlkPSJzdmc4IgogICBpbmtzY2FwZTp2ZXJzaW9uPSIxLjAgKDQwMzVhNGYsIDIwMjAtMDUtMDEpIgogICBzb2RpcG9kaTpkb2NuYW1lPSJyZXNlYXJjaGdhdGUuc3ZnIj4KICA8dGl0bGUKICAgICBpZD0idGl0bGU4MzMiPlJlc2VhcmNoR2F0ZTwvdGl0bGU+CiAgPGRlZnMKICAgICBpZD0iZGVmczIiIC8+CiAgPHNvZGlwb2RpOm5hbWVkdmlldwogICAgIGlkPSJiYXNlIgogICAgIHBhZ2Vjb2xvcj0iI2ZmZmZmZiIKICAgICBib3JkZXJjb2xvcj0iIzY2NjY2NiIKICAgICBib3JkZXJvcGFjaXR5PSIxLjAiCiAgICAgaW5rc2NhcGU6cGFnZW9wYWNpdHk9IjAuMCIKICAgICBpbmtzY2FwZTpwYWdlc2hhZG93PSIyIgogICAgIGlua3NjYXBlOnpvb209IjAuNzUiCiAgICAgaW5rc2NhcGU6Y3g9IjE5MiIKICAgICBpbmtzY2FwZTpjeT0iMjU2IgogICAgIGlua3NjYXBlOmRvY3VtZW50LXVuaXRzPSJweCIKICAgICBpbmtzY2FwZTpjdXJyZW50LWxheWVyPSJsYXllcjYiCiAgICAgaW5rc2NhcGU6ZG9jdW1lbnQtcm90YXRpb249IjAiCiAgICAgc2hvd2dyaWQ9InRydWUiCiAgICAgdW5pdHM9InB4IgogICAgIHNob3dndWlkZXM9InRydWUiCiAgICAgaW5rc2NhcGU6Z3VpZGUtYmJveD0idHJ1ZSIKICAgICBpbmtzY2FwZTp3aW5kb3ctd2lkdGg9IjE2ODAiCiAgICAgaW5rc2NhcGU6d2luZG93LWhlaWdodD0iOTU4IgogICAgIGlua3NjYXBlOndpbmRvdy14PSIwIgogICAgIGlua3NjYXBlOndpbmRvdy15PSIyMyIKICAgICBpbmtzY2FwZTp3aW5kb3ctbWF4aW1pemVkPSIxIj4KICAgIDxpbmtzY2FwZTpncmlkCiAgICAgICB0eXBlPSJ4eWdyaWQiCiAgICAgICBpZD0iZ3JpZDEwIgogICAgICAgZW1wc3BhY2luZz0iNCIgLz4KICAgIDxzb2RpcG9kaTpndWlkZQogICAgICAgcG9zaXRpb249IjM4NC4wMDAwMSw2NCIKICAgICAgIG9yaWVudGF0aW9uPSIwLDEiCiAgICAgICBpZD0iZ3VpZGUyMyIKICAgICAgIGlua3NjYXBlOmxhYmVsPSJiYXNlbGluZSIKICAgICAgIGlua3NjYXBlOmxvY2tlZD0idHJ1ZSIKICAgICAgIGlua3NjYXBlOmNvbG9yPSJyZ2IoMCwwLDI1NSkiIC8+CiAgPC9zb2RpcG9kaTpuYW1lZHZpZXc+CiAgPG1ldGFkYXRhCiAgICAgaWQ9Im1ldGFkYXRhNSI+CiAgICA8cmRmOlJERj4KICAgICAgPGNjOldvcmsKICAgICAgICAgcmRmOmFib3V0PSIiPgogICAgICAgIDxkYzpmb3JtYXQ+aW1hZ2Uvc3ZnK3htbDwvZGM6Zm9ybWF0PgogICAgICAgIDxkYzp0eXBlCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vcHVybC5vcmcvZGMvZGNtaXR5cGUvU3RpbGxJbWFnZSIgLz4KICAgICAgICA8ZGM6dGl0bGU+UmVzZWFyY2hHYXRlPC9kYzp0aXRsZT4KICAgICAgICA8Y2M6bGljZW5zZQogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9PRkwiIC8+CiAgICAgICAgPGRjOnN1YmplY3Q+CiAgICAgICAgICA8cmRmOkJhZz4KICAgICAgICAgICAgPHJkZjpsaT5SZXNlYXJjaEdhdGU8L3JkZjpsaT4KICAgICAgICAgIDwvcmRmOkJhZz4KICAgICAgICA8L2RjOnN1YmplY3Q+CiAgICAgIDwvY2M6V29yaz4KICAgICAgPGNjOkxpY2Vuc2UKICAgICAgICAgcmRmOmFib3V0PSJodHRwOi8vc2NyaXB0cy5zaWwub3JnL09GTCI+CiAgICAgICAgPGNjOnBlcm1pdHMKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvcHViL09GTC9SZXByb2R1Y3Rpb24iIC8+CiAgICAgICAgPGNjOnBlcm1pdHMKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvcHViL09GTC9EaXN0cmlidXRpb24iIC8+CiAgICAgICAgPGNjOnBlcm1pdHMKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvcHViL09GTC9FbWJlZGRpbmciIC8+CiAgICAgICAgPGNjOnBlcm1pdHMKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9zY3JpcHRzLnNpbC5vcmcvcHViL09GTC9EZXJpdmF0aXZlV29ya3MiIC8+CiAgICAgICAgPGNjOnJlcXVpcmVzCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vc2NyaXB0cy5zaWwub3JnL3B1Yi9PRkwvTm90aWNlIiAvPgogICAgICAgIDxjYzpyZXF1aXJlcwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML0F0dHJpYnV0aW9uIiAvPgogICAgICAgIDxjYzpyZXF1aXJlcwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML1NoYXJlQWxpa2UiIC8+CiAgICAgICAgPGNjOnJlcXVpcmVzCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vc2NyaXB0cy5zaWwub3JnL3B1Yi9PRkwvRGVyaXZhdGl2ZVJlbmFtaW5nIiAvPgogICAgICAgIDxjYzpyZXF1aXJlcwogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3NjcmlwdHMuc2lsLm9yZy9wdWIvT0ZML0J1bmRsaW5nV2hlblNlbGxpbmciIC8+CiAgICAgIDwvY2M6TGljZW5zZT4KICAgIDwvcmRmOlJERj4KICA8L21ldGFkYXRhPgogIDxnCiAgICAgaW5rc2NhcGU6Z3JvdXBtb2RlPSJsYXllciIKICAgICBpZD0ibGF5ZXI2IgogICAgIGlua3NjYXBlOmxhYmVsPSJpY29uIj4KICAgIDxwYXRoCiAgICAgICBzdHlsZT0iZmlsbDojOGFiMGJmO2ZpbGwtb3BhY2l0eToxO3N0cm9rZS13aWR0aDowLjA2NjMxNDgiCiAgICAgICBkPSJNIDIyOC42NTk2OCw0MDguMTAxNDkgQyAyMDguNDY2MTYsMzg2LjEyNzg0IDE4MS40NDIyLDM1MC44MzMyNCAxNTguOTYyLDMxMC45MzUxNiBjIDM3LjEzNDQ0LC04LjY4NDg3IDY0LjYzMzE4LC00My41NTAyNyA2NC42MzMxOCwtNzguMzgzODYgMCwtNTEuMzQ1NyAtMzkuODgwODksLTc1LjE3Njc4IC05Mi4xNjMyOSwtNzUuMTc2NzggLTI3LjAyMjY2LDAgLTQ4LjU4MzM2NywxLjM2NTQ2IC02OC43MTU0NDEsMS4zNjU0NiAtMTguMzY5MTU1LDAgLTM2LjcyMjQwMSwwIC00OC4xNTQxNTMsLTAuNDQ0NyB2IDEzLjMwNDY2IGwgMTcuNDAxMDk3LDMuMTkxNTkgYyAxMS45Njg4NDUsMi4zMDIxNCAxOC43ODI0OSw3Ljc2NTc5IDE4Ljc4MjQ5LDM2LjIxNDk0IHYgMTgwLjU5OTE5IGMgMCwyOC40MzU0MyAtNi44MTMyMjIsMzMuOTI4NzEgLTE4Ljc4MjQ5LDM2LjE5OTUzIGwgLTE3LjQwMTA5NywzLjI1MjU4IHYgMTMuMjU5MDggYyAxMi4zNjc5OTgsLTAuNDQ0NyAzMy45MTI4MjgsLTEuMzUxMyA1NS40NzM1MzMsLTEuMzUxMyAyMC42MjQwNCwwIDQ3LjIxNzQ1MSwwLjkwNjYgNTguNjgwNTYxLDEuMzUxMyB2IC0xMy4yNTkwOCBsIC0yMy44NDY5OCwtMy4yNTI1OCBjIC0xMi4zNjY2OTYsLTEuNzk2IC0xOC44MTM4NDgsLTcuNzY0MSAtMTguODEzODQ4LC0zNi4xOTk1MyB2IC03Ni41NDIyNSBjIDExLjAwMjUwNywwLjkyMDc4IDIwLjYyNDAzOCwwLjkyMDc4IDM1LjMyNTU0OCwwLjkyMDc4IDI3Ljk1OTMyLDQ5Ljk0ODg3IDU0LjU1MTAyLDg3LjU1ODk2IDY5LjY1MjE0LDEwNC45NjE4IDEzLjc3OSwxNi41MjU4NyAzNC44NDk0OSwyNy4wNTQwMSA2MS40NDE2MiwyNy4wNTQwMSA3Ljc4MTI4LDAgMTYuMDIzMTYsLTEuMzY3MjEgMjEuMDU0NTgsLTMuNjgzMTUgdiAtMTEuODkzNjEgYyAtMTYuNDk2MjUsMCAtMzIuOTkyMDQsLTExLjQ3NzI1IC00NC44NzAyLC0yNC4zMjEzIHogTSAxMTkuMDYzOTIsMjk1LjM0NDI0IGMgLTE1LjU5MDkzLDAgLTIyLjQzNDI1NywtMC40MTQ2MyAtMzMuMDA3OTM5LC0xLjQxMTAzIFYgMTc4LjkxOTM1IGMgMTAuNTczNjgyLC0wLjkyMDc5IDI0Ljc1MTg2OSwtMC45MjA3OSAzNy4xMzU3MzksLTAuOTIwNzkgMzguNTMxMzIsMCA2MS40Mjc0NSwyMC4xOTUyNSA2MS40Mjc0NSw1Ni44MzkwMiAwLDM2LjIxNDk1IC0yNC43MzY0MSw2MC41MDY2NiAtNjUuNTU1MjUsNjAuNTA2NjYgeiBtIDEzOS45MzQxNiwtMTE1LjcwMzUgYyAtMC40NjA2LC0yLjQwODcgLTAuODc1MjQsLTUuMjE3MDEgLTEuMjQzMDMsLTguNDU1NDQgLTAuMzgzMjgsLTMuMjY4NDYgLTAuNjMwMzMsLTcuMTA0MTEgLTAuNzgyNDMsLTExLjYzMTA4IC0wLjE1Mzg2LC00LjQ5NTYyIC0wLjIxNTA5LC05Ljk4ODkxIC0wLjIxNTA5LC0xNi4yODA5OCAwLC02LjMyMjk3IDAuMDYxNiwtMTEuNzM5OCAwLjIxNTA5LC0xNi4yNTA4NyAwLjE1MjU4LC00LjUyNjk4IDAuMzk5MTUsLTguNDA5OSAwLjc4MjQzLC0xMS42NDcwMyAwLjM2Nzc4LC0zLjIyMjkxIDAuNzgyNDMsLTYuMDQ2NyAxLjI0MzAzLC04LjQ1NTQxIDAuNDQ0NywtMi40MjQ2IDEuMDI3NzYsLTQuNzI2NzkgMS43ODA1NSwtNi45ODE2NzggMy45NDI2MiwtMTEuOTIzMjY0IDEwLjQwMzk2LC0yMC44ODQ4NjEgMTkuNDQyNDYsLTI2LjkwMTUwMyA5LjAyMjYsLTYuMDE1MzQyIDIwLjA4NjUzLC05LjAzNjc1MiAzMy4yMjE0NywtOS4wMzY3NTIgNi43NTIyNCwwIDEyLjg3NTg3LDAuNzgyNDI5IDE4LjMwNzczLDIuMzc3MzQ1IDUuMzcwODIsMS41ODA3NjEgMTAuMTg4NywzLjc2MDQ1NyAxNC40Mzg5Nyw2LjU1MjQzNSA0LjE4OTI2LDIuNzYyMzM5IDcuODExMzcsNi4wMTUzNDkgMTAuODAzMTcsOS42OTg0NzYgMy4wNjg2OCwzLjY5NzI4NSA1LjUzODgyLDcuNjQxNjI0IDcuNDczMTksMTEuODc3NzE4IDAuNzUxMDYsMS4zNDk1ODMgMC41Mzc1MiwyLjQyNDYwOCAtMC42NTk5OCwzLjE3NTY3NiBsIC0xNi42ODAxMiw2Ljg1ODc3MyBjIC0xLjQxMjc1LDAuNzUyNzggLTIuNDA4NjksMC4zNjkwNyAtMy4xOTI4NywtMS4xMjA1OCAtMy43NDI4MywtNi45MzYxMyAtNi45MzYxMSwtMTEuNDkyNzUgLTEyLjE4MjgxLC0xNC44MDY3NTcgLTUuMzA5NDEsLTMuMjk5ODUzIC0xMC4yMjAwNywtNC44NjUxMjkgLTE4LjMwNzcsLTQuODY1MTI5IC04Ljc5MzE0LDAgLTEyLjcyMTYsMS43NDg3MzQgLTE4LjIyOTUzLDUuNjkzMTAyIC01LjQ2MzY0LDMuODY3MDA4IC05LjQ2OTAxLDguNzMyMTM3IC0xMS44MzIxNywxNi4xNDI2MTQgLTAuNDc0MzcsMS4zMzU4MSAtMC45MDUzLDIuOTkzNDggLTEuNDExMDQsNC45NDIwMyAtMC40MTQ2MywxLjk3OTkxIC0wLjc1Mjc4LDQuNDA0NTIgLTAuOTY2Myw3LjI0MjQ3IC0wLjIxNTQ0LDIuODM5MjUgLTAuNDE0NjQsNi4zNTMwOCAtMC41OTg1NSwxMC40OTY3OCAtMC4xMjI0Niw0LjE0MzczIC0wLjE4NDAzLDkuMTc2ODYgLTAuMTg0MDMsMTUuMDA3ODggMCw1Ljg2MTA3IDAuMDYxNiwxMC44OTQyIDAuMTg0MDMsMTUuMDM3OTQgMC4xODQwMiw0LjEyNzc5IDAuMzgzMjgsNy42NDE2MiAwLjU5ODU1LDEwLjQ4MDg2IDAuMjEzNSwyLjg2OTMxIDAuNTUxNjcsNS4yNzgwNSAwLjk2NjMsNy4yNzM4NiAwLjUwNjE2LDEuOTE4NDYgMC45MzY2NywzLjU3NDgzIDEuNDExMDQsNC45MjYxMyAyLjM2MzE2LDcuMzgwODQgNS43NDAzNywxMS40MTU4NSAxMC43MTIwNSwxNC42NTQ2NyA0LjkxMDY4LDMuMjgzOTQgMTAuNTU2NTEsNS42NDc1NSAxOS4zNDk2NSw1LjY0NzU1IDcuODExMzQsMCAxNC45NjIzMiwtMi4yMjUyNCAxOS42MjYzNCwtNS42MTc0NiA0LjYxOTc5LC0zLjM5MDk0IDguNDU1NDQsLTcuODcxMDggMTAuMTc0NTIsLTEzLjk5NDcyIDAuNzUyODEsLTIuNTc4NDMgMS43MTkxNiwtNS43ODU4OCAyLjM3OTEsLTkuNzEzOTMgMC41OTg1MywtMy45Mjg0MyAwLjU5ODUzLC04LjA4NjMyIDAuNTk4NTMsLTEzLjgyNDk3IDAsLTAuOTA2NiAtMC41MDc0LC0xLjM2NzIxIC0xLjM1MTMsLTEuMzY3MjEgaCAtMjYuNzE2MzEgYyAtMS41MDM4MywwIC0yLjI0MDczLC0wLjczNTU4IC0yLjI0MDczLC0yLjIzODk4IHYgLTE1LjMxNDY2IGMgMCwtMS41MTkzMSAwLjczNjksLTIuMjU2NjEgMi4yNDA3MywtMi4yNTY2MSBoIDQ5LjAyODA3IGMgMS41MzUyMSwwIDIuMjU2NjEsMC43MzY4OCAyLjI1NjYxLDIuMjU2NjEgdiAxMy4wODkzOCBjIDAsNi45MzYxNCAwLDEzLjM2NjEyIC0wLjcyMTQsMTkuMzIwMDEgLTAuNjkxMzEsNS45NTM5MiAtMS42MjYyNSwxMS4xMDk1MSAtMi44MDc4OSwxNC44NjgyMiAtMy43NDQ1MywxMS43NzEyIC05LjY4MjU0LDIwLjE0OTY5IC0xOC43ODI0NywyNi4zOTQwNSAtOS4xMzEzLDYuMjkxNjUgLTIwLjkwMDc5LDkuNjgyNTcgLTMzLjY4MzgsOS42ODI1NyAtMTMuMTM0OTQsMCAtMjQuMTk4OTEsLTMuMDIxODMgLTMzLjIyMTQ3LC05LjAyMjYgLTkuMDM4NSwtNi4wNzYzNSAtMTUuNDk5ODQsLTE0Ljk5MjM5IC0xOS40NDI0NiwtMjYuOTE1NjYgLTAuNzUyNzksLTIuMjU0ODkgLTEuMzM1ODUsLTQuNTcyNTMgLTEuNzgwNTUsLTYuOTk3MTQgeiIKICAgICAgIGlkPSJwYXRoMiIgLz4KICA8L2c+Cjwvc3ZnPgo=" /></a> <a href="https://twitter.com/TamayoLeiva_J"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgaW5rc2NhcGU6dmVyc2lvbj0iMS4wICg0MDM1YTRmLCAyMDIwLTA1LTAxKSIKICAgc29kaXBvZGk6ZG9jbmFtZT0idHdpdHRlci5zdmciCiAgIGlkPSJzdmc4NTMiCiAgIHZlcnNpb249IjEuMSIKICAgdmlld0JveD0iMCAwIDI0IDI0IgogICByb2xlPSJpbWciPgogIDxtZXRhZGF0YQogICAgIGlkPSJtZXRhZGF0YTg1OSI+CiAgICA8cmRmOlJERj4KICAgICAgPGNjOldvcmsKICAgICAgICAgcmRmOmFib3V0PSIiPgogICAgICAgIDxkYzpmb3JtYXQ+aW1hZ2Uvc3ZnK3htbDwvZGM6Zm9ybWF0PgogICAgICAgIDxkYzp0eXBlCiAgICAgICAgICAgcmRmOnJlc291cmNlPSJodHRwOi8vcHVybC5vcmcvZGMvZGNtaXR5cGUvU3RpbGxJbWFnZSIgLz4KICAgICAgPC9jYzpXb3JrPgogICAgPC9yZGY6UkRGPgogIDwvbWV0YWRhdGE+CiAgPGRlZnMKICAgICBpZD0iZGVmczg1NyIgLz4KICA8c29kaXBvZGk6bmFtZWR2aWV3CiAgICAgaW5rc2NhcGU6Y3VycmVudC1sYXllcj0ic3ZnODUzIgogICAgIGlua3NjYXBlOndpbmRvdy1tYXhpbWl6ZWQ9IjAiCiAgICAgaW5rc2NhcGU6d2luZG93LXk9IjIzIgogICAgIGlua3NjYXBlOndpbmRvdy14PSIwIgogICAgIGlua3NjYXBlOmN5PSIxMiIKICAgICBpbmtzY2FwZTpjeD0iMTIiCiAgICAgaW5rc2NhcGU6em9vbT0iMzAuNjI1IgogICAgIHNob3dncmlkPSJmYWxzZSIKICAgICBpZD0ibmFtZWR2aWV3ODU1IgogICAgIGlua3NjYXBlOndpbmRvdy1oZWlnaHQ9Ijk1NyIKICAgICBpbmtzY2FwZTp3aW5kb3ctd2lkdGg9IjEzNzUiCiAgICAgaW5rc2NhcGU6cGFnZXNoYWRvdz0iMiIKICAgICBpbmtzY2FwZTpwYWdlb3BhY2l0eT0iMCIKICAgICBndWlkZXRvbGVyYW5jZT0iMTAiCiAgICAgZ3JpZHRvbGVyYW5jZT0iMTAiCiAgICAgb2JqZWN0dG9sZXJhbmNlPSIxMCIKICAgICBib3JkZXJvcGFjaXR5PSIxIgogICAgIGJvcmRlcmNvbG9yPSIjNjY2NjY2IgogICAgIHBhZ2Vjb2xvcj0iI2ZmZmZmZiIKICAgICBpbmtzY2FwZTpkb2N1bWVudC1yb3RhdGlvbj0iMCIgLz4KICA8dGl0bGUKICAgICBpZD0idGl0bGU4NDkiPlR3aXR0ZXIgaWNvbjwvdGl0bGU+CiAgPHBhdGgKICAgICBzdHlsZT0iZmlsbDojOGFiMGJmO2ZpbGwtb3BhY2l0eToxIgogICAgIGlkPSJwYXRoODUxIgogICAgIGQ9Ik0yMy45NTMgNC41N2ExMCAxMCAwIDAxLTIuODI1Ljc3NSA0Ljk1OCA0Ljk1OCAwIDAwMi4xNjMtMi43MjNjLS45NTEuNTU1LTIuMDA1Ljk1OS0zLjEyNyAxLjE4NGE0LjkyIDQuOTIgMCAwMC04LjM4NCA0LjQ4MkM3LjY5IDguMDk1IDQuMDY3IDYuMTMgMS42NCAzLjE2MmE0LjgyMiA0LjgyMiAwIDAwLS42NjYgMi40NzVjMCAxLjcxLjg3IDMuMjEzIDIuMTg4IDQuMDk2YTQuOTA0IDQuOTA0IDAgMDEtMi4yMjgtLjYxNnYuMDZhNC45MjMgNC45MjMgMCAwMDMuOTQ2IDQuODI3IDQuOTk2IDQuOTk2IDAgMDEtMi4yMTIuMDg1IDQuOTM2IDQuOTM2IDAgMDA0LjYwNCAzLjQxNyA5Ljg2NyA5Ljg2NyAwIDAxLTYuMTAyIDIuMTA1Yy0uMzkgMC0uNzc5LS4wMjMtMS4xNy0uMDY3YTEzLjk5NSAxMy45OTUgMCAwMDcuNTU3IDIuMjA5YzkuMDUzIDAgMTMuOTk4LTcuNDk2IDEzLjk5OC0xMy45ODUgMC0uMjEgMC0uNDItLjAxNS0uNjNBOS45MzUgOS45MzUgMCAwMDI0IDQuNTl6IiAvPgo8L3N2Zz4K" /></a> <a href="https://github.com/TamayoLeivaJ"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgcm9sZT0iaW1nIgogICB2aWV3Qm94PSIwIDAgMjQgMjQiCiAgIHZlcnNpb249IjEuMSIKICAgaWQ9InN2ZzgzNyIKICAgc29kaXBvZGk6ZG9jbmFtZT0iZ2l0aHViLnN2ZyIKICAgaW5rc2NhcGU6dmVyc2lvbj0iMS4wICg0MDM1YTRmLCAyMDIwLTA1LTAxKSI+CiAgPG1ldGFkYXRhCiAgICAgaWQ9Im1ldGFkYXRhODQzIj4KICAgIDxyZGY6UkRGPgogICAgICA8Y2M6V29yawogICAgICAgICByZGY6YWJvdXQ9IiI+CiAgICAgICAgPGRjOmZvcm1hdD5pbWFnZS9zdmcreG1sPC9kYzpmb3JtYXQ+CiAgICAgICAgPGRjOnR5cGUKICAgICAgICAgICByZGY6cmVzb3VyY2U9Imh0dHA6Ly9wdXJsLm9yZy9kYy9kY21pdHlwZS9TdGlsbEltYWdlIiAvPgogICAgICA8L2NjOldvcms+CiAgICA8L3JkZjpSREY+CiAgPC9tZXRhZGF0YT4KICA8ZGVmcwogICAgIGlkPSJkZWZzODQxIj4KICAgIDxmaWx0ZXIKICAgICAgIGhlaWdodD0iMy40Mjg5OTUxIgogICAgICAgeT0iLTEuMjE0NDk3NiIKICAgICAgIHdpZHRoPSIzLjM3MTY4OSIKICAgICAgIHg9Ii0xLjE4NTg0NDUiCiAgICAgICBpZD0iZmlsdGVyODQyIgogICAgICAgc3R5bGU9ImNvbG9yLWludGVycG9sYXRpb24tZmlsdGVyczpzUkdCIgogICAgICAgaW5rc2NhcGU6Y29sbGVjdD0iYWx3YXlzIj4KICAgICAgPGZlR2F1c3NpYW5CbHVyCiAgICAgICAgIGlkPSJmZUdhdXNzaWFuQmx1cjg0NCIKICAgICAgICAgc3RkRGV2aWF0aW9uPSIxMS43OTAyNzciCiAgICAgICAgIGlua3NjYXBlOmNvbGxlY3Q9ImFsd2F5cyIgLz4KICAgIDwvZmlsdGVyPgogIDwvZGVmcz4KICA8c29kaXBvZGk6bmFtZWR2aWV3CiAgICAgaW5rc2NhcGU6ZG9jdW1lbnQtcm90YXRpb249IjAiCiAgICAgcGFnZWNvbG9yPSIjZmZmZmZmIgogICAgIGJvcmRlcmNvbG9yPSIjNjY2NjY2IgogICAgIGJvcmRlcm9wYWNpdHk9IjEiCiAgICAgb2JqZWN0dG9sZXJhbmNlPSIxMCIKICAgICBncmlkdG9sZXJhbmNlPSIxMCIKICAgICBndWlkZXRvbGVyYW5jZT0iMTAiCiAgICAgaW5rc2NhcGU6cGFnZW9wYWNpdHk9IjAiCiAgICAgaW5rc2NhcGU6cGFnZXNoYWRvdz0iMiIKICAgICBpbmtzY2FwZTp3aW5kb3ctd2lkdGg9IjE2NzUiCiAgICAgaW5rc2NhcGU6d2luZG93LWhlaWdodD0iODc4IgogICAgIGlkPSJuYW1lZHZpZXc4MzkiCiAgICAgc2hvd2dyaWQ9ImZhbHNlIgogICAgIGlua3NjYXBlOnpvb209IjI5LjY2NjY2NyIKICAgICBpbmtzY2FwZTpjeD0iMi41ODMxNTE3IgogICAgIGlua3NjYXBlOmN5PSIxMiIKICAgICBpbmtzY2FwZTp3aW5kb3cteD0iNSIKICAgICBpbmtzY2FwZTp3aW5kb3cteT0iMjMiCiAgICAgaW5rc2NhcGU6d2luZG93LW1heGltaXplZD0iMCIKICAgICBpbmtzY2FwZTpjdXJyZW50LWxheWVyPSJzdmc4MzciIC8+CiAgPHRpdGxlCiAgICAgaWQ9InRpdGxlODMzIj5HaXRIdWIgaWNvbjwvdGl0bGU+CiAgPHBhdGgKICAgICBkPSJtIDEyLDAuMjk3IGMgLTYuNjMsMCAtMTIsNS4zNzMgLTEyLDEyIDAsNS4zMDMgMy40MzgsOS44IDguMjA1LDExLjM4NSAwLjYsMC4xMTMgMC44MiwtMC4yNTggMC44MiwtMC41NzcgMCwtMC4yODUgLTAuMDEsLTEuMDQgLTAuMDE1LC0yLjA0IEMgNS42NzIsMjEuNzg5IDQuOTY4LDE5LjQ1NSA0Ljk2OCwxOS40NTUgNC40MjIsMTguMDcgMy42MzMsMTcuNyAzLjYzMywxNy43IDIuNTQ2LDE2Ljk1NiAzLjcxNywxNi45NzEgMy43MTcsMTYuOTcxIGMgMS4yMDUsMC4wODQgMS44MzgsMS4yMzYgMS44MzgsMS4yMzYgMS4wNywxLjgzNSAyLjgwOSwxLjMwNSAzLjQ5NSwwLjk5OCBDIDkuMTU4LDE4LjQyOSA5LjQ2NywxNy45IDkuODEsMTcuNiA3LjE0NSwxNy4zIDQuMzQ0LDE2LjI2OCA0LjM0NCwxMS42NyA0LjM0NCwxMC4zNiA0LjgwOSw5LjI5IDUuNTc5LDguNDUgNS40NDQsOC4xNDcgNS4wMzksNi45MjcgNS42ODQsNS4yNzQgYyAwLDAgMS4wMDUsLTAuMzIyIDMuMywxLjIzIDAuOTYsLTAuMjY3IDEuOTgsLTAuMzk5IDMsLTAuNDA1IDEuMDIsMC4wMDYgMi4wNCwwLjEzOCAzLDAuNDA1IDIuMjgsLTEuNTUyIDMuMjg1LC0xLjIzIDMuMjg1LC0xLjIzIDAuNjQ1LDEuNjUzIDAuMjQsMi44NzMgMC4xMiwzLjE3NiAwLjc2NSwwLjg0IDEuMjMsMS45MSAxLjIzLDMuMjIgMCw0LjYxIC0yLjgwNSw1LjYyNSAtNS40NzUsNS45MiAwLjQyLDAuMzYgMC44MSwxLjA5NiAwLjgxLDIuMjIgMCwxLjYwNiAtMC4wMTUsMi44OTYgLTAuMDE1LDMuMjg2IDAsMC4zMTUgMC4yMSwwLjY5IDAuODI1LDAuNTcgQyAyMC41NjUsMjIuMDkyIDI0LDE3LjU5MiAyNCwxMi4yOTcgYyAwLC02LjYyNyAtNS4zNzMsLTEyIC0xMiwtMTIiCiAgICAgaWQ9InBhdGg4MzUiCiAgICAgc3R5bGU9ImZpbGw6IzAwMDAwMDtmaWxsLW9wYWNpdHk6MTtzdHJva2Utd2lkdGg6MS4wMDE1NztzdHJva2UtbWl0ZXJsaW1pdDo0O3N0cm9rZS1kYXNoYXJyYXk6bm9uZTtmaWxsLXJ1bGU6bm9uemVybztzdHJva2U6bm9uZTtzdHJva2Utb3BhY2l0eToxIiAvPgogIDxwYXRoCiAgICAgc3R5bGU9ImZpbGw6IzhhYjBiZjtmaWxsLW9wYWNpdHk6MTtmaWxsLXJ1bGU6bm9uemVybztzdHJva2U6bm9uZTtzdHJva2Utd2lkdGg6MC4wMzM3NjA5O3N0cm9rZS1taXRlcmxpbWl0OjQ7c3Ryb2tlLWRhc2hhcnJheTpub25lO3N0cm9rZS1vcGFjaXR5OjEiCiAgICAgZD0iTSA3Ljk5MzI2NCwyMy41ODUxMzMgQyA2LjkxODU4MDEsMjMuMjAxMjAyIDUuODAyMDc3OSwyMi42MDc4NjIgNC45MDg5OTQ0LDIxLjk0NjA2NyA0LjM0ODQwNDYsMjEuNTMwNjU3IDQuMDIyNDM1MSwyMS4yNDgwODUgMy41MDk2MjQ4LDIwLjczMyAxLjc4MTM2MTEsMTguOTk3MDcxIDAuNjUyOTY4MTUsMTYuODM4ODQzIDAuMjA4MzM4NzQsMTQuNDE4Nzc1IC0wLjA0MTg0ODk4LDEzLjA1NzAzMSAtMC4wMzYwNzIxOSwxMS41MDg2MzQgMC4yMjQ0NjgwNywxMC4wOTU1MDYgMC41MzY1MjgyMSw4LjQwMjk0MjEgMS4yMjY0MTE5LDYuNzg1NDk5NSAyLjI1ODk5MzgsNS4zMjU1MjkzIDIuODAxNDIyMSw0LjU1ODU4ODQgMy43NDUyMDY4LDMuNTYyMjExOSA0LjUwNDUsMi45NTQ4OTMxIDcuNzA2MzE1NCwwLjM5MzkyOTY2IDExLjk1MDIzNCwtMC4zNjA3OTQ0MiAxNS44MTY1NjQsMC45NDMxOTQyMyBjIDEuNjk3MzMsMC41NzI0NTQ0NyAzLjE5OTU1OCwxLjQ4NTQ5MDI3IDQuNDg5ODk5LDIuNzI4OTA0NjcgMS45MTk4MiwxLjg1MDAwMTUgMy4xNDIxOTYsNC4yMjczODY1IDMuNTU0MDIxLDYuOTEyMTcxMSAwLjA4NDAxLDAuNTQ3NjkgMC4wOTA1NSwwLjY2ODI1NCAwLjA5MTM1LDEuNjg1MzkzIDkuNGUtNCwxLjE5NTgwOSAtMC4wMjQ3NCwxLjQ4MTYzNCAtMC4yMTg1ODgsMi40MzI0MDYgLTAuNzA5NDk3LDMuNDc5OTY2IC0yLjk1NDA2Miw2LjQ2NjYzMiAtNi4xMzMyMzgsOC4xNjEwMTYgLTAuNzY2MTIyLDAuNDA4MzE1IC0xLjc3NDg3OSwwLjc5OTgzNiAtMi4wNjA3ODQsMC43OTk4MzYgLTAuMjI3OTg4LDAgLTAuNDE0NTQ1LC0wLjEwMDE3MSAtMC41MDA1NiwtMC4yNjg3NzMgLTAuMDY1NjEsLTAuMTI4NjExIC0wLjA2Nzg4LC0wLjE5NTU3NSAtMC4wNjgwNywtMi4wMTM1OTIgLTEuNjhlLTQsLTEuNDgyMTE4IC0wLjAxMDcsLTEuOTUwMTk5IC0wLjA0OTc0LC0yLjIwOTI2IC0wLjA4Nzk0LC0wLjU4MzQyNCAtMC4yOTUzNywtMS4wODA4MTMgLTAuNTkwNjE2LC0xLjQxNjIzMSAtMC4wNzcyNSwtMC4wODc3NiAtMC4xMjIxNDcsLTAuMTU5NTU5IC0wLjA5OTc4LC0wLjE1OTU1OSAwLjEwNTE4NywwIDAuNzk3MDUzLC0wLjExNDAzNyAxLjE0Mjc0OCwtMC4xODgzNTQgMC44NDUxMDQsLTAuMTgxNjc5IDEuNjk2MzQxLC0wLjUzNjcyNiAyLjI2MDIzOSwtMC45NDI3MzUgMC4zNDg4OTgsLTAuMjUxMjA5IDAuODQ3NTczLC0wLjc2NjgwMiAxLjA3NDA0NiwtMS4xMTA0ODQgMC40MzcyODcsLTAuNjYzNjAxIDAuNjk3NjY0LC0xLjQwODA5NCAwLjg1NDY5NSwtMi40NDM4MjEgMC4wNjk1MywtMC40NTg1NzggMC4wNTg0NiwtMS44Njg1MzIgLTAuMDE3NDMsLTIuMjIwNDk3IEMgMTkuMzcyNDExLDkuODkwMjg5OCAxOS4wODQ5NCw5LjI3ODA1NjIgMTguNTkxMzc2LDguNjU5MTc2MSBMIDE4LjQxMzY1NCw4LjQzNjMyOTcgMTguNDU5NjMyLDguMzAyODc3NiBDIDE4LjYzNTA5Myw3Ljc5MzU5NjIgMTguNjgzNzIsNi44NTg0MTg0IDE4LjU2NjIzLDYuMjUyODA5MSAxOC40NzUzNTIsNS43ODQzNzczIDE4LjMxNzA3Myw1LjI5OTQzMTYgMTguMjQxOTM2LDUuMjU5MjE5NSAxOC4wODk4LDUuMTc3Nzk5MiAxNy40NzU3NTksNS4yNjE4OTU5IDE2Ljk3NjQxLDUuNDMyNTQyMyAxNi40Mzg4NzIsNS42MTYyMzkxIDE1Ljc2MDg5Nyw1Ljk2Mzk2MiAxNS4xOTMxMDYsNi4zNDcxNzI1IGwgLTAuMTk5MDM0LDAuMTM0MzMxMyAtMC4zNzQsLTAuMDg3NDQ3IEMgMTMuNzE3MTUsNi4xODI5Mzg3IDEyLjk2MTY2Nyw2LjA5ODA1OTMgMTEuOTg1NDk0LDYuMDk4MDU5MyBjIC0wLjkyMzE5OSwwIC0xLjY2OTU4MywwLjA3ODM5MyAtMi40NTg1MjIxLDAuMjU4MjIwMSBDIDkuMzIzMDM5Myw2LjQwMjc2MjUgOS4xMTU5ODU3LDYuNDQ5Mzg4OSA5LjA2Njg1MzEsNi40NTk4OTMgOC45OTkyODMsNi40NzQzMzg4IDguOTI5NDg0MSw2LjQ0NjgzNzMgOC43ODAzMzYyLDYuMzQ3MDAxNiA4LjIzMDA4Niw1Ljk3ODY3OTUgNy4zOTQ1MTEyLDUuNTYyNjAzMSA2Ljg3MTM4NTcsNS4zOTY0MzUxIDYuNDk5ODA0NCw1LjI3ODQwNCA1Ljk4MjU4NTIsNS4yMDE4MjIxIDUuNzk1ODE1OSw1LjIzNzE4MTMgNS42ODI4MDcsNS4yNTg1NzYgNS42NjUxNTI4LDUuMjc3NDM2MyA1LjYwMzIzODksNS40NDI5MTYzIDUuMzMzMDE0Nyw2LjE2NTE1NTUgNS4yNDgxMzYsNy4xNTYwNDcgNS4zOTQ5Mzk0LDcuODc0NjQ5NSA1LjQzMzI0NzQsOC4wNjIxNjczIDUuNDg1NzQyNCw4LjI2NTg5NjMgNS41MTE1OTQ2LDguMzI3MzgwOCA1LjU1ODA2NjMsOC40Mzc5MDQ1IDUuNTU2NjY2NCw4LjQ0MTQ4ODQgNS4zODgwMzk4LDguNjQzNzQyNiA0Ljg5OTc3MjksOS4yMjkzNzkgNC41NzM3ODY5LDkuOTI2NzAxNyA0LjQxNDEwODMsMTAuNzI3MDk1IGMgLTAuMTMyMTA3NSwwLjY2MjE5MyAtMC4wOTg5NywxLjgwNjgzOCAwLjA3ODMyMywyLjcwNTQ4OSAwLjMzODQwNDQsMS43MTUyNzYgMS4zMDE1ODg3LDIuOTM2NjEgMi44MjY2NzUxLDMuNTg0MjczIDAuNjA4MjM0MiwwLjI1ODMgMS41Mzg3Njk4LDAuNDkzNTA3IDIuMjM0NDAwMywwLjU2NDc3NiBsIDAuMjI0Njc5MSwwLjAyMzAyIC0wLjE1OTMxNDgsMC4xNzQ5NSBDIDkuMzc1OTIwMiwxOC4wNDYzOTYgOS4xMjkzMjM5LDE4LjU5NzQwNSA5LjA1NjAyOTUsMTkuMDM3MjQ3IDkuMDI3NTAxOSwxOS4yMDg0NCA4Ljk3ODg5NzgsMTkuMjM5NzczIDguNTQ4Mjc2NSwxOS4zNjQ1NjIgNy4yOTkwMDIzLDE5LjcyNjU4NyA2LjMzMzA3ODIsMTkuMzYzODc0IDUuNjE2NzU2NCwxOC4yNjM3NTIgNS4xNjIzMjE4LDE3LjU2NTgzMyA0LjYzNzMzMDksMTcuMTQ2MzAzIDQuMDQzMDk3LDE3LjAwNjIxIDMuNjYwNzk1NSwxNi45MTYwOCAzLjI4MTg5NzcsMTYuOTYxNTkgMy4xODUxMzgzLDE3LjEwOTI2OCBjIC0wLjEwMTA4MzEsMC4xNTQyNzIgMC4wNTQ5NjQsMC4zNDg1OTQgMC41NDUxMzMyLDAuNjc4ODQxIDAuNDgyNTIzNywwLjMyNTA5NiAwLjg5MTc3MDUsMC44OTI1MiAxLjIyNzg3MiwxLjcwMjQ1NSAwLjM2MDQ1LDAuODY4NjEgMS4wNzg2NTgxLDEuNDQ4OTA1IDIuMDQwMTQwNSwxLjY0ODM4NiAwLjM4NDIzNjYsMC4wNzk3MiAxLjE4MjEyNTQsMC4wODY1NiAxLjYxODU3NTUsMC4wMTM4NyAwLjE2Njg1MzksLTAuMDI3NzkgMC4zMjIzMzE1LC0wLjA1MDc5IDAuMzQ1NTA1NiwtMC4wNTExMSAwLjAzMjgzOCwtNC41OGUtNCAwLjA0MjEzNSwwLjIzNjI0IDAuMDQyMTM1LDEuMDcyNzYyIDAsMS4wNTM2NDYgLTAuMDAxNDIsMS4wNzYwMjggLTAuMDc3MTgsMS4yMTkyOTkgLTAuMDg4MDcxLDAuMTY2NTQ1IC0wLjIzNzU1NzksMC4yNTkxMzIgLTAuNDYzNjA3OSwwLjI4NzE0NSAtMC4xMTc3MTcsMC4wMTQ1OSAtMC4yMTk1MTY0LC0wLjAwNjEgLTAuNDcwNDQ4MiwtMC4wOTU3OCB6IgogICAgIGlkPSJwYXRoODU3IiAvPgogIDxwYXRoCiAgICAgc3R5bGU9ImZpbGw6bm9uZTtmaWxsLW9wYWNpdHk6MTtmaWxsLXJ1bGU6bm9uemVybztzdHJva2U6bm9uZTtzdHJva2Utd2lkdGg6MC4wMzQwMTU3NTtzdHJva2UtbWl0ZXJsaW1pdDo0O3N0cm9rZS1kYXNoYXJyYXk6bm9uZTtzdHJva2Utb3BhY2l0eToxO2ZpbHRlcjp1cmwoI2ZpbHRlcjg0MikiCiAgICAgZD0iTSA4LjA5NDM4NzYsMjMuNjAwMjEzIEMgNy40NzIyNjQxLDIzLjM4OTgzNyA2LjQ1ODg1NCwyMi45MTE2MzUgNS44MzU5NjA3LDIyLjUzNDUxOSA0LjYxMDc3MDEsMjEuNzkyNzU3IDMuNTY4MzY3NSwyMC44Nzc3OTcgMi42NjI1OTg3LDE5Ljc0OTEyNyAxLjI0NDc2NDYsMTcuOTgyMzc3IDAuMzc5NTUzMTEsMTUuODk0MjggMC4xMDExODIxMiwxMy41Njc0MTYgMC4wNDQyMTM3NCwxMy4wOTEyMjUgMC4wNDMzMjc5NSwxMS41NjQyNDMgMC4wOTk3MzUxNCwxMS4wNzMwMzQgMC4zNzY2NzE1Miw4LjY2MTM5NTggMS4yOTA0NCw2LjQ4MzIzNjcgMi43NzIzNzAxLDQuNzAyMjQ3MyAzLjE1Mjk5OTcsNC4yNDQ4MDQ5IDQuMDk1NzEyOSwzLjMwNzc5MjUgNC40ODc2NDYxLDIuOTk3MzQzOSA2LjMwNTE0ODIsMS41NTc3MDc3IDguMjc0ODE5MSwwLjcyMTkxMzQ5IDEwLjU3MTkxNiwwLjQxNTU5NDEyIGMgMC42MjkwNjgsLTAuMDgzODg2NyAyLjI1MjQyOCwtMC4wODMwODgyIDIuODY1MTY4LDAuMDAxNDA5IDMuMzM4ODI5LDAuNDYwNDI4MDUgNi4xMTkzMzcsMi4wNzE0MDE1OCA4LjEyMTAzMyw0LjcwNTE3MjA4IDEuMTUzNjIyLDEuNTE3OTAwOSAxLjg5OTkyOSwzLjI3OTU3MzYgMi4yNjgxMzMsNS4zNTM5NzI4IDAuMDkxOCwwLjUxNzE2OSAwLjA5Mzk3LDAuNTU4ODM3IDAuMDkzNzEsMS43OTM1MTUgLTMuMDNlLTQsMS4zNzg5NjMgLTAuMDE2LDEuNTU1MzYyIC0wLjIyMzA1NywyLjUwNTgwOCAtMC43NTkzNSwzLjQ4NTU3MyAtMy4wMTk1MjYsNi40NTY2MzEgLTYuMTUwMDM2LDguMDg0Mzc2IC0wLjU0OTgyLDAuMjg1ODg2IC0xLjA0MDQyLDAuNTAxNTEyIC0xLjQ5NDA3MywwLjY1NjY2OCAtMC4yOTM2NzYsMC4xMDA0NDIgLTAuMzg1Mjg2LDAuMTE3Njk2IC0wLjU2NDM2OCwwLjEwNjI5MyAtMC4xODgzMDEsLTAuMDExOTkgLTAuMjI3ODMyLC0wLjAyNzIzIC0wLjMzMTY0OSwtMC4xMjc4NTQgLTAuMDY0ODMsLTAuMDYyODMgLTAuMTI1OTY1LC0wLjE2MDcxNiAtMC4xMzU4NjQsLTAuMjE3NTE5IC0wLjAwOTksLTAuMDU2OCAtMC4wMzE0OSwtMS4wMTMzOSAtMC4wNDc5NywtMi4xMjU3NSAtMC4wMzI2MSwtMi4yMDA0ODcgLTAuMDMzNTgsLTIuMjExNTkgLTAuMjM5ODExLC0yLjc0NjM1MiAtMC4wNTcyNiwtMC4xNDg0NzMgLTAuMTg2NDc1LC0wLjM4NDc2IC0wLjI4NzE1MywtMC41MjUwODMgLTAuMTAwNjc3LC0wLjE0MDMyMyAtMC4xODMwNDksLTAuMjU5ODM2IC0wLjE4MzA0OSwtMC4yNjU1ODUgMCwtMC4wMDU3IDAuMDk0OCwtMC4wMjEwNCAwLjIxMDY3NCwtMC4wMzM5OCAwLjM4OTAyMiwtMC4wNDM0NSAxLjE5OTUyOSwtMC4yMTg5NDMgMS41ODMxNjksLTAuMzQyNzg2IDAuNzk3NTYsLTAuMjU3NDYxIDEuNDAyNzMyLC0wLjU3OTc5MiAxLjg4MTcyNiwtMS4wMDIyNiAwLjg2MTAwNSwtMC43NTkzOTYgMS4zMzI0NDgsLTEuNjUxNDI5IDEuNTg3MTM5LC0zLjAwMzA3OSAwLjA4NjA2LC0wLjQ1NjczNCAwLjEyMzg0MiwtMS42OTgzNDkgMC4wNjgyNywtMi4yNDM3OTMgQyAxOS41MTM4MDEsMTAuMjAyNDQyIDE5LjIwMjMyNyw5LjQ0MTUyODIgMTguNjc0MDYsOC43NDE2MjEyIEwgMTguNDQwODgyLDguNDMyNjgwNSAxOC41MDgyNDgsOC4xODk1ODQ0IEMgMTguNzEyOTYxLDcuNDUwODUxNiAxOC42Nzk1NTYsNi40MzcxMjMzIDE4LjQyMzkxLDUuNjMwMTk1MSAxOC4zMDE4NjgsNS4yNDQ5NzY5IDE4LjI4MTAwNCw1LjIyNDcxOTIgMTguMDA2MzEyLDUuMjI0NzE5MiBjIC0wLjY3OTYyNSwwIC0xLjY4MjkxNywwLjM3NTA2NzEgLTIuNjMxNTg1LDAuOTgzNzgyMiBDIDE1LjA0MjkzOCw2LjQyMTM5MzkgMTQuOTcxOTEyLDYuNDU0Mzg3OCAxNC44ODU5NjIsNi40MzU1NDg1IDEzLjc4MTA4NCw2LjE5MzM3MDMgMTMuNDIyODU3LDYuMTQwOTAxNiAxMi42MTYwNDUsNi4xMDMwNzc3IDExLjUzMjMxMiw2LjA1MjI3MTUgMTAuNTIzNzY1LDYuMTI2NTgyNSA5LjYxMTI0MTUsNi4zMjQ0NzQ3IDguOTUxNTg4Nyw2LjQ2NzUyODUgOS4wMjE1Njg5LDYuNDczNzMzMSA4LjY0NzIyNjIsNi4yMzkwMDIgNy42NDA4NzczLDUuNjA3OTcwOSA2LjYyMzcyMDksNS4yMjQ3MTkyIDUuOTU1MzAzNCw1LjIyNDcxOTIgYyAtMC4yMzQ3NTMxLDAgLTAuMjU2Njc1LDAuMDA2MTkgLTAuMzA5OTUzLDAuMDg3NTAxIC0wLjAzMTUzMywwLjA0ODEyNiAtMC4wOTMxMiwwLjIxMTA1ODEgLTAuMTM2ODYsMC4zNjIwNzIgQyA1LjI2NDY4NjIsNi41MTYwMzg4IDUuMjUyNjc0LDcuNjQyMDc2IDUuNDgwNjA3LDguMjg4MTA3NSBMIDUuNTI4NDk5LDguNDIzODQ4NyA1LjI2ODE5OTUsOC43NzE3NzgyIGMgLTAuNDI5NjYsMC41NzQzMDUyIC0wLjcxMzgzODIsMS4yMDk1NDM0IC0wLjg3MTE0ODEsMS45NDczMjI4IC0wLjA5Mjc0NCwwLjQzNDk2NiAtMC4wODMwNjgsMS42OTUyMTUgMC4wMTgwNSwyLjM1MTAzOCAwLjI4MzA3MTIsMS44MzU5MjUgMS4xNDM2NTU5LDMuMDkxNTAzIDIuNjI3MjIzLDMuODMzMDczIDAuNTcwMDQ1NSwwLjI4NDk0MSAxLjU0MTY1NDgsMC41NjU5MjIgMi4yODI0MDAzLDAuNjYwMDUgMC4xNjY4NTM5LDAuMDIxMiAwLjMyMjg5MDYsMC4wNDQ0NSAwLjM0Njc0ODcsMC4wNTE2NyAwLjAyODQ4NCwwLjAwODYgLTAuMDA0MzksMC4wNzM1NSAtMC4wOTU3MjksMC4xODkxMDMgLTAuMjI3NDgzMywwLjI4Nzc4OCAtMC40NjkyNTI5LDAuODI3NDE5IC0wLjUyMzk0NTksMS4xNjk0NTEgLTAuMDI4MTksMC4xNzYyODkgLTAuMTE2NjM1LDAuMjUyMzI4IC0wLjM5MTEwODQsMC4zMzYyNDggQyA3Ljg1OTAzNDIsMTkuNTU0ODM4IDcuMjczMDc0NSwxOS41MzUwMjEgNi42ODAxOTk0LDE5LjI0Mjc1NCA2LjQ1NDc4MDQsMTkuMTMxNjMgNi4zMzQzNzAyLDE5LjA0MTY4NSA2LjEwNDAzOTYsMTguODEyMzcxIDUuOTQzODg4NSwxOC42NTI5MjcgNS43NDc0NzY1LDE4LjQyMzg3NiA1LjY2NzU2ODYsMTguMzAzMzcxIDUuMDY1NjIxNiwxNy4zOTU1OTkgNC40MDg2Njg4LDE2Ljk1NTA1NiAzLjY1NjkxNiwxNi45NTUwNTYgYyAtMC4yNjA1MTQ2LDAgLTAuMzcyODIsMC4wMjg1OSAtMC40NjE3OTE3LDAuMTE3NTYyIC0wLjE0MzM2NTcsMC4xNDMzNjYgLTAuMDAxNDEsMC4zNTM1MDkgMC40Njc0NTQ3LDAuNjkxOTg2IDAuNTIyNjY5NywwLjM3NzMyMiAwLjgxMDA2NCwwLjc1NDA3MiAxLjIwMjk0MzcsMS41NzY5NjUgMC4yNzQyNTA5LDAuNTc0NDIzIDAuNDE3ODg1OCwwLjc5NTI5MyAwLjY5MDYxOTUsMS4wNjE5ODIgMC4zNDE0NjE3LDAuMzMzODkzIDAuNjg5Mzk5LDAuNTM2MzgzIDEuMTc5OTMzNywwLjY4NjY5MSAwLjI4MjQyNTQsMC4wODY1NCAwLjMzMTc1MDEsMC4wOTA5NiAxLjAzODA4NzMsMC4wOTMgMC40MDc4NjUyLDAuMDAxMiAwLjg0MjMzNjYsLTAuMDEwODIgMC45NjU0OTIsLTAuMDI2NjcgbCAwLjIyMzkxOSwtMC4wMjg4MiAwLjAxMjAzNiwwLjUzNDQzOSBjIDAuMDA2NjIsMC4yOTM5NDIgMC4wMDQ0NSwwLjc4NzQ4MiAtMC4wMDQ4MiwxLjA5Njc1NiAtMC4wMTgyOSwwLjYxMDIyNCAtMC4wMzE5NDMsMC42NTU2MTMgLTAuMjQ1MDY3LDAuODE0NzM0IC0wLjEyNTIyMjMsMC4wOTM0OSAtMC4zOTkyNjc2LDAuMTA1MDA5IC0wLjYzMTMzNzUsMC4wMjY1MyB6IgogICAgIGlkPSJwYXRoODU5IiAvPgogIDxwYXRoCiAgICAgc3R5bGU9ImZpbGw6IzAwMDAwMDtmaWxsLW9wYWNpdHk6MDtmaWxsLXJ1bGU6bm9uemVybztzdHJva2U6bm9uZTtzdHJva2Utd2lkdGg6MC4wMzM3NjA5O3N0cm9rZS1taXRlcmxpbWl0OjQ7c3Ryb2tlLWRhc2hhcnJheTpub25lO3N0cm9rZS1vcGFjaXR5OjEiCiAgICAgZD0iTSA4LjAzMTk4MjUsMjMuNTc4NjE0IEMgNy4zNzk1MzcyLDIzLjM1MTU4MyA2LjQ5NTA1NjYsMjIuOTMwODIxIDUuODM3NzI3LDIyLjUzNDc3MiA0LjIzNDg1MzMsMjEuNTY5MDIxIDIuODU5Njc1NiwyMC4yMTY0OTUgMS44NTA4OTg5LDE4LjYxMzYxMyAxLjUyNTkyOSwxOC4wOTcyNTcgMS4wMTA5MTIxLDE3LjA1Njk5NCAwLjgwMjA1NTQxLDE2LjQ5NTA5NSAwLjQ4NzI3NDMxLDE1LjY0ODIyMyAwLjIzMzg5MjU1LDE0LjYxMzUyMyAwLjExODEzMDcsMTMuNzAyMjQ3IDAuMDQ0NzgzMjIsMTMuMTI0ODU3IDAuMDM0Mjg5NzgsMTEuNjM5MTIxIDAuMDk5NTYxMTMsMTEuMDczMDM0IDAuMzIyODQ4NDYsOS4xMzY1MDA3IDAuOTA3MjI2OTgsNy40NzE0NTMyIDEuOTMyMDE5Myw1Ljg1MTg4MDMgMi42NTA3OTEsNC43MTU5Mzk0IDMuNzQ0MzEyOCwzLjU0NDI1MjcgNC44NzUyODY1LDIuNjk4MjI3NSA2LjU0MTcyMDgsMS40NTE2NTExIDguNTQxMDc3OCwwLjY2MzQwMzc4IDEwLjcwNjc0NywwLjM5OTE3MjEgYyAwLjU2OTUyLC0wLjA2OTQ4NjYgMi4xNDIxOTksLTAuMDU5ODI0MyAyLjcxMzQ4MywwLjAxNjY3MTMgMS4yOTg3MSwwLjE3Mzg5ODY4IDIuNDM3MjE1LDAuNDk5MjM3MzIgMy41ODk4ODgsMS4wMjU4NDUgMi4zNTgyNjIsMS4wNzczOTA3IDQuMzgyNTE3LDMuMDA1MTI3OCA1LjU4MjkzNiw1LjMxNjczODYgMC41OTgxOTgsMS4xNTE5MzI5IDAuOTg0NDY2LDIuMzE2NDE0OSAxLjIzMzE5NiwzLjcxNzcyMSAwLjA5MTgsMC41MTcxNjkgMC4wOTM5NywwLjU1ODgzNyAwLjA5MzcxLDEuNzkzNTE1IC0zLjAzZS00LDEuMzc4OTYzIC0wLjAxNiwxLjU1NTM2MiAtMC4yMjMwNTcsMi41MDU4MDggLTAuMzE5MjU2LDEuNDY1NDUgLTAuOTAyNjI1LDIuODUxOTk0IC0xLjcxNDA2NSw0LjA3Mzk2NiAtMC40OTQwMDIsMC43NDM5MzMgLTAuODU5NDQyLDEuMTkwNDgxIC0xLjQ5OTkwOCwxLjgzMjgwOSAtMS4xNTgxMjMsMS4xNjE0ODcgLTIuMzYzODk4LDEuOTgxMDM5IC0zLjgzNDE4OCwyLjYwNjA1NiAtMC42MDU3OTMsMC4yNTc1MiAtMC45MjUzNTgsMC4zNDk0NjcgLTEuMTYyMTgzLDAuMzM0Mzg3IC0wLjE4NjIyOSwtMC4wMTE4NiAtMC4yMjYyMDIsLTAuMDI3MzQgLTAuMzI5NzgyLC0wLjEyNzczNSAtMC4wNjQ4MywtMC4wNjI4MyAtMC4xMjYxNjksLTAuMTYwNzE2IC0wLjEzNjMxNiwtMC4yMTc1MTkgLTAuMDEwMTUsLTAuMDU2OCAtMC4wMzE0NywtMS4wMTMzOSAtMC4wNDczOSwtMi4xMjU3NSAtMC4wMjg3OSwtMi4wMTIxMjUgLTAuMDI5MzgsLTIuMDI0MTQ4IC0wLjExNDU0MiwtMi4zNTAxNDEgLTAuMTAyNTA3LC0wLjM5MjM3NyAtMC4yMzc5MjMsLTAuNjg2NjQ2IC0wLjQ0NjEyMywtMC45Njk0NjYgLTAuMDg0NzksLTAuMTE1MTc5IC0wLjE0OTMxNywtMC4yMTI5NzkgLTAuMTQzMzk0LC0wLjIxNzMzMiAwLjAwNTksLTAuMDA0NCAwLjEzOTcwMSwtMC4wMjY1MSAwLjI5NzI4NiwtMC4wNDkyNCAxLjIzMDQxOCwtMC4xNzc0NiAyLjI4NDg1NCwtMC41NDUzMTUgMy4wMjI0MDQsLTEuMDU0NDEgMC4zMjc3NjcsLTAuMjI2MjQyIDAuODU1NTQ0LC0wLjc1MzgzMyAxLjA4MjQ2MSwtMS4wODIwODEgMC40ODIxNTYsLTAuNjk3NDYyIDAuODE2NjI2LC0xLjY4MTgxIDAuOTIyMzQ3LC0yLjcxNDQ3MSAwLjA1MDg2LC0wLjQ5Njc5MSAwLjAxNjc1LC0xLjgzOTUwMSAtMC4wNTQxMiwtMi4xMzAyNzUgQyAxOS4zNTg0MDMsOS44NDE2NTg1IDE5LjA2Njc4NSw5LjI0MTM5NjUgMTguNTg4NDYsOC42MjY4ODU3IEwgMTguNDQyMzQsOC40MzkxNjQyIDE4LjUxMDg1OSw4LjE3MTgyOTUgQyAxOC42NjA2NzEsNy41ODczMTYyIDE4LjY4MTg5LDYuNzI5MDAxNyAxOC41NjA3NzYsNi4xNTI2NjcgMTguNDg2Niw1Ljc5OTY4NyAxOC4zNTc3MDYsNS4zOTY1NjYzIDE4LjI4Njg2MSw1LjI5NTk4MjMgYyAtMC4wNTUwNSwtMC4wNzgxNTkgLTAuMDY2MjIsLTAuMDc5ODc2IC0wLjM5NzAzMSwtMC4wNjEwMDIgLTAuNDEyMjQzLDAuMDIzNTIgLTAuODM5ODA3LDAuMTM1MzgzMyAtMS4zNTA0MzUsMC4zNTMzMTI0IC0wLjM3MzQyNCwwLjE1OTM3MjggLTEuMDA5OTgyLDAuNTA3NTY4IC0xLjM2MjcxNSwwLjc0NTQwMDIgTCAxNC45ODc2NDYsNi40NjExNTEgMTQuNTY2Mjk4LDYuMzY2NTg0NiBDIDEzLjU1MTMxNCw2LjEzODc4MzggMTIuNDg4Mjg4LDYuMDQ5OTU4MiAxMS4zOTc3NTgsNi4xMDE4MjQ4IDEwLjY2NzYzOSw2LjEzNjU0OTYgMTAuMzA2MzQ2LDYuMTgwNTA4NyA5LjY5ODMxNzcsNi4zMDg1OTg2IDguOTMxNzE2Niw2LjQ3MDA5MzcgOS4wNDMwODc0LDYuNDgyMjUwMSA4LjU0NjczNCw2LjE4MjkwMSA3LjQ5MDYzMDksNS41NDU5Njc2IDYuNTQ0NTUxNiw1LjIwNDA3NDggNS45MTQxODA4LDUuMjMxNTUxOCA1LjY4NjEzMjYsNS4yNDE0OTIyIDUuNjgzNzg0MSw1LjI0MjU3NzYgNS42MjM0OTksNS4zNjU5MDc5IDUuMzQ3NzAzOSw1LjkzMDEyNiA1LjIzMjQ2MDEsNy4xOTM1NTY1IDUuMzkwNDA1MSw3LjkyMTM0ODQgNS40Mjg2MjcxLDguMDk3NDcyIDUuNDc1NDkyOCw4LjI4MjI2ODYgNS40OTQ1NTA2LDguMzMyMDA3NiA1LjUyNjA1OCw4LjQxNDIzODIgNS41MDgyMTUsOC40NTA2NTc5IDUuMjk3ODU4OCw4LjczMzQ4MSA0Ljg4NDA4OSw5LjI4OTc5MzMgNC42NDE3NzI5LDkuNzgzNDgzNCA0LjQ3MzA1ODMsMTAuNDEzOTE1IGMgLTAuMTIzNTA4NiwwLjQ2MTUxMiAtMC4xNjU5Nzg5LDAuOTM3MjgxIC0wLjE0MzA4ODIsMS42MDI5MzkgMC4wNjI0NjUsMS44MTY0ODIgMC41MDk1NjgyLDMuMDI5ODIzIDEuNDcyMjgyNywzLjk5NTQ1NCAwLjM4Mzc4MDksMC4zODQ5NDQgMC43MzkxNTM4LDAuNjM5NDQzIDEuMjQ3MTkxLDAuODkzMTcgMC41NjYxMjEyLDAuMjgyNzM3IDEuNTQwNDAxMiwwLjU2NDQwMSAyLjI3NTI4MDksMC42NTc3ODQgMC4xNjY4NTM5LDAuMDIxMiAwLjMyMjg5MDYsMC4wNDQ0NSAwLjM0Njc0ODcsMC4wNTE2NyAwLjAyODQ4NCwwLjAwODYgLTAuMDA0MzksMC4wNzM1NSAtMC4wOTU3MjksMC4xODkxMDMgLTAuMDc2NTA5LDAuMDk2NzkgLTAuMTg5MDM5NCwwLjI3NTI4MyAtMC4yNTAwNjc1LDAuMzk2NjQ4IC0wLjExNDAwMDQsMC4yMjY3MSAtMC4yODc0Njg4LDAuNzMzNjcgLTAuMjg3NDY4OCwwLjg0MDEyNSAwLDAuMTQ4MTIzIC0wLjI5NjMxLDAuMjc5NzI5IC0wLjg3ODQzMjMsMC4zOTAxNTYgLTAuNTQ4NzEyOCwwLjEwNDA4OSAtMS4wMDQxNjEsMC4wNDYxNSAtMS40Nzk1NzY0LC0wLjE4ODIxIEMgNi4zMjY3NzkyLDE5LjA2ODUyOSA1LjkxODUxNTMsMTguNjc3MDU1IDUuNTk4NDI0NywxOC4yMDU0NjYgNS4yMTczNzQ4LDE3LjY0NDA2NyA0Ljc5MzM5MTIsMTcuMjc2OTIyIDQuMjc3OTQyNCwxNy4wNjIwMDYgYyAtMC4xODQ2NTQxLC0wLjA3Njk5IC0wLjI1OTY1MjEsLTAuMDg4MzYgLTAuNTg4MDkzMSwtMC4wODkxNyAtMC4zNTI2NTAyLC04LjY2ZS00IC0wLjM4MTk1MDIsMC4wMDQyIC0wLjQ2MzQ4MywwLjA4MDc5IC0wLjA0Nzg2OSwwLjA0NDk1IC0wLjA4NzAzNSwwLjExNjk1MyAtMC4wODcwMzUsMC4xNjAwMTYgMCwwLjEwODQwMyAwLjE5NTc5OTQsMC4zMTQ1NzMgMC41MjMyNDcyLDAuNTUwOTYgMC41MTY3NzE1LDAuMzczMDYzIDAuODExMTgwNywwLjc1NjIwNyAxLjE4ODIwMDUsMS41NDYzMiAwLjI2NzAxMzIsMC41NTk1NzYgMC4zOTYyMjg1LDAuNzY1NDk5IDAuNjQzMjY1MSwxLjAyNTEzNiAwLjM1NjE5NiwwLjM3NDM2NiAwLjc2NzYwNzYsMC42MTY5OTUgMS4zMDk3MTk3LDAuNzcyNDAzIDAuMjMzNjc4MSwwLjA2Njk5IDAuMzQwMzM2NSwwLjA3NTM1IDAuOTcwMzk4OSwwLjA3NjEgMC4zODkzMjU4LDQuNThlLTQgMC44MTU3MDEzLC0wLjAxMTg3IDAuOTQ3NTAwOCwtMC4wMjczOSBsIDAuMjM5NjM1OSwtMC4wMjgyMiAwLjAxMzk5OSwwLjU2NzU1MSBjIDAuMDA3NywwLjMxMjE1MyAwLjAwNjA0LDAuODA1NDI3IC0wLjAwMzY4LDEuMDk2MTY0IC0wLjAxNTU1LDAuNDY0OTQzIC0wLjAyNjA1MSwwLjUzOTgwMyAtMC4wODcxNzgsMC42MjE1MTMgLTAuMTkxODg3NCwwLjI1NjQ5OCAtMC40NDY2NjczLDAuMzA1NjQ2IC0wLjg1MjQ1NzEsMC4xNjQ0NDMgeiIKICAgICBpZD0icGF0aDg2MSIgLz4KPC9zdmc+Cg==" /></a> <a href="https://www.linkedin.com/in/javier-ignacio-tamayo-leiva-94613267/"><img height="30" width="30" src="data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiIHN0YW5kYWxvbmU9Im5vIj8+CjxzdmcKICAgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIgogICB4bWxuczpjYz0iaHR0cDovL2NyZWF0aXZlY29tbW9ucy5vcmcvbnMjIgogICB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiCiAgIHhtbG5zOnN2Zz0iaHR0cDovL3d3dy53My5vcmcvMjAwMC9zdmciCiAgIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIKICAgeG1sbnM6c29kaXBvZGk9Imh0dHA6Ly9zb2RpcG9kaS5zb3VyY2Vmb3JnZS5uZXQvRFREL3NvZGlwb2RpLTAuZHRkIgogICB4bWxuczppbmtzY2FwZT0iaHR0cDovL3d3dy5pbmtzY2FwZS5vcmcvbmFtZXNwYWNlcy9pbmtzY2FwZSIKICAgcm9sZT0iaW1nIgogICB2aWV3Qm94PSIwIDAgMjQgMjQiCiAgIHZlcnNpb249IjEuMSIKICAgaWQ9InN2ZzgzNyIKICAgc29kaXBvZGk6ZG9jbmFtZT0ibGlua2VkaW4uc3ZnIgogICBpbmtzY2FwZTp2ZXJzaW9uPSIxLjAgKDQwMzVhNGYsIDIwMjAtMDUtMDEpIj4KICA8bWV0YWRhdGEKICAgICBpZD0ibWV0YWRhdGE4NDMiPgogICAgPHJkZjpSREY+CiAgICAgIDxjYzpXb3JrCiAgICAgICAgIHJkZjphYm91dD0iIj4KICAgICAgICA8ZGM6Zm9ybWF0PmltYWdlL3N2Zyt4bWw8L2RjOmZvcm1hdD4KICAgICAgICA8ZGM6dHlwZQogICAgICAgICAgIHJkZjpyZXNvdXJjZT0iaHR0cDovL3B1cmwub3JnL2RjL2RjbWl0eXBlL1N0aWxsSW1hZ2UiIC8+CiAgICAgIDwvY2M6V29yaz4KICAgIDwvcmRmOlJERj4KICA8L21ldGFkYXRhPgogIDxkZWZzCiAgICAgaWQ9ImRlZnM4NDEiIC8+CiAgPHNvZGlwb2RpOm5hbWVkdmlldwogICAgIGlua3NjYXBlOmRvY3VtZW50LXJvdGF0aW9uPSIwIgogICAgIHBhZ2Vjb2xvcj0iI2ZmZmZmZiIKICAgICBib3JkZXJjb2xvcj0iIzY2NjY2NiIKICAgICBib3JkZXJvcGFjaXR5PSIxIgogICAgIG9iamVjdHRvbGVyYW5jZT0iMTAiCiAgICAgZ3JpZHRvbGVyYW5jZT0iMTAiCiAgICAgZ3VpZGV0b2xlcmFuY2U9IjEwIgogICAgIGlua3NjYXBlOnBhZ2VvcGFjaXR5PSIwIgogICAgIGlua3NjYXBlOnBhZ2VzaGFkb3c9IjIiCiAgICAgaW5rc2NhcGU6d2luZG93LXdpZHRoPSIxNTkzIgogICAgIGlua3NjYXBlOndpbmRvdy1oZWlnaHQ9IjgwNyIKICAgICBpZD0ibmFtZWR2aWV3ODM5IgogICAgIHNob3dncmlkPSJmYWxzZSIKICAgICBpbmtzY2FwZTp6b29tPSIyOS42NjY2NjciCiAgICAgaW5rc2NhcGU6Y3g9IjExLjA4OTg4OCIKICAgICBpbmtzY2FwZTpjeT0iMC40MzgyMDIyNSIKICAgICBpbmtzY2FwZTp3aW5kb3cteD0iMCIKICAgICBpbmtzY2FwZTp3aW5kb3cteT0iMjMiCiAgICAgaW5rc2NhcGU6d2luZG93LW1heGltaXplZD0iMCIKICAgICBpbmtzY2FwZTpjdXJyZW50LWxheWVyPSJzdmc4MzciIC8+CiAgPHRpdGxlCiAgICAgaWQ9InRpdGxlODMzIj5MaW5rZWRJbiBpY29uPC90aXRsZT4KICA8cGF0aAogICAgIGQ9Ik0yMC40NDcgMjAuNDUyaC0zLjU1NHYtNS41NjljMC0xLjMyOC0uMDI3LTMuMDM3LTEuODUyLTMuMDM3LTEuODUzIDAtMi4xMzYgMS40NDUtMi4xMzYgMi45Mzl2NS42NjdIOS4zNTFWOWgzLjQxNHYxLjU2MWguMDQ2Yy40NzctLjkgMS42MzctMS44NSAzLjM3LTEuODUgMy42MDEgMCA0LjI2NyAyLjM3IDQuMjY3IDUuNDU1djYuMjg2ek01LjMzNyA3LjQzM2MtMS4xNDQgMC0yLjA2My0uOTI2LTIuMDYzLTIuMDY1IDAtMS4xMzguOTItMi4wNjMgMi4wNjMtMi4wNjMgMS4xNCAwIDIuMDY0LjkyNSAyLjA2NCAyLjA2MyAwIDEuMTM5LS45MjUgMi4wNjUtMi4wNjQgMi4wNjV6bTEuNzgyIDEzLjAxOUgzLjU1NVY5aDMuNTY0djExLjQ1MnpNMjIuMjI1IDBIMS43NzFDLjc5MiAwIDAgLjc3NCAwIDEuNzI5djIwLjU0MkMwIDIzLjIyNy43OTIgMjQgMS43NzEgMjRoMjAuNDUxQzIzLjIgMjQgMjQgMjMuMjI3IDI0IDIyLjI3MVYxLjcyOUMyNCAuNzc0IDIzLjIgMCAyMi4yMjIgMGguMDAzeiIKICAgICBpZD0icGF0aDgzNSIKICAgICBzdHlsZT0iZmlsbDojOGFiMGJmO2ZpbGwtb3BhY2l0eToxIiAvPgo8L3N2Zz4K" /></a>
### Pontificia Universidad Católica de Valparaíso
### Junio 22, 2021

---

## Análisis exploratorio de datos (*Exploratory data analysis*)
.pull-left[
### Estadística descriptiva
- {gtsummary} package 
- {ggpubr} package 
- Visualizando distribuciones 
 - Histograms, Density plot, Heatmaps 
 - Medidas de Tendencia 
 - q-q plots, ECDF plots 
 - Boxplot, Violin plots 
- Análisis de normalidad 
 - Shapiro–Wilk test 
 - Kolmogorov-Smirnov test 
- Análisis de Homocedasticidad 
 - Levene test 
 - Bartlett’s test 
 - Fligner-Killeen test 
]

.pull-right[
### Inferencia Estadística
- Análisis de correlación 
 - Pearson's *r* 
 - Spearman's *ρ* (rho) 
 - Kendall's *τ* (tau) 
- Reducción de dimensión 
 - Análisis de componentes principales (*PCA*) 
- Modelos Lineales 
 - `lm()` 
- Análisis de varianza 
 - F-test, T-test, chi-square 
 - ANOVA test 
 - Kruskal-Wallis test 
 - Scheirer-Ray-Hare test 
- {report} package 
]

---
## Dependencias de esta clase

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages(c("tidyverse","ggplot2","gtsummary","car","flextable","RColorBrewer","ggpubr",
                   "Hmisc","corrplot","broom","scales","viridis","ggrepel","ggbiplot",
                   "rcompanion","report","FSA"))
```

```r
# Cargar desde la librería
sapply(c("tidyverse","ggplot2","gtsummary","car","flextable","RColorBrewer","ggpubr",
 "Hmisc","corrplot","broom","scales","viridis","ggrepel","ggbiplot","rcompanion",
 "report","FSA"), 
 require, character.only = TRUE, quietly = TRUE)
```
 
.bg-washed-green.b--dark-green.ba.bw2.br3.shadow-5.ph4.mt4[
Si bien los paquetes más importantes de esta sesión se irán mostrando durante la presentación, te recomiendo instalar todos los paquetes y cargar las librerías para poder ir siguiendo en paralelo en tú sesión de R los ejercicios.

.tr[
— _Consejos_
]]

---
class: inverse center middle

# Análisis exploratorio de datos

---
## Análisis exploratorio de datos

El análisis exploratorio de datos o “EDA” (por sus siglas en inglés **E**xploratory **D**ata **A**nalysis), es como se conoce -en estadística- el proceso por el cual un investigador inspecciona un set de datos con la finalidad de generar preguntas, procesar y adquirir conocimiento (procesa datos, genera resultados), y refina o genera nuevas preguntas.

### Proceso 
1.	Generar preguntas basandose en los datos. 
1.	Obtener resultados (procesa datos, gráficos, modelos, etc.). 
1.	Refinar preguntas y/o generar nuevas preguntas.

Si bien el nombre puede generar la impresión de un proceso estandarizado, en la práctica no hay reglas que limiten los análisis o técnicas utilizadas para generar el proceso iterativo. Sin embargo, existen análisis que se vuelven recurrentes durante el proceso.

.footnote[Análisis exploratorio de datos [(EDA) R4DS-ES](https://es.r4ds.hadley.nz/análisis-exploratorio-de-datos-eda.html)]

---
class: inverse center middle

## Estadística descriptiva

---
## Los datos (ggplot2::diamonds)

<div id="htmlwidget-0a17acc57d6032340e1d" style="width:100%;height:auto;" class="datatables html-widget"></div>
<script type="application/json" data-for="htmlwidget-0a17acc57d6032340e1d">{"x":{"filter":"none","fillContainer":false,"data":[["1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20"],[0.23,0.21,0.23,0.29,0.31,0.24,0.24,0.26,0.22,0.23,0.3,0.23,0.22,0.31,0.2,0.32,0.3,0.3,0.3,0.3],["Ideal","Premium","Good","Premium","Good","Very Good","Very Good","Very Good","Fair","Very Good","Good","Ideal","Premium","Ideal","Premium","Premium","Ideal","Good","Good","Very Good"],["E","E","E","I","J","J","I","H","E","H","J","J","F","J","E","E","I","J","J","J"],["SI2","SI1","VS1","VS2","SI2","VVS2","VVS1","SI1","VS2","VS1","SI1","VS1","SI1","SI2","SI2","I1","SI2","SI1","SI1","SI1"],[61.5,59.8,56.9,62.4,63.3,62.8,62.3,61.9,65.1,59.4,64,62.8,60.4,62.2,60.2,60.9,62,63.4,63.8,62.7],[55,61,65,58,58,57,57,55,61,61,55,56,61,54,62,58,54,54,56,59],[326,326,327,334,335,336,336,337,337,338,339,340,342,344,345,345,348,351,351,351],[3.95,3.89,4.05,4.2,4.34,3.94,3.95,4.07,3.87,4,4.25,3.93,3.88,4.35,3.79,4.38,4.31,4.23,4.23,4.21],[3.98,3.84,4.07,4.23,4.35,3.96,3.98,4.11,3.78,4.05,4.28,3.9,3.84,4.37,3.75,4.42,4.34,4.29,4.26,4.27],[2.43,2.31,2.31,2.63,2.75,2.48,2.47,2.53,2.49,2.39,2.73,2.46,2.33,2.71,2.27,2.68,2.68,2.7,2.71,2.66]],"container":"<table class=\"display\">\n <thead>\n <tr>\n <th> <\/th>\n <th>carat<\/th>\n <th>cut<\/th>\n <th>color<\/th>\n <th>clarity<\/th>\n <th>depth<\/th>\n <th>table<\/th>\n <th>price<\/th>\n <th>x<\/th>\n <th>y<\/th>\n <th>z<\/th>\n <\/tr>\n <\/thead>\n<\/table>","options":{"pageLength":8,"columnDefs":[{"className":"dt-right","targets":[1,5,6,7,8,9,10]},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false,"lengthMenu":[8,10,25,50,100]}},"evals":[],"jsHooks":[]}</script>

Precio de mas de 50,000 diamantes cortados

---
## paquete-**gtsummary**

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("gtsummary")

# Instalar la versión beta desde GitHub
install.packages("remotes")
remotes::install_github("ddsjoberg/gtsummary")

# Para guardar tablas en formato PDF, Word o PowerPoint 
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("flextable")

# Instalar la versión beta desde GitHub
devtools::install_github("davidgohel/flextable")
```

```r
# Cargar desde la librería
library(gtsummary)
library(flextable)
```

.footnote[[{gtsummary}](http://www.danieldsjoberg.com/gtsummary/index.html) package 
 [{gtsummary}](http://www.danieldsjoberg.com/gtsummary/reference/index.html) package reference]

---
## Resumir data con **gtsummary**

```r
dmnd2 <- diamonds %>% select(cut) # Seleccionar solo columnas "class" desde mpg. Se pueden seleccionar 
 # varias columnas en "select()". Ej select(class, year, ...)
 # También se pueden usar todas las columnas sí se saltan este paso.
tbl_summary(dmnd2) %>% modify_header(label = "**Samples**") # Resumir data
```

<div id="kxmfasuojh" style="overflow-x:auto;overflow-y:auto;width:auto;height:auto;">
<style>html {
 font-family: -apple-system, BlinkMacSystemFont, 'Segoe UI', Roboto, Oxygen, Ubuntu, Cantarell, 'Helvetica Neue', 'Fira Sans', 'Droid Sans', Arial, sans-serif;
}

#kxmfasuojh .gt_table {
  display: table;
  border-collapse: collapse;
  margin-left: auto;
  margin-right: auto;
  color: #333333;
  font-size: 16px;
  font-weight: normal;
  font-style: normal;
  background-color: #FFFFFF;
  width: auto;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #A8A8A8;
  border-right-style: none;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #A8A8A8;
  border-left-style: none;
  border-left-width: 2px;
  border-left-color: #D3D3D3;
}

#kxmfasuojh .gt_heading {
  background-color: #FFFFFF;
  text-align: center;
  border-bottom-color: #FFFFFF;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
}

#kxmfasuojh .gt_title {
  color: #333333;
  font-size: 125%;
  font-weight: initial;
  padding-top: 4px;
  padding-bottom: 4px;
  border-bottom-color: #FFFFFF;
  border-bottom-width: 0;
}

#kxmfasuojh .gt_subtitle {
  color: #333333;
  font-size: 85%;
  font-weight: initial;
  padding-top: 0;
  padding-bottom: 4px;
  border-top-color: #FFFFFF;
  border-top-width: 0;
}

#kxmfasuojh .gt_bottom_border {
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
}

#kxmfasuojh .gt_col_headings {
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
}

#kxmfasuojh .gt_col_heading {
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: normal;
  text-transform: inherit;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
  vertical-align: bottom;
  padding-top: 5px;
  padding-bottom: 6px;
  padding-left: 5px;
  padding-right: 5px;
  overflow-x: hidden;
}

#kxmfasuojh .gt_column_spanner_outer {
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: normal;
  text-transform: inherit;
  padding-top: 0;
  padding-bottom: 0;
  padding-left: 4px;
  padding-right: 4px;
}

#kxmfasuojh .gt_column_spanner_outer:first-child {
  padding-left: 0;
}

#kxmfasuojh .gt_column_spanner_outer:last-child {
  padding-right: 0;
}

#kxmfasuojh .gt_column_spanner {
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  vertical-align: bottom;
  padding-top: 5px;
  padding-bottom: 6px;
  overflow-x: hidden;
  display: inline-block;
  width: 100%;
}

#kxmfasuojh .gt_group_heading {
  padding: 8px;
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: initial;
  text-transform: inherit;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
  vertical-align: middle;
}

#kxmfasuojh .gt_empty_group_heading {
  padding: 0.5px;
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: initial;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  vertical-align: middle;
}

#kxmfasuojh .gt_from_md > :first-child {
  margin-top: 0;
}

#kxmfasuojh .gt_from_md > :last-child {
  margin-bottom: 0;
}

#kxmfasuojh .gt_row {
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
  margin: 10px;
  border-top-style: solid;
  border-top-width: 1px;
  border-top-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
  vertical-align: middle;
  overflow-x: hidden;
}

#kxmfasuojh .gt_stub {
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: initial;
  text-transform: inherit;
  border-right-style: solid;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
  padding-left: 12px;
}

#kxmfasuojh .gt_summary_row {
  color: #333333;
  background-color: #FFFFFF;
  text-transform: inherit;
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
}

#kxmfasuojh .gt_first_summary_row {
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
}

#kxmfasuojh .gt_grand_summary_row {
  color: #333333;
  background-color: #FFFFFF;
  text-transform: inherit;
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
}

#kxmfasuojh .gt_first_grand_summary_row {
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
  border-top-style: double;
  border-top-width: 6px;
  border-top-color: #D3D3D3;
}

#kxmfasuojh .gt_striped {
  background-color: rgba(128, 128, 128, 0.05);
}

#kxmfasuojh .gt_table_body {
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
}

#kxmfasuojh .gt_footnotes {
  color: #333333;
  background-color: #FFFFFF;
  border-bottom-style: none;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 2px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
}

#kxmfasuojh .gt_footnote {
  margin: 0px;
  font-size: 90%;
  padding: 4px;
}

#kxmfasuojh .gt_sourcenotes {
  color: #333333;
  background-color: #FFFFFF;
  border-bottom-style: none;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 2px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
}

#kxmfasuojh .gt_sourcenote {
  font-size: 90%;
  padding: 4px;
}

#kxmfasuojh .gt_left {
  text-align: left;
}

#kxmfasuojh .gt_center {
  text-align: center;
}

#kxmfasuojh .gt_right {
  text-align: right;
  font-variant-numeric: tabular-nums;
}

#kxmfasuojh .gt_font_normal {
  font-weight: normal;
}

#kxmfasuojh .gt_font_bold {
  font-weight: bold;
}

#kxmfasuojh .gt_font_italic {
  font-style: italic;
}

#kxmfasuojh .gt_super {
  font-size: 65%;
}

#kxmfasuojh .gt_footnote_marks {
 font-style: italic;
 font-weight: normal;
 font-size: 65%;
}
</style>
<table class="gt_table">
 
 <thead class="gt_col_headings">
 <tr>
 <th class="gt_col_heading gt_columns_bottom_border gt_left" rowspan="1" colspan="1">Samples</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">N = 53,9401</th>
 </tr>
 </thead>
 <tbody class="gt_table_body">
 <tr><td class="gt_row gt_left">cut</td>
<td class="gt_row gt_center"></td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Fair</td>
<td class="gt_row gt_center">1,610 (3.0%)</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Good</td>
<td class="gt_row gt_center">4,906 (9.1%)</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Very Good</td>
<td class="gt_row gt_center">12,082 (22%)</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Premium</td>
<td class="gt_row gt_center">13,791 (26%)</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Ideal</td>
<td class="gt_row gt_center">21,551 (40%)</td></tr>
 </tbody>
 
 <tfoot>
 <tr class="gt_footnotes">
 <td colspan="2">
 
 
 1
 
 
 n (%)
 
 
 </td>
 </tr>
 </tfoot>
</table>
</div>

---
## Resumir data con **gtsummary**

```r
# Crear tabla resumen 
diamonds %>% tbl_cross(row = cut, col = color) %>% bold_labels() 
```
.medium[
<div id="klyrhnejfm" style="overflow-x:auto;overflow-y:auto;width:auto;height:auto;">
<style>html {
 font-family: -apple-system, BlinkMacSystemFont, 'Segoe UI', Roboto, Oxygen, Ubuntu, Cantarell, 'Helvetica Neue', 'Fira Sans', 'Droid Sans', Arial, sans-serif;
}

#klyrhnejfm .gt_table {
  display: table;
  border-collapse: collapse;
  margin-left: auto;
  margin-right: auto;
  color: #333333;
  font-size: 16px;
  font-weight: normal;
  font-style: normal;
  background-color: #FFFFFF;
  width: auto;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #A8A8A8;
  border-right-style: none;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #A8A8A8;
  border-left-style: none;
  border-left-width: 2px;
  border-left-color: #D3D3D3;
}

#klyrhnejfm .gt_heading {
  background-color: #FFFFFF;
  text-align: center;
  border-bottom-color: #FFFFFF;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
}

#klyrhnejfm .gt_title {
  color: #333333;
  font-size: 125%;
  font-weight: initial;
  padding-top: 4px;
  padding-bottom: 4px;
  border-bottom-color: #FFFFFF;
  border-bottom-width: 0;
}

#klyrhnejfm .gt_subtitle {
  color: #333333;
  font-size: 85%;
  font-weight: initial;
  padding-top: 0;
  padding-bottom: 4px;
  border-top-color: #FFFFFF;
  border-top-width: 0;
}

#klyrhnejfm .gt_bottom_border {
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
}

#klyrhnejfm .gt_col_headings {
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
}

#klyrhnejfm .gt_col_heading {
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: normal;
  text-transform: inherit;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
  vertical-align: bottom;
  padding-top: 5px;
  padding-bottom: 6px;
  padding-left: 5px;
  padding-right: 5px;
  overflow-x: hidden;
}

#klyrhnejfm .gt_column_spanner_outer {
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: normal;
  text-transform: inherit;
  padding-top: 0;
  padding-bottom: 0;
  padding-left: 4px;
  padding-right: 4px;
}

#klyrhnejfm .gt_column_spanner_outer:first-child {
  padding-left: 0;
}

#klyrhnejfm .gt_column_spanner_outer:last-child {
  padding-right: 0;
}

#klyrhnejfm .gt_column_spanner {
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  vertical-align: bottom;
  padding-top: 5px;
  padding-bottom: 6px;
  overflow-x: hidden;
  display: inline-block;
  width: 100%;
}

#klyrhnejfm .gt_group_heading {
  padding: 8px;
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: initial;
  text-transform: inherit;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
  vertical-align: middle;
}

#klyrhnejfm .gt_empty_group_heading {
  padding: 0.5px;
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: initial;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  vertical-align: middle;
}

#klyrhnejfm .gt_from_md > :first-child {
  margin-top: 0;
}

#klyrhnejfm .gt_from_md > :last-child {
  margin-bottom: 0;
}

#klyrhnejfm .gt_row {
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
  margin: 10px;
  border-top-style: solid;
  border-top-width: 1px;
  border-top-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 1px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 1px;
  border-right-color: #D3D3D3;
  vertical-align: middle;
  overflow-x: hidden;
}

#klyrhnejfm .gt_stub {
  color: #333333;
  background-color: #FFFFFF;
  font-size: 100%;
  font-weight: initial;
  text-transform: inherit;
  border-right-style: solid;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
  padding-left: 12px;
}

#klyrhnejfm .gt_summary_row {
  color: #333333;
  background-color: #FFFFFF;
  text-transform: inherit;
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
}

#klyrhnejfm .gt_first_summary_row {
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
}

#klyrhnejfm .gt_grand_summary_row {
  color: #333333;
  background-color: #FFFFFF;
  text-transform: inherit;
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
}

#klyrhnejfm .gt_first_grand_summary_row {
  padding-top: 8px;
  padding-bottom: 8px;
  padding-left: 5px;
  padding-right: 5px;
  border-top-style: double;
  border-top-width: 6px;
  border-top-color: #D3D3D3;
}

#klyrhnejfm .gt_striped {
  background-color: rgba(128, 128, 128, 0.05);
}

#klyrhnejfm .gt_table_body {
  border-top-style: solid;
  border-top-width: 2px;
  border-top-color: #D3D3D3;
  border-bottom-style: solid;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
}

#klyrhnejfm .gt_footnotes {
  color: #333333;
  background-color: #FFFFFF;
  border-bottom-style: none;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 2px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
}

#klyrhnejfm .gt_footnote {
  margin: 0px;
  font-size: 90%;
  padding: 4px;
}

#klyrhnejfm .gt_sourcenotes {
  color: #333333;
  background-color: #FFFFFF;
  border-bottom-style: none;
  border-bottom-width: 2px;
  border-bottom-color: #D3D3D3;
  border-left-style: none;
  border-left-width: 2px;
  border-left-color: #D3D3D3;
  border-right-style: none;
  border-right-width: 2px;
  border-right-color: #D3D3D3;
}

#klyrhnejfm .gt_sourcenote {
  font-size: 90%;
  padding: 4px;
}

#klyrhnejfm .gt_left {
  text-align: left;
}

#klyrhnejfm .gt_center {
  text-align: center;
}

#klyrhnejfm .gt_right {
  text-align: right;
  font-variant-numeric: tabular-nums;
}

#klyrhnejfm .gt_font_normal {
  font-weight: normal;
}

#klyrhnejfm .gt_font_bold {
  font-weight: bold;
}

#klyrhnejfm .gt_font_italic {
  font-style: italic;
}

#klyrhnejfm .gt_super {
  font-size: 65%;
}

#klyrhnejfm .gt_footnote_marks {
 font-style: italic;
 font-weight: normal;
 font-size: 65%;
}
</style>
<table class="gt_table">
 
 <thead class="gt_col_headings">
 <tr>
 <th class="gt_col_heading gt_columns_bottom_border gt_left" rowspan="2" colspan="1">Characteristic</th>
 <th class="gt_center gt_columns_top_border gt_column_spanner_outer" rowspan="1" colspan="7">
 color
 </th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="2" colspan="1">Total</th>
 </tr>
 <tr>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">D</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">E</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">F</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">G</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">H</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">I</th>
 <th class="gt_col_heading gt_columns_bottom_border gt_center" rowspan="1" colspan="1">J</th>
 </tr>
 </thead>
 <tbody class="gt_table_body">
 <tr><td class="gt_row gt_left" style="font-weight: bold;">cut</td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td>
<td class="gt_row gt_center"></td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Fair</td>
<td class="gt_row gt_center">163</td>
<td class="gt_row gt_center">224</td>
<td class="gt_row gt_center">312</td>
<td class="gt_row gt_center">314</td>
<td class="gt_row gt_center">303</td>
<td class="gt_row gt_center">175</td>
<td class="gt_row gt_center">119</td>
<td class="gt_row gt_center">1,610</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Good</td>
<td class="gt_row gt_center">662</td>
<td class="gt_row gt_center">933</td>
<td class="gt_row gt_center">909</td>
<td class="gt_row gt_center">871</td>
<td class="gt_row gt_center">702</td>
<td class="gt_row gt_center">522</td>
<td class="gt_row gt_center">307</td>
<td class="gt_row gt_center">4,906</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Very Good</td>
<td class="gt_row gt_center">1,513</td>
<td class="gt_row gt_center">2,400</td>
<td class="gt_row gt_center">2,164</td>
<td class="gt_row gt_center">2,299</td>
<td class="gt_row gt_center">1,824</td>
<td class="gt_row gt_center">1,204</td>
<td class="gt_row gt_center">678</td>
<td class="gt_row gt_center">12,082</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Premium</td>
<td class="gt_row gt_center">1,603</td>
<td class="gt_row gt_center">2,337</td>
<td class="gt_row gt_center">2,331</td>
<td class="gt_row gt_center">2,924</td>
<td class="gt_row gt_center">2,360</td>
<td class="gt_row gt_center">1,428</td>
<td class="gt_row gt_center">808</td>
<td class="gt_row gt_center">13,791</td></tr>
 <tr><td class="gt_row gt_left" style="text-align: left; text-indent: 10px;">Ideal</td>
<td class="gt_row gt_center">2,834</td>
<td class="gt_row gt_center">3,903</td>
<td class="gt_row gt_center">3,826</td>
<td class="gt_row gt_center">4,884</td>
<td class="gt_row gt_center">3,115</td>
<td class="gt_row gt_center">2,093</td>
<td class="gt_row gt_center">896</td>
<td class="gt_row gt_center">21,551</td></tr>
 <tr><td class="gt_row gt_left" style="font-weight: bold;">Total</td>
<td class="gt_row gt_center">6,775</td>
<td class="gt_row gt_center">9,797</td>
<td class="gt_row gt_center">9,542</td>
<td class="gt_row gt_center">11,292</td>
<td class="gt_row gt_center">8,304</td>
<td class="gt_row gt_center">5,422</td>
<td class="gt_row gt_center">2,808</td>
<td class="gt_row gt_center">53,940</td></tr>
 </tbody>
 
 
</table>
</div>
]

```r
# Guardar en formato Word
diamonds %>% tbl_cross(row = color, col = cut) %>% bold_labels()  %>% 
        as_flex_table() %>% # Transforma el elemento para su compatibilidad con {flextable}
        flextable::save_as_docx(path = "./diamonds.docx") # Nombre y ruta al documento Word a crear
```

---
## paquete-**ggpubr**

El paquete "**ggpubr**" proporciona algunas funciones que facilitan el uso de **ggplot2** para crear y personalizar gráficos con una estética orientada y apropiada para su publicación. Los gráficos y objetos creados son complementarios con las funciones y elementos en el universo **tidyverse** y **ggplot2**.

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("ggpubr")

# Instalar la versión beta desde GitHub
install.packages("devtools")
devtools::install_github("kassambara/ggpubr")
```

```r
# Cargar desde la librería
library(ggpubr)
```

.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/index.html) package]

---
class: inverse center middle

## Estadística descriptiva
### Visualizando distribuciones

---
## Histograms 
Gráfica la distribución de una variable continua dividiendola en `bins` y calculando sus frecuencías. 
.left-code[
Versión {ggpubr} package

```r
gghistogram(diamonds, x = "price", 
            fill = "royalblue", 
            bins = 100)
```
Versión {ggplot2} package

```r
ggplot(diamonds, aes(x=price), 
       fill = "royalblue")+
  geom_histogram(bins = 100)
```
> El paramétro `bins` permite seleccionar el número de observaciones agrupadas por cada barra.

]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-hp-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/gghistogram.html) histogram [{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_histogram.html) histogram]

---
## Density plot
Calcula y gráfica la estimación de la densidad de los datos, es una versión suavizada del histograma. El area bajo la curva suma 1. 
.left-code[
Versión {ggpubr} package

```r
ggdensity(diamonds, x = "price")
```
Versión {ggplot2} package

```r
ggplot(data = diamonds, aes(price)) +
 geom_density() +
 theme_classic()
```
]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-dp-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggdensity.html) density plot [{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_density.html) density plot]

---
## Density plot
Calcula y gráfica la estimación de la densidad de los datos, es una versión suavizada del histograma. El area bajo la curva suma 1. 
.left-code[
Versión {ggpubr} package

```r
ggdensity(diamonds, 
          x = "price", 
          fill = "color")
```
Versión {ggplot2} package

```r
ggplot(data = diamonds, 
       aes(price, fill = color)) +
  geom_density() +
  theme_classic()
```
> Podemos separar los grupos al incorporar una tercera variable

]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-dp2-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggdensity.html) density plot [{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_density.html) density plot]

---
## Heatmaps
Heatmaps plot muestran un mapa de calor donde se observa la distribución de los datos (muy útiles para evitar overplot).
.left-code[
Versión {ggplot2} package

```r
ggplot(diamonds %>% 
 slice_sample(n=100), 
 aes(x = price, y = carat)) + 
 geom_density_2d_filled() +
 theme_classic()
```
> Para facilitar la visualización y a modo de ejemplo hemos realizado un submuestreo de solo 100 diamantes con `slice_sample(n=100)`

]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-hm-out-1.png" width="576" />
> Versión {ggplot2} package

]
.footnote[[{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_density_2d.html) density plot 2d]

---
## Heatmaps
Heatmaps plot muestran un mapa de calor donde se observa la distribución de los datos (muy útiles para evitar overplot).
.left-code[
Versión {ggplot2} package

```r
ggplot(diamonds %>% 
 slice_sample(n=100), 
 aes(x = price, y = carat)) + 
 geom_density_2d_filled(contour_var = "count") +
 theme_classic()
```
> También podemos cambiar la escala desde probabilidad a cuentas con `contour_var = "count"`.
]

.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-hm2-out-1.png" width="576" />
> Versión {ggplot2} package

]
.footnote[[{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_density_2d.html) density plot 2d]

---
## Medidas de Tendencia central
.left-code-wide[
Versión {ggplot2} package

```r
ggplot(data = diamonds, aes(price)) +
  geom_density(fill = "#00CC99") +
  geom_vline(aes(xintercept=mean(price)), color="black") +
  geom_vline(aes(xintercept=median(price)), color="red") +
  geom_vline(aes(xintercept=diamonds %>% 
                            add_count(price,sort=TRUE) %>% 
                            slice_head(n=1) %>% 
                            select(price) %>% 
                            as.numeric()), color="blue") +
  theme_classic()
```
> Podemos adicionar los valores de tendencia central Promedio (negro), Media (rojo), Moda (azul)

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-mt-out-1.png" width="432" />
> Versión {ggplot2} package

]
.footnote[[{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_abline.html) reference line]

---
## Q-Q plot
Q-Q plot (*quantile-quantile* plot) dibuja la correlación entre la muestra de interés y una distribución normal.
.left-code[
Versión {ggpubr} package

```r
ggqqplot(diamonds, x="price")
```

Versión {ggplot2} package

```r
ggplot(diamonds, aes(sample=price))+
 stat_qq() + stat_qq_line()
```
]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-qq-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggqqplot.html) qq plot [{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_qq.html) qq plot]

---
## ECDF plot
ECDF plot (*Empirical Cumulative Density Function* plot) permite observar la distribución empírica acumulativa de los datos.
.left-code[
Versión {ggpubr} package

```r
ggecdf(diamonds, x="price")
```

Versión {ggplot2} package

```r
ggplot(diamonds, aes(x=price))+
 stat_ecdf()
```
]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-ecdf-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggecdf.html) ecdf plot [{ggplot2}](https://ggplot2.tidyverse.org/reference/stat_ecdf.html) ecdf plot]

---
## ECDF plot
ECDF plot (*Empirical Cumulative Density Function* plot) permite observar la distribución empírica acumulativa de los datos.
.left-code[
Versión {ggpubr} package

```r
ggecdf(diamonds, 
       x="price", 
       color = "cut")
```

Versión {ggplot2} package

```r
ggplot(diamonds, 
 aes(x=price,color=cut))+
 stat_ecdf()
```
]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-ecdf2-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggecdf.html) ecdf plot [{ggplot2}](https://ggplot2.tidyverse.org/reference/stat_ecdf.html) ecdf plot]

---
## Boxplot 
El boxplot muestra de forma compacta la distribución de una variable continua (Mediana, 1ºQ-4ºQ) y "outliers". 
.left-code[
Versión {ggpubr} package

```r
ggboxplot(diamonds, 
          x = "cut", 
          y = "price" , 
          fill = "cut", 
          outlier.shape = NA)
```
Versión {ggplot2} package

```r
ggplot(data = diamonds, 
 aes(x=cut, 
 y=price, 
 fill=cut)) +
 geom_boxplot(outlier.shape = NA) +
 theme_classic()
```
]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-bp-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggboxplot.html) boxplot [{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_boxplot.html) boxplot]

---
## Violin plot 
Un violin plot es un diagrama de densidad que se muestra con la estructura de un boxplot. 
.left-code[
Versión {ggpubr} package

```r
ggviolin(diamonds, 
          x = "color", 
          y = "price" , 
          fill = "color")
```
Versión {ggplot2} package

```r
ggplot(data = diamonds, 
 aes(x=color, 
 y=price, 
 fill=color)) +
 geom_violin(trim = FALSE) +
 theme_classic()
```
]
.right-plot[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-vp-out-1.png" width="576" />
> Versión {ggpubr} package

]
.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/ggviolin.html) violin plot [{ggplot2}](https://ggplot2.tidyverse.org/reference/geom_violin.html) violin plot]

---
class: inverse center middle

## Estadística descriptiva
### Análisis de normalidad

---
## Pueba de normalidad Shapiro-Wilk
### Shapiro-Wilk (SW) test
.pull-left[
H0 = La variable muestra una distribución normal 
H1 = La variable no muestra una distribución normal 
> Para grupos pequeños N < 5000

### Interpretación 
p-value > alfa: No rechazar H0 (normal) 
p-value < alfa: Rechazar H0 (no normal) 
> alfa hipotético 5% (0,05)
]

.pull-right[
.left[

```r
dmnds<- diamonds %>% slice_sample(n = 5000)
shapiro.test(dmnds$price)
```
.medium[

```
## 
## 	Shapiro-Wilk normality test
## 
## data: dmnds$price
## W = 0.79564, p-value < 2.2e-16
```
]
]
.left[

```r
# 1) Crear serie entre -50 y 50, cada 5 "x"
# 2) Crear distribución normal "y"
x = seq(-50, 50, by=5) 
y = dnorm(x, mean(x), sd(x)) 
shapiro.test(dnorm(y, mean(y), sd(y)))
```
.medium[

```
## 
## 	Shapiro-Wilk normality test
## 
## data:  dnorm(y, mean(y), sd(y))
## W = 0.91786, p-value = 0.07866
```
]
]
]

---
## Pueba de normalidad Kolmogorov-Smirnov
### Kolmogorov-Smirnov (KS) test
.pull-left[
H0 = La variable muestra una distribución normal 
H1 = La variable no muestra una distribución normal 
> Menos robusto que SW, pero no limita el N

### Interpretación 
p-value > alfa: No rechazar H0 (normal) 
p-value < alfa: Rechazar H0 (no normal) 
> alfa hipotético 5% (0,05)
]

.pull-right[
.left[

```r
# Múltiples variables independientes
ks.test(diamonds$price, "pnorm", mean=mean(diamonds$price),sd=sd(diamonds$price))
```
.medium[

```
## 
## 	One-sample Kolmogorov-Smirnov test
## 
## data: diamonds$price
## D = 0.18467, p-value < 2.2e-16
## alternative hypothesis: two-sided
```
]
]
.left[

```r
x = seq(-50, 50, by=5) 
y = dnorm(x, mean(x), sd(x)) 
ks.test(y, "pnorm", mean=mean(y), sd=sd(y))
```
.medium[

```
## 
## 	One-sample Kolmogorov-Smirnov test
## 
## data:  y
## D = 0.13645, p-value = 0.8291
## alternative hypothesis: two-sided
```
]
]
]

---
## Distribución de referencia 
`stat_overlay_normal_density()`de **ggpubr** permite dibujar una distribución normal sobre la distribución de la muestra de interés.

```r
ggdensity(data=diamonds, x = "price", fill = "lightgreen") + 
  stat_overlay_normal_density(color = "darkblue", linetype = "dashed")
```

.center[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-sond-out-1.png" width="360" />
]

.footnote[[{ggpubr}](http://rpkgs.datanovia.com/ggpubr/reference/stat_overlay_normal_density.html) package - function]

---
## Transfromación de datos
En la situaciones donde no se cumple el supuesto de normalidad, se puede considerar transformar los datos. 
- Asimetria a la derecha (o positivo): Generalmente Moda < Mediana < Media 
- Asimetria a la izquierda (o negativo): Generalmente Moda > Mediana > Media

.pull-left[
- square-root para asimetría moderada: 
 - sqrt(x) Asimetria positiva, 
 - sqrt(max(x+1) - x) Asimetria negativa 
- log para asimetría mayor: 
 - log10(x) Asimetria positiva, 
 - log10(max(x+1) - x) Asimetria negativa 
- inverso para asimetría severa: 
 - 1/x Asimetria positiva, 
 - 1/(max(x+1) - x) Asimetria negativa 
]
.pull-right[

```r
# Log transform
## log10(x)
dmnds<- diamonds %>% 
 mutate(price_log = log10(price))

## log10(max(x+1) - x)
dmnds<- diamonds %>% 
mutate(price_log = log10(max(price+1)-price))
```
]

--
.bg-washed-red.b--dark-red.ba.bw2.br3.shadow-5.ph3.mt3[
La transformación no siempre será exitosa, frente a lo que se puede optar por métodos no paramétricos.
.tr[
— _Advertencia_
]]

---
class: inverse center middle

## Estadística descriptiva
### Análisis de Homocedasticidad

---
## Pueba de homocedasticidad Levene test
### Levene test
.pull-left[
Se usa para probar la homogeneidad de varianza (homocedasticidad) en k grupos de muestras, donde k puede ser mayor a dos. Está adaptado para datos distribuidos normalmente. 
H0 = Los grupos muestra homocedasticidad 
H1 = Los grupos (al menos 2) no muestra homocedasticidad

### Interpretación 
p-value > alfa: No rechazar H0 (homocedasticidad) 
p-value < alfa: Rechazar H0 (no homocedasticidad) 
> alfa hipotético 5% (0,05)
]

.pull-right[

```r
library(car)
# Una variable independiente
leveneTest(y = diamonds$price, 
           group = diamonds$cut, 
           center = "median")
```
.small[

```
## Levene's Test for Homogeneity of Variance (center = "median")
## Df F value Pr(>F) 
## group 4 123.6 < 2.2e-16 ***
## 53935 
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
```
]
]

---
## Pueba de homocedasticidad Bartlett’s test
### Bartlett’s test
.pull-left[
Se usa para probar la homogeneidad de varianza (homocedasticidad) en k grupos de muestras, donde k puede ser mayor a dos. Está adaptado para datos distribuidos normalmente. 
H0 = Los grupos muestra homocedasticidad 
H1 = Los grupos (al menos 2) no muestra homocedasticidad

### Interpretación 
p-value > alfa: No rechazar H0 (homocedasticidad) 
p-value < alfa: Rechazar H0 (no homocedasticidad) 
> alfa hipotético 5% (0,05)
]

.pull-right[

```r
# Una variable independiente
bartlett.test(price ~ cut, data = diamonds)
```
.medium[

```
## 
## 	Bartlett test of homogeneity of variances
## 
## data: price by cut
## Bartlett's K-squared = 406.7, df = 4, p-value < 2.2e-16
```
]

```r
# Múltiples variables independientes
bartlett.test(price ~ interaction(cut,color), data=diamonds)
```
.medium[

```
## 
## 	Bartlett test of homogeneity of variances
## 
## data: price by interaction(cut, color)
## Bartlett's K-squared = 1965.3, df = 34, p-value < 2.2e-16
```
]
]

---
## Pueba de homocedasticidad Fligner-Killeen test
### Fligner-Killeen test
.pull-left[
Para probar homocedasticidad en k grupos de muestras, donde k puede ser mayor a dos. Más robusto contra las desviaciones de la normalidad o cuando hay problemas relacionados con valores atípicos (outliers). 
H0 = Los grupos muestra homocedasticidad 
H1 = Los grupos (al menos 2) no muestra homocedasticidad

### Interpretación 
p-value > alfa: No rechazar H0 (homocedasticidad) 
p-value < alfa: Rechazar H0 (no homocedasticidad) 
> alfa hipotético 5% (0,05)
]

.pull-right[

```r
# Una variable independiente
fligner.test(price ~ cut, data = diamonds)
```
.small[

```
## 
## 	Fligner-Killeen test of homogeneity of variances
## 
## data: price by cut
## Fligner-Killeen:med chi-squared = 1676.4, df = 4, p-value < 2.2e-16
```
]

```r
# Múltiples variables independientes
fligner.test(price ~ interaction(cut,color), data=diamonds)
```
.small[

```
## 
## 	Fligner-Killeen test of homogeneity of variances
## 
## data: price by interaction(cut, color)
## Fligner-Killeen:med chi-squared = 3332.1, df = 34, p-value < 2.2e-16
```
]
]

---
class: inverse center middle

## Inferencia Estadística
### Análisis de correlación

---
## paquete-**Hmisc** y **corrplot**

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("Hmisc")
install.packages("corrplot")

# Instalar la versión beta desde GitHub
install.packages("devtools")
devtools::install_github("harrelfe/Hmisc")
devtools::install_github("taiyun/corrplot", build_vignettes = TRUE)
```

```r
# Cargar desde la librería
library(Hmisc)
library(corrplot)
```

.footnote[[{Hmisc}](https://hbiostat.org/R/Hmisc/) package [{corrplot}](https://github.com/taiyun/corrplot) package]

---
## Correlación
### ¿Qué es la correlación?

La correlación es una medida estadística que nos informa sobre la asociación entre dos variables. Describe cómo se comporta una variable frente a algún cambio en la variable asociada.

### Interpretación 
Si las variables aumentan o disminuyen juntas, entonces tendrán una correlación positiva. Si las variables son opuestas (una aumenta y la otra disminuye), entonces la correlación será negativa entre ellas. Si varían de forma independiente, la correlación entre ellos será cercana a cero. 
--
.bg-washed-red.b--dark-red.ba.bw2.br3.shadow-5.ph3.mt3[
1. Las correlaciones son la medida de dependencia lineal entre dos variables, solo se pueden aplicar entre pares.
1. Son sensibles al número de muestras. Si el número de muestras es pequeños (<20 ~ 30) el análisis pierde poder estadístico (óptimo ≥ 100) 
.tr[
— _Advertencia_
]]

---
### Métodos de Correlación
.left[
Pearson: Evalúa la relación lineal entre dos variables continuas.

```r
cor(diamonds$price, diamonds$carat, use ="pairwise.complete.obs", method = "pearson")
```

```
## [1] 0.9215913
```
]

.left[
Spearman: Evalúa la relación monótona. El coeficiente de correlación de Spearman se basa en los valores clasificados por ranking para cada variable en lugar de los datos brutos.

```r
dmnds <- diamonds %>% mutate(color=as.numeric(dplyr::recode(color, "D"="1", "E"="2", "F"="3", "G"="4", "H"="5", "I"="6", "J"="7")))
cor(dmnds$price, dmnds$color, use ="pairwise.complete.obs", method = "spearman")
```

```
## [1] 0.1501422
```
]

.left[
Kendall: Evalúa la correlación de rangos: la semejanza en el ordenamiento de los datos cuando se agrupan en rangos.

```r
dmnds <- dmnds %>% mutate(cut=as.numeric(dplyr::recode(cut,"Ideal"="1", "Premium"="2", "Very Good"="3", "Good"="4", "Fair"="5")))
cor(dmnds$color, dmnds$cut, use ="pairwise.complete.obs", method = "kendall")
```

```
## [1] -0.01399794
```
]

---
### ¿Cómo gráfico múltiples correlaciones?
.pull-left[

```r
# Seleccionar variables numéricas
dmnds <- diamonds %>% select_if(is.numeric) %>% as.matrix()
cor<-rcorr(dmnds,type="spearman") #"pearson" 
# Ajustar valor de p (múltiples comparaciones)
cor$P.adj<-p.adjust(cor$P, method = "fdr") 
dim(cor$P.adj) <- dim(cor$P)

# Gráfico
corrplot(cor$r, 
 method="circle", 
 type="upper", 
 col=brewer.pal(n=8, name="PuOr"), 
 tl.col="black", tl.srt=90, 
 p.mat = cor$P.adj, 
 sig.level = 0.05, 
 insig = "blank", 
 pch.cex=0.9, 
 pch.col = "white", 
 diag=TRUE, 
 title = NULL, 
 mar = c(.1, .1, .1, .1))
```
.footnote[[{corrplot}](https://github.com/taiyun/corrplot) package]
]
.pull-right[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-corrplot-out-1.png" width="504" />
]

---
class: inverse center middle

## Inferencia Estadística
### Reducción de dimensión

---
## paquete-**ggbiplot** y **broom**

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("ggbiplot")
install.packages("broom")
install.packages("ggrepel")

# Instalar la versión beta desde GitHub
install.packages("devtools")
devtools::install_github("vqv/ggbiplot")
devtools::install_github("tidymodels/broom")
devtools::install_github("slowkow/ggrepel")
```

```r
# Cargar desde la librería
library(ggbiplot)
library(broom)
library(ggrepel)
```

.footnote[[{ggbiplot}](https://github.com/vqv/ggbiplot) package [{broom}](https://broom.tidymodels.org) package [{ggrepel}](https://ggrepel.slowkow.com/index.html) package]

---
## Análisis de componentes principales (*PCA*)

El análisis de componentes principales o PCA por su nombre en inglés (**P**rincipal **C**omponent **A**nalysis) es una técnica útil durante el análisis exploratorio de datos. Permite visualizar de mejor manera la variación presente en un conjunto de datos con n > 2 variables. Es particularmente útil en el caso de conjuntos de datos "grandes", donde se tiene muchas variables para cada una de las muestras.

Los componentes principales son la extracción de características de la estructura subyacente en los datos. Con las variables "antiguas" creamos "nuevas" variables independientes, donde cada "nueva" variable independiente (PC1, ...) es una combinación de cada una de las "antiguas" variables independientes. Luego ordenamos estas nuevas variables según qué tan bien predicen nuestra variable dependiente.

--
### Interpretación
En el análisis se busca encontrar la línea recta que mejor distribuye los datos cuando se proyectan en ella. Este es conocido como el primer componente principal o PC1, y se entiende como el eje que muestra la mayor varianza en los datos. Además, tenemos "eigenvectors" que representan direcciones y "eigenvalues" que representan magnitud o importancia de las "antiguas" variables independientes. Eigenvalues más grandes se correlacionan con variables más importantes.

.footnote[[PCA in a nutshell](https://towardsdatascience.com/a-one-stop-shop-for-principal-component-analysis-5582fb7e0a9c)]

---
## Estandarización de datos *z-score*

La estandarización, también conocida como z-score, consiste en restar la media y dividir por la desviación estándar. Al hacerlo, cada valor refleja la distancia desde la media en unidades de desviación estándar.

```r
# z-score por variable
diamonds %>% 
  mutate(price_zscore = (price - mean(price))/sd(price))

# z-score por set de datos
dmnds_scaled <- diamonds %>% select_if(is.numeric) %>% 
 scale() # Función para escalar datos presente en R base
```

--
.bg-washed-red.b--dark-red.ba.bw2.br3.shadow-5.ph3.mt3[
La estandarización es necesaria en algoritmos que asumen datos centrados en cero como el análisis de componentes principales (PCA).
.tr[
— _Advertencia_
]]

---
## Análisis de componentes principales (*PCA*)

```r
dmnds.scaled.pca <- diamonds %>% select_if(is.numeric) %>% 
 scale() %>% prcomp()
summary(dmnds.scaled.pca)
```

```
## Importance of components:
##                           PC1    PC2     PC3     PC4     PC5     PC6     PC7
## Standard deviation     2.1826 1.1340 0.83115 0.41684 0.20077 0.18151 0.11135
## Proportion of Variance 0.6806 0.1837 0.09869 0.02482 0.00576 0.00471 0.00177
## Cumulative Proportion  0.6806 0.8642 0.96294 0.98776 0.99352 0.99823 1.00000
```

.left-code-wide[

```r
dmnds.scaled.pca %>%
 broom::tidy(matrix = "eigenvalues") %>%
 ggplot(aes(PC, percent)) +
 geom_col(fill = "#208C8C", alpha = 1) +
 scale_x_continuous(breaks = 1:7) +
 scale_y_continuous(labels = scales::percent_format(accuracy =1), 
 expand = expansion(mult = c(0, 0.01))) +
 coord_cartesian(xlim = c(0, 7), ylim = c(0, 0.7)) +
 labs(x="Principal Component",
 y="Percentage of Contribution") +
 theme_classic()
```
]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-pca-out-1.png" width="360" />
]

---
## Análisis de componentes principales (*PCA*)
### Gráfico basico
.left-code-wide[

```r
ggbiplot(dmnds.scaled.pca, obs.scale = 1, var.scale = 1)
```
> El biplot es una forma típica de representar un PCA. En resumén permite la representación conjunta de las muestras y los eigenvectors/eigenvalues. Sin embargo, si se trabaja con muchas muestras y/o variables, la visualización puede tornarse compleja.
]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-pca2-out-1.png" width="396" />
]

---
## Análisis de componentes principales (*PCA*)
### Gráfico **ggplot2**
.left-code-wide[

```r
dmnds.scaled.pca %>% 
 broom::augment(diamonds) %>% # Addicionar el set de datos original
 ggplot(aes(.fittedPC1, .fittedPC2, fill = color)) + 
 geom_point(size = 1, colour="white", pch=21, alpha=1) +
 guides(fill = guide_legend(title.position = "top", 
 title.hjust = .5)) +
 labs(fill="Color",
 x="PC1 (68.06% variance explained)",
 y="PC2 (18.37% variance explained)") +
 theme_classic()
```
]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-pca3-out-1.png" width="396" />
]

---
## Análisis de componentes principales (*PCA*)
### Gráfico **ggplot2**
.left-code-wide[

```r
dmnds.scaled.pca %>% 
  broom::augment(diamonds) %>% # Addicionar el set de datos original
  ggplot(aes(.fittedPC1, .fittedPC2, fill = color)) + 
  geom_point(size = 1, colour="white", pch=21, alpha=1) +
  guides(fill = guide_legend(title.position = "top", 
                             title.hjust = .5)) +
  coord_cartesian(xlim = c(-10,10), ylim = c(-10,10)) +
  labs(fill="Color",
       x="PC1 (68.06% variance explained)",
       y="PC2 (18.37% variance explained)") +
  theme_classic()
```
> Es muy importante centrar el gráfico (0,0) para poder correctamente interpretar la distribución.

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-pca4-out-1.png" width="396" />
]

---
## Análisis de componentes principales (*PCA*)
### Gráfico **ggplot2** Eigenvectors/Eigenvalues
.left-code-wide[

```r
dmnds.scaled.pca %>% 
 broom::tidy(matrix = "rotation") %>% # extract rotation matrix
 pivot_wider(names_from = "PC", names_prefix = "PC", values_from = "value") %>% # Crea una tabla larga fusionando columnas
 ggplot(aes(PC1, PC2)) +
 geom_segment(xend = 0, yend = 0, 
 arrow = arrow(angle = 20, ends = "first",
 type = "closed", 
 length = grid::unit(8,"pt"))) +
 geom_text_repel(aes(PC1, PC2, label = column),
 hjust =0.5, nudge_x = 0.05, 
 nudge_y = 0.02, color = "#208C8C") +
 coord_cartesian(xlim = c(-0.5, 0.5), 
 ylim = c(-0.8, 0.8)) +
 theme_classic()
```
]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-pca5-out-1.png" width="396" />
]

---
## Análisis de componentes principales (*PCA*)
### Eigenvectors/Eigenvalues

```r
dmnds.s.pca.eing <- dmnds.scaled.pca %>% 
 broom::tidy(matrix = "rotation") %>% # extract rotation matrix
 pivot_wider(names_from = "PC",names_prefix = "PC", values_from = "value")
```

<table>
 <thead>
 <tr>
 <th style="text-align:left;"> column </th>
 <th style="text-align:right;"> PC1 </th>
 <th style="text-align:right;"> PC2 </th>
 <th style="text-align:right;"> PC3 </th>
 <th style="text-align:right;"> PC4 </th>
 <th style="text-align:right;"> PC5 </th>
 <th style="text-align:right;"> PC6 </th>
 <th style="text-align:right;"> PC7 </th>
 </tr>
 </thead>
<tbody>
 <tr>
 <td style="text-align:left;"> carat </td>
 <td style="text-align:right;"> 0.4524455 </td>
 <td style="text-align:right;"> -0.0346960 </td>
 <td style="text-align:right;"> 0.0054948 </td>
 <td style="text-align:right;"> -0.0683594 </td>
 <td style="text-align:right;"> 0.1339995 </td>
 <td style="text-align:right;"> 0.7681511 </td>
 <td style="text-align:right;"> 0.4258803 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> depth </td>
 <td style="text-align:right;"> -0.0009161 </td>
 <td style="text-align:right;"> -0.7306797 </td>
 <td style="text-align:right;"> -0.6728293 </td>
 <td style="text-align:right;"> -0.0472480 </td>
 <td style="text-align:right;"> -0.0887383 </td>
 <td style="text-align:right;"> 0.0144503 </td>
 <td style="text-align:right;"> -0.0556003 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> table </td>
 <td style="text-align:right;"> 0.0995161 </td>
 <td style="text-align:right;"> 0.6750674 </td>
 <td style="text-align:right;"> -0.7280695 </td>
 <td style="text-align:right;"> -0.0595406 </td>
 <td style="text-align:right;"> -0.0103761 </td>
 <td style="text-align:right;"> -0.0252683 </td>
 <td style="text-align:right;"> -0.0020493 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> price </td>
 <td style="text-align:right;"> 0.4255193 </td>
 <td style="text-align:right;"> -0.0352579 </td>
 <td style="text-align:right;"> 0.1054495 </td>
 <td style="text-align:right;"> -0.8497782 </td>
 <td style="text-align:right;"> -0.0537721 </td>
 <td style="text-align:right;"> -0.2733095 </td>
 <td style="text-align:right;"> -0.0828143 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> x </td>
 <td style="text-align:right;"> 0.4532125 </td>
 <td style="text-align:right;"> 0.0035126 </td>
 <td style="text-align:right;"> 0.0395088 </td>
 <td style="text-align:right;"> 0.2429951 </td>
 <td style="text-align:right;"> 0.0889802 </td>
 <td style="text-align:right;"> 0.1984606 </td>
 <td style="text-align:right;"> -0.8286582 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> y </td>
 <td style="text-align:right;"> 0.4472649 </td>
 <td style="text-align:right;"> 0.0021579 </td>
 <td style="text-align:right;"> 0.0541888 </td>
 <td style="text-align:right;"> 0.3284606 </td>
 <td style="text-align:right;"> -0.7740579 </td>
 <td style="text-align:right;"> -0.2152666 </td>
 <td style="text-align:right;"> 0.2088571 </td>
 </tr>
 <tr>
 <td style="text-align:left;"> z </td>
 <td style="text-align:right;"> 0.4459537 </td>
 <td style="text-align:right;"> -0.0890352 </td>
 <td style="text-align:right;"> -0.0396034 </td>
 <td style="text-align:right;"> 0.3170073 </td>
 <td style="text-align:right;"> 0.6033966 </td>
 <td style="text-align:right;"> -0.4986704 </td>
 <td style="text-align:right;"> 0.2799579 </td>
 </tr>
</tbody>
</table>

---
## Análisis de componentes principales (*PCA*)
### Extraer Eigenvalues por muestra

```r
dmnds.s.pca.eing.sample <- predict(dmnds.scaled.pca, newdata = diamonds) %>% 
 as.data.frame()
```

---
class: inverse center middle

## Inferencia Estadística
### Modelos lineales

---
## Modelos lineales

###¿Qué es una regresión lineal?

Una regresión lineal es un modelo (ecuación) estadístico que se usa para analizar la relación entre una variable de respuesta (comunmente llamada **Y**) y una o más variables y sus interacciones (comunmente llamada **X** o variables explicativas). La regresión lineal asume una relación lineal entre la variable de respuesta y las variables explicativas. Esto significa que puede ajustar una línea entre las dos (o más variables).

Un modelo básico es: 
**Y** = *A* + *B* \* **X**

Donde *A* corresponde al intersecto (donde empiezas a medir cuando **X** es cero) y *B* a la pendiente (el cambio de **Y** con respecto a **X**).

La idea del modelamiento es generar un **modelo ajustado** de la mejor forma a nustras variables, encontrando el *A* y *B* que mejor se adapten a nuestro set de datos.

---
## Modelos lineales

###¿Cómo calcular una regresión lineal en **R**?

Cómo mencionamos anteriormente, **R** es una _Lingua Franca_ para estadística y ciencia de datos. En **R** base (conjunto de funciones y paquetes que se cargan por defecto al cargar tú sesión de **R**) existe la función `lm` que te permite ajustar un modelo lineal a tús variablas de interés.

La estructura básica del comando es la siguiente: 
lm([variable de respuesta o **Y**] ~ [variable explicativa o **X**], data = [nombre del set de datos])

En nuestro caso intentaremos predecir el precio de los diamantes (variable de respuesta o **Y**), en base a el tamaño o carat (variable explicativa o **X**) con el siguiente código.

```r
# Primero obtener regresión lineal y guardarla en objeto llamado "lm.pc"
lm.pc <- lm(price ~ carat, data = diamonds)
# Segundo obtener resumén del resultado
summary(lm.pc) 
```

---
## Modelos lineales

###¿Cómo interpretar una regresión lineal?

```r
lm.pc <- lm(price ~ carat, data = diamonds)
summary(lm.pc) 
```

```
## 
## Call:
## lm(formula = price ~ carat, data = diamonds)
## 
## Residuals:
## Min 1Q Median 3Q Max 
## -18585.3 -804.8 -18.9 537.4 12731.7 
## 
## Coefficients:
## Estimate Std. Error t value Pr(>|t|) 
## (Intercept) -2256.36 13.06 -172.8 <2e-16 ***
## carat 7756.43 14.07 551.4 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1549 on 53938 degrees of freedom
## Multiple R-squared: 0.8493,	Adjusted R-squared: 0.8493 
## F-statistic: 3.041e+05 on 1 and 53938 DF, p-value: < 2.2e-16
```

---
## Modelos lineales

###¿Cómo gráficar una regresión lineal básica?
.left-code-wide[

```r
ggplot(data=diamonds, aes(x=carat, y=price)) +
  geom_smooth(method = lm) +
  theme_classic()
```
> `geom_smooth(method = "...")` nos permite ajustar la función a el modelo que nos interesa gráficar, en base a funciones de modelamiento presentes en **R** como `lm` o en base a un modelo que nosotros generemos.

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-lm-out-1.png" width="396" />
]
.footnote[[geom_smooth()](https://ggplot2.tidyverse.org/reference/geom_smooth.html) reference]

---
## Modelos lineales

###¿Cómo gráficar una regresión lineal para publicación?
.left-code-wide[

```r
ggplot(data=diamonds, aes(x=carat, y=price)) +
 geom_point(alpha=1/100, color="#56B4E9") +
 geom_smooth(method = lm) +
 scale_y_continuous(labels = scales::dollar) +
 ggpubr::stat_regline_equation(label.x = 1, label.y = 31000, aes(label = paste(..eq.label.., ..adj.rr.label.., sep = "~~~~"))) +
 theme_classic()
```
> La función `geom_point(alpha=1/100, color="#56B4E9")` nos permite mostrar los datos originales, sin saturar la imagen al aplicar una transparencia de 1/100 en color "#56B4E9" (Malibu Blue). 
> La función `stat_regline_equation(...)` de {**ggpubr**} nos permite incorporar en el gráfico valores importantes como el adjust R^2 y la ecuación.

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-lm2-out-1.png" width="396" />
]
.footnote[[stat_regline_equation()](https://rpkgs.datanovia.com/ggpubr/reference/stat_regline_equation.html) reference]

---
## Modelos lineales

###¿Cómo gráficar una regresión lineal para publicación?
.left-code-wide[

```r
ggplot(data=diamonds, aes(x=carat, y=price, color=cut)) +
  geom_smooth(aes(fill = cut, color = cut), method = lm) +
  scale_y_continuous(labels = scales::dollar) +
  ggpubr::stat_regline_equation(label.y=seq(22500, 30000, 1500), aes(label =  paste(..eq.label.., ..adj.rr.label.., sep = "~~~~"))) +
  scale_color_viridis(begin =0, end = 0.6, option = "D", discrete = TRUE) +
  scale_fill_viridis(begin =0, end = 0.6, option = "D", discrete = TRUE) +
  theme_classic()
```
> Al asignar color a los grupos dentro de la variable "cut", generamos A) que los puntos se distribuyan según el color de su grupo "cut" y B) que la regresión se aplique a cada grupo de forma independiente. Sin embargo, deben prestar atención a de el factor "cut" se asigna de forma general en `ggplot(...)` y además en `geom_smooth()`

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-lm3-out-1.png" width="468" />
]
.footnote[[stat_regline_equation()](https://rpkgs.datanovia.com/ggpubr/reference/stat_regline_equation.html) reference]

---
class: inverse center middle

## Inferencia Estadística
### Análisis de varianza

---
## Análisis de varianza

Son una colección de análisis que permiten examinar si las medias de los grupos (poblaciones, muestras) difieren entre sí. Existen métodos que difieren en el número de grupos que pueden comparar. Además existen métodos parametricos (requieren que se cumplan los supuestos) y no parámetricos (no requieren que se cumplan todos los supuestos) que son más o menos robustos frente al incumplimiento de algunos de los tres supuestos en los datos.

### Supuestos
1. Independencia (los elementos de un grupo no están relacionados con los del otro grupo) 
1. Normalidad (los elementos de una muestra tienen distribución normal) 
1. Homocedasticidad (las variaciones de los grupos son iguales)

--
.bg-washed-red.b--dark-red.ba.bw2.br3.shadow-5.ph3.mt3[
La independencia de las muestras corresponde en gran medida al correcto diseño experimental, y es uno de los parámetros que debe ser cumplido para la mayor parte de los métodos no parámetricos.
.tr[
— _Advertencia_
]]

---
## paquete-**rcompanion** y **FSA**

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("rcompanion")
install.packages("FSA")
```

```r
# Cargar desde la librería
library(rcompanion)
library(FSA)
```

.footnote[[{rcompanion}](https://rdrr.io/cran/rcompanion/man/) package [{rcompanion}](https://rcompanion.org/handbook/) handbook [FSA](https://github.com/droglenc/FSA) package]

---
## T-test
T-test se utiliza para comprobar la igualdad de las medias entre dos grupos. 
> T-test requiere que los dos grupos cumplan los 3 supuestos.

.left-code-wide[

```r
dmnds <- diamonds %>% filter(cut %in% c("Fair","Ideal")) 
t.test(price ~ cut, dmnds)
```

```
## 
## 	Welch Two Sample t-test
## 
## data: price by cut
## t = 9.7484, df = 1894.8, p-value < 2.2e-16
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## 719.9065 1082.5251
## sample estimates:
## mean in group Fair mean in group Ideal 
## 4358.758 3457.542
```
]

--
.right-plot-narrow[
### Interpretación 
p-value > alfa: No rechazar H0 (igualdad) 
p-value < alfa: Rechazar H0 (no igualdad) 
> alfa hipotético 5% (0,05)
]

---
## T-test
### Gráfico
.left-code-wide[

T-test se utiliza para comprobar la igualdad de las medias entre dos grupos. 
> ¿Tienen "Fair" y "Ideal" la misma media en precio?

```r
ggplot(data = dmnds, aes(x=cut, y=price, fill=cut)) +
 geom_boxplot(outlier.shape = NA) +
 stat_compare_means(method = "t.test") +
 theme_classic()
```
]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-ttest-out-1.png" width="432" />
]

---
## F-test
F-test se utiliza para comprobar la igualdad de las varianzas entre dos grupos. 
> F-test requiere que los dos grupos cumplan los 3 supuestos. 
> ¿Tienen "Fair" y "Ideal" la misma varianza en precio?

.left-code-wide[

```r
var.test(price ~ cut, dmnds)
```

```
## 
## 	F test to compare two variances
## 
## data:  price by cut
## F = 0.87399, num df = 1609, denom df = 21550, p-value = 0.0003133
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##  0.8144715 0.9399720
## sample estimates:
## ratio of variances 
##           0.873995
```
]

--
.right-plot-narrow[
### Interpretación 
p-value > alfa: No rechazar H0 (igualdad) 
p-value < alfa: Rechazar H0 (no igualdad) 
> alfa hipotético 5% (0,05)

]

---
## chi-square test
chi-square test evalúa si existe una asociación significativa entre las categorías de dos grupos. 
> ¿Son las mediciones de los pacientes "A" y "B" independientes?

.left-code-wide[

```r
# Crearemos grupos a comparar 
df <- tibble(PacienteA = c(89,37,30), 
 PacienteB = c(20,15,50))
chisq.test(df) # Evaluar p-value
```

```
## 
## 	Pearson's Chi-squared test
## 
## data:  df
## X-squared = 40.593, df = 2, p-value = 1.533e-09
```
]

--
.right-plot-narrow[
### Interpretación 
p-value > alfa: No rechazar H0 (Independientes) 
p-value < alfa: Rechazar H0 (Dependientes) 
> alfa hipotético 5% (0,05)

]

---
## ANOVA test (One-way)
ANOVA es una técnica estadística que se utiliza para comparar las medias de más de dos grupos. 
> ANOVA requiere que los grupos cumplan los 3 supuestos.

```r
dmnds.aov <- aov(price ~ cut, data = diamonds)
summary(dmnds.aov)
```

```
## Df Sum Sq Mean Sq F value Pr(>F) 
## cut 4 1.104e+10 2.760e+09 175.7 <2e-16 ***
## Residuals 53935 8.474e+11 1.571e+07 
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
```

--
.left[
#### Post-Hoc test

```r
# Tukey's HSD (honestly significant difference) test
TukeyHSD(dmnds.aov)
```
]

---
## ANOVA test (One-way)
### Gráfico
.left-code-wide[

```r
ggplot(data = diamonds, 
       aes(x=cut, y=price, fill=cut)) +
  geom_boxplot(outlier.shape = NA) +
  stat_compare_means(method = "anova") +
  stat_compare_means(aes(label = ..p.signif..), 
                     method = "t.test", 
                     ref.group = ".all.", 
                     label.y = 18000) +
  theme_classic()
```
> La segunda comparación incluye el Post-Hoc test "t.test" para comparar entre grupos.

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-ANOVA-out-1.png" width="432" />
]
.footnote[[ggpubr](http://rpkgs.datanovia.com/ggpubr/reference/stat_compare_means.html) `stat_compare_means()` function]

---
## ANOVA test (Two-way)

```r
dmnds.aov2 <- aov(price ~ cut*color, data = diamonds)
summary(dmnds.aov2)
```

```
## Df Sum Sq Mean Sq F value Pr(>F) 
## cut 4 1.104e+10 2.760e+09 181.405 <2e-16 ***
## color 6 2.551e+10 4.251e+09 279.371 <2e-16 ***
## cut:color 24 1.653e+09 6.889e+07 4.527 1e-12 ***
## Residuals 53905 8.203e+11 1.522e+07 
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
```
--
.left[
#### Post-Hoc test

```r
# Tukey's HSD (honestly significant difference) test
TukeyHSD(dmnds.aov2)
```
]

---
## Kruskal-Wallis test

```r
kruskal.test(price ~ cut, data = diamonds)
```

```
## 
## 	Kruskal-Wallis rank sum test
## 
## data: price by cut
## Kruskal-Wallis chi-squared = 978.62, df = 4, p-value < 2.2e-16
```

--
.left[
#### Post-Hoc test

```r
# Mann–Whitney–Wilcoxon o Wilcoxon rank-sum test
pairwise.wilcox.test(diamonds$price, diamonds$cut, p.adjust.method="fdr")
```
.medium[

```
## 
## 	Pairwise comparisons using Wilcoxon rank sum test 
## 
## data: diamonds$price and diamonds$cut 
## 
## Fair Good Very Good Premium
## Good 7.9e-14 - - - 
## Very Good < 2e-16 0.023 - - 
## Premium 1.5e-05 1.6e-12 < 2e-16 - 
## Ideal < 2e-16 < 2e-16 < 2e-16 < 2e-16
## 
## P value adjustment method: fdr
```
]
]

---
## Kruskal-Wallis test
### Gráfico
.left-code-wide[

```r
ggplot(data = diamonds, aes(x=cut, y=price, fill=cut)) +
  geom_boxplot(outlier.shape = NA) +
  stat_compare_means() +
  stat_compare_means(aes(label = ..p.signif..), 
                     method = "wilcox.test", 
                     ref.group = ".all.", 
                     label.y = 18000) +
  theme_classic()
```
> La segunda comparación incluye el Post-Hoc test "wilcox.test".

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-KW-out-1.png" width="468" />
]
.footnote[[ggpubr](http://rpkgs.datanovia.com/ggpubr/reference/stat_compare_means.html) `stat_compare_means()` function]

---
## Kruskal-Wallis test
### Gráfico
.left-code-wide[

```r
my_comparisons <- list(c("Fair", "Ideal"), 
 c("Fair", "Premium"), 
 c("Fair", "Very Good"))
ggplot(data = diamonds, aes(x=cut, y=price, fill=cut)) +
 geom_boxplot(outlier.shape = NA) +
 stat_compare_means() +
 stat_compare_means(aes(label = ..p.signif..), 
 method = "wilcox.test", 
 comparisons = my_comparisons, 
 label.y = c(15000, 16000, 17000)) +
 theme_classic()
```
> También podemos adicionar comparaciones específicas entre los grupos.

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-KW2-out-1.png" width="468" />
]
.footnote[[ggpubr](http://rpkgs.datanovia.com/ggpubr/reference/stat_compare_means.html) `stat_compare_means()` function]

---
## Kruskal-Wallis test
### Gráfico
.left-code-wide[

```r
my_comparisons <- list(c("Fair", "Ideal"), 
 c("Fair", "Premium"), 
 c("Fair", "Very Good"))
ggplot(data = diamonds, aes(x=cut, y=price, fill=cut)) +
 geom_boxplot(outlier.shape = NA) +
 stat_compare_means() +
 stat_compare_means(aes(label = ..p.value..), 
 method = "wilcox.test", 
 comparisons = my_comparisons, 
 label.y = c(15000, 16000, 17000)) +
 theme_classic()
```
> También podemos cambiar desde marcas a p-values cambiando por `..p.value..`

]
.right-plot-narrow[
<img src="data:image/png;base64,#Sesion_002_files/figure-html/plot-KW3-out-1.png" width="468" />
]
.footnote[[ggpubr](http://rpkgs.datanovia.com/ggpubr/reference/stat_compare_means.html) `stat_compare_means()` function]

---
## Scheirer-Ray-Hare test
La prueba de Scheirer-Ray-Hare es una prueba no paramétrica que se utiliza para un diseño bidireccional (two-way).
.left[

```r
scheirerRayHare(price ~ cut * color, data = diamonds)
```
.medium[

```
## 
## DV:  price 
## Observations:  53940 
## D:  0.9999996 
## MS total:  242464795
```

```
##              Df     Sum Sq       H    p.value
## cut           4 2.3728e+11  978.62 0.0000e+00
## color         6 2.9932e+11 1234.50 0.0000e+00
## cut:color    24 2.1919e+10   90.40 1.2373e-09
## Residuals 53905 1.2520e+13
```
]
]

--
.left[
#### Post-Hoc test
> Luego se analizan los grupos significativos.

```r
dunnTest(price ~ cut, data = diamonds, method="bh")
dunnTest(price ~ color, data = diamonds, method="bh")
```
]

---
## paquete-**report**

```r
# Instalar desde CRAN (The Comprehensive R Archive Network)
install.packages("report")

# Instalar la versión beta desde GitHub
install.packages("remotes")
remotes::install_github("easystats/report")
```

```r
# Cargar desde la librería
library(report)
```

.footnote[[{report}](https://easystats.github.io/report/index.html) package]

---
## paquete-**report**

.bg-washed-yellow.b--gold.ba.bw2.br3.shadow-5.ph4.mt4[
“*From R to your manuscript*”

El objetivo principal **report** es reducir la brecha entre la salida de datos de R y el formato de resultados contenidos en un manuscrito. **report** Genera automáticamente informes de modelos y marcos de datos de acuerdo con las pautas de las mejores prácticas (por ejemplo, el estilo de la APA), lo que garantiza la estandarización y la calidad en los informes de resultados.

.tr[
— _report package_
]]

--
.bg-washed-red.b--dark-red.ba.bw2.br3.shadow-5.ph3.mt3[
**report** es un paquete en desarrollo, no todos los análisis tienen aún su función de reporte. Sin embargo, es probable que en un tiempo se adicionen nuevas funciones en base a contribuidores (una de las muchas ventajas del _open-source_ & _open-science_).
.tr[
— _Advertencia_
]]

.footnote[[{report}](https://easystats.github.io/report/index.html) package]

---
## Funciones útiles

`report(model)` ayuda a la interpretación de los resultados obtenidos y te permite generar un reporte del modelo generado y sus resultados ya estandarizados con un estilo apropiado para su publicación (**es una guía**). Algunas funciones y resultados actualmente soportados corresponden a: Correlations, t-tests, Regression models, ANOVAs, Mixed models, Bayesian models, Bayes factors, Structural Equation Models (SEM), y Model comparison.

### Ejemplo

```r
lm.pc.rep <- report(lm.pc)
as.report_text(lm.pc.rep, summary=TRUE)
```

We fitted a linear model to predict price with carat. The model's explanatory power is substantial (R2 = 0.85, adj. R2 = 0.85). The model's intercept is at -2256.36 (95% CI [-2281.95, -2230.77]). Within this model:

- The effect of carat is statistically significant and positive (beta = 7756.43, 95% CI [7728.86, 7784.00], t(53938) = 551.41, p < .001, Std. beta = 0.92)

.footnote[[{report}](https://easystats.github.io/report/reference/report.html) reference]

---
class: center, middle

# Muchas Gracias!!

Presentación creada en R con los paquetes:

[xaringan](https://github.com/yihui/xaringan) 
[xaringanthemer](https://github.com/gadenbuie/xaringanthemer) 
[R Markdown](https://rmarkdown.rstudio.com).