skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Multi-step ahead predictive model for blood glucose concentrations of type-1 diabetic patients
Abstract Continuous monitoring of blood glucose (BG) levels is a key aspect of diabetes management. Patients with Type-1 diabetes (T1D) require an effective tool to monitor these levels in order to make appropriate decisions regarding insulin administration and food intake to keep BG levels in target range. Effectively and accurately predicting future BG levels at multi-time steps ahead benefits a patient with diabetes by helping them decrease the risks of extremes in BG including hypo- and hyperglycemia. In this study, we present a novel multi-component deep learning model that predicts the BG levels in a multi-step look ahead fashion. The model is evaluated both quantitatively and qualitatively on actual blood glucose data for 97 patients. For the prediction horizon (PH) of 30 mins, the average values forroot mean squared error(RMSE),mean absolute error(MAE),mean absolute percentage error(MAPE), andnormalized mean squared error(NRMSE) are$$23.22 \pm 6.39$$ 23.22 ± 6.39 mg/dL, 16.77 ± 4.87 mg/dL,$$12.84 \pm 3.68$$ 12.84 ± 3.68 and$$0.08 \pm 0.01$$ 0.08 ± 0.01 respectively. When Clarke and Parkes error grid analyses were performed comparing predicted BG with actual BG, the results showed average percentage of points in Zone A of$$80.17 \pm 9.20$$ 80.17 ± 9.20 and$$84.81 \pm 6.11,$$ 84.81 ± 6.11 , respectively. We offer this tool as a mechanism to enhance the predictive capabilities of algorithms for patients with T1D.  more » « less
Award ID(s):
1910539
PAR ID:
10360998
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
11
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H$$_2$$ 2 O and D$$_2$$ 2 O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density$$\rho (T)$$ ρ ( T ) , isothermal compressibility$$\kappa _T(T)$$ κ T ( T ) , and self-diffusion coefficientsD(T) of H$$_2$$ 2 O and D$$_2$$ 2 O are in excellent agreement with available experimental data; the isobaric heat capacity$$C_P(T)$$ C P ( T ) obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H$$_2$$ 2 O and D$$_2$$ 2 O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H$$_2$$ 2 O and D$$_2$$ 2 O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H$$_2$$ 2 O, from PIMD simulations, is located at$$P_c = 167 \pm 9$$ P c = 167 ± 9  MPa,$$T_c = 159 \pm 6$$ T c = 159 ± 6  K, and$$\rho _c = 1.02 \pm 0.01$$ ρ c = 1.02 ± 0.01  g/cm$$^3$$ 3 . Isotope substitution effects are important; the LLCP location in q-TIP4P/F D$$_2$$ 2 O is estimated to be$$P_c = 176 \pm 4$$ P c = 176 ± 4  MPa,$$T_c = 177 \pm 2$$ T c = 177 ± 2  K, and$$\rho _c = 1.13 \pm 0.01$$ ρ c = 1.13 ± 0.01  g/cm$$^3$$ 3 . Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water,$$P_c = 203 \pm 4$$ P c = 203 ± 4  MPa,$$T_c = 175 \pm 2$$ T c = 175 ± 2  K, and$$\rho _c = 1.03 \pm 0.01$$ ρ c = 1.03 ± 0.01  g/cm$$^3$$ 3 ). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of$$T_c$$ T c for D$$_2$$ 2 O and, particularly, H$$_2$$ 2 O suggest that improved water models are needed for the study of supercooled water. 
    more » « less
  2. A<sc>bstract</sc> We present a study of$$ {\Xi}_c^0\to {\Xi}^0{\pi}^0 $$ Ξ c 0 Ξ 0 π 0 ,$$ {\Xi}_c^0\to {\Xi}^0\eta $$ Ξ c 0 Ξ 0 η , and$$ {\Xi}_c^0\to {\Xi}^0{\eta}^{\prime } $$ Ξ c 0 Ξ 0 η decays using the Belle and Belle II data samples, which have integrated luminosities of 980 fb−1and 426 fb−1, respectively. We measure the following relative branching fractions$$ {\displaystyle \begin{array}{c}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0{\pi}^0\right)/\mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right)=0.48\pm 0.02\left(\textrm{stat}\right)\pm 0.03\left(\textrm{syst}\right),\\ {}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0\eta \right)/\mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right)=0.11\pm 0.01\left(\textrm{stat}\right)\pm 0.01\left(\textrm{syst}\right),\\ {}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0{\eta}^{\prime}\right)/\mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right)=0.08\pm 0.02\left(\textrm{stat}\right)\pm 0.01\left(\textrm{syst}\right)\end{array}} $$ B Ξ c 0 Ξ 0 π 0 / B Ξ c 0 Ξ π + = 0.48 ± 0.02 stat ± 0.03 syst , B Ξ c 0 Ξ 0 η / B Ξ c 0 Ξ π + = 0.11 ± 0.01 stat ± 0.01 syst , B Ξ c 0 Ξ 0 η / B Ξ c 0 Ξ π + = 0.08 ± 0.02 stat ± 0.01 syst for the first time, where the uncertainties are statistical (stat) and systematic (syst). By multiplying by the branching fraction of the normalization mode,$$ \mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right) $$ B Ξ c 0 Ξ π + , we obtain the following absolute branching fraction results$$ {\displaystyle \begin{array}{c}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0{\pi}^0\right)=\left(6.9\pm 0.3\left(\textrm{stat}\right)\pm 0.5\left(\textrm{syst}\right)\pm 1.3\left(\operatorname{norm}\right)\right)\times {10}^{-3},\\ {}\mathcal{B}\left({\Xi}_c^0\to {\Xi}^0\eta \right)=\left(1.6\pm 0.2\left(\textrm{stat}\right)\pm 0.2\left(\textrm{syst}\right)\pm 0.3\left(\operatorname{norm}\right)\right)\times {10}^{-3},\\ {}\mathcal{B}\left({\varXi}_c^0\to {\Xi}^0{\eta}^{\prime}\right)=\left(1.2\pm 0.3\left(\textrm{stat}\right)\pm 0.1\left(\textrm{syst}\right)\pm 0.2\left(\operatorname{norm}\right)\right)\times {10}^{-3},\end{array}} $$ B Ξ c 0 Ξ 0 π 0 = 6.9 ± 0.3 stat ± 0.5 syst ± 1.3 norm × 10 3 , B Ξ c 0 Ξ 0 η = 1.6 ± 0.2 stat ± 0.2 syst ± 0.3 norm × 10 3 , B Ξ c 0 Ξ 0 η = 1.2 ± 0.3 stat ± 0.1 syst ± 0.2 norm × 10 3 , where the third uncertainties are from$$ \mathcal{B}\left({\Xi}_c^0\to {\Xi}^{-}{\pi}^{+}\right) $$ B Ξ c 0 Ξ π + . The asymmetry parameter for$$ {\Xi}_c^0\to {\Xi}^0{\pi}^0 $$ Ξ c 0 Ξ 0 π 0 is measured to be$$ \alpha \left({\Xi}_c^0\to {\Xi}^0{\pi}^0\right)=-0.90\pm 0.15\left(\textrm{stat}\right)\pm 0.23\left(\textrm{syst}\right) $$ α Ξ c 0 Ξ 0 π 0 = 0.90 ± 0.15 stat ± 0.23 syst
    more » « less
  3. A<sc>bstract</sc> In this paper we explorepp→W±(ℓ±ν)γto$$ \mathcal{O}\left(1/{\Lambda}^4\right) $$ O 1 / Λ 4 in the SMEFT expansion. Calculations to this order are necessary to properly capture SMEFT contributions that grow with energy, as the interference between energy-enhanced SMEFT effects at$$ \mathcal{O}\left(1/{\Lambda}^2\right) $$ O 1 / Λ 2 and the Standard Model is suppressed. We find that there are several dimension eight operators that interfere with the Standard Model and lead to the same energy growth, ~$$ \mathcal{O}\left({E}^4/{\Lambda}^4\right) $$ O E 4 / Λ 4 , as dimension six squared. While energy-enhanced SMEFT contributions are a main focus, our calculation includes the complete set of$$ \mathcal{O}\left(1/{\Lambda}^4\right) $$ O 1 / Λ 4 SMEFT effects consistent with U(3)5flavor symmetry. Additionally, we include the decay of theW±→ ℓ±ν, making the calculation actually$$ \overline{q}{q}^{\prime}\to {\ell}^{\pm}\nu \gamma $$ q ¯ q ± νγ . As such, we are able to study the impact of non-resonant SMEFT operators, such as$$ \left({L}^{\dagger }{\overline{\sigma}}^{\mu }{\tau}^IL\right)\left({Q}^{\dagger }{\overline{\sigma}}^{\nu }{\tau}^IQ\right) $$ L σ ¯ μ τ I L Q σ ¯ ν τ I Q Bμν, which contribute to$$ \overline{q}{q}^{\prime}\to {\ell}^{\pm}\nu \gamma $$ q ¯ q ± νγ directly and not to$$ \overline{q}{q}^{\prime}\to {W}^{\pm}\gamma $$ q ¯ q W ± γ . We show several distributions to illustrate the shape differences of the different contributions. 
    more » « less
  4. Abstract The total charm-quark production cross section per unit of rapidity$$\textrm{d}\sigma ({{\textrm{c}}\overline{\textrm{c}}})/\textrm{d}y$$ d σ ( c c ¯ ) / d y , and the fragmentation fractions of charm quarks to different charm-hadron species$$f(\textrm{c}\rightarrow {\textrm{h}}_{\textrm{c}})$$ f ( c h c ) , are measured for the first time in p–Pb collisions at$$\sqrt{s_\textrm{NN}} = 5.02~\text {Te}\hspace{-1.00006pt}\textrm{V} $$ s NN = 5.02 Te V at midrapidity ($$-0.96<0.04$$ - 0.96 < y < 0.04 in the centre-of-mass frame) using data collected by ALICE at the CERN LHC. The results are obtained based on all the available measurements of prompt production of ground-state charm-hadron species:$$\textrm{D}^{0}$$ D 0 ,$$\textrm{D}^{+}$$ D + ,$$\textrm{D}_\textrm{s}^{+}$$ D s + , and$$\mathrm {J/\psi }$$ J / ψ mesons, and$$\Lambda _\textrm{c}^{+}$$ Λ c + and$$\Xi _\textrm{c}^{0}$$ Ξ c 0 baryons. The resulting cross section is$$ \textrm{d}\sigma ({{\textrm{c}}\overline{\textrm{c}}})/\textrm{d}y =219.6 \pm 6.3\;(\mathrm {stat.}) {\;}_{-11.8}^{+10.5}\;(\mathrm {syst.}) {\;}_{-2.9}^{+8.3}\;(\mathrm {extr.})\pm 5.4\;(\textrm{BR})\pm 4.6\;(\mathrm {lumi.}) \pm 19.5\;(\text {rapidity shape})+15.0\;(\Omega _\textrm{c}^{0})\;\textrm{mb} $$ d σ ( c c ¯ ) / d y = 219.6 ± 6.3 ( stat . ) - 11.8 + 10.5 ( syst . ) - 2.9 + 8.3 ( extr . ) ± 5.4 ( BR ) ± 4.6 ( lumi . ) ± 19.5 ( rapidity shape ) + 15.0 ( Ω c 0 ) mb , which is consistent with a binary scaling of pQCD calculations from pp collisions. The measured fragmentation fractions are compatible with those measured in pp collisions at$$\sqrt{s} = 5.02$$ s = 5.02 and 13 TeV, showing an increase in the relative production rates of charm baryons with respect to charm mesons in pp and p–Pb collisions compared with$$\mathrm {e^{+}e^{-}}$$ e + e - and$$\mathrm {e^{-}p}$$ e - p collisions. The$$p_\textrm{T}$$ p T -integrated nuclear modification factor of charm quarks,$$R_\textrm{pPb}({\textrm{c}}\overline{\textrm{c}})= 0.91 \pm 0.04\;\mathrm{(stat.)} ^{+0.08}_{-0.09}\;\mathrm{(syst.)} ^{+0.05}_{-0.03}\;\mathrm{(extr.)} \pm 0.03\;\mathrm{(lumi.)}$$ R pPb ( c c ¯ ) = 0.91 ± 0.04 ( stat . ) - 0.09 + 0.08 ( syst . ) - 0.03 + 0.05 ( extr . ) ± 0.03 ( lumi . ) , is found to be consistent with unity and with theoretical predictions including nuclear modifications of the parton distribution functions. 
    more » « less
  5. A<sc>bstract</sc> Using data samples of 983.0 fb−1and 427.9 fb−1accumulated with the Belle and Belle II detectors operating at the KEKB and SuperKEKB asymmetric-energye+ecolliders, singly Cabibbo-suppressed decays$$ {\Xi}_c^{+}\to p{K}_S^0 $$ Ξ c + p K S 0 ,$$ {\Xi}_c^{+}\to \Lambda {\pi}^{+} $$ Ξ c + Λ π + , and$$ {\Xi}_c^{+}\to {\Sigma}^0{\pi}^{+} $$ Ξ c + Σ 0 π + are observed for the first time. The ratios of branching fractions of$$ {\Xi}_c^{+}\to p{K}_S^0 $$ Ξ c + p K S 0 ,$$ {\Xi}_c^{+}\to \Lambda {\pi}^{+} $$ Ξ c + Λ π + , and$$ {\Xi}_c^{+}\to {\Sigma}^0{\pi}^{+} $$ Ξ c + Σ 0 π + relative to that of$$ {\Xi}_c^{+}\to {\Xi}^{-}{\pi}^{+}{\pi}^{+} $$ Ξ c + Ξ π + π + are measured to be$$ {\displaystyle \begin{array}{c}\frac{\mathcal{B}\left({\Xi}_c^{+}\to p{K}_S^0\right)}{\mathcal{B}\left({\Xi}_c^{+}\to {\Xi}^{-}{\pi}^{+}{\pi}^{+}\right)}=\left(2.47\pm 0.16\pm 0.07\right)\%,\\ {}\frac{\mathcal{B}\left({\Xi}_c^{+}\to \Lambda {\pi}^{+}\right)}{\mathcal{B}\left({\Xi}_c^{+}\to {\Xi}^{-}{\pi}^{+}{\pi}^{+}\right)}=\left(1.56\pm 0.14\pm 0.09\right)\%,\\ {}\frac{\mathcal{B}\left({\Xi}_c^{+}\to {\Sigma}^0{\pi}^{+}\right)}{\mathcal{B}\left({\Xi}_c^{+}\to {\Xi}^{-}{\pi}^{+}{\pi}^{+}\right)}=\left(4.13\pm 0.26\pm 0.22\right)\%.\end{array}} $$ B Ξ c + p K S 0 B Ξ c + Ξ π + π + = 2.47 ± 0.16 ± 0.07 % , B Ξ c + Λ π + B Ξ c + Ξ π + π + = 1.56 ± 0.14 ± 0.09 % , B Ξ c + Σ 0 π + B Ξ c + Ξ π + π + = 4.13 ± 0.26 ± 0.22 % . Multiplying these values by the branching fraction of the normalization channel,$$ \mathcal{B}\left({\Xi}_c^{+}\to {\Xi}^{-}{\pi}^{+}{\pi}^{+}\right)=\left(2.9\pm 1.3\right)\% $$ B Ξ c + Ξ π + π + = 2.9 ± 1.3 % , the absolute branching fractions are determined to be$$ {\displaystyle \begin{array}{c}\mathcal{B}\left({\Xi}_c^{+}\to p{K}_S^0\right)=\left(7.16\pm 0.46\pm 0.20\pm 3.21\right)\times {10}^{-4},\\ {}\mathcal{B}\left({\Xi}_c^{+}\to \Lambda {\pi}^{+}\right)=\left(4.52\pm 0.41\pm 0.26\pm 2.03\right)\times {10}^{-4},\\ {}\mathcal{B}\left({\Xi}_c^{+}\to {\Sigma}^0{\pi}^{+}\right)=\left(1.20\pm 0.08\pm 0.07\pm 0.54\right)\times {10}^{-3}.\end{array}} $$ B Ξ c + p K S 0 = 7.16 ± 0.46 ± 0.20 ± 3.21 × 10 4 , B Ξ c + Λ π + = 4.52 ± 0.41 ± 0.26 ± 2.03 × 10 4 , B Ξ c + Σ 0 π + = 1.20 ± 0.08 ± 0.07 ± 0.54 × 10 3 . The first and second uncertainties above are statistical and systematic, respectively, while the third ones arise from the uncertainty in$$ \mathcal{B}\left({\Xi}_c^{+}\to {\Xi}^{-}{\pi}^{+}{\pi}^{+}\right) $$ B Ξ c + Ξ π + π +
    more » « less