# Multicriteria Optimization Model for the Study of the Efficacy of Skin Antiaging Therapy

Authors: Maria Crişan; Guido Cappare; Luciana Neamţiu; Ioana Chiorean; Liana Lupşa; Diana Crişan; Radu Badea
Journal: Computational and Mathematical Methods in Medicine
 (2012)

---

## Abstract

The evolution of the cutaneous structure after topical treatment with P63 antiaging complex, assessed with high frequency ultrasound, is studied by means of multicriteria optimization model. Due to the fact that the impact of the treatment may influence the quality of life, a medical index which measures, from this point of view, the efficacy of the treatment is given, also taking into account medical and economical aspects.


---

## Body

## 1. Introduction

The basic idea of Pharmaco-Economics studies is to gain a physical and psychical comfort state for as long as possible, with the smallest amount of money. Therefore, according to these studies (see [1]) which consider a treatment in terms of results related to costs, one of the following five types of analysis is used: and cost-effects (CEA), cost-minimization (CMA), cost-utility (CUA), cost-efficiency (CEAC), and cost-benefit (CBA). Then, the data obtained after each of this analysis are used to compare two or more treatments. Cost-utility analysis was developed to help decision makers compare the value of alternative interventions that have very different health benefits, and it facilitates these comparisons without recourse to placing monetary values on different health states. The primary outcome of a cost-utility analysis is the cost per quality-adjusted life years (QALYs) or incremental cost-effectiveness ratio (ICER), which is calculated as the difference in the expected cost of two interventions, divided by the difference in the expected QALYs produced by the two interventions. QALYs measure health as a combination of the duration of life and the health-related quality of life. Also, there is another index, denoted by NB or INB, which means incremental net benefit, defined by(1)INB=λΔe-Δc,
where λ is the willingness to pay.These indexes are largely used in the literature (see, e.g., [2–5]). Unfortunately, they may not always reflect with sufficient accuracy all aspects of medical outcomes of treatment, perception and impact of treatment on the patient’s psyche, economic effects, and so forth.Using the multicriteria optimization technique, in [6] a new index, called medicoeconomic index of a treatment (denoted by MEI), is introduced. Its construction may use all the desired aspects. It permits a simultaneous comparison of two or more medical treatments. And, in addition, due to the fact that it emphasizes the importance of every aspect in a general context, it gains an increased flexibility (see [7–10]).There are papers in the literature which use multifactorial decisions to compare the medical treatments (see the book [11] and the studies [12, 13]).The purpose of this paper is to show how the MEI index may best quantify the effect of applying a skin rejuvenation treatment, taking into account multiple aspects: changes in skin parameters occurred after treatment application, side effects, how the patient perceives the treatment and outcome of application, cost of the treatment.Quantification, materialized by building the MEI index, permits comparisons among several types of treatments or comparing the results obtained with the same treatment applied to different groups of patients, groups that differ either by age or by phototype of the patients.


## 2. The Model on a Specific Sample

