# Cost-Utility Analysis of STN1013001, a Latanoprost Cationic Emulsion, versus Other Latanoprost Formulations (Latanoprost) in Open-Angle Glaucoma or Ocular Hypertension and Ocular Surface Disease in France

Carlo Lazzaro, Cecile van Steen, Florent Aptel, Cedric Schweitzer, Luigi Angelillo

► To cite this version:

Carlo Lazzaro, Cecile van Steen, Florent Aptel, Cedric Schweitzer, Luigi Angelillo. Cost-Utility Anal- ysis of STN1013001, a Latanoprost Cationic Emulsion, versus Other Latanoprost Formulations (La- tanoprost) in Open-Angle Glaucoma or Ocular Hypertension and Ocular Surface Disease in France. Journal of Ophthalmology, 2022, 2022, pp.3837471. 10.1155/2022/3837471. hal-03696350

HAL Id: hal-03696350

https://hal.science/hal-03696350

Submitted on 16 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research Article

# Cost-Utility Analysis of STN1013001, a Latanoprost Cationic Emulsion, versus Other Latanoprost Formulations (Latanoprost) in Open-Angle Glaucoma or Ocular Hypertension and Ocular Surface Disease in France

Carlo Lazzaro<sup>1</sup>,<sup>†</sup> Cécile van Steen<sup>1</sup>,<sup>‡</sup> Florent Aptel<sup>1</sup>,<sup>§</sup> Cédric Schweitzer<sup>1</sup>,<sup>¶</sup> and Luigi Angelillo<sup>2</sup>

Correspondence should be addressed to Carlo Lazzaro; carlo.lazzaro@tiscalinet.it

Received 4 December 2021; Accepted 5 February 2022; Published 29 April 2022

Academic Editor: Carlo Cagini

Copyright © 2022 Carlo Lazzaro et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

*Purpose*. To investigate the cost utility of STN1013001, a latanoprost cationic emulsion, versus Latanoprost in patients with open-angle glaucoma or ocular hypertension (OAG/OHT) and concomitant ocular surface disease (OSD) in France. *Methods*. An early Markov model, including 7 health states and a 1-year cycle length, was developed to estimate the cost utility of STN1013001 versus Latanoprost from the French health system perspective over a 5-year time horizon. The model was populated with pooled data (treatment adherence, quality of life, disease progression, and resource utilization) collected, via a questionnaire, from a convenience sample of 5 French glaucoma specialists. Remaining data were retrieved from published sources. Half-cycle correction and 2.5% real social discount rate were applied to costs (in €2020), life years saved (LYS), and quality-adjusted life years (QALYs). The incremental cost-utility ratio (ICUR) was contrasted against the informal willingness-to-pay (WTP) range for incremental LYS or QALY gained (€30,000–€50,000) suggested for France. One-way and probabilistic sensitivity analyses tested the robustness of the baseline ICUR. *Results*. Over a 5-year time horizon, STN1013001 resulted in an incremental 0.35 QALYs gained at an incremental cost of €7.39 compared to Latanoprost, resulting in an ICUR of €21.26. This is well below the lower limit of the unofficial WTP range proposed for France. Sensitivity analyses confirmed the robustness of the baseline results. *Conclusion*. Once on the market, STN1013001 will provide the French health system with a cost-effective treatment versus Latanoprost for OAG/OHT + OSD patients. These results should be confirmed by future economic evaluations carried out alongside empirical trials.

## 1. Introduction

Glaucoma is one of the most frequent causes of irreversible blindness, hence posing a substantial public health problem [1]. The clinical and economic burden of glaucoma are expected to become even more alarming in the future, as the number of patients suffering from glaucoma is projected to rise steeply, affecting approximately 112 million people

worldwide by 2040 compared with 63.4 million in 2013 [1–3]. Open-angle glaucoma (OAG) is the most prevalent form of glaucoma, accounting for approximately 86% of all glaucoma diagnoses worldwide [1].

In France, about 500,000 patients are diagnosed with open-angle glaucoma (OAG) or ocular hypertension (OHT). Similar numbers of OAG/OHT patients are expected to be undiagnosed [4].

<footnote n="1"><sup>1</sup>Studio di Economia Sanitaria, Milan, Italy</footnote>

<footnote n="2"><sup>2</sup>Santen GmbH, München, Germany</footnote>

<footnote n="3"><sup>3</sup>Department of Ophthalmology, University Hospital, CHU Grenoble-Alpes, Grenoble, France</footnote>

<footnote n="4"><sup>4</sup>CHU Bordeaux, Department of Ophthalmology, University of Bordeaux, 33000 Bordeaux, France</footnote>

<footnote n="5"><sup>5</sup>Inserm, Bordeaux Population Health Research Center, Team LEHA, UMR 1219, 33000 Bordeaux, France</footnote>
Glaucoma is a progressive and chronic disease and generally develops according to the following stages: ocular hypertension (OHT, stage 0), early glaucoma (stage 1), moderate glaucoma (stage 2), advanced glaucoma (stage 3), severe glaucoma (stage 4), and end-stage glaucoma/blindness (stage 5) [4-7] (Supplementary Material Table S1).

Patients generally only suffer from glaucoma symptoms once the disease has progressed to its advanced stages and the patient is experiencing irreversible vision loss [7]. Hence, early diagnosis and appropriate treatment are vital to maintain vision in OAG/OHT patients [3].

According to the most recent available data, during a 5-year time span, the annual cost per patient affected by OHT or OAG in France varies between Euros (€)7322 (OHT) and €8488 (OAG) (values expressed in €2000–€2004) [8].

Despite existing treatment options, low treatment adherence rates remain an important barrier to effective OAG/OHT treatment [9, 10].

In addition, a significant part of OAG/OHT patients (35%–65%, depending on the number of hypotensive drops administered) experience ocular surface disease (OSD) and related symptoms [11]. The effective management of OSD symptoms, including dry eye disease (DED), remains an unmet therapeutic need [12, 13].

The concomitant OSD symptoms can dramatically affect patients' activities of daily living, notably lower their adherence and persistence to OAG/OHT therapies, and reduce their health-related quality of life (HRQoL), also known as utility (that usually ranges between 0, i.e., death or HRQoL perceived worse than death and 1, i.e., perfect health, bounds included) [14-16]. Insufficient effective management of concomitant OSD can deteriorate IOP control and increase the risk of OAG/OHT progression [17]. Hence, there is a strong case for new therapeutic options that address these unmet needs and help avoid irreversible vision loss.

STN1013001 (Santen, Osaka, Japan), formerly DE-130A, is a latanoprost cationic emulsion for the treatment of OAG/OHT with concomitant OSD [18]. The cationic nano-emulsion (Novasorb® technology) possesses tear film stabilization and anti-inflammatory properties, representing a technical innovation for the effective management OAG/OHT with concomitant OSD [13].

In a 3-month Phase 2 trial, STN1013001 proved similar to Latanoprost at reducing IOP (-6.0% vs. -5.4%; *p* > 0.05) and significantly reduced both OSD-related signs and symptoms (-36.0% vs. -7.0%; *p* < 0.05, respectively) versus baseline in the per protocol study population [19].

