Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Chlorine resistant desalination membranes based on directly sulfonated poly(arylene ether sulfone) copolymers

The present invention provides a membrane, kit, and method of making a hydrophilic-hydrophobic random copolymer membrane. The hydrophilic-hydrophobic random copolymer membrane includes a hydrophilic-hydrophobic random copolymer. The hydrophilic-hydrophobic random copolymer includes one or more hydrophilic monomers having a sulfonated polyarylsulfone monomer and a second monomer and one or more hydrophobic monomers having a non-sulfonated third monomer and a fourth monomer. The sulfonated polyarylsulfone monomer introduces a sulfonate into the hydrophilic-hydrophobic random copolymer prior to polymerization.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Low-rank separated representation surrogates of high-dimensional stochastic functions: Application in Bayesian inference

This study introduces a non-intrusive approach in the context of low-rank separated representation to construct a surrogate of high-dimensional stochastic functions, e.g., PDEs/ODEs, in order to decrease the computational cost of Markov Chain Monte Carlo simulations in Bayesian inference. The surrogate model is constructed via a regularized alternative least-square regression with Tikhonov regularization using a roughening matrix computing the gradient of the solution, in conjunction with a perturbation-based error indicator to detect optimal model complexities. The model approximates a vector of a continuous solution at discrete values of a physical variable. The required number of random realizations to achieve a successful approximation linearly depends on the function dimensionality. The computational cost of the model construction is quadratic in the number of random inputs, which potentially tackles the curse of dimensionality in high-dimensional stochastic functions. Furthermore, this vector-valued separated representation-based model, in comparison to the available scalar-valued case, leads to a significant reduction in the cost of approximation by an order of magnitude equal to the vector size. The performance of the method is studied through its application to three numerical examples including a 41-dimensional elliptic PDE and a 21-dimensional cavity flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Adaptive Quantum Generative Training using an Unbounded Loss Function

We propose a generative quantum learning algorithm using the Adaptive Derivative-Assembled Problem Tailored ansatz (ADAPT) framework in which the loss function to be minimized is the maximal quantum Rényi divergence of order two, an unbounded function that mitigates barren plateaus which inhibit training variational circuits. We benchmark this method against other state-of-the-art adaptive algorithms by learning random two-local thermal states. We perform numerical experiments of up to 12 qubits comparing our method learning algorithms that use linear objective functions and show that Rényi-ADAPT is capable of constructing shallow quantum circuits competitive with existing methods, while the gradients remain favorable resulting from the maximal Rényi divergence loss function.

quantum algorithms, quantum machine learning, quan↗
DOE OSTI2011

Randomized Trial of Postoperative Adjuvant Therapy in Stage II and III Rectal Cancer to Define the Optimal Sequence of Chemotherapy and Radiotherapy: 10-Year Follow-Up

Purpose: To determine the optimal sequence of postoperative adjuvant chemotherapy and radiotherapy in patients with Stage II or III rectal cancer. Methods and Materials: A total of 308 patients were randomized to early (n = 155) or late (n = 153) radiotherapy (RT). Treatment included eight cycles of chemotherapy, consisting of fluorouracil 375 mg/m{sup 2}/day and leucovorin 20 mg/m{sup 2}/day, at 4-week intervals, and pelvic radiotherapy of 45 Gy in 25 fractions. Radiotherapy started on Day 1 of the first chemotherapy cycle in the early RT arm and on Day 1 of the third chemotherapy cycle in the late RT arm. Results: At a median follow-up of 121 months for surviving patients, disease-free survival (DFS) at 10 years was not statistically significantly different between the early and late RT arms (71% vs. 63%; p = 0.162). A total of 36 patients (26.7%) in the early RT arm and 49 (35.3%) in the late RT arm experienced recurrence (p = 0.151). Overall survival did not differ significantly between the two treatment groups. However, in patients who underwent abdominoperineal resection, the DFS rate at 10 years was significantly greater in the early RT arm than in the late RT arm (63% vs. 40%; p = 0.043). Conclusions: After the long-term follow-up duration, this study failed to show a statistically significant DFS advantage for early radiotherapy with concurrent chemotherapy after resection of Stage II and III rectal cancer. Our results, however, suggest that if neoadjuvant chemoradiation is not given before surgery, then early postoperative chemoradiation should be considered for patients requiring an abdominoperineal resection.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Effectiveness of supervised exercise, home-based exercise, or walk advice strategies on walking performance and muscle endurance in patients with intermittent claudication (SUNFIT trial): a randomized clinical trial

