conteo de secuenciaciones por paciente

dinamyc.py

@@ -57,6 +57,17 @@ Total_pacientes = len(data_o.index)
 data = data_o.loc[data_o["incluido"]]
 Total_incluidos = len(data)
 
+secuencias = visitas_datos[0][["id", "patient_visit_id"]]
+secuencias.set_index("id", inplace=True)
+for visita in [1, 2, 3]:
+    otro = visitas_datos[visita][["id", "patient_visit_id"]]
+    otro.set_index("id", inplace=True)
+    t = secuencias.join(otro, lsuffix="l-"+str(visita), how='outer')
+    secuencias = t
+secuencias["NSecs"] = secuencias.count(axis=1)
+                          
+
+
 # grafica estática de visitas/secuenciaciones
 scat = px.scatter(data, x="NDias_glucosa", y="Comidas_Aisladas",size="AUCs_glucosa")
 vis_sec = {
@@ -69,7 +80,7 @@ vis_sec = {
                 "histnorm": "percent",                
                 #"orientation": "v",
             },
-            { "x": data["Secuenciaciones"], "type": "histogram", "name": "Secuenciaciones", "histnorm": "percent"},
+            { "x": secuencias["NSecs"], "type": "histogram", "name": "Secuenciaciones", "histnorm": "percent"},
         ],
         "layout": {
             #"title": {"text": title + " " + etapalab + str(etapaDet)},