As a head-to-head Phase III trial of STN1013001 versus Latanoprost is ongoing (NCT04133311) [20], STN1013001 is not yet available on the French market.

In this instance, early-stage health economic models help assess the economic value of drugs that are in development [21, 22].

This paper reports on an early Markov model-supported cost-utility analysis (CUA) [14-16, 23] comparing STN1013001 to other latanoprost formulations that are currently available in France (Latanoprost) in notional OAG/OHT + OSD patients, following the national health system perspective [15, 16].

The aim of this early health economic model is to anticipate evidence about the sustainability of STN1013001 for the French health system before it enters the market [21, 22].

## 2. Material and Methods

2.1. Decision Model. The Markov model includes 7 mutually exclusive health states (OAG/OHT stages 0–5 and gender and age-specific all-cause mortality) [4, 6, 23, 24]. It compares two hypothetical cohorts of patients on STN1013001 or Latanoprost (1000 notional patients each) over a 5-year time horizon, adopting a 1-year cycle length (365.25 days per cycle to take leap years into account) [25] (Figure 1).

The Markov model, which consisted of 904 parameters, was mainly populated with pooled data obtained from a convenience sample [26] of 5 French ophthalmologists with extensive experience in diagnosing and treating OAG/OHT (on average 902 patients followed up per year). The ophthalmologists practice in various settings (teaching hospital: 3; private office: 1; and private eye clinic: 1). Remaining data were obtained from the literature.

For both the hypothetical cohorts, OAG/OHT stage 0 was the starting Markov health state. At each Markov cycle, notional patients can remain in the same OAG/OHT stage, move to more severe OAG/OHT stages, or pass away (all cause-mortality) according to a transition probability matrix [14–16, 23] (Table S2). In line with the expected progression of the disease, backward transitions from more to less severe OAG/OHT stages were not allowed.

Costs, life years saved (LYS) and quality-adjusted life years (QALYs) were discounted at 2.50% real social discount rate [14–16, 27]. The real social discount rate was set at 0% and 4.5% in one-way sensitivity analysis (OWSA) [14–16, 27].

The half-cycle correction (i.e., the assumption that notional patients accrue 6-month costs, LYS and QALYs the year they die) was applied [14, 23].

The Markov model-supported CUA was created with Excel per Windows® 2010 (Microsoft, Redmond, WA, USA).

2.1.1. Data Collection. Between July and December 2020, a questionnaire and the target product profile of STN1013001 were sent out by e-mail to the 5 French glaucoma experts from high-volume care practice centres.

The glaucoma specialists were requested to complete the questionnaire based on the target product profile of STN1013001 as well as their previous experience with Latanoprost.

The following OAG/OHT stage-specific data were collected for both hypothetical cohorts of patients according to ophthalmologists' opinion [28]: expected number of patients on STN1013001 or Latanoprost per year, annual probabilities of remaining in the same OAG/OHT stage or moving to a more severe one, annual treatment adherence probabilities, HRQoL (utility values), and healthcare resource utilization associated with the diagnosis and management and follow-up of OAG/OHT and OSD.
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /><figcaption>FIGURE 1: Markov model. OAG/OHT =open-angle glaucoma or ocular hypertension.</figcaption></figure>

All questionnaires were returned complete. When necessary, follow-up teleconferences were scheduled with the ophthalmologists for clarifications or missing data management.

The annual probability of OSD (STN1013001 = 0.762; Latanoprost = 0.837; *p* > 0.05) was obtained from the Phase II trial comparing STN1013001 versus Latanoprost [19].

**2.1.2. Quality-Adjusted Life Years.** QALYs were obtained by multiplying LYS by OAG/OHT stage-specific utility and OSD-related disutility values. The latter was assumed to be the same for notional patients on STN1013001 or Latanoprost and was assumed to be equal to the disutility associated with severe DED (−0.120), in line with clinical practice [29]. Utility values were either obtained from literature (OAG/OHT stages 0 and 5) [30] and were assumed to be the same for both the hypothetical cohorts of patients. Utility values for OAG/OHT stages 1–4 were elicited from the experts, who were assumed to be good proxies given their extensive experience [28] (Table 1).

Utility value for death was set at 0 [15, 16].

**2.2. Cost of STN1013001 and Latanoprost.** In line with the French health system perspective adopted in this CUA, the following cost items were included in the Markov model: STN1013001, Latanoprost, add-on therapies in case of insufficient IOP control on STN1013001 or Latanoprost monotherapy, drugs for OSD management; tests and specialist visits for OAG/OHT diagnosis, patients' follow-up, and OSD management.

The treatment duration of one pack of STN1013001 was assumed to cover, on average, 30 days of treatment. One pack of Latanoprost was assumed to last 29.73 days, taking into account the availability of latanoprost formulations covering 30 days (30 unit dose containers) and 28 days (multidose bottles which are generally discarded 28 days after opening).

Following a similar approach, a 29-day and a 30-day mean treatment duration were estimated for add-on therapies and OSD-related drugs, respectively.

STN1013001 and Latanoprost therapies were assumed to last 365.25 days a year.

The unit cost per diem for STN1013001 was calculated on the grounds of its estimated ex-factory price provided by Santen (Table 1). The unit cost per diem for Latanoprost was based on the ex-factory price per month of treatment of all latanoprost formulations currently available on the French market in November 2020 [32], weighted for their current market share [31]. The market shares were calculated based on IQVIA MIDAS Sales data to moving annual total 3rd quarter of 2020. The same approach was adopted to calculate the unit costs per diem of OSD-related drugs and add-on therapies, such as timolol, in case of insufficient IOP control on STN1013001 or Latanoprost monotherapy.

Since all drugs are self-administered by the patient at home, drug administrations were not costed.

Tests and specialist visits for OAG/OHT diagnosis, patients' follow-up, and OSD management were valued at the current outpatient tariffs [33], which were assumed to be acceptable proxies for the real costs borne by the healthcare facilities to provide those healthcare services [35].

Cost were expressed in €2020 values.

**2.3. Cost-Utility Analysis.** In CUA, the difference in costs (incremental cost—ΔC) of STN1013001 and Latanoprost is divided by the difference in QALYs (incremental QALYs—ΔQALYs) [15, 16]. The CUA results are summarized in the incremental cost-utility ratio (ICUR) [15, 16].

The ICUR informs the third-party payer about the sustainability for the healthcare budget of the cost per QALY of a given healthcare programme versus comparator.

This study was conducted in line with the tenets of the Declaration of Helsinki [36].