LetT be a medical treatment, applied against skin aging, whose effect we want to quantify. In particular we considered treatment with P63 antiaging complex.From the treatment point of view, the patients were classified in two groups of study: the first group included persons to whom the new medication to be tested was applied, and the second group, the control group, to whom the old medication was given. It is important to mention that both groups were homogenous relatively to patient age, health state, and phenotype. From medical point of view, the obtained results refer to thickness of the epidermis (mm), thickness of the dermis (mm), number of low echogenity pixels (LEP), medium and high echogenity pixels (MEP and HEP), and the LEPs/LEPi ratio (low echogenic pixels in the upper dermis/low echogenic pixels in the lower dermis). All these imagistic parameters quantify both the skin aging process and the efficacy of the antiaging treatments. For each of these parameters, the mean values for every study group, age group, and phenotype group are given.Taking this into account, for each treatmentTi,i∈{1,2}, each age group j,j∈{1,2,3}, and each phenotype h,h∈{2,3}, we indicate the following:(i)
Eijh1,Eijh2, the mean values of the thickness of epidermis of the patients who take the treatment i, are in age group j and have phenotype h, before the treatment, and, respectively, after the treatment;(ii)
Dijh1,Dijh2, the mean values of the thickness of dermis of the patients who take the treatment i, are in age group j and have phototype h, before the treatment, respectively, after the treatment;(iii)
Lijh1,Lijh2, the mean values of the LEP of the patients who take the treatment i, are in age group j and have phototype h, before the treatment, respectively after the treatment;(iv)
Mijh1,Mijh2, the mean values of the MEP of the patients who take the treatment i, are in age group j and have phototype h, before the treatment, respectively, after the treatment;(v)
Hijh1,Hijk2, the mean values of the HEP of the patients who take the treatment i, are in age group j and have phototype h, before the treatment, respectively, after the treatment;(vi)
Fijh1,Fijk2, the mean values of LEPi/LEPu;(vii)
Rijh%, the percent of those with no adverse effects (erythema, pruritus, ocular disturbance, etc.);(viii)
Cijh, the mean value of the cost of treatment for the period of time when the study was performed, taking into account a person in group i, of age j and a phototype h.The patients were asked to fill in a questionnaire and give two marks, denoted byq and s∈{1,…,10}, concerning the convenience of treatment application and the consequential effect of the treatment, respectively. We denote by Qijh and Sijk, respectively, the mean values of the marks given by the patients in age group j for the treatment i, and which have phototype h.LetK be the maximum accepted value for the cost of the treatment for one patient during the study period of time.


### 2.1. Methods

The method used for comparison was that of introducing a medicoeconomic index, by means of multicriteria optimization, which permits a rigorous observation of the evolution of our patients.For a protocol treatmentTi, the percentage Pijh% of patients from treatment group i, age group j, and phototype h, is known.Our problem was the following: having the previous information and taking into account the cost of the new treatment and the opinions of the patients, we wanted to determine the efficacy of the new treatment, and consequently, if it was worth using it. The expected medical results are increase of the mean thickness of epidermis and dermis, decrease of the mean number of low echogenity pixels, the increase of medium and high echogenity pixels, and the increase of LEPs/LEPi ratio. From economical point of view, it is desirable that the mean value of the cost for the new treatment to be as low as possible.Finally, from the patients’ point of view, the effect of the treatment is expected to be maximum possible and its application as comfortable as it can be.In order to solve this problem, we compared the treatments at a general level (group 1 versus group 2 of treatments), but also on each age group and phenotype group. In every case, we had to solve a mathematical multicriteria optimization problem.We introduced the variablesx1 and x2, which take only the values 0 and 1. We have: x1=1, if the first treatment protocol is chosen, x1=0, if not; x2=1, if the second treatment protocol is chosen, x2=0, if not. Due to the fact that one patient has to follow one and only one treatment protocol, we have the obvious condition x1+x2=1.Then we consider the functions:fk:{0,1}2→R,k∈{1,…,60}, given for every j∈{1,2,3} and h∈{2,3}, respectively, by(2)f2(j-1)+(h-1)(x1,x2)=P1jh(E1jh2E1jh1)x1+P2jh(E2jh2E2jh1)x2,∀(x1,x2)∈{0,1}2,f6+2(j-1)+(h-1)(x1,x2)=P1jh(D1jh2D1jh1)x1+P2jh(D2jh2D2jh1)x2,∀(x1,x2)∈{0,1}2,f12+2(j-1)+(h-1)(x1,x2)=P1jh(L1jh1L1jh2)x1+P2jh(L2jh1L2jh2)x2,∀(x1,x2)∈{0,1}2,f18+2(j-1)+(h-1)(x1,x2)=P1jh(M1jh2M1jh1)x1+P2jh(M2jh2M2jh1)x2,∀(x1,x2)∈{0,1}2,f24+2(j-1)+(h-1)(x1,x2)=P1jh(H1jh2H1jh1)x1+P2jh(H2jh2H2jh1)x2,∀(x1,x2)∈{0,1}2,f30+2(j-1)+(h-1)(x1,x2)=P1jh(F1jh2F1jh1)x1+P2jh(F2jh2F2jh1)x2,∀(x1,x2)∈{0,1}2,f36+2(j-1)+(h-1)(x1,x2)=-(R1jhP1jh)x1-(R2jh2P2jh)x2,∀(x1,x2)∈{0,1}2,f42+2(j-1)+(h-1)(x1,x2)=P1jh(Q1jh10)x1+P2jh(Q2jh10)x2,∀(x1,x2)∈{0,1}2,f48+2(j-1)+(h-1)(x1,x2)=P1jh(S1jh10)x1+P2jh(S2jh10)x2,∀(x1,x2)∈{0,1}2,f54+2(j-1)+(h-1)(x1,x2)=-P1jh(C1jhK)x1-P2jh(C2jhK)x2,∀(x1,x2)∈{0,1}2.Obviously, our purpose is to determine, simultaneously, their maximum. Generally, this fact cannot be fulfilled; therefore, a compromised solution is to be accepted: a synthesis function is considered and each of the previous functions gets a weight, according to the importance of every criterion in the problem.So, firstly we attribute weights to the main elements of our final purpose:(i)
getting the best medical results (purpose indicated by the functionsf1,…,f36,f49,…,f54);(ii)
the best reception of the patient (functionsf43,…,f54);(iii)
the lowest possible cost of the treatment (indicated by the functionsf55,…,f60).Each of these purposes gets a nonnegative weight denoted bytu(u∈{1,2,3}), with the property that t1+t2+t3=1. Obviously, the three weights may be or may be not equal.Remark 1.
In some other cases, we may distinguish these weights even further. More precisely, the weightt2, corresponding to the way in which the patient is contained with the treatment may be considered as the sum between two others weights: tQ, for the convenience of treatment application and tS for the degree of content after the treatment. Also, the weight t1 may be seen as the sum of 7 weights, tE,tD,tL,tM,tH,tF,tR, connected with the importance given to the increase of the thickness of the epidermis, dermis, decreasing of low echogenity pixels, increasing of medium and high echogenity pixels and increasing of ratio LEPi/LEPu, respectively decreasing of the adverse reactions. We change the notation of t3 with tC.