Abstract Aims Supervised exercise is a guideline-recommended treatment in intermittent claudication (IC). Hospital-based supervised exercise programmes (SEPs) are underutilized, while home-based structured exercise programmes (HSEPs) have attracted interest. The results from HSEP in IC are inconsistent and may confer no benefit over walk advice (WA) and be less effective than SEP. The aim of the study was to compare the effectiveness of best medical treatment, including Nordic pole WA alone, or WA + SEP or WA + HSEP for patients with IC. Methods and results This three-armed, multicentre randomized clinical trial enrolled patients with IC; all patients received best medical treatment including walking poles and the advice of regular Nordic pole walking (WA). For HSEP and SEP, additional exercise programmes were provided. The primarily investigated hypothesis was a non-inferiority analysis of SEP vs. HSEP regarding the 6-min walk test (6MWT) maximum distance, with a pre-defined non-inferiority margin of 50 m. Supporting outcomes included muscle endurance tests and the walking impairment questionnaire. Outcomes were assessed at baseline, 3, 6, and 12 months by a blinded evaluator. Altogether 166 patients (mean age 72 years; 59% males) were randomized. In HSEP and SEP, 24 and 26% patients, respectively, were fully exercise adherent. All three groups improved pain-free walking distance over time, but there were no significant intergroup differences. The intergroup 6MWT difference between SEP and HSEP from 0 to 12 months was –11.6 m, 95% confidence interval: –36.4 to 13.0 m (i.e. within the pre-specified non-inferiority margin). Conclusion The HSEP was non-inferior to SEP in patients with IC. There were no significant differences observed between the three groups at 1 year. Registration ClinicialTrials.gov: NCT02341716.

Sandberg, Anna (ORCID:0000000258432784)↗
DOE OSTI2018

Variational Principle for Mass Transport

A variation principle for mass transport in solids is derived that recasts transport coefficients as minima of local thermodynamic average quantities. The result is independent of diffusion mechanisms and applies to amorphous and crystalline systems. This unifies different computational approaches for diffusion and provides a framework for the creation of new approximation methods with error estimation. It gives a different physical interpretation of the Green function. Lastly, the variational principle quantifies the accuracy of competing approaches for a nontrivial diffusion problem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Parametric models for samples of random functions

A new class of parametric models, referred to as sample parametric models, is developed for random elements that match sample rather than the first two moments and/or other global properties of these elements. The models can be used to characterize, e.g., material properties at small scale in which case their samples represent microstructures of material specimens selected at random from a population. The samples of the proposed models are elements of finite-dimensional vector spaces spanned by samples, eigenfunctions of Karhunen–Loève (KL) representations, or modes of singular value decompositions (SVDs). The implementation of sample parametric models requires knowledge of the probability laws of target random elements. Numerical examples including stochastic processes and random fields are used to demonstrate the construction of sample parametric models, assess their accuracy, and illustrate how these models can be used to solve efficiently stochastic equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Stochastic generation of explicit pore structures by thresholding Gaussian random fields

We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗
DOE OSTI2015

Bayesian inference on the sphere beyond statistical isotropy

We present a general method for Bayesian inference of the underlying covariance structure of random fields on a sphere. We employ the Bipolar Spherical Harmonic (BipoSH) representation of general covariance structure on the sphere. We illustrate the efficacy of the method as a principled approach to assess violation of statistical isotropy (SI) in the sky maps of Cosmic Microwave Background (CMB) fluctuations. SI violation in observed CMB maps arise due to known physical effects such as Doppler boost and weak lensing; yet unknown theoretical possibilities like cosmic topology and subtle violations of the cosmological principle, as well as, expected observational artefacts of scanning the sky with a non-circular beam, masking, foreground residuals, anisotropic noise, etc. We explicitly demonstrate the recovery of the input SI violation signals with their full statistics in simulated CMB maps. Our formalism easily adapts to exploring parametric physical models with non-SI covariance, as we illustrate for the inference of the parameters of a Doppler boosted sky map. Our approach promises to provide a robust quantitative evaluation of the evidence for SI violation related anomalies in the CMB sky by estimating the BipoSH spectra along with their complete posterior.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dimensional Interpolation for Random Walk

In this work, we employ a simple and accurate dimensional interpolation formula for the shapes of random walks at D = 3 and D = 2 based on the analytically known solutions at both limits D = ∞ and D = 1. The results obtained for the radius of gyration of an arbitrary shaped object have about 2% error compared with accurate numerical results at D = 3 and D = 2. We also calculated the asphericity for a three-dimensional random walk using the dimensional interpolation formula. The results agree very well with the numerically simulated results. The method is general and can be used to estimate other properties of random walks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Communication: Random phase approximation renormalized many-body perturbation theory

