• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Smart Healthcare Using Data-Driven Prediction of Immunization Defaulters in Expanded Program on Immunization(EPI)

    2021-12-14 03:50:04SadafQaziMuhammadUsmanAzharMahmoodAaqifAfzaalAbbasiMuhammadAttiqueandYunyoungNam
    Computers Materials&Continua 2021年1期

    Sadaf Qazi,Muhammad Usman,Azhar Mahmood,Aaqif Afzaal Abbasi,Muhammad Attique and Yunyoung Nam

    1Predicitive Analytics Laboratory,Shaheed Zulfikar Ali Bhutto Institute of Science and Technology,Islamabad,44000,Pakistan

    2Department of Software Engineering,Foundation University Islamabad,Islamabad,44000,Pakistan

    3Department of Software,Sejong University,Seoul,05006,Korea

    4Department of Computer Science and Engineering,Soonchunhyang University,Asan,31538,Korea

    Abstract:Immunization is a noteworthy and proven tool for eliminating lifethreating infectious diseases,child mortality and morbidity.Expanded Program on Immunization(EPI)is a nation-wide program in Pakistan to implement immunization activities,however the coverage is quite low despite the accessibility of free vaccination.This study proposes a defaulter prediction model for accurate identification of defaulters.Our proposed framework classifies defaulters at five different stages:defaulter,partially high,partially medium,partially low,and unvaccinated to reinforce targeted interventions by accurately predicting children at high risk of defaulting from the immunization schedule.Different machine learning algorithms are applied on Pakistan Demographic and Health Survey(2017–18)dataset.Multilayer Perceptron yielded 98.5% accuracy for correctly identifying children who are likely to default from immunization series at different risk stages of being defaulter.In this paper,the proposed defaulters’prediction framework is a step forward towards a data-driven approach and provides a set of machine learning techniques to take advantage of predictive analytics.Hence,predictive analytics can reinforce immunization programs by expediting targeted action to reduce dropouts.Specially,the accurate predictions support targeted messages sent to at-risk parents’ and caretakers’ consumer devices(e.g.,smartphones)to maximize healthcare outcomes.

    Keywords:Smart healthcare;routine immunization;predictive analytics;defaulters;vaccination;machine learning;targeted messaging

    1 Introduction

    World Health Organization(WHO)started a global vaccination program,Expanded Program on Immunization(EPI)in 1974.The aim of the program was to protect children from ten infectious and fatal childhood diseases such as Diphtheria,Pertussis,Tetanus,Polio,Tuberculosis,Measles,Haemophilus Influenza Meningitis,Yellow fever,Hepatitis B and Pneumococcal,by making the corresponding vaccines available to all.WHO has indicated that around two to three million lives are savedper annumdue to vaccination.Tab.1 provides a schedule of immunization recommended by WHO to prevent childhood diseases in early stages of life prior to exposure of deadly diseases.Reducing the child morbidity and mortality is the major objective of immunization services.In 1978,Pakistan initiated its EPI with similar motive to vaccinate children for fatal diseases and improve the immunization facilities.PDHS conducted in 2017–18 showed that high child mortality rate is due to low immunization coverage,i.e.,only 66% across the country.In year 2017,74 childhood deaths out of per 1000 live births were reported.Government of Pakistan and its partners are putting substantial efforts to achieve the projected benchmark rate of immunization as there were multiple reasons for incomplete vaccinations among children in Pakistan as reported in Riaz et al.[1].However,it is unrealistic to assume that just vaccinating every child will decrease the mortality rate;as this process has multidimensional structure along with the contribution of several other complex factors including keeping a record of all children and all vaccines.Within Pakistan,EPI has faced multiple challenges such as parental and female awareness,education and inefficient health worker attitudes[2].Pakistan is trying to catch up with the national targets for immunizations.It is indeed a concern for the public health policy makers of Pakistan as it has its global impact.

    Table 1:Routine immunization schedule

    The focus of previous studies was on the identification of hurdles and provocations in the process.Most of the researchers centered their work on issues like maternal illiteracy,inadequate data,lack of parental knowledge and awareness program of vaccination,erroneous reporting,distance from vaccination centers etc.Unfortunately,the available data is not utilized to analyze and to discover trend in data for forecasting precise future outcomes.Considering the gaps found in literature,utilization of data for analysis and pattern discovery to improve the immunization coverage is kept as the major focus of our study.Health care domain can significantly benefit from the combined power of Machine learning and Big data where predictive analytics technology can play an important role in forecasting immunization outcomes as well.

    Considering these factors,Chandir et al.[3]focused on developing a model for identifying defaulters,children who are more likely to miss a scheduled vaccine,in order to enhance the immunization coverage,by utilizing the strength of predictive analytics.However,their developed model did not consider those children who missed one or two vaccine doses in immunization schedule,and predicted defaulters and non-defaulters only.As a result,children missing few doses would be categorized at high risk of defaulting as well.

    Motivated by aforementioned limitation,we propose a model in which we categorize immunization series into five stages:Not-vaccinated,Fully-immunized,Partially-High,Partially-medium,and Partially-Low on the basis of immunization doses.These categories were beneficial in recognize recognizing the level of risk associated with defaulting,making it possible to efficiently target children with limited resources and time.Hence,to predict the immunization stage of a child,various machine learning algorithms were used by this study,including Decision Tree,Support Vector Machine,Multilayer Perceptron,and Na?ve Bayes.98% accuracy was achieved by using Multilayer Perceptron.PDHS(2017–18)was used to extract children’s vaccination data for the purpose of experimentation.Overall,predictive analytics can be beneficial for the population having high rate of incomplete or partial immunization by making targeted efforts to reach defaulters on time.The accurate predictions in turn support smart healthcare[4–6]by sending targeted alerts to at-risk children’s parents and caretakers via their consumer devices,e.g.,smartphones.

    2 Related Work

    It is indicated by several studies that coverage is not the only factor to be considered for EPI but it constitutes of areas such as operational cost and stock maintenance of vaccine,value of data,electronic registries,health information and recommendation systems to monitor vaccinators.Literature discusses the proposed theories and rationales to develop the understanding of these areas and their influence on the performance of immunization program.In the next few paragraphs,we divided the literature with respect to the use of technology and analytical techniques used so far on vaccination data.

    Socio-economic factors have been the focus of many researcher’s study to determine the factors of low vaccine completion coverage.Owais et al.[7]conducted a study on a low-income community in Karachi to evaluate the impact of educational intervention,for mothers with low literacy rate,on their child.It was evident from the results that the designed educational intervention for less literate mothers elevated the immunization completion rate of DPT-3/Hepatitis B vaccine by 39%,thus proving to have vital implications for achieving better immunization ratios.A similar work by Onsomu et al.[8]focused on the importance of maternal education’s role in a child immunization completion process in Kenya that remains an under vaccinated country.It was found that the likelihood of immunization completion rate of women having above primary level education is much higher.Even after considering the individual and community level aspects in model building,it was shown that maternal literacy level is statistically imperative for a child’s well-being.According to our knowledge,maternal education has been never considered as a feature to predict child immunization rate in any of the prior researches.

    Apart from literacy rate,several other socio-economic attributes play a vital role in immunization coverage.Rammohan et al.[9]considered health facilities in India and other contributing socioeconomic factors in addition to education level to assess the childhood immunization rate at district level.Results showed that an individual’s income is strongly associated with the coverage rate of immunization and makes an imperative contribution in prediction as well.It was highlighted by their study that maternal education holds a statistically substantial importance and has a positive impact on immunization outcomes.

    Crouch et al.[10]highlighted the importance of demographic and socioeconomic features as predictors of immunization.Alongside these features,data collected through immunization survey at national level was also considered for model building.Model was built using features like mother’s education level and age,culture,race,first or later born child and their geographic location being one of the most significant contributors.Demographic and socioeconomic elements were found to be the major predictors of immunization and it was suggested that if these characteristics are used in a right manner then it could benefit health officials and policy makers to devise policies which will result in improved outcomes.On the other hand,Ameen et al.[11]conducted a similar study by considering socio-economic features and immunization status as variables from a survey conducted in 2011.It was concluded that the variation in immunization coverage was mainly based on child’s birthplace,maternal literacy level and availability of immunization card.

    Moreover,several researchers have acknowledged the significance of timely vaccinating the child and emphasized on its importance to prevent disabilities and fatal diseases in children.Qazi et al.[12]highlighted that delayed vaccinated children in KPK,Pakistan were not included in current vaccination reports.Due to this,if a child gets vaccinated after the specified time,he/she will be considered as fully immunized individual in the national level report.Hence,it tends to overestimate the vaccination coverage results.Problems related to non-vaccinated or partially vaccinated children are overlooked and rarely examined.However,Walton et al.[13]realized the significance of timeliness of vaccination and measured the quality of immunization services through routine data.It was highlighted through this study that delayed vaccinations highly affected the immunity of children leaving them at a high risk of catching infectious diseases.Moreover,there were no studies found focusing on the risk associated with partial or no vaccination affecting the whole immunization program.Moreover,none of the previous studies have utilized this information to improve the immunization coverage or to prevent future delays.

    The proposed approach filled in a major gap of predictive analytics by identifying early and on time future outcomes.Chandir et al.[3]took the first step in this area to recognize the usage of predictive analytics and its capacity to enhance immunization coverage by assessing a limited scope of this approach.

    Nevertheless,the previous studies have several limitations such as model was able to predict the potential candidates to miss their upcoming vaccination.Such individuals were termed as defaulters.However,even though a child misses only one dose,the model puts that child at an elevated risk of defaulting.

    This motivated us to build a prediction model for accurate predictions of defaulters.This paper addresses the power of predictive analytics for decision making,targeted interventions and increasing immunization coverage,so far lacking in the literature.Accurate prediction is at the center of a consumer deviceenabled targeted messaging that can lead to desirable smart healthcare outcomes.In the next section,the overview of proposed framework is presented for accurate defaulter prediction.

    3 Proposed Framework for Defaulter Prediction

    Based upon the shortcomings found in previous work as discussed in previous section,we propose a framework to predict immunization defaulters.This section provides an overview of proposed framework,details of dataset and experimental setup of the study.

    3.1 Vaccine Defaulter Identification

    The development of this framework was motivated by the lack of accurate identification of defaulters in an EPI.Fig.1 represents our defaulter prediction framework.

    The main contribution of this data-driven framework is to accurately identify the defaulters at different levels of risk of dropping out based on their vaccination status and the number of doses they have taken so far and categorizes children into five stages:Fully-immunized,partially-high,partially-medium,partially-low and un-vaccinated.The framework consists of three main components:Data Preprocessing,Attribute Selection,Model Building and Evaluation.

    3.2 Pakistan Demographic and Health Survey

    Dataset was extracted from PDHS(2017–18).This Survey contains information about basic demographics and health indicators.It provides an overview of population,maternal and child health issues in the country.The intent of this survey is to help policymakers and program managers to design strategies to improve the health of the population.Originally,the dataset consisted of 50,486 data records and 1186 variables(including variables related to household,education,and poverty).

    Figure 1:Framework for prediction of vaccine defaulters

    Survey data contains demographic information of children such as parent’s age,region,residence,parental language,families’ wealth index,child’s vaccine information,i.e.,if a child is vaccinated or not.All attributes have multiple values except maternal age and number of children under 5 years.Tab.2 shows the multiple values of each attribute.

    Table 2:Variable values details

    3.3 Experimental Setup

    All experiments have been performed using Core i5 with 8 GB RAM.SPSS was used to extract the immunization and demographic related attributes from a dataset of Demographic and Health Survey.A number of machine learning algorithms were used for classification and association rule mining technique has been applied on the given dataset.

    4 Component of the Proposed Framework

    In this section,the components of the framework are explained in step by step manner.

    4.1 Data Preprocessing

    Data for the PDHS study was collected through conducting interviews from different people and through questionnaires.Since the survey data was collected manually,there are possibilities of anomalies and variation in data.In such scenario,data preprocessing becomes a vital step to prepare the data for an efficient and effective model building,as unclean data has a negative impact on model’s performance.The dataset contains numerous records with multiple missing values.The details of the main steps of the preprocessing process are given below,

    4.2 Removing Missing Values

    Dealing with the missing values is an important step of data preprocessing.Various missing values were removed or replaced including spaces,‘Null’,special character‘?’to attain good results.There are only two possible options in immunization i.e.,either a child vaccinated or not vaccinated for a specific disease.According to the nature data,it was not possible to replace the missing values with mean,min,max etc.Hence,the missing values were removed from the data and r the total number of instances were reduced from 50,486 to 7153 as shown in Fig.2.

    Figure 2:Removed missing values

    4.3 Data Decoding

    In this step,nominal values were generated using the numeric demographic values.Survey dataset contains the metadata to decode demographic variables’ numeric value.In order to gain better understanding of variables,information was extracted from the metadata and demographic variables were decoded to nominal values,leaving out mother’s age,the number of children under 5 in a household.

    4.4 Discretization of Data

    A subset of categorical values was created from the dataset of large number of numeric values.The process of discretization helps is data evaluation and management process.Mother’s age was present in the data as a categorical attribute ‘maternal age’ indicating the age of child’s mother.After discretization,the values of the ‘respondent age’ attribute are distributed in three groups namely Young,Middle-aged and Old.15 to 24 were grouped in Young category,25–34 as Middle-aged and the rest as Old.

    4.5 Attribute Selection

    Immunization data contains 1186 attributes.All these attributes do not hold the same importance in terms for contribution to the intended outcome,which makes it vital to perform dimensionality reduction on dataset and to exclude less important or irrelevant attributes from data.Selected features include child’s immunization status such as getting vaccinated according to the schedule,child’s basic demographics including his residence,district,maternal age,wealth index and number of children under 5 years of age in a household.We selected 19 features based on domain knowledge as shown in Fig.3.

    Figure 3:Attribute selection using domain knowledge

    4.6 Model Building

    Once the data is completely preprocessed,prediction model is built using various classifiers.In first step,training and testing set is generated by splitting the dataset into two parts.The training set constitutes of 70%of the data and rest of the 30%is taken as the testing set.To train the model,only training set is used,whereas,to measure the performance of the model on unseen data,testing set is utilized.For classification,Artificial Neural Network,Na?ve Bayes,Decision Tree and Support Vector Machine are applied to the training set.Tab.3 shows the sample dataset after data preprocessing and attribute selection steps.

    4.7 Class Label Assignment

    Before passing training set to model building,training data was labeled into five categories by classifying the vaccination status of children.Children who have not taken even a single vaccination dose were labeled as unvaccinated.Children who have taken 5 or less than 5 doses out of total 13 were labeled as partially-low and if a child has taken 6 out of 10 doses,he/she is labeled as partially-medium.In case a child has taken 11 or more than 11 doses out of 13 was labeled as partially-high.Another label assigned was fully immunized where children has completed the vaccination series by taking all the doses of vaccination schedule.In the next step,all these labels are then passed to the model to perform further processing.We applied four different classifiers(J48,SVM,Na?ve Bayes,and Multilayer Perceptron)using Weka and noted down accuracy,TP rate,FP rate,Precision,Recall and ROC area for each classifier.Multilayer perceptron performed best and identified defaulters with an accuracy of 98.50%.In the next few lines,the reason behind our choice of these classifiers has been explained.

    4.8 Multilayer Perceptron

    Artificial Intelligence is the heart of new technologies with the capability to deliver optimal and costeffective healthcare solutions in real-time.With the increasing amount of healthcare data and the increasing power of analytical techniques,it gave a paradigm shift to healthcare.We applied this classifier on our dataset because it has the ability to detect complex non-linear relationships between dependent and non-dependent variables.This is really important because many of the input-output interactions are non-linear and dynamic in real-life.

    Table 3:Sample dataset after data preprocessing

    4.9 Decision Tree

    Decision Tree helps in understanding the data relationship in a tree hierarchy structure.This model is simple,explainable and useful with the simple dataset and fewer features and generates better predictions.It also requires less effort for data preparation during data pre-processing step.Survey data mostly contains missing values whereas the decision tree building process is less affected by missing values to a significant extent.The other main reason for using this algorithm is that visualizations of a complex model of the Decision Tree can be simplified.Even a naive person can make sense of logic.The number of hyperparameters to be tuned is nearly zero and also supports classification problems with more than two classes.

    4.10 Na?ve Bayes

    Na?ve Bayes is considered as one of the simple and easiest algorithms to predict outcomes.We applied Na?ve Bayes because predicting the class of the test data set is easy and quick using this algorithm.In multiclass prediction,it performs really well.It also performs better as compared to other models such as logistic regression when assuming independence holds,and less training data is needed whereas,it can also work easily with the missing values.

    4.11 Support Vector Machine

    The support vector machine is a versatile and powerful supervised machine learning algorithm.We applied the support vector machine because it can solve the issue of unbalanced data as a cost-sensitive classifier.Using a technique called the kernel trick to transform the results;it finds an optimal boundary between the possible outputs based on these transformations.

    5 Results and Discussion

    In the previous step,different machine learning algorithms have been applied to the training dataset and the performance of those algorithms has been measured using different performance metrics.According to four outcomes produced,Multilayer Perceptron performed best and yielded 98.5%accuracy.Decision Tree performed well after Multilayer Perceptron and yielded 98.1% accuracy.The support vector machine produced 94.68% accuracy and Na?ve Bayes yielded 94.5% accuracy.Fig.4 shows a comparison of the accuracy of each machine learning algorithm.

    Figure 4:Prediction accuracy of machine learning algorithms

    Among other performance metrics used in the literature,AUC is one of the most widely used metric to test the accuracy of the model.It is the plot between Sensitivity,i.e.,True Positive rate on y-axis and Specificity,i.e.,False positive rate on x-axis.Fig.5 shows the AUC of Multilayer Perceptron that gave better results.It tells us about the measure of separability between two classes.If the value is closer to 1,it means that the model is performing well and is able to classify 1 s as 1 s and 0 s as 0 s.Performance of the model decreases as the value move towards zero.Hence,a value closer to zero indicates poor model performance with less separability.

    Figure 5:AUC of multilayer perceptron

    However,Multilayer Perceptron outperforms the other algorithms in terms of accuracy,but accuracy alone is not the only metric to measure the performance of the model.Algorithms were analyzed on some other metrics like Precision,TP rate,and FP rate.Tab.4 shows the results of all performance metrics against each algorithm.

    Table 4:Performance metrics of each algorithm

    Out of 7153 children,only 1955(27.33%)children were fully-immunized.Tab.5 shows the increasing dropout rate with each successive dose with the highest dropout for Measles 2,i.e.,72.6%.The same table also shows the inter-vaccination defaulting rate ranging from low(6.6%)between Penta 2 and Penta 3 to high defaulting rate(15.7%)between Measles 1 and Measles 2.

    Immunization trend according to PDHS clearly shows that the rate of vaccination is high in early scheduled vaccines but it is decreasing with each subsequent vaccine in immunization visits which should not be the case.Using machine learning algorithms on demographic and immunization data,we can identify that either the child is at high risk or at low risk of missing their next immunization visit.This study aims not only to identify the child at risk of defaulting from next immunization visits but also the rate of risk as high,medium and low of defaulting.

    Table 5:Defaulting rate

    Several predictive analytics algorithms have been used in research to identify the children who are most likely to skip their upcoming immunization visit.Policymakers and officials can benefit greatly from predictive analytics by performing analysis on EPI dataset to decrease the frequency of dropout.The results of predictive analytics provide on-time information extracted from data to make targeted interventions and alerts possible via consumer devices,thereby improving smart healthcare outcomes.

    Despite paramount efforts to plan vaccine products that will come to get the requirements from endusers,the industry might need intuition into how product traits would influence immunization programs in low and middle-income countries.Similarly,priorities may be ignored by national-level decisionmakers.National decision-makers,especially from the GAVI supported countries,the vaccine alliance,must pay attention to the advantages and disadvantages of peculiar product aspects along with budget limitations while making up decisions related to which vaccines to propose.A timely vaccination reduce the death rate in children up to 2.5 millionper annumand has a great potential to save more lives as there is a rapid development in new vaccines to prevent the main causes of child mortality such as pneumonia,diarrhea,and malaria.Despite of this progress,several discrepancies persist in immunization coverage between several nations and the basic vaccine dosage is unavailable to many children worldwide.

    Government has been doing substantial efforts to enhance immunization coverage and reduce the rate of missed opportunities.Unfortunately,there is no tool available for health providers to timely identify children who are likely to miss next vaccination dose.With the passage of time,digital immunization registries are expanding and crafting a distinctive chance to generate real-time insights from vaccination data through exploration and translation,as a result,enhancing immunization coverage and other related services.

    To achieve complete vaccination coverage along with timeliness and age-appropriate vaccination,modeling and visualization of large datasets could be performed to identify the risk factors associated with dropout rate or delayed vaccination sessions.In this era of big data where massive data is being generated and collected,predictive analytics is utilized in several ways on reliable data to make the analysis more efficient.

    Nowadays,socio-economic factors,awareness and knowledge development have more significant impact on the decision making about vaccination than decades ago.Elevated income,improved education and knowledge level,positive attitude towards immunization,and health insurance are associated with the increase in immunization rates.Furthermore,immunization rate can be improved further through targeted policy efforts such as educational interventions for less-educated mothers.

    In addition to this,regular examining and assessment of district-level immunization coverage are important to identify local contextual factors Which could be better adapted to customized approaches such as a mix of health workers,more intensive education campaigns,socio-economic improvements and user-friendly health infrastructure provision.

    There has been no research conducted in the past to identify the associations between demographic and socioeconomic factors with the defaulters.It is important to understand the factors due to which children fail to complete their vaccination schedule or miss their vaccination at all.More education,higher incomes,added knowledge,and affirmative immunization attitudes,and health insurance are all linked to higher immunization rates.Moreover,Targeted policy efforts could help in the improvement of immunization rates such as targeted educational interventions for less-educated mothers.

    In a recent study[14],authors have checked the impact of languages spoken in Pakistan on the immunization of children.Interestingly,the results revealed the same finding that Balochi language speaking families are less likely to make their children immunized.So in this case,the decision-makers or health care provider can take steps to deliver message in their language or can design campaigns on vaccination in Baluchi language so that people who only understand Baluchi language could know about vaccination and its schedule and bring their children to complete his/her vaccination on time.

    It is important to discuss the role of Lady Health Workers(LHWs)who are generally recruited from the local community to disseminate the information related to maternal and child health.Although,the LHWs have the ability to communicate with the families using the same language,still in majority of the Pakistani community mothers are not taken as a head of the family and they rely heavily on the male head of the family regarding the decisions related to child health.On the other hand,a typical LHW in Pakistan has to support not only EPI program but there are a large number of tasks associated with respect to different healthcare programs and it is very hard for a LHW to cover all the 22 different tasks as reported in Hafeez et al.[15].

    An international target of immunization coverage of 90% and 80% at national and district level respectively is set by WHO.To meet this international criterion and to upsurge immunization coverage nation-wide,Pakistan require effectual and focused strategies.Predictive analytics,as presented in this paper,can help in identification of children who are going to miss their next visit so that those areas where the ratio of incomplete or partial immunization is high can get benefits of predictive analytics to reduce the incomplete vaccination series and ultimately increasing rate of complete immunization coverage.

    6 Conclusion

    Lack of vaccination coverage in Pakistan is due to several reasons.Numerous numbers of children are not fully vaccinated,late for their age-appropriate vaccine dose,or withdraw from their vaccination schedule.In this paper,the defaulters’ prediction model has been proposed and several classifiers were applied on PDHS dataset including Decision Tree,Support Vector Machine,Na?ve Bayes and Multilayer Perceptron.98% accuracy for identification of children at different stages of risk of dropping out was achieved through proposed model.Our framework is advancement towards data-driven approach that also benefits from predictive analytics by providing a set of machine learning techniques.Therefore,if immunization programs are powered by predictive analytics,smart health outcomes could be improved by expediting targeted action via consumer devices.In the field of immunization,predictive modeling is a novel idea,and its potential in revolutionizing immunization service delivery is yet to be identified.

    Acknowledgement:This research was conducted by Predictive Analytics(PA)Lab SZABIST Islamabad,established under the National Center in Big Data and Cloud Computing(NCBC)Pakistan.PA Lab aims to assist in the healthcare domain of government institutions by providing predictive analytics tools and solutions for(NTP,EPI,LHWP,NSP,MNCH)programs.

    Funding Statement:This research was supported by Korea Institute for Advancement of Technology(KIAT)grant funded by the Korea Government(MOTIE)(P0012724,The Competency Development Program for Industry Specialist)and the Soonchunhyang University Research Fund.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    成年人黄色毛片网站| 成人国产一区最新在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲精华国产精华精| 操美女的视频在线观看| netflix在线观看网站| 久久久精品区二区三区| 亚洲精品久久成人aⅴ小说| 国产精品.久久久| 1024视频免费在线观看| 亚洲精品国产一区二区精华液| 女人爽到高潮嗷嗷叫在线视频| 侵犯人妻中文字幕一二三四区| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成77777在线视频| 久久精品亚洲av国产电影网| 欧美久久黑人一区二区| 一二三四在线观看免费中文在| 免费黄频网站在线观看国产| 日韩欧美免费精品| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品一区二区www | 少妇精品久久久久久久| 欧美激情高清一区二区三区| 少妇的丰满在线观看| 精品高清国产在线一区| 亚洲国产成人一精品久久久| 国产一区二区三区综合在线观看| 中文字幕人妻熟女乱码| 亚洲精品国产色婷婷电影| 精品第一国产精品| 窝窝影院91人妻| 超碰97精品在线观看| 国产一卡二卡三卡精品| 中文欧美无线码| 欧美精品av麻豆av| 国产亚洲av高清不卡| 久热这里只有精品99| 高清黄色对白视频在线免费看| 老司机午夜福利在线观看视频 | 国产一区二区三区综合在线观看| 18禁黄网站禁片午夜丰满| xxxhd国产人妻xxx| 丝袜在线中文字幕| 999久久久国产精品视频| 99精品欧美一区二区三区四区| a在线观看视频网站| 高清在线国产一区| 侵犯人妻中文字幕一二三四区| 国产av一区二区精品久久| 亚洲av欧美aⅴ国产| 国产激情久久老熟女| 欧美精品啪啪一区二区三区| 两个人免费观看高清视频| 欧美在线一区亚洲| 欧美精品av麻豆av| 亚洲国产欧美日韩在线播放| 亚洲av电影在线进入| 99久久国产精品久久久| 亚洲专区字幕在线| 精品一品国产午夜福利视频| 黄片大片在线免费观看| 国产1区2区3区精品| 每晚都被弄得嗷嗷叫到高潮| 岛国毛片在线播放| 国产精品美女特级片免费视频播放器 | 亚洲一区中文字幕在线| 一个人免费看片子| 黄色视频,在线免费观看| 亚洲黑人精品在线| 最近最新中文字幕大全电影3 | 女人高潮潮喷娇喘18禁视频| 少妇裸体淫交视频免费看高清 | 免费观看人在逋| 色婷婷久久久亚洲欧美| 国产精品久久久人人做人人爽| 久久免费观看电影| 在线播放国产精品三级| 9热在线视频观看99| 日韩 欧美 亚洲 中文字幕| 国产成人精品久久二区二区91| 国产亚洲欧美精品永久| 变态另类成人亚洲欧美熟女 | 少妇猛男粗大的猛烈进出视频| 久久婷婷成人综合色麻豆| 啦啦啦 在线观看视频| 日韩有码中文字幕| 高清欧美精品videossex| 国产成人系列免费观看| 这个男人来自地球电影免费观看| 国产在线精品亚洲第一网站| 欧美+亚洲+日韩+国产| 国产亚洲欧美在线一区二区| 十分钟在线观看高清视频www| 两性午夜刺激爽爽歪歪视频在线观看 | 五月天丁香电影| 少妇精品久久久久久久| 欧美精品一区二区大全| 青青草视频在线视频观看| 少妇粗大呻吟视频| 精品久久久久久久毛片微露脸| 夜夜夜夜夜久久久久| 又黄又粗又硬又大视频| 宅男免费午夜| 日本wwww免费看| 咕卡用的链子| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 欧美成人免费av一区二区三区 | 亚洲国产成人一精品久久久| 性色av乱码一区二区三区2| 一级片'在线观看视频| 午夜激情av网站| 国产欧美日韩一区二区三| 久久九九热精品免费| 啦啦啦在线免费观看视频4| 日韩欧美国产一区二区入口| 亚洲国产精品一区二区三区在线| 91老司机精品| 男人舔女人的私密视频| cao死你这个sao货| 久久人人爽av亚洲精品天堂| 涩涩av久久男人的天堂| 日韩大码丰满熟妇| 黄片播放在线免费| 女人高潮潮喷娇喘18禁视频| 交换朋友夫妻互换小说| 男女高潮啪啪啪动态图| 亚洲全国av大片| 黑丝袜美女国产一区| 精品午夜福利视频在线观看一区 | 国产免费福利视频在线观看| 久久人妻av系列| 纯流量卡能插随身wifi吗| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级| www.自偷自拍.com| 久久天躁狠狠躁夜夜2o2o| 人人澡人人妻人| 大型黄色视频在线免费观看| 色精品久久人妻99蜜桃| 一区二区三区国产精品乱码| 一区二区日韩欧美中文字幕| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看| 亚洲中文av在线| 狂野欧美激情性xxxx| av免费在线观看网站| 精品乱码久久久久久99久播| 桃花免费在线播放| 成人特级黄色片久久久久久久 | 色视频在线一区二区三区| 91字幕亚洲| 亚洲七黄色美女视频| 欧美国产精品一级二级三级| 亚洲精品在线观看二区| 男女边摸边吃奶| 高潮久久久久久久久久久不卡| 国产精品影院久久| 免费看十八禁软件| 香蕉久久夜色| 人妻一区二区av| 一区二区av电影网| 亚洲专区国产一区二区| 日韩视频一区二区在线观看| 色尼玛亚洲综合影院| 老汉色∧v一级毛片| 国产精品自产拍在线观看55亚洲 | 下体分泌物呈黄色| 精品少妇黑人巨大在线播放| 亚洲成a人片在线一区二区| 国产精品久久久久久精品电影小说| 他把我摸到了高潮在线观看 | 久久国产精品男人的天堂亚洲| 欧美亚洲 丝袜 人妻 在线| 中文字幕av电影在线播放| 国产一区有黄有色的免费视频| 国产黄色免费在线视频| 99精品在免费线老司机午夜| 精品一品国产午夜福利视频| 免费高清在线观看日韩| 极品教师在线免费播放| 成人18禁在线播放| 国产日韩欧美亚洲二区| 欧美成狂野欧美在线观看| 99久久99久久久精品蜜桃| 精品一区二区三卡| 视频区图区小说| www日本在线高清视频| 18禁国产床啪视频网站| 亚洲人成伊人成综合网2020| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频| 黄色毛片三级朝国网站| 精品视频人人做人人爽| 精品一区二区三区四区五区乱码| 丰满饥渴人妻一区二区三| 欧美黑人欧美精品刺激| 男女边摸边吃奶| 午夜福利视频精品| 丁香六月天网| 交换朋友夫妻互换小说| 国产精品电影一区二区三区 | 亚洲av欧美aⅴ国产| av国产精品久久久久影院| 搡老岳熟女国产| 桃红色精品国产亚洲av| 热99国产精品久久久久久7| 一本一本久久a久久精品综合妖精| 日本撒尿小便嘘嘘汇集6| 又紧又爽又黄一区二区| 狠狠狠狠99中文字幕| 日韩免费高清中文字幕av| 亚洲美女黄片视频| 69精品国产乱码久久久| 少妇 在线观看| 美女高潮到喷水免费观看| 高清av免费在线| 丝袜在线中文字幕| 一级片免费观看大全| 看免费av毛片| av视频免费观看在线观看| 18禁黄网站禁片午夜丰满| 九色亚洲精品在线播放| 免费观看人在逋| 91老司机精品| 免费少妇av软件| 两性夫妻黄色片| 午夜福利影视在线免费观看| 日韩三级视频一区二区三区| 国产高清videossex| 精品亚洲乱码少妇综合久久| 国产熟女午夜一区二区三区| 国产精品99久久99久久久不卡| 国产午夜精品久久久久久| 亚洲精品乱久久久久久| 不卡一级毛片| 亚洲av欧美aⅴ国产| 热re99久久国产66热| kizo精华| 大码成人一级视频| 极品少妇高潮喷水抽搐| 99热国产这里只有精品6| 九色亚洲精品在线播放| 狠狠婷婷综合久久久久久88av| 天堂8中文在线网| 男女午夜视频在线观看| 久久午夜亚洲精品久久| 九色亚洲精品在线播放| 老司机影院毛片| 天堂8中文在线网| 日韩视频在线欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 国产视频一区二区在线看| 免费高清在线观看日韩| 啦啦啦中文免费视频观看日本| 国产一区有黄有色的免费视频| 99香蕉大伊视频| 水蜜桃什么品种好| 久久99一区二区三区| 亚洲中文日韩欧美视频| 91精品三级在线观看| 熟女少妇亚洲综合色aaa.| 99久久人妻综合| 欧美精品一区二区免费开放| 在线观看66精品国产| 国内毛片毛片毛片毛片毛片| 日本av手机在线免费观看| 国产高清激情床上av| 99久久精品国产亚洲精品| 久久久久久免费高清国产稀缺| 国产福利在线免费观看视频| 香蕉久久夜色| 国产一区二区三区在线臀色熟女 | 日本vs欧美在线观看视频| 亚洲熟女毛片儿| 满18在线观看网站| 老司机深夜福利视频在线观看| 亚洲精品粉嫩美女一区| 亚洲欧美日韩高清在线视频 | 国产精品免费大片| tube8黄色片| 最黄视频免费看| 黑人巨大精品欧美一区二区mp4| 狠狠婷婷综合久久久久久88av| 国精品久久久久久国模美| 1024香蕉在线观看| 黄色丝袜av网址大全| 少妇精品久久久久久久| 成人国产av品久久久| 久久久国产成人免费| 变态另类成人亚洲欧美熟女 | 精品一区二区三区视频在线观看免费 | 天堂俺去俺来也www色官网| 男女之事视频高清在线观看| 久久精品亚洲精品国产色婷小说| 精品国产一区二区三区久久久樱花| 精品少妇久久久久久888优播| 80岁老熟妇乱子伦牲交| 夫妻午夜视频| av欧美777| 久久久久视频综合| 免费黄频网站在线观看国产| 成人国产av品久久久| 午夜老司机福利片| 男人舔女人的私密视频| bbb黄色大片| 超碰97精品在线观看| 亚洲一区二区三区欧美精品| 黄色成人免费大全| 久久人妻熟女aⅴ| 黄色怎么调成土黄色| 老熟女久久久| 国产精品久久久久久人妻精品电影 | avwww免费| 十分钟在线观看高清视频www| 伊人久久大香线蕉亚洲五| 波多野结衣一区麻豆| 欧美成人午夜精品| 中国美女看黄片| www日本在线高清视频| 久久中文字幕人妻熟女| 精品久久久精品久久久| 久久国产精品大桥未久av| 啦啦啦视频在线资源免费观看| 天堂动漫精品| 啦啦啦中文免费视频观看日本| 日韩欧美三级三区| 91老司机精品| 十八禁网站网址无遮挡| 国产99久久九九免费精品| 1024香蕉在线观看| 亚洲伊人久久精品综合| 色婷婷久久久亚洲欧美| 俄罗斯特黄特色一大片| 女人爽到高潮嗷嗷叫在线视频| 色94色欧美一区二区| 国产av国产精品国产| 中文字幕人妻丝袜一区二区| 午夜免费鲁丝| 黄色视频不卡| 男男h啪啪无遮挡| 男女床上黄色一级片免费看| 免费在线观看日本一区| 捣出白浆h1v1| 国产免费视频播放在线视频| 一个人免费在线观看的高清视频| 精品少妇内射三级| 国产精品一区二区在线观看99| 黄色毛片三级朝国网站| 90打野战视频偷拍视频| netflix在线观看网站| 午夜福利免费观看在线| 搡老熟女国产l中国老女人| 女人久久www免费人成看片| 99精品在免费线老司机午夜| 五月开心婷婷网| 日日夜夜操网爽| 在线天堂中文资源库| 久久青草综合色| 十八禁网站免费在线| 黄色a级毛片大全视频| 亚洲男人天堂网一区| 亚洲欧洲日产国产| 日韩欧美国产一区二区入口| 亚洲色图 男人天堂 中文字幕| 欧美成人午夜精品| 亚洲av日韩在线播放| 午夜免费鲁丝| 99在线人妻在线中文字幕 | 国产成人欧美| bbb黄色大片| 久久人人97超碰香蕉20202| 日韩有码中文字幕| av天堂在线播放| 99久久精品国产亚洲精品| 极品少妇高潮喷水抽搐| 成年女人毛片免费观看观看9 | 亚洲第一青青草原| 久久久久国产一级毛片高清牌| 好男人电影高清在线观看| bbb黄色大片| 日韩欧美国产一区二区入口| 十分钟在线观看高清视频www| 亚洲美女黄片视频| 大片电影免费在线观看免费| 人妻 亚洲 视频| 国产精品久久电影中文字幕 | 一二三四社区在线视频社区8| 黄色怎么调成土黄色| 日本黄色日本黄色录像| 色综合婷婷激情| 国产一区二区三区在线臀色熟女 | 日日摸夜夜添夜夜添小说| 女性生殖器流出的白浆| av电影中文网址| 中文字幕人妻丝袜制服| 成年版毛片免费区| 日韩大码丰满熟妇| 亚洲国产欧美日韩在线播放| 亚洲国产av新网站| 乱人伦中国视频| 操美女的视频在线观看| 国产真人三级小视频在线观看| 美女扒开内裤让男人捅视频| 国产精品免费一区二区三区在线 | 丝袜美足系列| 国产伦人伦偷精品视频| 成年版毛片免费区| 国产精品偷伦视频观看了| 老汉色∧v一级毛片| 精品国产乱码久久久久久小说| 日韩欧美国产一区二区入口| 亚洲第一欧美日韩一区二区三区 | 首页视频小说图片口味搜索| 男女无遮挡免费网站观看| 国产男女内射视频| 日韩欧美三级三区| 美女福利国产在线| 免费观看人在逋| 一边摸一边抽搐一进一出视频| 夜夜爽天天搞| 亚洲精品在线美女| 热99re8久久精品国产| videos熟女内射| 欧美一级毛片孕妇| 1024香蕉在线观看| 国产成人免费无遮挡视频| 国产成人精品在线电影| 在线十欧美十亚洲十日本专区| 亚洲国产欧美一区二区综合| 国产一区二区三区视频了| 国产欧美日韩一区二区三区在线| 无人区码免费观看不卡 | 黄色丝袜av网址大全| 美女高潮喷水抽搐中文字幕| 男男h啪啪无遮挡| 久久人妻福利社区极品人妻图片| 水蜜桃什么品种好| 成人黄色视频免费在线看| 久久久久久久国产电影| 女人精品久久久久毛片| 日韩欧美一区视频在线观看| 色视频在线一区二区三区| 性高湖久久久久久久久免费观看| 又黄又粗又硬又大视频| 视频在线观看一区二区三区| 老熟妇仑乱视频hdxx| 性少妇av在线| 国产又色又爽无遮挡免费看| 精品第一国产精品| 国产成人精品在线电影| 色老头精品视频在线观看| 汤姆久久久久久久影院中文字幕| 男人舔女人的私密视频| 婷婷丁香在线五月| 久久中文看片网| 国产成人免费无遮挡视频| 日韩欧美免费精品| 亚洲五月色婷婷综合| 国产熟女午夜一区二区三区| 91成人精品电影| 欧美日韩亚洲国产一区二区在线观看 | 黑人巨大精品欧美一区二区蜜桃| 777米奇影视久久| 他把我摸到了高潮在线观看 | 亚洲七黄色美女视频| 亚洲精品成人av观看孕妇| 别揉我奶头~嗯~啊~动态视频| av天堂在线播放| 国产亚洲欧美精品永久| 亚洲成av片中文字幕在线观看| 亚洲国产欧美网| 最近最新中文字幕大全免费视频| 日本黄色视频三级网站网址 | 在线观看舔阴道视频| 交换朋友夫妻互换小说| 黄色 视频免费看| 国产97色在线日韩免费| 久久天堂一区二区三区四区| 成人精品一区二区免费| 欧美精品一区二区大全| 每晚都被弄得嗷嗷叫到高潮| 日韩中文字幕视频在线看片| 国产精品免费一区二区三区在线 | 欧美日韩av久久| 成人手机av| 韩国精品一区二区三区| 亚洲avbb在线观看| 99在线人妻在线中文字幕 | 欧美在线黄色| 黄色怎么调成土黄色| 俄罗斯特黄特色一大片| 在线观看免费视频日本深夜| 精品久久久久久电影网| 80岁老熟妇乱子伦牲交| 午夜福利视频精品| av网站免费在线观看视频| 老司机靠b影院| 日韩大片免费观看网站| 久久中文看片网| 精品卡一卡二卡四卡免费| 天天操日日干夜夜撸| 999精品在线视频| 一本—道久久a久久精品蜜桃钙片| 久久久精品免费免费高清| 国产av国产精品国产| 日本a在线网址| 亚洲三区欧美一区| 国产成人免费观看mmmm| 欧美老熟妇乱子伦牲交| 9191精品国产免费久久| 国产黄色免费在线视频| 国产在线一区二区三区精| svipshipincom国产片| 大香蕉久久网| 老司机亚洲免费影院| 国产精品偷伦视频观看了| 国产日韩欧美亚洲二区| 91大片在线观看| 大陆偷拍与自拍| 黑人巨大精品欧美一区二区mp4| 99久久国产精品久久久| 考比视频在线观看| 午夜日韩欧美国产| 视频在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 久久午夜综合久久蜜桃| 国产精品av久久久久免费| 久久国产精品影院| aaaaa片日本免费| 久热这里只有精品99| 国产1区2区3区精品| 亚洲黑人精品在线| 午夜精品国产一区二区电影| 国产亚洲欧美精品永久| 国产三级黄色录像| 成年版毛片免费区| 国产成人av教育| 人妻 亚洲 视频| 精品少妇内射三级| 最新的欧美精品一区二区| 97在线人人人人妻| 久久热在线av| 狠狠婷婷综合久久久久久88av| 精品国产亚洲在线| 精品欧美一区二区三区在线| 在线天堂中文资源库| 激情在线观看视频在线高清 | 色播在线永久视频| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产一区二区入口| 国产高清国产精品国产三级| 国产有黄有色有爽视频| 蜜桃国产av成人99| 欧美另类亚洲清纯唯美| 欧美激情极品国产一区二区三区| 久久久久网色| 国产成人av教育| 一进一出好大好爽视频| 国产在视频线精品| 成人黄色视频免费在线看| 亚洲国产看品久久| 日日爽夜夜爽网站| 久久天躁狠狠躁夜夜2o2o| 亚洲伊人色综图| 国产日韩一区二区三区精品不卡| 免费观看av网站的网址| 激情视频va一区二区三区| 久久久久国产一级毛片高清牌| 19禁男女啪啪无遮挡网站| 国产精品一区二区精品视频观看| 亚洲国产av新网站| 精品乱码久久久久久99久播| av不卡在线播放| 最黄视频免费看| 午夜福利一区二区在线看| 久久精品国产综合久久久| 中文欧美无线码| 国产又色又爽无遮挡免费看| 中文字幕高清在线视频| 视频在线观看一区二区三区| 免费黄频网站在线观看国产| 精品少妇久久久久久888优播| 国产精品久久久久久人妻精品电影 | a在线观看视频网站| 亚洲国产毛片av蜜桃av| 中文字幕av电影在线播放| 视频在线观看一区二区三区| 日韩成人在线观看一区二区三区| 99精品在免费线老司机午夜| 中国美女看黄片| 国产一区有黄有色的免费视频| 啦啦啦在线免费观看视频4| 麻豆av在线久日| 欧美av亚洲av综合av国产av| 啦啦啦中文免费视频观看日本| 五月开心婷婷网| 国产精品1区2区在线观看. | 国产片内射在线| 中文字幕另类日韩欧美亚洲嫩草| 免费在线观看影片大全网站| 丰满人妻熟妇乱又伦精品不卡| 中文字幕另类日韩欧美亚洲嫩草| 老熟妇仑乱视频hdxx| 国产高清videossex| 性色av乱码一区二区三区2| 亚洲专区中文字幕在线| 搡老乐熟女国产| 天天添夜夜摸| 女同久久另类99精品国产91| 欧美一级毛片孕妇| 一个人免费在线观看的高清视频| 国产无遮挡羞羞视频在线观看| 夫妻午夜视频|