Since no patient was enrolled in the present CUA [14, 34], no Ethics Committee approval of the study protocol (questionnaire included) was needed according to the current French legislation [37].
<table><caption>TABLE 1: Base case analysis—methods—unit cost for healthcare resources, utility, and disutility values (95% CI)<sup>a,b</sup> (costs in €2020).</caption><thead><tr><th>Model main items</th><th>Point estimate (95% CI)</th><th>Source</th></tr></thead><tbody><tr><td><strong>OAG/OHT stages 0-5</strong></td><td></td><td></td></tr><tr><td><strong>Drugs</strong></td><td></td><td></td></tr><tr><td><strong>OAG/OHT medications<sup>c,d</sup></strong></td><td></td><td></td></tr><tr><td>STN1013001</td><td>€0.30</td><td>Santen GmbH2020</td></tr><tr><td>Latanoprost</td><td>€0.24</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td><strong>Add-on therapies<sup>d,e</sup></strong></td><td></td><td></td></tr><tr><td>Acetazolamide</td><td>€0.16</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Brimonidine</td><td>€0.12</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Brinzolamide</td><td>€0.12</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Brinzolamide + brimonidine</td><td>€0.28</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Dorzolamide</td><td>€0.14</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Timolol</td><td>€0.18</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Timolol + dorzolamide</td><td>€0.20</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td><strong>OSD therapies<sup>e</sup></strong></td><td></td><td></td></tr><tr><td>Cyclosporin</td><td>€3.67</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td>Preservative-free lubricant</td><td>€0.27</td><td>IQVIA [31], Ministère des Affaires Sociales et de la Santé [32]</td></tr><tr><td><strong>Healthcare procedures<sup>f</sup></strong></td><td></td><td></td></tr><tr><td>Breakup time test</td><td>€19.20<br/>(€15.44; €22.96)</td><td>Sécurité Sociale [33] (elaborated on: Tarif Secteur 1; code: 02.01.06.03 BBQP001)</td></tr><tr><td>Diurnal curve of intraocular pressure measurement</td><td>€41.66 (€33.49; €49.83)</td><td>Sécurité Sociale [33] (elaborated on: Tarif Secteur 1; code: 02.01.06.01 BHQP001)</td></tr><tr><td>Fluorescein test</td><td>€19.20<br/>(€15.44; €22.96)</td><td>Sécurité Sociale [33] (elaborated on: Tarif Secteur 1; code: 02.01.06.03 BBQP001)</td></tr><tr><td>Gonioscopy</td><td>€17.28 (€13.89; €20.67)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.06.01 BHQP002)</td></tr><tr><td>Lissamine test</td><td>€19.20<br/>(€15.44; €22.96)</td><td>Sécurité Sociale [33] (elaborated on: Tarif Secteur 1; code: 02.01.06.03 BBQP001)</td></tr><tr><td>Optical coherence tomography retinal nerve fiber layer</td><td>€56.54 (€45.46; €67.62)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.05 BZQK001)</td></tr><tr><td>Retinal nerve fiber thickness assessment</td><td>€26.78 (€21.53; €32.03)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.04 BGQP009)</td></tr><tr><td>Schirmer test</td><td>€19.20<br/>(€15.44; €22.96)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.06.03 BBQP001)</td></tr><tr><td>Slit lamp examination</td><td>€28.29 (€24.92; €37.06)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.06.03 BGQP002)</td></tr><tr><td>Tonometry</td><td>€41.66 (€11.21; € 16.67)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.06.01 BHQP001)</td></tr><tr><td>Visual field test</td><td>€39.43<br/>(€31.70; € 47.16)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.06.02 BLQP004)</td></tr><tr><td><strong>Specialist visits</strong></td><td></td><td></td></tr><tr><td>Ophthalmologist</td><td>€25.32 (€20.36; €30.28)</td><td>Sécurité Sociale [33] (Tarif Secteur 1; code: 02.01.06.02 BLQP010)</td></tr><tr><td><strong>Utility and disutility values<sup>g</sup></strong></td><td>STN1013001<sup>h</sup><br/>0.900 (0.885; 0.915)</td><td>Latanoprost<sup>l</sup><br/>0.900 (0.884; 0.915)</td></tr><tr><td>OAG/OHT stage 0</td><td></td><td>van Gestel et al. [29]</td></tr></tbody></table>
<table><caption>TABLE 1: Continued.</caption><thead><tr><th>Model main items</th><th>Point estimate (95% CI)</th><th>Source</th></tr></thead><tbody><tr><td>OAG/OHT stage 1</td><td>0.897 (0.880;<br>0.913)</td><td>0.890<br>(0.872; 0.908)</td><td>Experts' opinion</td></tr><tr><td>OAG/OHT stage 2</td><td>0.879 (0.860;<br>0.896)</td><td>0.866 (0.845;<br>0.885)</td><td>Experts' opinion</td></tr><tr><td>OAG/OHT stage 3</td><td>0.862 (0.840;<br>0.883)</td><td>0.849 (0.825;<br>0.871)</td><td>Experts' opinion</td></tr><tr><td>OAG/OHT stage 4</td><td>0.825 (0.795;<br>0.853)</td><td>0.808 (0.776;<br>0.838)</td><td>Experts' opinion</td></tr><tr><td>OAG/OHT stage 5</td><td>0.790 (0.750;<br>0.828)</td><td>0.790 (0.748;<br>0.829)</td><td>van Gestel et al. [29],<br>Canadian Agency for drugs and Technology in<br>Health [30]</td></tr><tr><td>OSD-related disutility</td><td>-0.120<br>(-0.231; -0.045)</td><td>-0.120<br>(-0.231; -0.045)</td><td></td></tr><tr><td>Death</td><td>0.000 (-)</td><td>0.000 (-)</td><td>Drummond et al. [15], Neumann et al. [16]</td></tr></tbody></table>

**2.4. Statistical Analysis.** The mean number (standard deviation—SD) of notional patients in each Markov health state was calculated.

Most of the parameters included in the Markov model were assigned a statistical distribution [14, 34].

The beta distribution was fitted to dichotomous events (i.e., OAG/OHT patient gender) and OAG/OHT stage-specific utility values, whereas the Dirichlet distribution was assigned to polytomous events (i.e., transition probabilities from less severe to more severe stages of OAG/OHT).

The gamma distribution was fitted to the volume of healthcare resources other than drugs as well as to OSD-related disutility value.

Finally, the normal distribution was fitted to unit cost of healthcare resources (if different from drugs).

The standard error of the mean was computed from the data collected or by imposing a coefficient of variation [38] on the sample estimates [14, 34].

The 95% confidence interval (95% CI) was calculated via the percentile method [34] for parameters that were assigned a theoretical probability distribution as well as for cost, LYS, QALYs, and a set of probabilities (treatment adherence, tests and specialist visits related to the diagnosis and management and follow-up of OAG/OHT and OSD, as well as OSD-related medications).

For parameters that were not given a statistical distribution, a range accompanied the mean.

**2.5. Sensitivity Analyses.** OWSA and probabilistic sensitivity analyses (PSA) were run to check the robustness of the baseline ICUR [14–16, 34].

**2.5.1. One-Way Sensitivity Analysis.** OWSA was performed on one parameter at a time by replacing its sample estimate with the lower and upper limits of its 95% CI or range [15, 16]. OWSA results were depicted on a Tornado diagram [15].

**2.5.2. Probabilistic Sensitivity Analysis.** PSA [14–16, 34, 39] investigated the uncertainty of the baseline ICUR via a 10,000-iteration Monte Carlo simulation [14–16, 34].