## 2.1. Methods

The method used for comparison was that of introducing a medicoeconomic index, by means of multicriteria optimization, which permits a rigorous observation of the evolution of our patients.For a protocol treatmentTi, the percentage Pijh% of patients from treatment group i, age group j, and phototype h, is known.Our problem was the following: having the previous information and taking into account the cost of the new treatment and the opinions of the patients, we wanted to determine the efficacy of the new treatment, and consequently, if it was worth using it. The expected medical results are increase of the mean thickness of epidermis and dermis, decrease of the mean number of low echogenity pixels, the increase of medium and high echogenity pixels, and the increase of LEPs/LEPi ratio. From economical point of view, it is desirable that the mean value of the cost for the new treatment to be as low as possible.Finally, from the patients’ point of view, the effect of the treatment is expected to be maximum possible and its application as comfortable as it can be.In order to solve this problem, we compared the treatments at a general level (group 1 versus group 2 of treatments), but also on each age group and phenotype group. In every case, we had to solve a mathematical multicriteria optimization problem.We introduced the variablesx1 and x2, which take only the values 0 and 1. We have: x1=1, if the first treatment protocol is chosen, x1=0, if not; x2=1, if the second treatment protocol is chosen, x2=0, if not. Due to the fact that one patient has to follow one and only one treatment protocol, we have the obvious condition x1+x2=1.Then we consider the functions:fk:{0,1}2→R,k∈{1,…,60}, given for every j∈{1,2,3} and h∈{2,3}, respectively, by(2)f2(j-1)+(h-1)(x1,x2)=P1jh(E1jh2E1jh1)x1+P2jh(E2jh2E2jh1)x2,∀(x1,x2)∈{0,1}2,f6+2(j-1)+(h-1)(x1,x2)=P1jh(D1jh2D1jh1)x1+P2jh(D2jh2D2jh1)x2,∀(x1,x2)∈{0,1}2,f12+2(j-1)+(h-1)(x1,x2)=P1jh(L1jh1L1jh2)x1+P2jh(L2jh1L2jh2)x2,∀(x1,x2)∈{0,1}2,f18+2(j-1)+(h-1)(x1,x2)=P1jh(M1jh2M1jh1)x1+P2jh(M2jh2M2jh1)x2,∀(x1,x2)∈{0,1}2,f24+2(j-1)+(h-1)(x1,x2)=P1jh(H1jh2H1jh1)x1+P2jh(H2jh2H2jh1)x2,∀(x1,x2)∈{0,1}2,f30+2(j-1)+(h-1)(x1,x2)=P1jh(F1jh2F1jh1)x1+P2jh(F2jh2F2jh1)x2,∀(x1,x2)∈{0,1}2,f36+2(j-1)+(h-1)(x1,x2)=-(R1jhP1jh)x1-(R2jh2P2jh)x2,∀(x1,x2)∈{0,1}2,f42+2(j-1)+(h-1)(x1,x2)=P1jh(Q1jh10)x1+P2jh(Q2jh10)x2,∀(x1,x2)∈{0,1}2,f48+2(j-1)+(h-1)(x1,x2)=P1jh(S1jh10)x1+P2jh(S2jh10)x2,∀(x1,x2)∈{0,1}2,f54+2(j-1)+(h-1)(x1,x2)=-P1jh(C1jhK)x1-P2jh(C2jhK)x2,∀(x1,x2)∈{0,1}2.Obviously, our purpose is to determine, simultaneously, their maximum. Generally, this fact cannot be fulfilled; therefore, a compromised solution is to be accepted: a synthesis function is considered and each of the previous functions gets a weight, according to the importance of every criterion in the problem.So, firstly we attribute weights to the main elements of our final purpose:(i)
getting the best medical results (purpose indicated by the functionsf1,…,f36,f49,…,f54);(ii)
the best reception of the patient (functionsf43,…,f54);(iii)
the lowest possible cost of the treatment (indicated by the functionsf55,…,f60).Each of these purposes gets a nonnegative weight denoted bytu(u∈{1,2,3}), with the property that t1+t2+t3=1. Obviously, the three weights may be or may be not equal.Remark 1.
In some other cases, we may distinguish these weights even further. More precisely, the weightt2, corresponding to the way in which the patient is contained with the treatment may be considered as the sum between two others weights: tQ, for the convenience of treatment application and tS for the degree of content after the treatment. Also, the weight t1 may be seen as the sum of 7 weights, tE,tD,tL,tM,tH,tF,tR, connected with the importance given to the increase of the thickness of the epidermis, dermis, decreasing of low echogenity pixels, increasing of medium and high echogenity pixels and increasing of ratio LEPi/LEPu, respectively decreasing of the adverse reactions. We change the notation of t3 with tC.


