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.
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 byINB=λΔ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
Let T 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 treatment Ti,i∈{1,2}, each age group j,j∈{1,2,3}, and each phenotype h,h∈{2,3}, we indicate the following:
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;
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;
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;
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;
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;
Fijh1,Fijk2, the mean values of LEPi/LEPu;
Rijh%, the percent of those with no adverse effects (erythema, pruritus, ocular disturbance, etc.);
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 by q 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.
Let K 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 treatment Ti, 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 variables x1 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, byf2(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:
getting the best medical results (purpose indicated by the functions f1,…,f36,f49,…,f54);
the best reception of the patient (functions f43,…,f54);
the lowest possible cost of the treatment (indicated by the functions f55,…,f60).
Each of these purposes gets a nonnegative weight denoted by tu(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 weight t2, 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 function F=(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
Λ(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:
ψ(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
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:
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
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:
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
MEIphenotiph=∑j=13[P1jh(tEE1jh2E1jh1+tDD1jh2D1jh1+tLL1jh2L1jh1+tMM1jh2M1jh1+tHH1jh2H1jh1)-tRR1jhP1jh]+∑j=13[P1jh(tFF1jh2F1jh1+tQQ1jh2Q1jh1+tSS1jh2S1jh1-tCC1jhK)].
Or, the easier formula is as follows:
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
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 take tE=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.
Group treated with P63 antiaging treatment.
Parameters
Mean value before treatment
Mean value after treatment
Thickness of epidermis (mm)
0.117
0.135
Thickness of dermis (mm)
1.537
1.719
LEPs/LEPi
1.149
1.574
Using the same weights, the data from Table 2 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.
Group treated with P63 complex, based on age groups.
Parameters
Mean b.t. 40–49
Mean a.t. 40–49
Mean b.t. 50–59
Mean a.t. 50–59
Mean b.t. 60–75
Mean a.t. 60–75
Thickness of epidermis (mm)
0.114
0.135333
0.11675
0.1345
0.12025
0.136125
Thickness of dermis (mm)
1.8585
1.901
1.386313
1.740563
1.367625
1.506625
LEPs/LEPi
1.204593
1.457353
1.096062
1.398047
1.213842
1.510156
The 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 Table 3 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.
Group treated with P63 antiaging treatment, based on phototype.
Parameters
Mean b.t. phototype 2
Mean a.t. phototype 2
Mean b.t. phototype 3
Mean a.t. phototype 3
Thickness of epidermis (mm)
0.114625
0.129938
0.12
0.141
Thickness of dermis (mm)
1.605313
1.768063
1.460571
1.644214
LEPs/LEPi
1.198661
1.580802
1.092623
1.567998
Placebo group.
Parameters
Mean value before treatment
Mean value after treatment
Thickness of epidermis (mm)
0.1358
0.13505
Thickness of dermis (mm)
1.49635
1.57155
LEPs/LEPi
0.924757
1.042922
Remark 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.
Acknowledgments
The authors gratefully acknowledge Cortex Technology for allowing the use of the ultrasound equipment and for their assistance with the application of the Dermascan to this project. The study is part of a national research grant CNCSIS-IDEI, no. 2624.
Pharmaceutical Research and Manufacturers of AmericaMethodological and conduct principles for pharmacoeconomic research1996822834BriggsA.ClaxtonK.SculpherM.2006Oxford, UKOxford University PressWillanA. R.Statistical analysis of cost-effectiveness data from randomized clinical trials2006633373462-s2.0-3374558727110.1586/14737167.6.3.337O'HaganA.StevensJ. W.A framework for cost-effectiveness analysis from clinical trial data20011043033152-s2.0-003495033910.1002/hec.617BriggsA.SculpherM.An introduction to Markov modelling for economic evaluation19981343974092-s2.0-003188607310.2165/00019053-199813040-00003LupsaL.A criterion for caracterizing a medical treatment that uses multicriteria programming in pharmacoeconomics1999Cluj-Napoca, RomaniaEditura Carpatica142146LupsaL.Construction of a medico-economic effectiveness index which characterizes a vaccine2001Cluj-Napoca, RomaniaSRIMA6165ChioreanI.LupsaL.NeamtiuL.SimosT. E.PsihoyiosG.TsitourasC.Markov models for the simulation of cancer screening process1048Proceedings of the International Conference on Numerical Analysis and Applied Mathematics2008American Institute of Physics143146NeamtiuL.2008Cluj-Napoca, RomaniaCasa Cărţii de ŞtiinţăGrosanC.AjithA.TiganS.ChangT. G.How to solve a multicriterion problem for which Pareto dominance relationship cannot be applied ? A case study from medicineProceedings of the 10th International Conference on Knowledge-Based Inteligent Information Engineering Systems (KES '06)2006SpringerLecture Notes in Computer ScienceTiganS.AchimasA.ComanI.DruganT.IacobE.2001Cluj-Napoca, RomaniaSRIMACimocaG.DollingerR.Specific quantitative and qualitative attribute for medical ranking/evaluation200181-2818TothI.GhermanL.ComanI.TiganS.On a multicriterial evaluation of the surgical techniques for anterior hypospadias20075245256CrisanM.CappareG.CrisanD.Ultrasonographic assessment of the cutaneous changes induced by topical antiaging INTERACTIVE P63 complex. Farmacia ISI-Science citation index expanded. CNCSIS A, In press