During each Monte Carlo iteration, a random value for each probability parameter was sampled [14, 34].

For methodological reasons, the reference literature on health economic modelling [34] recommends not to include parameters that are set by clinical guidelines and national regulatory agencies in the PSA, as these parameters are not subject to uncertainty. Consequently, the following parameters were not included in the PSA: drug posology and cost; real social discount rate for costs, LYS and QALYs (in line with the French guidelines on the economic evaluation of healthcare programmes) [27].

The conjoint density of ΔC and ΔQALYs provided by the MC simulation were reported on the cost-effectiveness plane (CEP) (SText definition S1) [40].

Once expressed in the Net Monetary Benefit (NMB) metric (SText definition S2) [14–16, 34, 39, 41, 42], PSA results were reported graphically via the nonparametric cost-effectiveness acceptability curve (CEAC) (SText definition S3) and cost-effectiveness acceptability frontier (CEAF) (SText definition S4). Given a set of threshold values decided by the third-party payer, the CEAC provides a graphical representation of the probability that a given healthcare programme is cost-effective, whereas the CEAF indicates

<footnote>a95% CI was calculated assuming a normal probability distribution [14,34]. b95% CI was not calculated for the unit cost of drugs since they are exogenous variables [34]. cMedications refer to STN1013001 and Latanoprost only. dCost per diem calculated on ex-factory price. eAdd-on therapies prescribed in addition to STN1013001 or Latanoprost due to poor IOP control. fBeing not funded by French health system, LipiView test was not included in cost calculation. g95% CI for utility and disutility values was calculated assuming a beta and a gamma probability distribution, respectively [14]. hNumber of observations per OAG/OHT stage (female %): 0 = 1560 (50.00%); 1 = 1280 (50.00%); 2 = 1280 (50.00%); 3 = 1000 (50.00%); 4 = 650 (50.00%); and 5 = 415 (50.00%). iNumber of observations per OAG/OHT stage: 0 = 1460 (50.00%); 1 = 1160 (57.43%); 2 = 1150 (50.00%); 3 = 930 (50.00%); 4 = 61 (50.00%); and 5 = 390 (50.00%). CI = confidence interval; IOP = intraocular pressure; OAG/OHT = open-angle glaucoma/ocular hypertension; and OSD = ocular surface disease. Bold fonts were used to define different sets of items that were used to populate the Markov model.</footnote>
which health technology is optimal given its expected highest NMB [14–16, 34, 43, 44].

## 3. Results

3.1. *Markov Model.* Both the hypothetical cohorts of patients were assumed to enter the Markov model in OAG/OHT—stage 0 at the age of 47.31 years (range: 45.00; 55.00 years) (Table S3).

During the 5-year time horizon, approximately 65% of notional patients on STN1013001 or Latanoprost remained in OAG/OHT stage 0 (STN1013001: mean: 672; SD: 215; Latanoprost: mean: 649; SD: 228) (Figures 2(a) and 2(b); Table S4). A small fraction of both the hypothetical cohorts of patients eventually moved to OAG/OHT stage 5 (STN1013001: mean: 5; SD: 7; Latanoprost: mean: 6; SD: 8) All-cause mortality was similar in both hypothetical cohorts (STN1013001: mean: 114 deaths; SD: 81; Latanoprost: mean: 115 deaths; SD: 82).

3.2. *Adherence to STN1013001 and Latanoprost.* Notional patients on STN1013001 reported higher treatment adherence probabilities compared to their Latanoprost counterparts across OAG/OHT stages 0–3 (Table S5). The probability of being treatment adherent is 0.080 (95% CI: 0.048; 0.113) higher in OAG/OHT stage 0 in year 2 and 0.088 (95% CI: 0.055; 0.122) higher in year 4 in STN1013001 notional patients compared to Latanoprost notional patients. In OAG/OHT stage 1 notional patients on STN1013001 report a 0.066 (95% CI: 0.027; 0.105) and 0.098 (95% CI: 0.060; 0.135) higher treatment adherence probability versus Latanoprost in years 5 and 4, respectively. A similar result was observed in OAG/OHT stage 3: adherence probabilities in patients in the STN1013001 hypothetical cohort were 0.047 (95% CI: 0.010; 0.083) and 0.081 (95% CI: 0.040; 0.123) higher for years 1 and 5, respectively, versus their Latanoprost counterparts.

For OAG/OG stages 4 and 5, the adherence probabilities were similar for both medications.

### 3.3. Base Case Analysis

#### 3.3.1. Healthcare Resources Consumption

*Diagnosis.* Healthcare resource utilization associated with the diagnosis of OAG/OHT was similar in both the STN1013001 and Latanoprost hypothetical cohorts (Table S6).

Slit lamp examination and tonometry were the most frequently performed tests for the diagnosis of OAG/OHT (STN1013001: 100.00%; mean: 1.18; 95% CI: 0.83; 1.59; Latanoprost: 100.00%; mean: 1.19; 95% CI: 0.97; 1.43 for both tests).

All notional patients undergo at least one ophthalmologist visit (STN1013001: mean: 1.24; 95% CI: 1.02; 1.49; Latanoprost: mean: 1.26; 95% CI: 1.05; 1.49).

3.3.2. *Add-On Therapies and Follow-Up.* According to experts' opinion, STN1013001 and Latanoprost notional patients have the same annual probability of requiring ≥1 add-

on therapies due to insufficient control of IOP on STN1013001 or Latanoprost monotherapy. This probability ranges from 40% (OAG/OHT stage 0) to 95% (OAG/OHT 5) (Table S7). Timolol (35% in OAG/OHT stage 0; 50% in OAG/OHT stages 1 and 3; 60% in OAG/OHT stages 2 and 4) and dorzolamide (15% in OAG/OHT stages 3 and 4; 20% in OAG/OHT stages 1 and 2; 45% in OAG/OHT stage 0) were the most frequently prescribed add-on therapies.

Healthcare resource utilization associated with the management and follow-up of OAG/OHT patients was similar for both STN1013001 and Latanoprost patients (Table S8).

Based on expert opinion, all patients in both hypothetical cohorts (stages 0–5) underwent at least 2 follow-up visits with an ophthalmologist. During the follow-up visit, a slit lamp examination and tonometry were performed in all patients.

3.3.3. *Management of OSD.* Across all disease stages, notional patients in the STN1013001 cohort were less likely to require medication for the management of OSD versus Latanoprost notional patients (Table S9). A statistically significantly lower proportion of notional patients in the STN1013001 cohort required preservative-free lubricants compared to the Latanoprost cohort, across all disease stages. Differences in prescription probabilities for preservative-free lubricants ranged from -5.88% (95% CI: -9.13%; -0.03%) to -11.34% (95% CI: -17.68%; -0.05%) in OAG/OHT stages 0 and 5, respectively. Additionally, ciclosporin was less likely to be required for OSD management in notional STN1013001 patients versus notional Latanoprost patients. Differences ranged from -4.31% (95% CI: -7.46%; -0.01%) to -5.39% (95% CI: -10.46%; -0.003%) in OAG/OHT stages 4 and 5, respectively.