We derive a renormalized many-body perturbation theory (MBPT) starting from the random phase approximation (RPA). This RPA-renormalized perturbation theory extends the scope of single-reference MBPT methods to small-gap systems without significantly increasing the computational cost. The leading correction to RPA, termed the approximate exchange kernel (AXK), substantially improves upon RPA atomization energies and ionization potentials without affecting other properties such as barrier heights where RPA is already accurate. Thus, AXK is more balanced than second-order screened exchange [A. Grüneis et al., J. Chem. Phys. 131, 154115 (2009)], which tends to overcorrect RPA for systems with stronger static correlation. Similarly, AXK avoids the divergence of second-order Møller-Plesset (MP2) theory for small gap systems and delivers a much more consistent performance than MP2 across the periodic table at comparable cost. RPA+AXK thus is an accurate, non-empirical, and robust tool to assess and improve semi-local density functional theory for a wide range of systems previously inaccessible to first-principles electronic structure calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Fractional charge and spin errors in self-consistent Green’s function theory

We examine fractional charge and spin errors in self-consistent Green’s function theory within a second-order approximation (GF2). For GF2, it is known that the summation of diagrams resulting from the self-consistent solution of the Dyson equation removes the divergences pathological to second-order Møller-Plesset (MP2) theory for strong correlations. In the language often used in density functional theory contexts, this means GF2 has a greatly reduced fractional spin error relative to MP2. The natural question then is what effect, if any, does the Dyson summation have on the fractional charge error in GF2? To this end, we generalize our previous implementation of GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like MP2, GF2 possesses only a very small fractional charge error, and consequently minimal many electron self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without altering the positive features that make it desirable. Furthermore, we find that GF2 has both less fractional charge and fractional spin errors than typical hybrid density functionals as well as random phase approximation with exchange.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

On numerical turbulence generation for test-particle simulations

A modified method is presented to generate artificial magnetic turbulence that is used for test-particle simulations. Such turbulent fields are obtained from the superposition of a set of wave modes with random polarizations and random directions of propagation. First, it is shown that the new method simultaneously fulfils requirements of isotropy, equal mean amplitude and variance for all field components, and vanishing divergence. Second, the number of wave modes required for a stochastic particle behavior is investigated by using a Lyapunov approach. For the special case of slab turbulence, it is shown that already for 16 wave modes the particle behavior agrees with that shown for considerably larger numbers of wave modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

A BAYESIAN APPROACH TO LOCATING THE RED GIANT BRANCH TIP MAGNITUDE. II. DISTANCES TO THE SATELLITES OF M31

In 'A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (Part I)', a new technique was introduced for obtaining distances using the tip of the red giant branch (TRGB) standard candle. Here we describe a useful complement to the technique with the potential to further reduce the uncertainty in our distance measurements by incorporating a matched-filter weighting scheme into the model likelihood calculations. In this scheme, stars are weighted according to their probability of being true object members. We then re-test our modified algorithm using random-realization artificial data to verify the validity of the generated posterior probability distributions (PPDs) and proceed to apply the algorithm to the satellite system of M31, culminating in a three-dimensional view of the system. Further to the distributions thus obtained, we apply a satellite-specific prior on the satellite distances to weight the resulting distance posterior distributions, based on the halo density profile. Thus in a single publication, using a single method, a comprehensive coverage of the distances to the companion galaxies of M31 is presented, encompassing the dwarf spheroidals Andromedas I-III, V, IX-XXVII, and XXX along with NGC 147, NGC 185, M33, and M31 itself. Of these, the distances to Andromedas XXIV-XXVII and Andromeda XXX have never before been derived using the TRGB. Object distances are determined from high-resolution tip magnitude posterior distributions generated using the Markov Chain Monte Carlo technique and associated sampling of these distributions to take into account uncertainties in foreground extinction and the absolute magnitude of the TRGB as well as photometric errors. The distance PPDs obtained for each object both with and without the aforementioned prior are made available to the reader in tabular form. The large object coverage takes advantage of the unprecedented size and photometric depth of the Pan-Andromeda Archaeological Survey. Finally, a preliminary investigation into the satellite density distribution within the halo is made using the obtained distance distributions. For simplicity, this investigation assumes a single power law for the density as a function of radius, with the slope of this power law examined for several subsets of the entire satellite sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A method to estimate the effect of deformable image registration uncertainties on daily dose mapping