## 3. Comparing the Results at a General Level and Constructing the MEI Index

In this case, the mathematical model of our medical problem is the following multicriterial optimization problem with the objective functionF=(f1,…,f60):{0,1}2→R60, previously given, and the set of possible solutions X={(x1,x2)∈{0,1}2∣x1+x2=1}.In order to solve this problem, we apply the weight method and we construct the synthesis functionΛ:{0,1}2→R,Λ(x1,x2)=∑i=160λifi(x1,x2), for every x=(x1,x2)∈{0,1}2, where λi is the importance (weight) given to the criteria generated by function fi. Taking into account the previous considerations, we set λ1=⋯=λ6=tE,
λ7=⋯=λ12=tD,
λ13=⋯=λ18=tL,
λ19=⋯=λ24=tM,
λ25=⋯=λ30=tH,
λ31=⋯=λ36=tF,
λ37=⋯=λ42=tR,
λ43=⋯=λ48=tQ,
λ49=⋯=λ54=tS,
λ55=⋯=λ60=tC.Any optimum solution of the problem(3)Λ(x1,x2)⟶max⁡,(x1,x2)∈X
is a max-efficient point (max-Pareto point), so is acceptable. Due to the fact that condition x1+x2=1 implies x1=1-x2, the solving of problem (3) reduces to the solving of the following problem:
(4)ψ(u)=∑i=160λifi(u,1-u)⟶max⁡,u∈{0,1}.
The simple form of the restrictions conduces immediately to the conclusion: if ψ(1)>ψ(0), then u=1 is the optimal solution of problem (4); if not, u=0 is the optimal solution. This mathematical result permits us to introduce the MEI index, called medicoeconomic index, attached to the protocol Ti, by
(5)MEITi=∑j=13∑h=12[pijh(tEEijh2Eijh1+tDDijh2Dijh1+tLLijh2Lijh1+tMMijh2Mijh1+tHHijh2Hijh1+tFFijh2Fijh1)-tRRijkpijk]+∑j=13∑h=12[pijh(tQQijh2Qijh1+tSSijh2Sijh1-tCCijhK)].
In case of small or homogenous samples, we may use the easier formula:
(6)MEITi=∑j=13∑h=12(tEEijh2Eijh1+tDDijh2Dijh1+tLLijh2Lijh1+tMMijh2Mijh1+tHHijh2Hijh1+tFFijh2Fijh1-tRRijh)+∑j=13∑h=12(tQQijh2Qijh1+tSSijh2Sijh1-tCCijhK).Remark 2.
Since for the construction of this index we consider only data related to the patient’s condition before and after treatment and data referring only this treatment, the index can be used to compare several treatments, unlike ICER and NB that allow comparison of only two treatments and of only two criteria related to cost and effectiveness.


