figuras de correlación en Macros y iAUCs

estadisticas.py

@@ -14,6 +14,40 @@ import flask
 
 server = flask.Flask(__name__)
 
+def update_correlacion_macro_iauc():
+
+    datos = comidas.sort_values(by="iAUC", ascending=True)
+    datos = datos[np.abs(stats.zscore(datos["iAUC"], nan_policy='omit'))<3]
+    datos.reset_index(inplace=True)
+    metricas = ["hc_total", "lipids_total","protein_total","fiber_total"]
+    pcts = [ metrica + "_pct" for metrica in metricas]
+    datos["total"] = datos[metricas].sum(axis=1)
+    datos = datos[datos["total"]>0]
+    datos[pcts] = datos[metricas].div(datos["total"], axis=0)
+    
+    fig = make_subplots(rows=2, cols=4, vertical_spacing=0.25)
+    lcol = 1
+    seq = 1
+    for metrica in metricas:
+        datosfil = datos[np.abs(stats.zscore(datos[metrica], nan_policy='omit')) < 3]     
+
+        fig.add_trace(
+            go.Scatter(x=datosfil["iAUC"], y=datosfil[metrica], name=metrica, mode='markers'),
+            row=1, col=lcol
+        )
+        fig.add_trace(
+            go.Scatter(x=datosfil["iAUC"], y=datosfil[metrica + "_pct"], name=metrica, mode='markers'),
+            row=2, col=lcol
+        )
+        fig["layout"]['xaxis'  + str(seq)]['title']="iAUC" 
+        fig["layout"]['yaxis'  + str(seq)]['title']= metrica
+        fig["layout"]['xaxis'  + str(seq + 4)]['title']="iAUC" 
+        fig["layout"]['yaxis'  + str(seq + 4)]['title']= "pct de " + metrica
+
+        lcol = lcol + 1
+        seq = seq + 1
+    return fig
+
 def datavis(datos):
     medidas = ["imc", "age", "glucose", "hba1c", "ct", "hdl"]    
     fig = make_subplots(rows=2, cols=3, vertical_spacing=0.25)
@@ -630,6 +664,12 @@ app.layout = html.Div(children=[
         html.Div(children= [html.H1("Distribución de %s de Macronutrientes de comidas", className = "header-title2")]),
         dcc.Graph(figure=comidas_porcentajes(comidas))
     ,]),
+    # relación  de  macronutrientes de comidas con sus iAUC
+    html.Div(children= [
+        html.Div(children= [html.H1("Relación de macros / iAUCs de las comidas ", className = "header-title2")]),
+        dcc.Graph(figure=update_correlacion_macro_iauc())
+    ,]),
+
     # # de comidas por componente dominante
     html.Div(children= [
         html.Div(children= [html.H1("Distribución de num. de comidas  por componente principal", className = "header-title2")]),