Purpose: To develop a statistical sampling procedure for spatially-correlated uncertainties in deformable image registration and then use it to demonstrate their effect on daily dose mapping. Methods: Sequential daily CT studies are acquired to map anatomical variations prior to fractionated external beam radiotherapy. The CTs are deformably registered to the planning CT to obtain displacement vector fields (DVFs). The DVFs are used to accumulate the dose delivered each day onto the planning CT. Each DVF has spatially-correlated uncertainties associated with it. Principal components analysis (PCA) is applied to measured DVF error maps to produce decorrelated principal component modes of the errors. The modes are sampled independently and reconstructed to produce synthetic registration error maps. The synthetic error maps are convolved with dose mapped via deformable registration to model the resulting uncertainty in the dose mapping. The results are compared to the dose mapping uncertainty that would result from uncorrelated DVF errors that vary randomly from voxel to voxel. Results: The error sampling method is shown to produce synthetic DVF error maps that are statistically indistinguishable from the observed error maps. Spatially-correlated DVF uncertainties modeled by our procedure produce patterns of dose mapping error that are different from that due to randomly distributed uncertainties. Conclusions: Deformable image registration uncertainties have complex spatial distributions. The authors have developed and tested a method to decorrelate the spatial uncertainties and make statistical samples of highly correlated error maps. The sample error maps can be used to investigate the effect of DVF uncertainties on daily dose mapping via deformable image registration. An initial demonstration of this methodology shows that dose mapping uncertainties can be sensitive to spatial patterns in the DVF uncertainties.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

PERIODOGRAMS FOR MULTIBAND ASTRONOMICAL TIME SERIES

This paper introduces the multiband periodogram, a general extension of the well-known Lomb–Scargle approach for detecting periodic signals in time-domain data. In addition to advantages of the Lomb–Scargle method such as treatment of non-uniform sampling and heteroscedastic errors, the multiband periodogram significantly improves period finding for randomly sampled multiband light curves (e.g., Pan-STARRS, DES, and LSST). The light curves in each band are modeled as arbitrary truncated Fourier series, with the period and phase shared across all bands. The key aspect is the use of Tikhonov regularization which drives most of the variability into the so-called base model common to all bands, while fits for individual bands describe residuals relative to the base model and typically require lower-order Fourier series. This decrease in the effective model complexity is the main reason for improved performance. After a pedagogical development of the formalism of least-squares spectral analysis, which motivates the essential features of the multiband model, we use simulated light curves and randomly subsampled SDSS Stripe 82 data to demonstrate the superiority of this method compared to other methods from the literature and find that this method will be able to efficiently determine the correct period in the majority of LSST’s bright RR Lyrae stars with as little as six months of LSST data, a vast improvement over the years of data reported to be required by previous studies. A Python implementation of this method, along with code to fully reproduce the results reported here, is available on GitHub.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Prophylaxis of Radiation-Induced Nausea and Vomiting Using 5-Hydroxytryptamine-3 Serotonin Receptor Antagonists: A Systematic Review of Randomized Trials

Purpose: To systematically review the effectiveness and safety of 5-hydroxytryptamine-3 receptor antagonists (5-HT3 RAs) compared with other antiemetic medication or placebo for prophylaxis of radiation-induced nausea and vomiting. Methods and Materials: We searched the following electronic databases: MEDLINE, Embase, the Cochrane Central Register of Controlled Clinical Trials, and Web of Science. We also hand-searched reference lists of included studies. Randomized, controlled trials that compared a 5-HT3 RA with another antiemetic medication or placebo for preventing radiation-induced nausea and vomiting were included. We excluded studies recruiting patients receiving concomitant chemotherapy. When appropriate, meta-analysis was conducted using Review Manager (v5) software. Relative risks were calculated using inverse variance as the statistical method under a random-effects model. We assessed the quality of evidence by outcome using the Grading of Recommendations Assessment, Development, and Evaluation approach. Results: Eligibility screening of 47 articles resulted in 9 included in the review. The overall methodologic quality was moderate. Meta-analysis of 5-HT3 RAs vs. placebo showed significant benefit for 5-HT3 RAs (relative risk [RR] 0.70; 95% confidence interval [CI] 0.57-0.86 for emesis; RR 0.84, 95% CI 0.73-0.96 for nausea). Meta-analysis comparing 5-HT3 RAs vs. metoclopramide showed a significant benefit of the 5-HT3 RAs for emetic control (RR 0.27, 95% CI 0.15-0.47). Conclusion: 5-Hydroxytryptamine-3 RAs are superior to placebo and other antiemetics for prevention of emesis, but little benefit was identified for nausea prevention. 5-Hydroxytryptamine-3 RAs are suggested for prevention of emesis. Limited evidence was found regarding delayed emesis, adverse events, quality of life, or need for rescue medication. Future randomized, controlled trials should evaluate different 5-HT3 antiemetics and new agents with novel mechanisms of action such at the NK{sub 1} receptor antagonists to determine the most effective drug. Delayed nausea and vomiting should be a focus of future study, perhaps concentrating on the palliative cancer population.

62 RADIOLOGY AND NUCLEAR MEDICINE↗