\documentclass[11pt,twoside]{article}\makeatletter \IfFileExists{xcolor.sty}% {\RequirePackage{xcolor}}% {\RequirePackage{color}} \usepackage{colortbl} \usepackage{wrapfig} \usepackage{ifxetex} \ifxetex \usepackage{fontspec} \usepackage{xunicode} \catcode`⃥=\active \def⃥{\textbackslash} \catcode`❴=\active \def❴{\{} \catcode`❵=\active \def❵{\}} \def\textJapanese{\fontspec{Noto Sans CJK JP}} \def\textChinese{\fontspec{Noto Sans CJK SC}} \def\textKorean{\fontspec{Noto Sans CJK KR}} \setmonofont{DejaVu Sans Mono} \else \IfFileExists{utf8x.def}% {\usepackage[utf8x]{inputenc} \PrerenderUnicode{–} }% {\usepackage[utf8]{inputenc}} \usepackage[english]{babel} \usepackage[T1]{fontenc} \usepackage{float} \usepackage[]{ucs} \uc@dclc{8421}{default}{\textbackslash } \uc@dclc{10100}{default}{\{} \uc@dclc{10101}{default}{\}} \uc@dclc{8491}{default}{\AA{}} \uc@dclc{8239}{default}{\,} \uc@dclc{20154}{default}{ } \uc@dclc{10148}{default}{>} \def\textschwa{\rotatebox{-90}{e}} \def\textJapanese{} \def\textChinese{} \IfFileExists{tipa.sty}{\usepackage{tipa}}{} \fi \def\exampleFont{\ttfamily\small} \DeclareTextSymbol{\textpi}{OML}{25} \usepackage{relsize} \RequirePackage{array} \def\@testpach{\@chclass \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else \ifnum \@lastchclass=7 5 \else \ifnum \@lastchclass=8 \tw@ \else \ifnum \@lastchclass=9 \thr@@ \else \z@ \ifnum \@lastchclass = 10 \else \edef\@nextchar{\expandafter\string\@nextchar}% \@chnum \if \@nextchar c\z@ \else \if \@nextchar l\@ne \else \if \@nextchar r\tw@ \else \z@ \@chclass \if\@nextchar |\@ne \else \if \@nextchar !6 \else \if \@nextchar @7 \else \if \@nextchar (8 \else \if \@nextchar )9 \else 10 \@chnum \if \@nextchar m\thr@@\else \if \@nextchar p4 \else \if \@nextchar b5 \else \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi \fi} \gdef\arraybackslash{\let\\=\@arraycr} \def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}} \def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}} \def\abbr{} \def\corr{} \def\expan{} \def\gap{} \def\orig{} \def\reg{} \def\ref{} \def\sic{} \def\persName{}\def\name{} \def\placeName{} \def\orgName{} \def\textcal#1{{\fontspec{Lucida Calligraphy}#1}} \def\textgothic#1{{\fontspec{Lucida Blackletter}#1}} \def\textlarge#1{{\large #1}} \def\textoverbar#1{\ensuremath{\overline{#1}}} \def\textquoted#1{‘#1’} \def\textsmall#1{{\small #1}} \def\textsubscript#1{\@textsubscript{\selectfont#1}} \def\textxi{\ensuremath{\xi}} \def\titlem{\itshape} \newenvironment{biblfree}{}{\ifvmode\par\fi } \newenvironment{bibl}{}{} \newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par } \newenvironment{citbibl}{}{\ifvmode\par\fi } \newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi } \newenvironment{docDate}{}{\ifvmode\par\fi } \newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi } \newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par } \newenvironment{msHead}{\vskip 6pt}{\par} \newenvironment{msItem}{\vskip 6pt}{\par} \newenvironment{rubric}{}{} \newenvironment{titlePart}{}{\par } \newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}} \newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}} \newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}} \newcolumntype{P}[1]{){\arraybackslash}p{#1}} \newcolumntype{B}[1]{){\arraybackslash}b{#1}} \newcolumntype{M}[1]{){\arraybackslash}m{#1}} \definecolor{label}{gray}{0.75} \def\unusedattribute#1{\sout{\textcolor{label}{#1}}} \DeclareRobustCommand*{\xref}{\hyper@normalise\xref@} \def\xref@#1#2{\hyper@linkurl{#2}{#1}} \begingroup \catcode`\_=\active \gdef_#1{\ensuremath{\sb{\mathrm{#1}}}} \endgroup \mathcode`\_=\string"8000 \catcode`\_=12\relax \usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry} \usepackage{framed} \definecolor{shadecolor}{gray}{0.95} \usepackage{longtable} \usepackage[normalem]{ulem} \usepackage{fancyvrb} \usepackage{fancyhdr} \usepackage{graphicx} \usepackage{marginnote} \renewcommand{\@cite}[1]{#1} \renewcommand*{\marginfont}{\itshape\footnotesize} \def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif} \pagestyle{fancy} \usepackage[pdftitle={Modeling Retention Indices of a Series Components Food and Pollutants of the Environment: Methods; OLS, LAD}, pdfauthor={}]{hyperref} \hyperbaseurl{} \paperwidth210mm \paperheight297mm \def\@pnumwidth{1.55em} \def\@tocrmarg {2.55em} \def\@dotsep{4.5} \setcounter{tocdepth}{3} \clubpenalty=8000 \emergencystretch 3em \hbadness=4000 \hyphenpenalty=400 \pretolerance=750 \tolerance=2000 \vbadness=4000 \widowpenalty=10000 \renewcommand\section{\@startsection {section}{1}{\z@}% {-1.75ex \@plus -0.5ex \@minus -.2ex}% {0.5ex \@plus .2ex}% {\reset@font\Large\bfseries}} \renewcommand\subsection{\@startsection{subsection}{2}{\z@}% {-1.75ex\@plus -0.5ex \@minus- .2ex}% {0.5ex \@plus .2ex}% {\reset@font\Large}} \renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}% {-1.5ex\@plus -0.35ex \@minus -.2ex}% {0.5ex \@plus .2ex}% {\reset@font\large}} \renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}% {-1ex \@plus-0.35ex \@minus -0.2ex}% {0.5ex \@plus .2ex}% {\reset@font\normalsize}} \renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}% {1.5ex \@plus1ex \@minus .2ex}% {-1em}% {\reset@font\normalsize\bfseries}} \def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt} \@tempdima 1.5em \begingroup \parindent \z@ \rightskip \@pnumwidth \parfillskip -\@pnumwidth \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par \endgroup} \def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}} \def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}} \def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}} \def\l@subparagraph{\@dottedtocline{5}{10em}{5em}} \@ifundefined{c@section}{\newcounter{section}}{} \@ifundefined{c@chapter}{\newcounter{chapter}}{} \newif\if@mainmatter \@mainmattertrue \def\chaptername{Chapter} \def\frontmatter{% \pagenumbering{roman} \def\thechapter{\@roman\c@chapter} \def\theHchapter{\roman{chapter}} \def\thesection{\@roman\c@section} \def\theHsection{\roman{section}} \def\@chapapp{}% } \def\mainmatter{% \cleardoublepage \def\thechapter{\@arabic\c@chapter} \setcounter{chapter}{0} \setcounter{section}{0} \pagenumbering{arabic} \setcounter{secnumdepth}{6} \def\@chapapp{\chaptername}% \def\theHchapter{\arabic{chapter}} \def\thesection{\@arabic\c@section} \def\theHsection{\arabic{section}} } \def\backmatter{% \cleardoublepage \setcounter{chapter}{0} \setcounter{section}{0} \setcounter{secnumdepth}{2} \def\@chapapp{\appendixname}% \def\thechapter{\@Alph\c@chapter} \def\theHchapter{\Alph{chapter}} \appendix } \newenvironment{bibitemlist}[1]{% \list{\@biblabel{\@arabic\c@enumiv}}% {\settowidth\labelwidth{\@biblabel{#1}}% \leftmargin\labelwidth \advance\leftmargin\labelsep \@openbib@code \usecounter{enumiv}% \let\p@enumiv\@empty \renewcommand\theenumiv{\@arabic\c@enumiv}% }% \sloppy \clubpenalty4000 \@clubpenalty \clubpenalty \widowpenalty4000% \sfcode`\.\@m}% {\def\@noitemerr {\@latex@warning{Empty `bibitemlist' environment}}% \endlist} \def\tableofcontents{\section*{\contentsname}\@starttoc{toc}} \parskip0pt \parindent1em \def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}} \newenvironment{reflist}{% \begin{raggedright}\begin{list}{} {% \setlength{\topsep}{0pt}% \setlength{\rightmargin}{0.25in}% \setlength{\itemsep}{0pt}% \setlength{\itemindent}{0pt}% \setlength{\parskip}{0pt}% \setlength{\parsep}{2pt}% \def\makelabel##1{\itshape ##1}}% } {\end{list}\end{raggedright}} \newenvironment{sansreflist}{% \begin{raggedright}\begin{list}{} {% \setlength{\topsep}{0pt}% \setlength{\rightmargin}{0.25in}% \setlength{\itemindent}{0pt}% \setlength{\parskip}{0pt}% \setlength{\itemsep}{0pt}% \setlength{\parsep}{2pt}% \def\makelabel##1{\upshape ##1}}% } {\end{list}\end{raggedright}} \newenvironment{specHead}[2]% {\vspace{20pt}\hrule\vspace{10pt}% \phantomsection\label{#1}\markright{#2}% \pdfbookmark[2]{#2}{#1}% \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}% }{} \def\TheFullDate{2016-01-15 (revised: 15 January 2016)} \def\TheID{\makeatother } \def\TheDate{2016-01-15} \title{Modeling Retention Indices of a Series Components Food and Pollutants of the Environment: Methods; OLS, LAD} \author{}\makeatletter \makeatletter \newcommand*{\cleartoleftpage}{% \clearpage \if@twoside \ifodd\c@page \hbox{}\newpage \if@twocolumn \hbox{}\newpage \fi \fi \fi } \makeatother \makeatletter \thispagestyle{empty} \markright{\@title}\markboth{\@title}{\@author} \renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@ \belowdisplayskip \abovedisplayskip \abovedisplayshortskip \z@ plus2\p@ \belowdisplayshortskip 4\p@ plus2\p@ minus2\p@ \def\@listi{\leftmargin\leftmargini \topsep 2\p@ plus1\p@ minus1\p@ \parsep 2\p@ plus\p@ minus\p@ \itemsep 1pt} } \makeatother \fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm} \fancyhf{} \setlength{\headheight}{14pt} \fancyhead[LE]{\bfseries\leftmark} \fancyhead[RO]{\bfseries\rightmark} \fancyfoot[RO]{} \fancyfoot[CO]{\thepage} \fancyfoot[LO]{\TheID} \fancyfoot[LE]{} \fancyfoot[CE]{\thepage} \fancyfoot[RE]{\TheID} \hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true} \fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}} \date{} \usepackage{authblk} \providecommand{\keywords}[1] { \footnotesize \textbf{\textit{Index terms---}} #1 } \usepackage{graphicx,xcolor} \definecolor{GJBlue}{HTML}{273B81} \definecolor{GJLightBlue}{HTML}{0A9DD9} \definecolor{GJMediumGrey}{HTML}{6D6E70} \definecolor{GJLightGrey}{HTML}{929497} \renewenvironment{abstract}{% \setlength{\parindent}{0pt}\raggedright \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}} \vskip16pt \textcolor{GJBlue}{\large\bfseries\abstractname\space} }{% \vskip8pt \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}} \vskip16pt } \usepackage[absolute,overlay]{textpos} \makeatother \usepackage{lineno} \linenumbers \begin{document} \author[1]{Fatiha Mebarki} \author[2]{Khadidja Amirat} \author[3]{Salima Ali Mokhnache} \author[4]{Djelloul Messadi} \author[5]{Khadidja Amirat} \affil[1]{ University Badji Mokhtar Annaba} \renewcommand\Authands{ and } \date{\small \em Received: 16 December 2015 Accepted: 1 January 2016 Published: 15 January 2016} \maketitle \begin{abstract} The gas chromatographic retention indices for (89 pyrazines of test and 25 of validation) on O V-101 and Carbowax -20M are successfuty modeled with the ald of a computer and the Software system. Structural descriptors are calculated and multiple linear regression analysis are used to generate model equations relating structural features to observed retention characteristics then was treated with two methods. The detection of influential observations for the standard least squares regression model is a problem which has been extensively studied. LAD regression diagnostics offers alternative dicapproaches whose main feature is the robustness. Here a nonparametric method for detecting influential observations is presented and compared with other classical diagnostics methods. Comparisons are between models generated for the two stationary was carried out with two methods, and descriptors that may encode differences in solute interactions with stationary phases of differing polarity are discussed and validated results in the state approached by the tests statistics: Test of Anderson-Darling, shapiro-wilk, Agostino, Jarque-Bera and the confidence interval thanks to the concept of robustness to check if the distribution of the errors is really approximate. \end{abstract} \keywords{LAD Regression, Robustness, Outliers, Leverage points, tests statistics.} \begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) \textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark} \end{textblock*} \begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) \textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}} \end{textblock*} \let\tabcellsep& \section[{Introduction}]{Introduction}\par ome compose food are volatile heterocyclic which are found in a natural way in our environment and the attraction which the men test for the flavours is ever contradicted during centuries and which have an interest in multiple fields, in particular in the food like flavour. Their presence in food results mainly, of process requiring a stage of cooking (partial or supplements), Egyptian civilization already used them for the kitchen.\par In the evaluation of the environmental risks, information on the fate in the environment, the properties, the behavior and the toxicity of a chemical substance is fundamental need.\par The volatile heterocyclic also constitutes a significant family of odorous molecules, particularly interesting in the field of the chemistry of the flavours. They represent more than one quarter of the 5000 volatile compounds insulated and characterized to date in our food.\par Pyrazines are heterocyclic very present in our food. More than 80 derived from pyrazines were identified in a great number of cooked food, like the bread, the meat, the torrefied coffee, the cocoa or the hazel nuts; they are very powerful aromatizing compounds. \hyperref[b24]{Mihara and Enomoto (1985)}, described a relation structure/retention for a unit of substituted pyrazines for which the increments of indices relating to various substituents on the cycle were given for a small series of substituents present. The method was then extended to integrate others substituents, by adding a term which takes account of the position on the cycle of a substituent compared to the others \hyperref[b25]{(Mihara \& Masuda, 1987)}. In a similar approach, \hyperref[b23]{Masuda and Mihara (1986)} describe the use of indices of connectivity modified to calculate in advance the indices of retention of a series of substituted pyrazines. The methods lead to good results, in so far as the increments of indices determined in experiments available for the unknown compounds are implied, which constitutes their principal defect. \hyperref[b0]{Stanton and Jurs (1989)}, used methodology QSRR to develop models connecting the structural characteristics of 107 variously substituted pyrazines, with their indices of retention obtained on two columns of very different polarities (OV-101 and Carbowax-20M). The equations were calculated using the multilinear regression, the choice of the explanatory variables (topological, electronic and physical properties) being realized by progressive elimination {\ref (Swall \& Jurs, 1983)}, among the 85 individual molecular descriptors obtained for each whole molecule. The indices of retention (IR) obtained on each column were treated separately, while drawing from the same sets of descriptors. The models calculated with 6 explanatory variables provide high standards errors (S = 23 units of index -u.i. -on OV-101 and S = 36.33 u.i. out of Carbowax -20 M) which do not predict good predictive capacities for these models, and which let suppose nonlinear relations between descriptors and property (IR) studied.\par The objective of this work aims at using methodology QSRR, the approach Method LAD /Least square (LAD/OLS), to model the indices of retention of (114) pyrazines reported from Davit T. Stanton and Peter C.Jurs (1989) and reported from Mihara and Enomoto pyrazines of test and 25 of validation) on O V-101 and Carbowax -20M are successfuty modeled with the ald of a computer and the Software system. Structural descriptors are calculated and multiple linear regression analysis are used to generate model equations relating structural features to observed retention characteristics then was treated with two methods. The detection of influential observations for the standard least squares regression model is a problem which has been extensively studied. LAD regression diagnostics offers alternative dicapproaches whose main feature is the robustness. Here a nonparametric method for detecting influential observations is presented and compared with other classical diagnostics methods. Comparisons are between models generated for the two stationary was carried out with two methods, and descriptors that may encode differences in solute interactions with stationary phases of differing polarity are discussed and validated results in the state approached by the tests statistics: Test of Anderson-Darling, shapiro-wilk, Agostino, Jarque-Bera and the confidence interval thanks to the concept of robustness to check if the distribution of the errors is really approximate. {\ref (1985)}, the molecular descriptors being only calculated starting from the chemical structure of the compounds.\par The linear statistical model for fixed purposes will be examined by two robust methods for the evaluation of the parameters of regression starting from estimates of the robust coefficients of regression most popular by the appendices. We based ourselves on the comparison between the two methods, the applicability (DA) will be discussed using the diagram of Williams who represents the residues of prediction standardized according to the values of the levers (hi) {\ref (Eriksson et al..2003;} {\ref Tropsha et al.2003)}. We present the tests statistics and graph of compatibility at the normal law for validated the results of the state approached between the two methods for a risk ?= 5\%. \section[{II.}]{II.} \section[{Methodology i. Descriptor Generation}]{Methodology i. Descriptor Generation}\par One used the molecular software of modeling Hyperchem 6.03, for to represent the molecules, then using semi-empirical method AM1 {\ref (Dewar et al.,. 1985};. Holder 1998) to obtain the final geometries. It is established \hyperref[b20]{(Levine, 2000)} that this Method gives good results when one treats small molecules (of less than one hundred atoms), like those considered in this work.\par The optimized geometries were transferred in the software dragon from data-processing software version 5. {\ref 4[19]}, for the calculation of 1320 descriptors while operating on 89 pyrazines of test; subsets of descriptors were chosen by genetic algorithm, these descriptors can be separate in four categories: topological descriptors of The topological, geometrical, physical, and electronic accounts of way and molecular indices of connectivity included. The geometrical descriptors included sectors of shade, the length with the reports/ratios of width, volumes of van der Waals, the surface, and principal moments of inertia. The calculated descriptors of physical property included the molecular refringency of polarizability and molar. The electronic descriptors included most positive and most negative described by {\ref Kaliszan.} By employing the software Mobydigs (Todeschini et al., 2009) \hyperref[b18]{[21]} and by maximizing the coefficient of prédiction Q 2 and minimal R 2 of S (the error). \section[{ii. Regression Analysis}]{ii. Regression Analysis}\par The analysis of the multiple linear regressions was carried out with two methods by software Matlab (R2009a) for (LAD) and Minitab 16 for (OLS).\par One considers the multiple model of regression given by \hyperref[b8]{[9]}:?? ?? = ?? 0 + ? ?? ?? ???1 ?? =1 ?? ???? + ?? ??\textbf{(1)}\par The detection of meaningless statements and `with action leverage according to the method of least squares is a problem which' was largely studied. The diagnosis by the regression LAD offers alternative approaches whose principal characteristic is the robustness. In our study a non-parametric method to detect the meaningless statements and the point's lever was applied and compared with the traditional method of diagnosis (least squares) \hyperref[b8]{[9]}.\par iii. Method of least squares OLS This one was carried out with the software Minitab 16 \hyperref[b28]{[33]}, method MLR applied to the multiple regression consists in defining the ? estimate which minimizes {\ref ([9, 17, 18]}:? ???? 2 = ? (yi-?0-? ????????) 2\textbf{(2)}\par iv. Least Absolute Deviations (LAD)\par The analysis of linear regression multiple was carried out with the software Matlab (R2009a) \hyperref[b26]{[31]}, by using the method of the least variations in absolute value, said method LAD (Least Absolute Deviations), is one of the principal alternatives to the method of least squares when it is a question of estimating the parameters of a model of regression, which minimizes the absolute values and not the values with the square of the term of erreur. La method stable-lad applied to the multiple regression consists in defining the ? estimates which minimize {\ref [9, 17, the 18]}:?|????| = ? |yi-?0-? ????????\textbf{(3)}\par III. \section[{The Data Set}]{The Data Set}\par One uses the molecular software Hyperchem The retention data for the114 compounds chromatographed on the OV-101 and CRW-20M stationary phases were taken from (113 taken from Davit T. Stanton and Peter C. Jurs (1) and 1 compound (2-VinylPyrazine) taken from Mihara and Enomoto \hyperref[b24]{[29]}) and are listed in table 1. \section[{IV.}]{IV.} \section[{Results and Discussion}]{Results and Discussion}\par An ideal model is one that has a high R value, allow standard error, and the fewest independent variables \hyperref[b0]{[1,}\hyperref[b8]{9]}. The best models found has 3 descriptors for each stationary phase by using the software Moby Digs \hyperref[b18]{[21]} are given below. The criterion for identifying a compound as an outlier was that compound being flagged by three or more of six standard statistical tests used to detect outliers in regression analysis. These tests were (1) residual, (2) standardized residual, (3) Studentized residual, (4) leverage, (5) DFFITS, (6) Cook's distance. The residual is the difference between the actual value and the value predicted by the regression equation. The standardized residual is the residual divided by the standard deviation of the regression equation. The Studentized residual is the residual of a prediction divided by its own standard deviation.\par Leverage allows for the determination of the influence of a point in determining the regression equation. DFFITS describes the difference in the fit of the equation caused by removal of a given observation, and Cook's distance describes the change in a model coefficient by the removal of a given point. The definition of each descriptor is given table 2: The coefficient of multiple determinations (R 2 ) indicates the amount of variance in the data set accounted for by the model. The standard error of the regression coefficient is given in each case, and n indicates the number of molecules involved in the regression analysis procedure \hyperref[b0]{[1,}\hyperref[b8]{9]}.\par a) The best models IR(OV-101): (MPC03, X1sol, GATS5e, AEigp, L3e,Qpos); -S=20.892, R 2 =99.30, n=89 compounds. IR(RWC): Se, Mp, X1sol, DP01, Mor06v, Tm; S=22.64, R 2 =99.22, n=89 compounds.\par Indeed Figure {\ref 1} reproduced the distributions of the standard residues di (ordinary residue report /root of the average square of the variations) according to the adjusted values, which seem random (without particular tendencies).That shows the constancy of variances ? 2 , it be-with saying their independence of the regresses and the adjusted dependent variable.\par The quasi-linearity (R = 0, 9951; OV-101 -R = 0, 9835; Carbowax-20M -critic = 0, 96048) of the diagram of the normal scores (Figure {\ref 2}) is an index of normality. Values of the statistics of Durbin-Watson (Durbin, \& Watson, 1951), [d= 1,33535; OV-101/D = 1,66161; Carbowax-20M] are the greater than higher values given by the tables, respectively for 3 regresses, and any reasonable risk ?, which establishes each time the independence of the residues. The diagnostic statistics joined together in Table {\ref 3} make it possible to make comparisons and to draw several conclusions \hyperref[b18]{[21]}. Values of R 2 and of R 2 (adj) show, each time, quality of adjustment, whereas the very weak differences between R 2 and Q 2 inform about the robustness of the models which are, moreover, very highly significant (high values of the statistics F of Fisher).\par Moreover, the similarity of SDEP and SDEC mean that the internal capacities of prediction models are not too dissimilar their capacities of adjustment.\par The validation by bootstrap (Q BOOT ) confirms all at the same time the capacity of internal prediction and the stability of the models. \section[{b) Robust Regression}]{b) Robust Regression}\par Any robust method must be reasonably effective once compared to the estimators of least squares; if the fundamental distribution of the errors is normal and primarily more effective independent than the estimators of least squares, when there are peripheral observations. There are various robust methods for the evaluation the parameters of regression. The principal goal of this section is the method LAD (nap of the absolute values of the errors) whose coefficient of regression qualifies the robustness among the additional data \hyperref[b15]{[16]}.\par i \section[{. Comparison Robust Regression of OLS and LAD}]{. Comparison Robust Regression of OLS and LAD}\par More particularly we will test 2 methods of estimate for the vector of the Parameters ((?? 0 * ,?? 1 * , ? , ?? ?? * ):\par -Method of least squares ordinary, more known and the most used.\par -The method LAD (Sum of the absolute values of the errors.)\par The large advantage of the method LAD is his robustness, i.e. that the estimators are not impact by the extreme values, (they are known as "robust"). It is thus particularly interesting to use the method LAD if one is in the presence of aberrant values in comparison with method OLS \hyperref[b7]{[8]}. \section[{ii. Comparison of hyperplanes of regression}]{ii. Comparison of hyperplanes of regression}\par Each equation on each column check the assumptions on the same linear statistical model for Fixes purposes for each method in comparison with the hyperplane calculated by LAD compared to the hyperplane calculated by the method of least squares.\par It is noticed that ??the calculated OLS are not very different for the regression with ?? the LAD on the two columns, except, ??1 the calculated OLS is almost the same ones as for the regression with ??1 the LAD on column CRW and ??4 the calculated OLS is almost the same ones as for the regression with ??4 the LAD on column OV-101.\par It is thus relevant to remake a checking of the presences of aberrant values by using the following stage (figure {\ref 3}):\par The hyperplane of regression can radically change, with the change of the coefficients of the hyperplane.\par iii \section[{. Graphical Comparisons of Alternative Regression Models}]{. Graphical Comparisons of Alternative Regression Models}\par The field of application was discussed using the diagram of Williams. The analysis of the residues shows that the observations (82,25) residues raised but it (48) point influence in the two estimates and the observation \hyperref[b11]{(12)} point influence with the LAD estimate and lever by least square also observation 4 residue raised with OLS and not lever with LAD in the whole of validation on column OV -101 and on column CRW -20M the observations (45) not influence in the two estimates and observation 16 point influence in the two estimates in the whole of validation.\par After elimination of the aberrant points collective between the two methods and after the secondary treatment one has the observation (12) point influence and the observations \hyperref[b0]{(1,}\hyperref[b20]{24)} residues raised in the two estimates but it (25) observation 4 residue raised with OLS and not lever with LAD also the observation 4 residue raised in the whole of validation in the two estimates on column OV -101 and on column CRW -20M the observations (45) not influence in the two estimates and observation 16 point influence in the two estimates in the whole of validation and on column CRW -20M the observations (24 {\ref 25 35)} residues raised but it (84)point influence in the two estimates and observation 8 point influence in the two estimates in the whole of validation.\par Thus finally the models in which the meaningless statements were removed become after elimination of the aberrant points collective [OV-101: test - \hyperref[b0]{(1,}\hyperref[b11]{12,}\hyperref[b20]{24)}, validation (4), CRW-20M: test - {\ref (24, 25, 35 84)}, validation \hyperref[b7]{(8)} \par It is noticed besides that ?? the OLS calculate more to approach which for the regression with ?? the LAD on the two columns into precise (??1 , ??3 ?????? ??4) the OLS calculate are almost the same ones as for the regression with (??1 , ??3 ?????? ??4) the LAD and on the same order with (??0 , ??5 ?????? ??6) on OV 101 and ??1the OLS calculate are almost the same ones as for the regression with ??1 the LAD on CRW -20M.\par The analysis of the residues shows that in this case All the point of lad method between (-2, 2), but it the analysis of the residues of OLS method shows that the observations [OV-101: test -(6,42), CRW-20M: test -(22, 24, 67 ,78), validation {\ref (7 ,13,14)}] the LAD estimate given good result On the other hand estimate OLS figure (4): iv. Graphical Comparisons of Alternative Regression Models 0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0,0 Column OV -101 Columns RW -20M Method LAD and OLS (test, validation) Fig. 4 : Diagram of Williams of the residues of prediction standardized according to the lever Lastly, it is noted that LAD is a robust estimator but loses stability in the presence of points aberrant.\par We note however the observation that the estimate the least square is near to the LAD estimate to which removed the aberrant values.\par To conform the approach between the two methods and to deduce the robust method between them, There is a package of tests of normality (of the standard errors or residues?) indeed, thanks to the concept of robustness, we can used simple techniques (descriptive e.g. statistics, technical graphs) to check if the distribution of the data is really approximate.\par Any test is associated a ? risk known as of first species years works us, we will adopt it risk ? = 5\%. \section[{c) Comparisons of the Tests of normality of the errors between the method LAD and OLS in the approached state}]{c) Comparisons of the Tests of normality of the errors between the method LAD and OLS in the approached state}\par The software Minitab 16 carries out automatically the estimate of the two principal parameters of the normal law (? the Mean (OV-101:0, CRW-20M:0), ? the variation-type(OV-101:13.26, CRW-20M:18.53) for OLS one applying the same principle with the method LAD but one used (it median (OV-101:-0.96, CRW-20M:0.01) ? variation-type (OV-101:13.84, CRW-20M:18.66) and with the number principal in the state approached to the two columns n=32. \section[{i. Test statistical a. Test of Anderson-Darling}]{i. Test statistical a. Test of Anderson-Darling}\par In our work, one finds us that AD [OV -101: (lad) = 0.250 with value of p>0.250, (OLS) = p=0.938 with value of p = 0.747, n=82]-RCW-20M: (lad) = 0.547 with value of p > 0.250, (OLS) = 0.165 with value of p=0.572 n=84] 0.1. To 5\%, the assumption of normality is compatible with the method LAD and OLS \hyperref[b28]{[33,}\hyperref[b29]{34,}\hyperref[b30]{35]}. \section[{b. test of Shapiro-Wilk}]{b. test of Shapiro-Wilk}\par It is particularly powerful for small manpower (n<50) for this that one using for valid the results of the validation.\par For a risk ? = 0. 05, the critical points read in the table of Shapiro-Wilk for n = 23 is W crit = 0. 914 and for n=24 and W crit = 0. 916. In our works, on (OV) [W LAD =0.9969, W MLR = 0.9877, n=24] and on CRW [W LAD = 0, 0.997, W MLR = 0,9227, n=23 ] W> W crit , with the risk of 5\%, the assumption of normality compatible with us is given (normal law) \hyperref[b29]{[34,}\hyperref[b30]{35]}. \section[{c. Test of D'Agostino}]{c. Test of D'Agostino}\par For ? = 0.05, the threshold critic is ?2 0:95(2) = 5.99.In our works, on (OV) [:(W LAD = 0,0072 with value of p = 0,99, W OLS = 0.042 with value of p = 0.97, n=82),: ] and on CRW [ (W LAD = 0,1202 with value of p = 0.94, W OLS = 0,00116 with value of p = 0.99, n=84), ] W 0.1 with the risk of 5\%, the assumption of normality compatible with us is given (normal law) \hyperref[b28]{[33,}\hyperref[b29]{34,}\hyperref[b30]{35]}. \section[{d. Test of Jarque-Bera}]{d. Test of Jarque-Bera}\par As the Test of Agostino It becomes particularly effective starting from N>20 for this that one using for valid the results.\par For ? = 0.05, the critical point is ?2 0:95(2) = 5.99. In our works, on (OV)[ (W LAD = 0,0971with value of p = 0.95, W OLS = 0.0949 with value of p = 0.95, n=82), ] and on CRW [ (W LAD = 0.1059 with value of p = 0.94, W OLS = 0,0979 with value of p = 0.95, n=84), ] W 0.1 than the risk of 5\%, the assumption of normality compatible with us is given (normal law). \hyperref[b28]{[33,}\hyperref[b29]{34,}\hyperref[b30]{35]} Completely all the statistical tests is accepted the data of the state approached between the two methods especially the test of Shapiro-Wilk the value of the method LAD closer to method OLS and the other tests the values of the method LAD is higher has the method MLR which explains than give them method LAD is effective and robust para for give method OLS.\par Completely all the statistical tests is accepted the data of the state approached between the two methods especially the test of Shapiro-Wilk the value of the method LAD closer to method MLR and the other tests the values of the method LAD is higher has the method OLS which explains than give them method LAD is effective and robust para for give method OLS.\par e. \section[{Interval of confidence}]{Interval of confidence}\par The confidence interval and the risk ?? = 0.05 constitute a complementary approach thus (an approach of estimate) the most used confidence interval is the confidence interval has 100 (1 -?) = 95 \%.\par The Column OV-101: LAD :(-28.11, 26.17), OLS (-25.9, 25.99).\par The Column CRW-20M: LAD (-36.56, 36.58), OLS {\ref (-36.34, 36.34}).\par These result is formed L approximate of two method.\par You can be 95\% confident that the 50th percentile for the population is between OV-101 (LAD:-3.96 and 2.027,-OLS:-2.87 and 2.87, CRW-20M (LAD:-3.98 and 4.00, OLS:-3.96 and 3.96) \hyperref[b28]{[33,}\hyperref[b29]{34,}\hyperref[b30]{35]}.\par V. \section[{Conclusion}]{Conclusion}\par The modeling of the indices of retention of 114 pyrazines (89 tests and 25 validations) eluted out of two columns various OV -101 and CRW-20M by two methods LAD and OLS are based on the following comparisons:\par a) The comparison of the equations of the hyperplanes L equations of OLS is closer to LAD after elimination of the aberrant points for the ?2 (LAD) ??2(OLS) and the other coefficient remaining with the same order for column OV-101 Pour the column Crw-20m the ?1 (LAD) ??1(OLS) and the other coefficient remaining with the same order after the secondary treatments for the checking of the presence of aberrant values (82, 48, 26, 25, 24,12, 1) on column OV -101 and item {\ref (45,} {\ref 82,}\hyperref[b30]{35,}\hyperref[b20]{24,}\hyperref[b21]{25)} for the column CRW-20M, and to be able to compare them By employing the following stage. \section[{b) Graphic comparison: The applicability was discussed using the diagram of Williams in dependence}]{b) Graphic comparison: The applicability was discussed using the diagram of Williams in dependence}\par Lastly, it is noted that LAD is a robust estimator but loses his stability in the presence of aberrant points.\par Used test of normality's of the errors by statistical test. One applied compatibility with the normal law, but to differing degrees using p-been worth. One notes that the touts test to accept the assumption of normality is that of Anderson-Darling, the test of Shapiro-Wilk His power is recognized in the literature.\par Lastly, the tests of Agostino and Jarque-Bera, based on the coefficients of asymmetry and flatness accepts readily the assumption of normality with one p-been worth sup 0.1 on the columns, Too one confirmed approached graphically by histogram of frequency in finished by the confidence interval.\par It general this study is shown that results by the two estimates theoretical (equation) and graph give good results expressed by the models.\begin{figure}[htbp] \noindent\textbf{5}\includegraphics[]{image-2.png} \caption{\label{fig_1}54 5 -}\end{figure} \begin{figure}[htbp] \noindent\textbf{1} \par \begin{longtable}{P{0.03636560069144339\textwidth}P{0.35557476231633534\textwidth}P{0.06722126188418323\textwidth}P{0.35557476231633534\textwidth}P{0.004407951598962834\textwidth}P{0.030855661192739848\textwidth}} n°\tabcellsep Compounds\tabcellsep ov-101\tabcellsep Compounds\tabcellsep IR(cw)\tabcellsep \\ 1\tabcellsep Pyrazine\tabcellsep 710\tabcellsep Pyrazine\tabcellsep \tabcellsep \\ 2 3 4\tabcellsep Methylpyrazine 2,3-dimethylpyrazine 2,5-dimethylpyrazine\tabcellsep 801 897 889\tabcellsep Methylpyrazine 2,3-dimethylpyrazine 2,5-dimethylpyrazine\tabcellsep \tabcellsep Year 2016\\ 5\tabcellsep 2,6-dimethylpyrazine\tabcellsep 889\tabcellsep 2,6-dimethylpyrazine\tabcellsep \tabcellsep \\ 6\tabcellsep Trimethylpyrazine\tabcellsep 981\tabcellsep Trimethylpyrazine\tabcellsep \tabcellsep 19\\ 7 8 9 10 11 n° 12 13 14 15 16\tabcellsep Trimethylpyrazine Ethylpyrazine 2-ethyl-5-methylpyrazine 2-ethyl-6-methylpyrazine 2,5-dimethyl-3-ethylpyrazine Compounds 2,6-dimethyl-6-ethylpyrazine 2,3-dimethyl-5-ethylpyrazine 2,3-diethylpyrazine 2,3-diethyl-5-methylpyrazine Propylpyrazine\tabcellsep 1067 894 980 977 1059 ov-101 1064 1066 1065 1137 986\tabcellsep Trimethylpyrazine Ethylpyrazine 2-ethyl-5-methylpyrazine 2-ethyl-6-methylpyrazine 2,5-dimethyl-3-ethylpyrazine Compounds 2,6-dimethyl-6-ethylpyrazine 2,3-dimethyl-5-ethylpyrazine 2,3-diethylpyrazine 2,3-diethyl-5-methylpyrazine Propylpyrazine\tabcellsep IR(cw)\tabcellsep Volume XVI Issue I Version I\\ 17 18\tabcellsep 2-methyl-3-propylpyrazine 2,3-dimethyl-5-propylpyrazine\tabcellsep 1072 1154\tabcellsep 2-methyl-3-propylpyrazine 2,3-dimethyl-5-propylpyrazine\tabcellsep \tabcellsep ( B )\\ 19 20 21 22 23 24 25 26 27 28 29 30 31 32\tabcellsep 2,5-dimethyl-3-propylpyrazine 2,6-methyl-3-propylpyrazine Isopropylpyrazine 2,3-dimethyl-5-isopropylpyrazine Butylpyrazine 2-butyl-3-methylpyrazine 3-butyl-3,5-dimethylpyrazine 3-butyl-3,6-dimethylpyrazine 5-butyl-2,3-dimethylpyrazine Isobutylpyrazine 2,3-dimethyl-5-isobutylpyrazine 2-isobutyl-3,5,6-trimethylpyrazine sec-butylpyrazine 5-sec-butyl-2,3-dimethylpyrazine\tabcellsep 1142 1151 949 1112 1088 1121 1184 1196 1254 1043 1200 1263 1040 1194\tabcellsep 2,5-dimethyl-3-propylpyrazine 2,6-methyl-3-propylpyrazine Isopropylpyrazine 2,3-dimethyl-5-isopropylpyrazine Butylpyrazine 2-butyl-3-methylpyrazine 3-butyl-3,5-dimethylpyrazine 3-butyl-3,6-dimethylpyrazine 5-butyl-2,3-dimethylpyrazine Isobutylpyrazine 2,3-dimethyl-5-isobutylpyrazine 2-isobutyl-3,5,6-trimethylpyrazine sec-butylpyrazine 5-sec-butyl-2,3-dimethylpyrazine\tabcellsep \tabcellsep Global Journal of Human Social Science -\\ 33\tabcellsep Pentylpyrazine\tabcellsep 1192\tabcellsep Pentylpyrazine\tabcellsep \tabcellsep \\ 34\tabcellsep 2,3-dimetyl-5-pentylpyrazine\tabcellsep 1352\tabcellsep 2,3-dimetyl-5-pentylpyrazine\tabcellsep \tabcellsep \\ 35\tabcellsep Isopentylpyrazine\tabcellsep 1157\tabcellsep Isopentylpyrazine\tabcellsep \tabcellsep \\ 36\tabcellsep 2,3-dimetyl-5-isopentylpyrazine\tabcellsep 1317\tabcellsep 2,3-dimetyl-5-isopentylpyrazine\tabcellsep \tabcellsep \\ 37\tabcellsep (2-methylbutyl)pyrazine\tabcellsep 1151\tabcellsep (2-methylbutyl)pyrazine\tabcellsep \tabcellsep \\ 38\tabcellsep 2,3-dimethyl-5-(2-methylbutyl)pyrazine\tabcellsep 1306\tabcellsep 2,3-dimethyl-5-(2-methylbutyl)pyrazine\tabcellsep \tabcellsep \\ 39\tabcellsep 2-(2-methylbutyl)-2,5,6-\tabcellsep 1363\tabcellsep 2-(2-methylbutyl)-2,5,6-\tabcellsep \tabcellsep \\ \tabcellsep trimethylpyrazine\tabcellsep \tabcellsep trimethylpyrazine\tabcellsep \tabcellsep \end{longtable} \par {\small\itshape [Note: © 2016 Global Journals Inc. (US)]} \caption{\label{tab_1}Table 1 :}\end{figure} \begin{figure}[htbp] \noindent\textbf{2} \par \begin{longtable}{P{0.14194499017681728\textwidth}P{0.7080550098231827\textwidth}} Name\tabcellsep Definition\\ MPC03\tabcellsep Molecular path count of order 03\\ \multicolumn{2}{l}{GATS5e Geary autocorrelation-lag 5/weighted by}\\ \tabcellsep atomic Sanderson electronegativityies\\ AEigp\tabcellsep Eigen value distance matrix sum from Polson\\ \tabcellsep arizability weight (Barysz matrix)\\ Qpos\tabcellsep total positive charge\\ Se\tabcellsep sum of atomic Sanderson electronegativityies\\ Mp\tabcellsep mean atomic polarizability (scaledon Carbon\\ \tabcellsep atom)\\ X1sol\tabcellsep salvation connectivity index chi-1\\ DP01\tabcellsep molecular profile no.01\\ Mor06v\tabcellsep (3D-MORSE-signal 06/weighted by atomic\\ \tabcellsep Vander Waals volumes\\ Tm\tabcellsep T (Total size index/weight atomic masses\end{longtable} \par \caption{\label{tab_2}Table 2 :}\end{figure} \footnote{Modeling Retention Indices of a Series Components Food and Pollutants of the Environment: Methods; OLS, LAD © 2016 Global Journals Inc. (US)} \footnote{Modeling Retention Indices of a Series Components Food and Pollutants of the Environment: Methods; OLS, LAD} \backmatter \begin{bibitemlist}{1} \bibitem[Ra ()]{b27}\label{b27} \textit{}, Matlab Ra . 2009a. \bibitem[Saporta ()]{b12}\label{b12} \textit{}, Gilbert Saporta . 2012. (Régression robuste) \bibitem[Pynnönen ()]{b6}\label{b6} ‘A Report on Least Absolute Deviation Regression with Ordinary Linear Programming’. Seppo Pynnönen , Timosalmi . \textit{Finnish Journal of Business Economics} 1994. 43 (1) p. . \bibitem[Dodge and Jureckova ()]{b17}\label{b17} \textit{Adaptive Regression}, Y Dodge , J Jureckova . 2000. New York: Springer-Verlag. \bibitem[Dodge and Valentin Rousson ()]{b8}\label{b8} \textit{Analyses de régression appliquée}, Y Dodge , Valentin Rousson . 2004. (paris) \bibitem[Chen et al.]{b9}\label{b9} \textit{Analysis of least absolute déviation}, Kani Chen , Zhiliang Ying , Hong Zhang , Lincheng Zhao . \bibitem[Damodar et al. ()]{b30}\label{b30} \textit{Basic Econometrics Fifth Edition}, N Damodar , Dawn C Gujarati , Porter . 2009. \bibitem[Mihara and Enomoto ()]{b24}\label{b24} ‘Calculation of retention indices of pyrazines on the basis of molecular structure’. S Mihara , N Enomoto . \textit{J. Chromatogr} 1985. 324 p. . \bibitem[Stanton and Jurs ()]{b0}\label{b0} ‘Computer-assisted prediction of gas chromatographic retention indexes of pyrazines’. D T Stanton , P C Jurs . \textit{Anal. Chem} 1989. 61 p. . \bibitem[Mihara and Masuda ()]{b25}\label{b25} ‘Correlation between molecular structures and retention indices of pyrazines’. S Mihara , H Masuda . \textit{J. Chromatogr} 1987. 402 p. . \bibitem[Ciuperca ()]{b11}\label{b11} \textit{Estimation robuste dans un modèle paramétrique avec rupture}, Gabriela Ciuperca . 2009. Bordeaux. \bibitem[Lee et al. ()]{b19}\label{b19} ‘Evaluation of predictive retention factors for phenolic compounds with QSPR equations’. Seung Ki Lee , Yulia Polyakova , Kyung Row , Ho . \textit{J. Liq. Chromatogr and Rel. Tech} 2004. 27 (4) p. . \bibitem[Soumaya]{b15}\label{b15} \textit{Indicateurs de la sensibilité de l'estimateur Least Absolute Déviation Assas Paris}, Rekaia Soumaya . \bibitem[Faria and Melfi ()]{b10}\label{b10} ‘Lad regression and nonparametric methods for detecting outliers and leverage points’. S Faria , G Melfi . \textit{Student} 2006. 5 p. . \bibitem[Levine ()]{b20}\label{b20} I N Levine . \textit{Quantum chemistry. 5}, (New Jersey) 2000. Prentice-Hall. (th ed.) \bibitem[Touhami et al. ()]{b2}\label{b2} ‘Modèle QSRR Hybrides Algorithme Génétique Régression Linéaire Multiple des indices de rétentions de pyrazines en chromatographie gazeuse’. Imen Touhami , Karima Mokrani , Djelloul Messadi . \textit{Lebanese Science Journal} 2012. 13 (1) . \bibitem[Machabert ()]{b7}\label{b7} \textit{Modèles en très grande dimension avec des outliers. Théorie, simulations, applications" paris}, Tiffany Machabert . 2014. \bibitem[Dr and Al -Noor ()]{b14}\label{b14} ‘Mohammad. Model of Robust Regression with Parametric and Nonparametric Methods’. Nadia H Dr , Asmaa A Al -Noor . \textit{Journal of Mathematical Theory and Modeling} 2013. 3 (5) . \bibitem[Parliment and Epstein ()]{b3}\label{b3} ‘Organoleptic properties of some alkyl-substituted Alkoxy-and alkylthiopyrazines’. T H Parliment , M F Epstein . \textit{J. Agric. Food Chem} 1973. 21 p. . \bibitem[Nornadiah Mohd ()]{b29}\label{b29} ‘Power Comparaisons of shapiro-wilk, Kolmogorovsmornov, lillieffors and Anderson-Darling tests’. Razali Nornadiah Mohd . \textit{Yab Bee Yah}, 2011. 2 p. . \bibitem[Kaliszan ()]{b4}\label{b4} ‘Quantitative relationships between molecular structure and Chromatographic retention’. R Kaliszan . \textit{CRC Crit. Rev. Anal. Chem} 1986. 16 p. . \bibitem[Kaliszan ()]{b5}\label{b5} \textit{Quantitative structure-chromatographic retention relationships}, R Kaliszan . 1987. New York: J. Wiley. \bibitem[Kaliszan ()]{b18}\label{b18} \textit{Quantitative structure-chromatographic retention relationships}, R Kaliszan . \url{http://www.disat.unimib.it22} 1987. New York: J. Wiley. (Moby Digs 1.1) \bibitem[Release 16.1, Statistical Software ()]{b28}\label{b28} \textit{Release 16.1, Statistical Software}, 2003. \bibitem[Saporta ()]{b13}\label{b13} \textit{Régression robuste Régression non-paramétrique}, Ndèyeniang-Gilbert Saporta . 2014. \bibitem[Dodge ()]{b16}\label{b16} \textit{Statistique: Dictionnaire encyclopédique}, Y Dodge . 2004. Springer-Verlag France Paris. \bibitem[Masuda et al. ()]{b22}\label{b22} ‘Synthesis of new pyrazines for flavor use’. H Masuda , Y Misaku , T Shibamoto . \textit{J. Agric. Food Chem} 1981. 29 p. . \bibitem[Berlin ()]{b1}\label{b1} \textit{The Pyrazine}, G B Berlin . 1982. New York: Wiley-Interscience. \bibitem[Buchbauer ()]{b26}\label{b26} \textit{Threshold-based structureactivity relationships of pyrazines with bellpepper Flavor}, G Buchbauer . 2000. \bibitem[Magnuson et al. ()]{b21}\label{b21} ‘Topological indices based on neighbor. 26. Symmetry: chemical and biological applications’. V R Magnuson , D K Harriss , S C Basak . \textit{Chemical Applications of Topology and Graph Theory. R.B. King}, (Amsterdam) 1983. Elsevier. p. . \bibitem[Masuda and Mihara ()]{b23}\label{b23} ‘Use of modified molecular connectivity indices to predict retention indices of monosubstituted alkyl, alkoxy, alkylthio, phenoxy and (phenylthio) pyrazines’. H Masuda , S Mihara . \textit{J. Chromatogr} 1986. 366 p. . \end{bibitemlist} \end{document}