Across all the disease stages, notional patients on STN1013001 were projected to have a lower probability of receiving a Schirmer test versus those on Latanoprost, ranging from -9.57% (95% CI: -15.93%; -3.11%) in OAG/OHT stage 5 to -7.29% (95% CI: -10.35%; -4.23%) in OAG/OHT stage 2, respectively (Table S10).

All notional patients across all disease stages received at least one slit lamp examination and two ophthalmologist visits for the management of OSD.

3.3.4. *Cost of STN1013001 and Latanoprost.* After 5 years, the mean costs per notional patient on STN1013001 or Latanoprost were similar (€2347.74 vs. €2340.34; ΔC: €7.39; 95% CI: -€159.39; €166.60) (Table 2).

The cost drivers were the healthcare resources consumed due to add-on therapies and follow-up (STN101300: 56.97%; Latanoprost: 57.23%) and OSD management (STN101300: 20.33%; Latanoprost: 24.20%).

Higher cost for OAG/OHT monotherapy totaled by notional patients on STN1013001 versus notional Latanoprost patients (€341.28 vs. €241.81; difference: €99.47; 95% CI: €95.04; €104.53) is partially offset by lower costs in terms of OSD management (€477.27 vs. €566.38; difference: -€89.11; 95% CI: -€160.16; -€19.73).
<figure><img src="image_4.png" /><figcaption>FIGURE 2: (a) Base case analysis-results-Markov trace-hypothetical cohort of patients on STN1013001. (b). Base case analysis-results-Markov trace-hypothetical cohort of patients on Latanoprost. OAG/OHT=open-angle glaucoma or ocular hypertension.</figcaption></figure>
<table><caption>TABLE 2: Base case analysis—results—costs per patient and cost-utility analysis (€2020).</caption><thead><tr><th>Items</th><th>STN1013001 (%)</th><th>Latanoprost (%)</th><th>Difference (%) [95% CI]<sup>a,b</sup></th></tr></thead><tbody><tr><td>French health system viewpoint</td><td></td><td></td><td></td></tr><tr><td>Cost</td><td></td><td></td><td></td></tr><tr><td>Diagnosis</td><td>€191.60 (8.16)</td><td>€192.75 (8.24)</td><td>-€1.15 (14.63) [-€32.62; €30.80]</td></tr><tr><td>Medications</td><td>€341.28 (14.54)</td><td>€241.81 (10.33)</td><td>€99.47 (-1070.48) [€95.04; €104.53]</td></tr><tr><td>Add-on therapies and follow-up</td><td>€1337.59 (56.97)</td><td>€1339.40 (57.23)</td><td>-€1.81 (22.14) [-€146.76; €143.26]</td></tr><tr><td>OSD management</td><td>€477.27 (20.33)</td><td>€566.38 (24.20)</td><td>-€89.11 (1133.72) [-€160.16; -€19.73]</td></tr><tr><td><strong>Overall</strong></td><td><strong>€2347.74 (100.00)</strong></td><td><strong>€2340.34 (100.00)</strong></td><td><strong>€7.39 (100) [-€159.39; €166.60]</strong></td></tr><tr><td>LYS and QALYs</td><td></td><td></td><td></td></tr><tr><td>LYS</td><td>4.228</td><td>4.223</td><td>0.005 [-0.013; 0.023]</td></tr><tr><td>QALYs</td><td>2.539</td><td>2.191</td><td>0.348 [0.285; 0.408]</td></tr><tr><td>Cost-utility analysis</td><td></td><td></td><td></td></tr><tr><td>Incremental costs (ΔC)</td><td>€7.39</td><td></td><td></td></tr><tr><td>Incremental QALYs (ΔQALYs)</td><td>0.348</td><td></td><td></td></tr><tr><td>ICUR (ΔC/ΔQALYs)</td><td>€21.26</td><td></td><td></td></tr></tbody></table>

Additionally, the mean 5-year cost for notional patients on STN1013001 is lower versus Latanoprost in OAG/OHT stages 4 (STN1013001: €2712.87; Latanoprost: €2912.55; difference: -€199.68; 95% CI: -€361.20; -€38.71) and 5 (STN1013001: €2895.87; Latanoprost: €3144.02; difference: -€248.16; 95% CI: -€471.90; -€30.39) (Figure S1).

### 3.3.5. Cost-Utility Analysis.
The ICUR of €21.26 per incremental QALY gained, resulting from an incremental QALY gain of 0.348 at an incremental cost of €7.39 for STN1013001, is well below the lower limit of the unofficial WTP range (€30,000–€50,000) proposed for France [45] (Table 2).

Notional patients on STN1013001 were estimated to gain statistically significantly more QALYs versus Latanoprost over the 5-year time horizon (2.539 vs. 2.191; ΔQALYs: 0.348; 95% CI: 0.285; 0.408) (Table 2).

Across all disease stages, higher QALYs were observed in the STN1013001 hypothetical cohort vs. the Latanoprost hypothetical cohort. This was most apparent in OAG/OHT stage 1 (ΔQALYs: 0.803; 95% CI: 0.716; 0.892) (Figure S2).

Over the 5-year time horizon, both hypothetical cohorts of patients totaled similar average LYS (STN1013001: 4.228; Latanoprost: 4.223; and ΔLYS: 0.005; 95% CI: -0.013; 0.023).

## 3.4. Sensitivity Analyses

### 3.4.1. One-Way Sensitivity Analysis.
OWSA results confirm the robustness of the base case findings. For all the hypotheses tested, STN1013001 remains highly cost-effective or strongly dominant (i.e., being less costly and producing more QALYs) versus Latanoprost, as the highest ICUR does not exceed €235 per QALY gained (Figure 3).

OWSA proves the baseline ICUR to be most sensitive to changes in the volume of follow-up tonometry (-786.85% to +889.91% vs. baseline ICUR) and visual field test (from -691.95% to +785.46% vs. baseline ICUR) performed in OAG/OHT stage 0 STN1013001 notional patients. Additionally, OWSA proves the baseline ICUR to be robust to

variations in the social discount rate for costs, LYS and QALYs.

### 3.4.2. Probabilistic Sensitivity Analysis.
Probabilistic sensitivity analysis shows that STN1013001 has a 49.12% probability of being strongly dominant (i.e., less costly and more effective) versus Latanoprost (Figure 4). The probability that STN1013001 is cost-effective versus the lower bound of the last available unofficial acceptability range for QALY gained (€30,000–€50,000) suggested for the French health system equals 100.00% [45] (Figure 5).

The CEAF highlights that the probability for STN1013001 to be the optimal alternative (i.e., the healthcare programme with the highest average NMB) for OAG/OHT + OSD patients ranges from 79.01% at a willingness to pay (WTP) for ΔQALY gained = €200.20 to 100.00% from a WTP of €1000 onward (Figure 6).

# 4. Discussion

The present study aimed to obtain a first cost-utility outline of STN1013001 versus Latanoprost for the French setting in OAG/OHT + OSD patients via a Markov model-supported CUA.

