this randomized clinical trail study was performed on 150 patients in ophthalmic research center of labbafinejad medical center ( tehran , iran ) during two years ending in 2007 . considering 90% power for detection of a 0.2-logmar difference ( equal to 2 snellen lines ) in visual acuity and being significant at the two - sided 5% level with an assumed standard deviation ( sd ) of 0.33 , 50 eyes in each group were required.15 we selected 103 patients eye with complete information from these 150 eyes .
the informed written consent was obtained from all patients and the protocol was approved by the ethical committee of ophthalmic research center of labbafinejad medical center ( tehran , iran ) . in this study , patients were randomly allocated into three following groups : 1 ) the ivb group , patients who received 1.25 mg ivb ( 32 subjects ) ; 2 ) the ivb / ivt group , patients who received 1.25 mg of ivb and 2 mg of ivt ( 32 subjects ) ; and 3 ) the mpc group , patients who underwent focal or modified grid laser ( 33 subjects ) . to decrease the side effect of ivb
, half doze was used instead of normal 4 mg injection . in the case of two eyes examination ,
each eye was studied separately and the best corrected visual acuity ( bcva ) and central macular thickness ( cmt ) of each patient were measured during the survey .
bcva was measured by the snellen chart and was recorded in logarithm of the minimum angle of resolution ( logmar ) scale .
retinal thickness was measured in a circle ( 3.5 mm in diameter ) centered on the fixation point .
mean thickness on the 1 mm circle centered on the fovea ( cmt ) was considered for statistical evaluation .
exclusion criteria were previous panretinal or focal laser photocoagulation , prior intraocular surgery or injection , history of glaucoma or ocular hypertension , va of 20/40 or better or worse than 20/300 , presence of iris neovascularization , high - risk proliferative diabetic retinopathy , and significant media opacity .
monocularity , pregnancy , serum creatinine  3 mg / dl , and uncontrolled diabetes mellitus were also among the exclusion criteria.15 the experiment was conducted in a triple - blind manner that patients , experimenters and data analyzers were not aware of the specific diagnostic methods .
for this reason , a non - therapeutic light was shed on the patients macula for 20 sec in ivb and ivb / ivt groups .
the injection was also simulated for mpc group.16 statistical modeling was based on the data recorded in weeks 6 , 12 , 24 , 36 , and 48 which the patient 's data was complete at these times .
bcva as a continuous variable and cmt as an ordinal variable were the response variables and the duration of diabetes ( > 10 years ) ( dd ) as the number of years having diabetes , sex and the baseline values of bcva and cmt were considered as confounder variables .
it should be noted that the cmt variable was ordered in three following groups:16 < 300 m , 300 - 399 m and > 400 m .
since the bivariate distribution of the recorded data was shown to be asymmetric or heavy - tailed based on the bivariate mardia 's skewness and kurtosis coefficients16 of 0.88 and 6.7 , asymmetric modeling was crucial .
after fitting different symmetric and skew - symmetric distribution and running the model , dics ( goodness of fit criteria of the models ) were compared in different models .
the results of comparing these criteria showed that the skew - symmetric models fitted better than the symmetric models . in addition , among asymmetric models , the distribution with more kurtosis was the best one .
the bivariate mixed linear model is shown in eq.1 and eq.2 . in this model ,
the constant random variable ( crv ) was used to compensate intra - subject variability while the bivariate model was chosen because of the relative correlation between cmt and bcva .
skew - symmetric distributions used in this model were either those that had more kurtosis in comparison with the normal one or those that had right skewness1517 such as skew - normal1 or skew - t.2 although the density functions of the skewed distributions are complicated , they are of interest for researchers because of the lower error variance and higher precision.17 in these models , traditional methods are not used to compute the coefficients . instead , bayesian methods are used . statistical modeling and inference
14.18 selection of the best fitted asymmetric model among different distributions was based on deviance information criteria ( dic ) which is a goodness of fit measure in bayesian methods .
data were reported as frequency ( percent ) for qualitative variables and mean  standard deviation ( sd ) for quantitative variables .
odds ratios ( or ) and their 95% confidence intervals ( cis ) were presented as the measure of association and p - values < 0.05 were considered as statistically significant .
their average diabetic duration was 10.24 1.2 years ( mean  sd ) . among the patients ,
twelve percent of the whole data were asymmetric or outlier based on some primary statistical method like z - score ( table 1 ) .
the summary statistics of best correlated visual acuity and central macular thickness at the beginning and during the follow - ups ( n=103 ) all confounders were significant in our model ( p < 0.05 ) ( table 2 ) .
these results are presented based on the best fitted model using dic . using this criterion ,
it was identified that the mixed asymmetric model of skew - t was the best model among the fitted models .
this distribution is one of the important asymmetric distributions that can cover more outliers than the normal distribution .
there was no significant difference between ivb and ivb / ivt groups for the treatment of dme ( p > 0.05 ) while these two groups showed significantly better therapeutic effect in comparison with the pmc group ( p < 0.05 ) ( table 2 ) .
considering the regression coefficients and the p - values , patients with longer diabetic duration showed better response to the treatment . in other words ,
the rate of decreasing cmt and increasing bcva was higher for patients whose diabetic duration was higher in comparison with those with lower diabetic duration .
it was also shown women showed better response to the treatment . taking into account the baseline values , patients with lower bcva and higher cmt showed better response to the treatment . it should be noted that since the cmt variable is ordinal , the odds ratio and its confidence interval was reported ( table 2 ) .
the regression results of the bivariate mixed linear asymmetric model for predictors of correlated visual acuity and central macular thickness in diabetic patients
using the specific design method in the analysis of clinical trials , it is possible not only to control the confounding variables but also to investigate the possible variation of the outcomes during time.19 non - normal data usually exists in medical records and it is rare to have normal data all the time . by removing asymmetric data that are non - normal data or outliers , it is possible to reduce the sample size so much that makes it difficult to analyze the data and correctly interpret the parameters.20 other proposed statistical methods used to analyze non - normal data , for example non - parametric analysis , have lower power in comparison with the parametric methods.15 therefore , in recent years , incorporating statistical methods that not only take into account the data correlation in repeated measurements but also tolerate the non - normality , has been of great interest among investigators .
longitudinal data analysis is recently used to more accurately analyze the clinical trials . in these methods , since there is a correlation between repeated measures , it is important to use approaches that take into account the data correlation as well as time - dependent effects of the treatment methods .
according to the results of the current study , intravitreal bevacizumab injection alone or in combination with intravitreal triamcinolone acetonide improved the bcva and significantly reduced the cmt that the association was not strictly proven to be significant in the previous researches using primary statistical methods.16 these improvements are much clearer up to the week 24 .
similar studies showed that the effect of bevacizumab injection is retained up to two weeks.21 in the literature , other investigators showed that intravitreal bevacizumab injection improves the vision in diabetic patients but intravitreal bevacizumab injection in combination with intravitreal triamcinolone acetonide does not improve the patients vision or reduce the cmt.2223 statistical analysis using repeated measurement on the same data set also showed intravitreal bevacizumab injection alone or in combination with intravitreal triamcinolone acetonide improved only the patients vision.15 the significancy of the patient 's sex and diabetic duration ( > 10 years ) on the patient 's responses to the treatment is in line with the previous studies.156 however , unlike the previous researches in which the effect of sex was not directly declared , we proved that women showed better responses to the treatments .
we also showed that the patient 's response to the treatment is better if the diabetic duration is higher which is in agreement with the previous researches.524 therefore , according to our statistical analysis , intravitreal bevacizumab injection alone or in combination with intravitreal triamcinolone acetonide was a suitable substitute to macular laser photocoagulation as a primary treatment for diabetic macular edema .
in addition , the aforementioned treatment methods had more positive effect on the patients that had lower bcva and higher cmt as the baseline values . in this paper ,
we not only modeled clinical outcomes in the presence of confounding parameters including disease duration and sex , but we also did not remove any asymmetric data or outliers that enabled to increase the power of tests in comparison with other similar studies . considering the bivariate analysis , two outcomes of bcva and cmt were analyzed simultaneously that indeed gave the opportunity to analyze their relationship that increased the accuracy of the inference in comparison with the previous investigations . in spite of all aforementioned advantages
, this updated statistical analysis has a lot of computational difficulties and requires programming software which is hard to use .
no significant difference was observed in this research between intravitreal bevacizumab injections alone or in combination with intravitreal triamcinolone acetonide .
however , both of these methods significantly showed better improvement in comparison with that of macular laser photocoagulation as a traditional treatment of diabetic macular edema .
this improvement was shown in increasing the bcva and decreasing the cmt during one year clinical follow - up .
this research showed that when we have non - normal data in clinical trials , skew - symmetric mixed effect models may lead to different results compared to symmetric mixed models . with regard to the advantages of this updated statistical method
mm carried out the design and coordinated the study , and run all modelling studies and prepared the manuscript .
ak and ik provided assistance in the design of the study , technical advice on modelling and revised the manuscript .
fz provided assistance for all statistical results . ms coordinated and carried out all the experiments and provided the technical advice on medical discussion .
 background : diabetic macular edema ( dme ) is one of the major causes of visual loss and increase in central macular thickness ( cmt ) . 
 the aim of this study was to determine the efficacy of a single intravitreal injection of bevacizumab ( ivb ) alone or in combination with intravitreal triamcinolone acetonide ( ivb / ivt ) versus macular laser photocoagulation ( mpc ) as primary treatment for dme when confounders were considered.methods:skew-symmetric bivariate mixed modeling according to best corrected visual acuity ( bcva ) and cmt was done on the data of 103 diabetic patients from ophthalmic research center of labbafinejad medical center ( tehran , iran ) to determine the best dme treatment by adjusting the effect of confounders.results:although there was no significant difference between ivb / ivt ( p > 0.05 ) , these two treatments increased bcva and decreased cmt better than mpc ( p < 0.05 ) . 
 the following three groups showed better treatment responses : 1 ) women , 2 ) patients with more diabetes duration , 3 ) patients whose cmt were higher and va were lower at the beginning of the clinical trial.conclusions:using skew - symmetric mixed effect model as updated statistical method in presence of asymmetric or outlier data , we received different results compared to the same investigation on this study by analyzing bcva and cmt simultaneously . 
 this research demonstrated the effect of ivb alone or in combination with intravitreal ivb / ivt on visual power and decreasing cmt during follow up . 
Methods
Results
Discussion
Conclusions
Authors Contributions