## 4. Special Cases

If we want to compare the overall effectiveness of treatment for two different age groups, we may follow the same steps, except that we consider only the data relative to an age group.For example, if we want to compare the efficiency of the treatment for each age group, i.e. 40–49, 50–59 and over 60, we construct the MEI index by taking(7)MEIgrupj=∑h=12[P1jh(tEE1jh2E1jh1+tDD1jh2D1jh1+tLL1jh2L1jh1+tMM1jh2M1jh1+tHH1jh2H1jh1)-tRR1jhP1jh]+∑h=12[P1jh(tFF1jh2F1jh1+tQQ1jh2Q1jh1+tSS1jh2S1jh1-tCC1jhK)].
In case of small or homogenous samples, an easier formula for computing the MEI index is the following:
(8)MEIgrupj=∑h=12(tEE1jh2E1jh1+tDD1jh2D1jh1+tLL1jh2L1jh1+tMM1jh2M1jh1+tHH1jh2H1jh1+tFF1jh2F1jh1-tRR1jh)+∑h=12(tQQ1jh2Q1jh1+tSS1jh2S1jh1-tCC1jhK).
Another special case may be the one connected with the phototype. In this case, we have
(9)MEIphenotiph=∑j=13[P1jh(tEE1jh2E1jh1+tDD1jh2D1jh1+tLL1jh2L1jh1+tMM1jh2M1jh1+tHH1jh2H1jh1)-tRR1jhP1jh]+∑j=13[P1jh(tFF1jh2F1jh1+tQQ1jh2Q1jh1+tSS1jh2S1jh1-tCC1jhK)].
Or, the easier formula is as follows:
(10)MEIphenotiph=∑j=13tEE1jh2E1jh1+tDD1jh2D1jh1+tLL1jh2L1jh1+tMM1jh2M1jh1+tHH1jh2H1jh1+tFF1jh2F1jh1+tQQ1jh2Q1jh1+tSS1jh2S1jh1-tCC1jhK-trR1jk.
Also, if we want a result based both on the age group and phototype, we have
(11)MEIj,h=tEE1jh2E1jh1+tDD1jh2D1jh1+tLL1jh2L1jh1+tMM1jh2M1jh1+tHH1jh2H1jh1+tFF1jh2F1jh1tQQ1jh2Q1jh1+tSS1jh2S1jh1-tCC1jh2K-tRR1jh.


## 5. Application of the P63 Antiaging Treatment