Its results demonstrated that, across a 5-year time horizon, STN1013001 is a highly cost-effective therapeutic alternative for this patient population versus Latanoprost from the French health system perspective. This was mainly driven by STN1013001's projected favorable safety profile versus Latanoprost in terms of OSD and consequent improvements in patients' HRQoL and treatment adherence. Despite currently available therapeutic options for the treatment of OAG/OHT, treatment adherence rates remain low and scant evidence is available on the effectiveness of interventions that can improve it [46]. Improving low treatment adherence represents a substantial therapeutic challenge in this patient population. Additionally, the effective management of concomitant OSD is a therapeutic

<footnote><sup>a</sup>(STN1013001-Latanoprost). <sup>b</sup>95% CI was calculated via the percentile method [34]. CI = confidence interval; ICUR = incremental cost-utility ratio; LYS = life year saved; OAG/OHT = open-angle glaucoma/ocular hypertension; OSD = ocular surface disease; and QALYs = quality-adjusted life years. Bold fonts were used to highlight some results.</footnote>
<figure><img src="image_5.png" /><figcaption>FIGURE 3: One-way sensitivity analysis—results concerning the first 10 parameters of the Markov model that causes the widest variations in base case ICUR STN1013001 (€21.26; NE sector of the cost-effectiveness plane) (€2020)ª. ªY- and X-axes intersect at the baseline ICUR. ICUR = incremental cost-utility ratio; LL 95% CI = lower limit 95% confidence interval; NE = north-east; OAG/OHT = open-angle glaucoma or ocular hypertension; and UL 95% CI = upper limit 95% confidence interval.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>FIGURE 4: Probabilistic sensitivity analysis-I-cost-effectiveness plane (10,000 out of 10,000 Monte Carlo iterations reported) (€2020). Base case ICUR STN1013001: €21.26 (NE sector of the cost-effectiveness plane). Number of Monte Carlo iterations (%) for each sector of the cost-effectiveness plane: NE = 5088 (50.88%); SE = 4912 (49.12%). ΔC = incremental cost; ΔQALYs = incremental quality-adjusted life years; ICUR = incremental cost-utility ratio; NE = north-east; and SE = south-east.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>FIGURE 5: Probabilistic sensitivity analysis-II-cost-effectiveness acceptability curve (1000 out of 1000 threshold values reported) (€2020). Base case ICUR STN1013001: €21.26 (NE sector of the cost-effectiveness plane). CEAC = cost-effectiveness acceptability curve; ICUR = incremental cost-utility ratio; and NE = north-east.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>FIGURE 6: Probabilistic sensitivity analysis-III-cost-effectiveness acceptability frontier (1000 out of 1000 threshold values reported) (€2020). Base case ICUR STN1013001: €21.26 (NE sector of the cost-effectiveness plane). CEAF: STN1013001 is the optimal strategy from a threshold value of €200.20 onwards. CEAF = cost-effectiveness acceptability frontier; ICUR = incremental cost-utility ratio; and NE = north-east.</figcaption></figure>

goal that is yet to be met in OAG/OHT patients. This research has shown that STN1013001 can potentially play an important role in addressing both these issues.

Patients in the STN1013001 hypothetical cohort were less likely to develop OSD and, consequently, incurred less costs associated with OSD management (e.g., specialist visits and prescription of preservative-free lubricants) versus Latanoprost notional patients. In line with previous research [7], effectively managing OSD in OAG/OHT patients was expected to significantly reduce the well-known negative impact of OSD on patients' treatment adherence and HRQoL. This was confirmed in the present study, in which higher treatment adherence probabilities were estimated in the STN1013001 hypothetical cohort versus Latanoprost based on experts' opinion. The importance of this is underlined by the fact that high treatment adherence is a reliable predictor of effectiveness

and cost-effectiveness of OAG/OHT drugs [47]. The improved treatment adherence, as well as the lower likelihood of developing OSD estimated in the hypothetical cohort of patients on STN1013001, resulted in substantial improvements in HRQoL versus notional patients on Latanoprost.

Another factor that contributed to the higher number of QALYs totaled in the STN1013001 hypothetical cohort was the higher number of patients in less advanced disease stages (stages 0 and 1) compared to the Latanoprost hypothetical cohort across the 5-year time horizon. In accordance with previous research [29], patients in less severe disease stages were assigned a higher utility than notional patients who were partitioned among more severe disease stages, hence contributing to the higher utility values observed in the STN1013001 hypothetical cohort.
Sensitivity analyses confirmed the base case results.

Regardless of the WTP-threshold, the CEP highlights that STN1013001 always produces higher QALYs than Latanoprost.

CEAC and CEAF show the probability for STN1013001 to be cost-effective or optimal at the entire informal acceptability range of €30,000–€50,000 per QALY gained [45] to be 100.00%.

Consequently, the probability of healthcare resources misallocation for the French health system due to STN1013001 funding is zero, even at very low threshold values per QALY gained.

The present research has the following limitations.

First, as STN1013001 is currently not yet available on the French market, an empirical CUA versus Latanoprost was unfeasible. In this scenario, early health economic models are a valuable methodological option for supporting OAG/ OHT healthcare decision-making [21, 22, 48–50].

In the literature, two Markov model-based health economic evaluations investigating the cost-effectiveness of Latanoprost in the French setting were identified. One compared Travoprost versus Latanoprost [51], whereas the second compared Travoprost to Latanoprost + Timolol [52]. Similar to the present study, both cost-effectiveness analyses adopted a 5-year time horizon. However, given important differences in the chosen patient population, comparators and outcomes (i.e., IOP reduction was used as an outcome in one of the health economic evaluations [51]), a meaningful comparison with the present study was unfeasible.

A second limitation is related to OAG/OHT stage-specific utility values. The utility values for OAG/OHT stages 0 and 5, that were taken from literature [29], represented the lower and the upper bounds of the scale and might have influenced the remaining OAG/OHT stage-specific utility values provided by experts. However, the adopted approach is similar to the usual layout of the visual analogue scale (VAS), which requests subjects to express their current health state in between the provided bounds of the VAS [15]. Furthermore, given their extensive experience in managing OAG/OHT, the glaucoma specialists were assumed to be reliable proxies for patients' utility. Additionally, comprehensive sensitivity analyses demonstrated that potential uncertainty surrounding the estimated utility values did not impact the base case results.

Lastly, in the present study, surgical or laser treatment was not considered as an alternative to OAG/OHT medications [48]. The rationale supporting this exclusion rests on the fact that laser therapy and surgical procedures are generally performed as second-line treatment options, after topical IOP-lowering eye drops have failed to sufficiently control IOP.

In conclusion, our results demonstrate that STN1013001 will provide the French health system with a cost-effective therapy versus Latanoprost for OAG/OHT + OSD patients.

These findings should be confirmed empirically when STN1013001 is available to patients in France.

**Data Availability**

Data are provided in the Supplementary Information files submitted alongside the manuscript.

**Conflicts of Interest**

CvS and LA are employees of Santen GmbH, München, Germany.

**Acknowledgments**

This study was funded by Santen GmbH, München, Germany. CL has received an unconditional research grant from Santen GmbH, München, Germany. Outside this research, in the past three years, CL has received research grants, speaker or consultancy fees from AstraZeneca S.p.A, Boehringer Ingelheim Italia S.p.A., CSL Behring S.p.A., Ferring S.p.A., Ipsen S.p.A., Roche S.p.A., Sanofi s.r.l., Santen GmbH, Shire S.p.A, and Sobi S.p.A. FA and CS received research fee for data provision from Santen GmbH, München, Germany. The authors express their gratitude for data provision to Prof. Antoine Labbé, Hospital Ophthalmology Unit, Teaching Hospital Department at Quinze-Vingts National Ophthalmology Hospital, Paris; Prof. Jean Paul Renard, Private Office of Ophthalmology, Paris, and Dr. Eric Sellem, Private Eye Clinic Centre Ophtalmologique Kléber, Lyon, Paris.

**Supplementary Materials**

SText. Probabilistic sensitivity analysis: essential glossary *Figure S1*. Base case analysis-results-mean cost per patient per OAG/OHT stage<sup>a,b</sup>. *Figure S2*. Base case analysis-results-mean QALYs per patient per OAG/OHT stage<sup>a,b</sup>. *Table S1*. Base case analysis-methods-OAG/OHT staging<sup>a</sup>. *Table S2*. Base case analysis-methods-transition probability matrix (95% CI)<sup>a</sup>. *Table S3*. Base case analysis-results-OAG/OHT patients' age (range). *Table S4*. Base case analysis-results-mean number (SD) of OAG/OHT notional patients in each Markov state during a 5-year time horizon. *Table S5*. Base case analysis-results-adherence probabilities to OAG/OHT medications (95% CI)<sup>a,b</sup>. *Table S6*. Base case analysis-results-healthcare resource average consumption (95% CI)<sup>a</sup>-diagnosis. *Table S7*. Base case analysis-results-healthcare resource average consumption-management and follow-up-I-add-on therapies and drugs (range)a. *Table S8*. Base case analysis-results-healthcare resource average consumption (95% CI)<sup>a</sup>-management and follow-up-II-healthcare procedures and specialist visits. *Table S9*. Base case analysis-results-healthcare resource average consumption-OSD management-I-drugs<sup>a,b</sup>. *Table S10*. Base case analysis-results-healthcare resource average consumption (95% CI)<sup>a,b</sup>-OSD management-II-healthcare procedures and specialist visits. (*Supplementary Materials*)

**References**

[1] Y. C. Tham, X. Li, T. Y. Wong, H. A. Quigley, T. Aung, and C. Y. Cheng, "Global prevalence of glaucoma and projections of glaucoma burden through 2040: a systematic review and meta-analysis," *Ophthalmology*, vol. 121, pp. 2081-2090, 2014.

[2] R. M. Feldman, G. A. Cioffi, J. M. Liebmann, and R. N. Weinreb, "Current knowledge and attitudes concerning cost-effectiveness in glaucoma pharmacotherapy: a glaucoma
specialists focus group study,” *Clinical Ophthalmology*, vol. 14, pp. 729–739, 2020.

[3] D. Kreft, G. Doblhammer, R. F. Guthoff, and S. Frech, "Prevalence, incidence, and risk factors of primary open-angle glaucoma - a cohort study based on longitudinal data from a German public health insurance," *BMC Public Health*, vol. 19, p. 851, 2019.

[4] C. E. Traverso, J. G. Walt, S. P. Kelly et al., "Direct costs of glaucoma and severity of the disease: a multinational long term study of resource utilisation in Europe," *British Journal of Ophthalmology*, vol. 89, pp. 1245–1249, 2005.

[5] European Glaucoma Society, *Terminology and Guidelines for Glaucoma*, PubliComm, Savona, Italy, 5th edition, 2020.

[6] R. P. Mills, D. L. Budenz, P. P. Lee et al., "Categorizing the stage of glaucoma from pre-diagnosis to end-stage disease," *American Journal of Ophthalmology*, vol. 141, pp. 24–30, 2006.

[7] A. van Gestel, C. A. Webers, J. L. Severens et al., "The long-term outcomes of four alternative treatment strategies for primary open-angle glaucoma," *Acta Ophthalmologica*, vol. 90, pp. 20–31, 2012.

[8] J. P. Nordmann, A. Lafuma, and G. Berdeaux, "Modelling the lifetime economic consequences of glaucoma in France," *Journal of Medical Economics*, vol. 12, pp. 9–16, 2009.

[9] C. Y. Kim, K. H. Park, J. Ahn et al., "Treatment patterns and medication adherence of patients with glaucoma in South Korea," *British Journal of Ophthalmology*, vol. 101, pp. 801–807, 2017.

[10] P. A. Newman-Casey, L. M. Niziol, B. W. Gillespie, N. K. Janz, P. R. Lichter, and D. C. Musch, "The association between medication adherence and visual field progression in the collaborative initial glaucoma treatment study," *Ophthalmology*, vol. 127, pp. 477–483, 2020.

[11] P. J. Pisella, P. Pouliquen, and C. Baudouin, "Prevalence of ocular symptoms and signs with preserved and preservative free glaucoma medication," *British Journal of Ophthalmology*, vol. 86, pp. 418–423, 2002.

[12] J. Garcia-Feijoo and J. R. Sampaolesi, "A multicenter evaluation of ocular surface disease prevalence in patients with glaucoma," *Clinical Ophthalmology*, vol. 6, pp. 441–446, 2012.

[13] F. Lallemand, P. Daull, S. Benita, R. Buggage, and J. S. Garrigue, "Successfully improving ocular drug delivery using the cationic nanoemulsion, novasorb," *Journal of Drug Delivery*, vol. 2012, Article ID 604204, 2012.

[14] A. Briggs, M. Sculpher, and K. Claxton, *Decision Modelling for Health Economic Evaluation*, Oxford University Press, Oxford, UK, 2006.

[15] M. F. Drummond, M. J. Sculpher, K. Claxton, G. L. Stoddart, and G. W. Torrance, *Methods for the Economic Evaluation of Health Care Programmes*, Oxford University Press, Oxford, UK, 4th edition, 2015.

[16] P. J. Neumann, T. G. Ganiats, L. B. Russell, G. D. Sanders, and J. E. Siegel, Eds., *Cost-Effectiveness in Health and Medicine*, Oxford University Press, Oxford, UK, 2nd edition, 2016.

[17] A. L. Mylla Boso, E. Gasperi, L. Fernandes, V. P. Costa, and M. Alves, "Impact of ocular surface disease treatment in patients with glaucoma," *Clinical Ophthalmology*, vol. 14, pp. 103–111, 2020.

[18] P. Daull, M. Amrane, and J. S. Garrigue, "Novasorb® cationic nanoemulsion and latanoprost: the ideal combination for glaucoma management?" *Journal of Eye Diseases and Disorders*, vol. 2, p. 14, 2017, https://www.longdom.org/open-access/novasorb-cationic-nanoemulsion-and-latanoprost-the-ideal-combination-forglaucoma-management.pdf.