We exemplify the application of this index for the characterization of the treatment based on P63 antiaging treatment (a metabolic dynamiser composed of alpha hydroxyacids, retinoids, a biomimetic peptidic complex, and gluconolactone incapsulated in liposomes), in the particular case of data presented in the following table and taken from [14]. Our mathematical model was applied on a sample of 30 patients with ages between 40 and 75 years old. Due to our purpose, we divided them in three age interval groups: the first group included patients aged 40–50 years, the second group included the aged ones 50–60, and the last group contained those over 60. All patients included in the study belonged to phototype class II or III, and, therefore, following this criteria, we grouped them into two classes, denoted by 2 and 3. The treatment cost equals 250  monetary units (m.u); the K constant is taken equal to 400. For those patients who did not follow the treatment, the cost was considered 100 m.u.The medical data concerning the group treated with P63 antiaging treatment, as well as the placebo group, are given in the following tables.In order to compute the MEI index in case of P63 antiaging complex treatment, and placebo one, formula (6) was used. It is important to mention that no side effects were noticed.We considered equally the weights given for the age groups. We taketE=1/8,tD=1/8,tF=1/2,tR=tQ=tS=0, and tC=1/4. Under these circumstances and using the data from Table 1, the MEI index for the treatment with P63 antiaging complex was TP63complex=0.812143. Using the data of Table 4 to build the MEI index for the placebo treatment, we obtained T0=0.569481. This means that the treatment is worth doing.Table 1
Group treated with P63 antiaging treatment.
ParametersMean value before treatmentMean value after treatmentThickness of epidermis (mm)0.1170.135Thickness of dermis (mm)1.5371.719LEPs/LEPi1.1491.574Using the same weights, the data from Table2 and denoting by “Mean” the mean value, b.t. = before treatment and a.t. = after treatment and computing MEI for every age group, we also got the following results: MEI40–49=0.724915,MEI50–59=0.897944, and MEI60–75=0.745013.Table 2
Group treated with P63 complex, based on age groups.
ParametersMean b.t. 40–49Mean a.t. 40–49Mean b.t. 50–59Mean a.t. 50–59Mean b.t. 60–75Mean a.t. 60–75Thickness of epidermis (mm)0.1140.1353330.116750.13450.120250.136125Thickness of dermis (mm)1.85851.9011.3863131.7405631.3676251.506625LEPs/LEPi1.2045931.4573531.0960621.3980471.2138421.510156The highest value was obtained for the 50–59 group, then for the group 60–75 and then for 40–49 group. It means that the age group 50–59 had the best response to the treatment.Using the same weights, the data from Table3 and computing MEI, according with the formula for phototypes 2 and 3, we get MEI2 = 0.782524, respectively, MEI3 = 0.84888. The greatest value is obtained for phototype 3. This result indicates that the treatment was best perceived for phototype 3.Table 3
Group treated with P63 antiaging treatment, based on phototype.
ParametersMean b.t. phototype 2Mean a.t. phototype 2Mean b.t. phototype 3Mean a.t. phototype 3Thickness of epidermis (mm)0.1146250.1299380.120.141Thickness of dermis (mm)1.6053131.7680631.4605711.644214LEPs/LEPi1.1986611.5808021.0926231.567998Table 4
Placebo group.
ParametersMean value before treatmentMean value after treatmentThickness of epidermis (mm)0.13580.13505Thickness of dermis (mm)1.496351.57155LEPs/LEPi0.9247571.042922Remark 3.
We mention that a change in the weights leads to a change of the index.


## 6. Conclusions

As mentioned before, starting with effective data concerning the treatment with P63 complex, we studied the possibility to quantify the outcome of antiaging skin treatment, from multiple points of view: medical parameters, side effects, cost, convenience of treatment application, the way in which the patient perceive the treatment.Using multicriteria optimization, we introduced a medicoeconomic index (MEI) capable to perform the required quantification. In its construction, only data connected with a specific treatment were used; therefore, it has the advantage of being able to compare more than 2 treatments.Further, we gave an application in which we computed the MEI index, in order to compare the P63 antiaging complex treatment with a placebo one (which provided only hydration of the skin). We also constructed specific indexes, according to age groups. This permitted us to compare the effects of P63 complex treatment for different age groups and to establish, from this point of view, a hierarchy of these groups. The same constructions and comparisons were made for the groups of phototypes 2 and 3.The application presented documents the importance and the facility of use of our index.



---

*Source: 103919-2012-02-28.xml*