[19] D. Ismail, M. Amrane, J. S. Garrigue, and R. Buggage, "A phase 2, randomized study evaluating the safety and efficacy of Catioprost® (unpreserved latanoprost 0.005% emulsion) compared to Travatan Z® in subjects with glaucoma and ocular surface disease," *Acta Ophthalmologica*, vol. 89, no. s248, 2011.

[20] U.S. National Library of Medicine, https://clinicaltrials.gov/ct2/show/NCT04133311?term=DE-130A&cond=glaucoma&draw=2&rank=1, 2019.

[21] L. Annemans, B. Genesté, and B. Jolain, "Early modelling for assessing health and economic outcomes of drug therapy," *Value in Health*, vol. 3, pp. 427–434, 2020.

[22] M. J. IJzerman, H. Koffijberg, E. Fenwick, and M. Krahn, "Emerging use of early health technology assessment in medical product development: a scoping review of the literature," *PharmacoEconomics*, vol. 35, no. 7, pp. 727–740, 2017.

[23] F. A. Sonnenberg and J. R. Beck, "Markov models in medical decision making: a practical guide," *Medical Decision Making*, vol. 13, pp. 322–339, 1993.

[24] Institut national de la statistique et des études économiques, https://www.insee.fr/fr/statistiques/3311422?sommaire=3311425(French), 2018.

[25] G. A. Wilkins, https://www.iau.org/static/publications/stylemanual1989.pdf, 1990.

[26] S. L. Lohr, *Sampling: Design and Analysis*, Brooks/Cole, Boston, MA, USA, 2nd edition, 2010.

[27] Haute Autorité de Santé (Has), https://www.has-sante.fr/upload/docs/application/pdf/2020-11/methodological_guidance_2020-choices_in_methods_for_economic_evalution.pdf, 2020.

[28] A. O'Hagan, C. E. Buck, A. Daneshkhah et al., *Uncertain Judgements: Eliciting Experts' Probabilities*, Wiley, Hoboken, NJ, USA, 2006.

[29] Canadian Agency for Drugs and Technology in Health (Cadth), https://cadth.ca/sites/default/files/cdr/pharmacoeconomic/sr0615-verkazia-pharmacoeconomic-review-report.pdf, 2020.

[30] A. van Gestel, C. A. Webers, H. J. Beckers et al., "The relationship between visual field loss in glaucoma and health-related quality-of-life," *Eye*, vol. 24, pp. 1759–1769, 2010.

[31] Iqvia, https://www.iqvia.com/solutions/commercialization/brand-strategy-and-management/market-measurement/midas, 2020.

[32] Ministère des Affaires sociales et de la Santé, https://medicprix.sante.gouv.fr/medicprix/welcome.do, 2020.

[33] Sécurité Sociale, https://www.ameli.fr/accueil-de-la-ccam/telechargement/index.php(French), 2020.

[34] A. H. Briggs, "Handling uncertainty in economic evaluation and presenting the results," in *Economic Evaluation in Health Care. Merging Theory with Practice*, M. Drummond and A. McGuire, Eds., pp. 172–214, Oxford University Press, Oxford, UK, 2001.

[35] W. Brouwer, F. Rutten, and M. Koopmanschap, "Costing in economic evaluations," in *Economic Evaluation in Health Care. Merging Theory with Practice*, M. Drummond and A. McGuire, Eds., pp. 68–93, Oxford University Press, Oxford, UK, 2001.

[36] World Health Organization, https://www.who.int/bulletin/archives/79%284%29373.pdf, 2001.

[37] European Network of Research Ethics Committees, https://www.eurecnet.org/legislation/france.html, 2020.

[38] M. Pagano and K. Gauvreau, *Principles of Biostatistics*, Brooks/Cole, Boston, MA, USA, 2nd edition, 2000.
[39] A. A. Stinnett and A. D. Paltiel, "Estimating CE ratios under second-order uncertainty: the mean ratio versus the ratio of means," *Medical Decision Making*, vol. 17, pp. 483-489, 1997.

[40] W. C. Black, "The CE plane: a graphic representation of cost-effectiveness," *Medical Decision Making*, vol. 10, no. 3, pp. 212-214, 1990.

[41] H. A. Glick, J. A. Doshi, S. A. Sonnad, and D. Polsky, *Economic Evaluation in Clinical Trials*, Oxford University Press, Oxford, UK, 2nd edition, 2014.

[42] A. A. Stinnett and J. Mullahy, "Net health benefits: a new framework for the analysis of uncertainty in cost-effectiveness analysis," *Medical Decision Making*, vol. 18, no. 2 Suppl, pp. S68-S80, 1998.

[43] E. Fenwick, K. Claxton, and M. Sculpher, "Representing uncertainty: the role of cost-effectiveness acceptability curves," *Health Economics*, vol. 10, no. 8, pp. 779-787, 2001.

[44] E. Fenwick, B. J. O'Brien, and A. Briggs, "Cost-effectiveness acceptability curves- facts, fallacies and frequently asked questions," *Health Economics*, vol. 13, no. 5, pp. 405-415, 2004.

[45] C. Cartier-Bechu, A. Gherardi, M. Sivignon et al., "Is there a threshold in France?: first exhaustive review of published health-economic appraisals by the haute autorite de sante (HAS), (French national authority for health)," *Value in Health*, vol. 19, p. PA490, 2016.

[46] F. Aptel and A. Toren, "Societal and economic impact of poor glaucoma medication adherence," *Ophthalmology*, vol. 127, no. 5, pp. 599-600, 2020.

[47] U. Thelen, D. Schnober, S. Schölzel et al., "Long-term cost and efficacy analysis of latanoprost versus timolol in glaucoma patients in Germany," *International Journal of Ophthalmology*, vol. 6, pp. 155-159, 2013.

[48] J. Bartelt-Hofer, L. Ben-Debba, and S. Flessa, "Systematic review of economic evaluations in primary open-angle glaucoma: decision analytic modeling insights," *Pharma-coEconomics - Open*, vol. 4, no. 1, pp. 5-12, 2020.

[49] M. J. Buxton, M. F. Drummond, B. A. Van Hout et al., "Modelling in economic evaluation: an unavoidable fact of life," *Health Economics*, vol. 6, no. 3, pp. 217-227, 1997.

[50] M. C. Weinstein, B. O'Brien, J. Hornberger et al., "principles of good practice for decision analytic modeling in health-care evaluation: report of the ISPOR task force on good research practices--modeling studies," *Value in Health*, vol. 6, pp. 9-17, 2003.

[51] S. Payet, P. Denis, G. Berdeaux, and R. Launois, "Assessment of the cost effectiveness of travoprost versus latanoprost as single agents for treatment of glaucoma in France," *Clinical Drug Investigation*, vol. 28, pp. 183-198, 2008.

[52] C. Le Pen, M. Ligier, and G. Berdeaux, "Cost-effectiveness and cost-utility analysis of travoprost versus latanoprost and timolol in the treatment of advanced glaucoma in five European countries: Austria, France, Germany, The Netherlands and the United Kingdom," *Journal of Medical Economics*, vol. 8, pp. 67-84, 2005.