Resumo
Este trabalho apresenta um método sequencial, envolvendo aplicação de Análise em Componentes Principais (ACP) e logit binomial para previsão de escolha por modo de transporte motorizado. A aplicação da ACP permite reduzir o banco de dados multicolinear a componentes não correlacionadas entre si. Tais components são utlizadas, posteriormente, como variáveis explicativas em modelos logísticos binomiais. Os dados utilizados para o desenvolvimento deste trabalho são provenientes da entrevista domiciliar da Pesquisa Origem-Destino de 2007, realizada na Região Metropolitana de São Paulo. Os modelos obtidos apresentaram bom poder preditivo e valores coerentes e significativos de parâmetros calibrados. Na etapa de validação, foram obtidas taxas de acertos variando entre 69% a 92%. Finalmente, o método proposto é razoável, sendo uma boa alternativa para o caso de dados multicolineares em métodos de regressão.
Palavras-chave:
Escolha modal; Análise em Componentes; Principais; Logit Binomial
Abstract
This paper presents a sequential method involving application of Principal Component Analysis (PCA) and binomial logit for forecasting of motorized travel mode choice. The application of ACP reduces the multicolinear database into uncorrelated components. Such components are used later as explanatory variables in binomial logit models. The data used are from the Origin-Destination Survey carried out in 2007 in São Paulo Metropolitan Area. The gotten models showed good accuracy and consistent and significant calibrated parameters. In the validation step, the hit rates were obtained ranging from 69% to 92%. Finally, the proposed method is reasonable to be a good alternative for the case of multicollinear data used in regression methods.
Key words:
travel mode choice; Principal Component Analysis; Logistc regression methods
Introdução
Este trabalho visa apresentar um método sequencial, envolvendo aplicação de Análise em Componentes Principais (ACP) e logit binomial para previsão de escolha por modo de transporte motorizado. A aplicação da ACP permite reduzir o banco de dados multicolinear a componentes não correlacionadas entre si. Tais componentes, extraídas pela ACP, são utilizadas como variáveis explicativas no modelo logit binomial. Desta forma, trata-se de uma abordagem exploratória-confirmatória que permite investigar preferências por modo de transporte motorizado. O método proposto é eficiente, sobretudo para o caso de banco de dados com multicolinearidade.
Técnicas de regressão múltipla são versáteis e poderosas. São aplicadas em uma infinidade de casos, onde se deseja encontrar uma relação entre uma única variável dependente e diversas variáveis independentes, com estimação de parâmetros a partir de diferentes critérios. Um estimador de um parâmetro ou um vetor de parâmetros desconhecidos é uma variável aleatória cujo valor pode ser calculado a partir de uma amostra. O vetor de parâmetros pode ser estimado por vários métodos, tais como Método dos Mínimos Quadrados (Regressão Linear Múltipla) e o Método da Máxima Verossimilhança (Regressão Logística – logit), que são os mais utilizados.
Multicolinearidade ocorre quando duas ou mais variáveis explicativas são muito correlacionadas entre si. Utilizando-se apenas modelos de regressão, torna-se difícil distinguir suas influências separadamente. Outra suposição do modelo de regressão é que nenhuma relação linear exata pode existir entre quaisquer covariáveis ou combinações lineares destas. Quando se viola esta hipótese têm-se o problema de multicolinearidade perfeita. Por outro lado, se as variáveis não estão correlacionadas entre si, denomina-se, este caso, ausência de multicolinearidade, sendo chamada de ortogonal a regressão com estas variáveis. O caso intermediário, muito comum em problemas reais, ocorre quando a correlação entre duas ou mais variáveis é alta, sendo esta situação chamada de alto grau de multicolinearidade.
Geralmente, a multicolinearidade não aumenta o poder preditivo de modelos de regressão, sendo uma tarefa usualmente difícil a seleção de variáveis explicativas multicolineares (Camminatiello e Lucadamo, 2010Camminatiello, I.; Lucadamo, A. (2010) Estimating Multinomial Logit Model with Multicollinear Data. Asian Journal of Mathematics & Statistics, vol 3,2, 93-101. DOI: 10.3923/ajms.2010.93.101.
https://doi.org/10.3923/ajms.2010.93.101...
). Um número alto de variáveis explicativas e correlacionadas pode tornar os modelos de regressão mais redundantes do que realmente bons. Para evitar os problemas provocados pela multicolinearidade o método mais simples é a eliminação, do modelo completo, das variáveis com os coeficientes estatisticamente não significativos para encontrar o melhor subconjunto de variáveis independentes. Outra alternativa, proposta por Hoerl e Kennard (1970)Hoerl, A.; Kennard, R.W. (1970) Biased Estimation for Nonorthogonal Problems. Technometrics, Vol. 12, No. 1 (Feb., 1970), pp. 55-67., é o método de regressão denominado de “Ridge”, que tem o objetivo de melhorar a precisão dos parâmetros estimados, sem o termo constante, por padronizar as variáveis independentes. Neter et al. (1989)Neter, J.; Wasserman, W.; Kutner, M.H. (1989). Applied Linear regression Models. 2nd Edition. Irwin Homewood IL., no entanto, colocam como principal limitação do modelo anterior, a impossibilidade de fazer inferências sobre seus parâmetros.
Outra forma de obter um modelo adequado, quando algumas variáveis independentes são muito correlacionadas, é a partir da técnica de Análise em Componentes Principais, que tem a vantagem de não descartar nenhuma variável explicativa. Alguns autores propuseram métodos diferentes, com perda de pouca informação a respeito da variância dos dados e redução significativa de variáveis independentes (Wold, 1985Wold, H. (1985) Partial Least Squares, Encyclopedia of statistical sciences, vol 6. Pp 581-591.; Frank et al., 1993Frank, I.E.; Friedman, J.H.; Wold, S.; Hastie, T.; Mallows, C. (1993) A statistical view of some chemometrics regression tools. Technometrics, c, pp 109-148.; Aguilera et al., 2006Aguilera, A.M.; Escabias, M.; Valderrama, M.J. (2006) Using principal components for estimating logistic regression with high-dimensional multicollinear data. Computational Statistics & Data Analysis 50, pp. 1905 – 1924. DOI:10.1016/j.csda.2005.03.011
https://doi.org/10.1016/j.csda.2005.03.0...
).
Camminatiello e Lucadamo (2010)Camminatiello, I.; Lucadamo, A. (2010) Estimating Multinomial Logit Model with Multicollinear Data. Asian Journal of Mathematics & Statistics, vol 3,2, 93-101. DOI: 10.3923/ajms.2010.93.101.
https://doi.org/10.3923/ajms.2010.93.101...
propuseram o modelo de regressão logística multinomial para dados multicolineares, desenvolvendo, a partir de dados simulados, uma extensão do modelo Principal Components Logistc Regression (PCLR). Aguilera et al. (2006)Aguilera, A.M.; Escabias, M.; Valderrama, M.J. (2006) Using principal components for estimating logistic regression with high-dimensional multicollinear data. Computational Statistics & Data Analysis 50, pp. 1905 – 1924. DOI:10.1016/j.csda.2005.03.011
https://doi.org/10.1016/j.csda.2005.03.0...
também desenvolveram método de extração de componentes e posterior uso em regressão logística com dados simulados. Chen et al. (2010)Chen, Z. J., Cheng, L., Deng, H. N., Zhang, J. K (2010) Analyzing Residential Travel Mode Choice Based on Principal Component Analysis. Proceedings of the 10th International Conference of Chinese Transportation Professionals. 2010. aplicaram metodologia similar para estudo de comportamento relativo a viagens. Na literatura de modelagem de transportes, por exemplo, o uso de tal método para aprimorar estimativas de demanda, resolvendo o problema de multicolinearidade, não é tão trivial.
Este trabalho visa melhorar estimativas, relacionadas a escolhas discretas de transportes (modo de transporte) através do uso de componentes principais em modelos logísticos tradicionais. Na Seção 1, serão apresentados conceitos básicos relativos às técnicas utilizadas; na Seção 2, são descritos materiais e método; na Seção 3 são apresentados os resultados obtidos e discussões. Finalmente, as conclusões são apresentadas e, em seguida, as referências bibliográficas utilizadas.
1. Técnicas abordadas
1.1 Análise em Componentes Principais (ACP)
A análise em componentes principais (ACP) é uma técnica de análise fatorial, que analisa inter-relações em uma estrutura multivariada. A ACP permite a transformação linear ótima de um conjunto de variáveis originais intercorrelacionadas em um novo conjunto de variáveis não correlacionadas. Esse procedimento se dá por uma rotação espacial dos dados, de forma que sejam salientadas as feições através da mudança no sistema de coordenadas e redução da massa de dados, com menor perda possível da informação (Jolliffe, 2002Jolliffe, I. T. (2002) Principal Component Analysis. 2a ed. Springer 2002. 518p.).
O entendimento da ACP está atrelado à aplicação da álgebra linear básica, em que são avaliados os autovetores e autovalores de uma matriz de covariância ou uma matriz de correlações. Dada uma matriz (m x n), em que m é o número de observações e n é o número de variáveis, a matriz de correlações é dada por:
Onde, x1 são os m valores observados para a variável 1 e xn são os m valores observados para a variável n.
A determinação das componentes principais é dada pelos autovetores (v). Os autovetores são calculados em função dos autovalores da matriz R. Sendo I a matriz identidade, os autovalores (λ) da matriz S são os escalares que satisfazem a equação característica:
Para cada autovalor, é possível calcular o respectivo autovetor v, conforme Equação 3.
No caso geral, a matriz de autovalores é uma matriz diagonal, onde o número de autovalores é igual à dimensão da matriz R (n). O resultado para os autovetores é equivalente à matriz quadrada R (n x n).
Um novo conjunto de variáveis (componentes W) pode ser obtido pela multiplicação dos coeficientes dos autovetores com o vetor de valores originais. Para isso, uma matriz quadrada A é composta usando os autovetores como colunas da matriz. As novas variáveis são combinações lineares das variáveis originais, calculadas conforme a Equação 4.
onde, A é a matriz composta pelos autovetores e X é o vetor de observações original.
A componente principal é o arranjo que melhor representa a distribuição dos dados e a componente secundária é perpendicular a componente principal, e assim por diante. A importância de cada componente é avaliada por meio de sua contribuição, dada em porcentagem e calculada pela variância dividida pela variância total (autovalores). A fim de obter, de forma preliminar, o número de componentes a extrair e de reduzir a massa de dados, a matriz inicial não rotacionada (A) é pertinente. Contudo, na maioria dos casos, esta solução não fornecerá a interpretação mais adequada das variáveis. Em geral, a rotação é desejável, porque simplifica a estrutura, melhorando a interpretação e reduzindo ambiguidades que, frequentemente, acompanham componentes não rotacionadas (Hair et al., 2010Hair Jr, J. F.; Black, W. C.; Babin, B. J.; Anderson, R. E. (2010) Multivariate Data Analysis. Prentice Hall. 7a ed. 785 p.). As rotações são caracterizadas em dois tipos: ortogonais e oblíquas. Neste trabalho a rotação ortogonal (Varimax) se mostrou conveniente.
1.2. Regressão Logística Múltipla (RLM)
Para prever escolhas discretas de modo de transporte, por exemplo, um determinado modo pode ser contrastado com uma opção alternativa, de maneira a ser transformado em uma probabilidade entre “0” e “1”. Um dos modelos de escolha discreta mais utilizados na demanda por transportes é o logístico (Ortúzar e Willumsen, 2011Ortúzar, J. D.; Willumsen, L. G. (2011) Modelling Transport. Londres: Wiley. 4ª ed. 586p.). A Regressão Logística se diferencia da Regressão Linear, pois seus cálculos envolvem variáveis dependentes qualitativas, ao invés de quantitativas.
A regressão logística é caracterizada pela quantidade de valores a serem discretizados. Caso haja apenas a duas opções, a regressão logística é dita binomial, caso haja mais opções, a regressão é generalizada e nomeada multinomial. Estas regressões são também denominadas de Modelos Logit Binomial e Logit Multinomial. A Regressão Logística permite o uso de um modelo (curva em S) para prever a probabilidade π de um evento categórico. A modelagem da curva S é dada por uma transformação logística da probabilidade π, conforme a Equação 5, da função logística g(x).
Através da Equação 5, deriva-se a equação de calibração (Equação 6).
Sendo x0=1, β0, β1,..., βn os coeficientes da equação de calibração e n o número de variáveis independentes. Em posse da equação de calibração e respectivos coeficientes, é simples calcular o valor de probabilidade π (Equação 7). A qualidade do ajuste pode ser mensurado por medidas estatísticas apropriadas, tais como: testes de regressão de Cox & Snell e Nagelkerke (Hair et al., 2010Hair Jr, J. F.; Black, W. C.; Babin, B. J.; Anderson, R. E. (2010) Multivariate Data Analysis. Prentice Hall. 7a ed. 785 p.).
2. Materiais e Método
2.1. Banco de dados e área de estudo
Os dados utilizados neste trabalho são referentes à Pesquisa Origem/Destino (O/D) realizada na Região Metropolitana de São Paulo (RMSP) em 2007 pela Companhia do Metropolitano de São Paulo. Na Pesquisa O/D da RMSP de 2007 foram levantadas informações de 30 mil domicílios, escolhidos aleatoriamente e distribuídos em 39 municípios, das 460 Zonas de Tráfego, na RMSP.
Os dados da Pesquisa O/D foram manipulados de forma que cada registro estava relacionado a informações por domicílio. Ao final foram utilizados 18.733 domicílios na amostra final. Originalmente, a Pesquisa O/D indica os modos realizados para cada viagem. Neste trabalho, os modos de transporte foram classificados como: Transporte Individual Motorizado, Transporte não Motorizado e Transporte Público. O Transporte Individual Motorizado inclui viagens em que o indivíduo dirige automóvel ou motocicleta. O Transporte não Motorizado refere-se ao uso de bicicleta e modo a pé. O Transporte Público inclui viagens por ônibus do município de São Paulo, ônibus de outros municípios, ônibus metropolitano, microônibus do município de São Paulo, microônibus de outros municípios, microônibus metropolitano, metrô e trem.
As variáveis originais, utilizadas para ACP foram: (1) coordenada geográfica x (UTM); (2) coordenada geográfica y (UTM); (3) Número de pessoas no domicílio; (4) Critério de Renda (variando de 1 a 8 – sendo 1 para maior renda e 8 para menor renda); (5) Quantidade de automóveis no domicílio; (6) Quantidade de motocicletas no domicílio; (7) Número total de viagens no domicílio; (8) Distância total de viagens realizadas no domicílio; (9) Distância de viagens realizadas por transporte individual motorizado no domicílio; (10) Distância de viagens realizadas por transporte não motorizado no domicílio; (11) Distância de viagens realizadas por transporte público no domicílio; (12) Quantidade de bicicletas no domicílio. A variável objeto do estudo (dicotômica) representa a preferência da maior parte das viagens motorizadas no domicílio: 0 – Transporte individual motorizado; 1 – Transporte Público.
2.2. Software
Este trabalho utilizou o software IBM SPSS Statistics 22 para a ACP e RLM, respectivamente.
2.3. Método
O método proposto para este trabalho é sequencial, formado por duas etapas principais: Aplicação da ACP e Calibração e Validação dos modelos logit binomial. A primeira etapa metodológica (ACP) envolve a transformação linear ótima do conjunto de dados das variáveis independentes em componentes. Para a sua aplicação, seguiu-se as seguintes etapas (Pitombo e Martins, 2014Pitombo, C. S., Gomes, M. M. (2014). Study of Work-Travel Related Behavior Using Principal Component Analysis. Open Journal of Statistics,4(11), 889.): (1) Cálculo da matriz de correlação das variáveis em estudo; (2) Extração e rotação das componentes, observando a variabilidade dos dados explicada; (3) Interpretação e nomenclatura para cada uma das componentes observando a contribuição (autovetores) de cada uma das variáveis.
A etapa de regressão logística utilizou as componentes, como variáveis independentes, e a variável preferência por modo motorizado como variável dependente. A calibração do modelo logit foi realizada com 70% da amostra, enquanto a validação foi realizada com 30% da amostra final. Vale ressaltar que foram calibradas equações logísticas de uma a n componentes principais, sendo n igual ao número de variáveis originais.
3. Resultados e discussões
3.1. Etapa 1: Aplicação da ACP
Como esperado, existem diversas variáveis correlacionadas entre si, fato que justifica a aplicação da ACP. Os valores de correlação mais altos encontrados foram entre (1) Quantidade de automóveis e critério de renda; (2) Total de viagens domiciliares e número de pessoas no domicílio; (3) Distância total de viagem e Distância de viagens realizadas por transporte público no domicílio.
Para extração das componentes mais significativas foram considerados os valores dos autovalores iguais ou superiores a um. Assim, foram extraídas cinco componentes com 68,7% da variabilidade dos dados explicada. Para finalidade de testes em modelos logísticos posteriormente, foram consideradas até 12 componentes, com 100% da variância explicada. A Tabela 1 descreve a variância explicada por cada uma das componentes.
Para se interpretar adequadamente e também dar nome às componentes, é necessário analisar os valores das cargas fatoriais (autovetores) das variáveis, ou seja, deve ser analisada a contribuição (negativa ou positiva) de cada uma das variáveis em relação a cada uma das doze componentes. Após o método de rotação Varimax, a Tabela 2 apresenta as variáveis e os valores dos autovetores em cada componente. Em seguida é apresentada a nomenclatura proposta.
(1) Componente 1: Distância de viagens realizadas por transporte público no domicílio; (2) Componente 2: Distâncias de viagens realizadas por modo individual motorizado no domicílio; (3) Renda; (4) Distâncias de viagens realizadas por modo não motorizado no domicílio; (5) Quantidade de bicicletas no domicílio; (6) Número de pessoas no domicílio; (7) Quantidade de motocicletas no domicílio; (8) Latitude; (9) Longitude; (10) Total de viagens domiciliares; (11) Quantidade de automóveis no domicílio; (12) Distâncias totais.
3.2. Etapa 2: Modelos logit binomial
Com 70% da amostra final de domicílios, foram gerados doze modelos logísticos com a inclusão de componentes por modelo. Os valores de parâmetros calibrados para cada variável independente foram considerados significativos. Além disso, a magnitude dos valores calibrados, bem como o sinal, em todos os casos, são coerentes. O poder preditivo do modelo (pseudo R2) melhora bastante quando são utilizadas 5 componentes apenas, para o caso de se considerar as componentes com autovalores maiores ou iguais a uma unidade. Quanto maior a variância explicada utilizada para calibração dos modelos, maior incremento da sua acurácia. Desta forma, para o caso do uso de doze componentes (100% da variância explicada), o pseudo R2 chega a um valor equivalente a 0,87. Uma das grandes vantagens do método é a garantia de que as componentes utilizadas (novas variáveis explicativas) não são correlacionadas entre si. A Tabela 3 apresenta principais resultados dos modelos logísticos.
Com 30% da amostra de domicílios restante, foi feita a validação dos doze modelos calibrados anteriormente. Após a obtenção de probabilidades e categorias modais estimadas, foram calculadas as taxas de acertos para todos os modelos. A Tabela 4 indica os resultados de taxa de acertos para a validação. O modelo que utiliza apenas a primeira componente (aproximadamente 23% da variância explicada) acerta 69,4%. Já o modelo que utiliza toda a variância explicada, através de doze componentes, tem 92,3% de acertos. A validação mostrou alta taxa de acertos para todos os modelos.
Conclusões
Este trabalho apresentou um método sequencial, envolvendo aplicação de Análise em Componentes Principais (ACP) e logit binomial para previsão de escolha por modo de transporte motorizado. Na primeira etapa foi realizada a extração das componentes a partir das variáveis originais socioeconômicas domiciliares, variáveis de viagens e locacionais. Em seguida, foi utilizado o método de rotação varimax para facilitar a interpretação das componentes. A partir das componentes extraídas na primeira etapa, são calibrados e validados n modelos (sendo n igual ao número total de variáveis originais e 100% da variância total explicada). Os modelos calibrados tiveram bom poder preditivo, com alta taxa de acertos, além disso, os parâmetros obtidos foram significativos e coerentes.
As componentes relativas às distâncias percorridas por cada modo motorizado foram aquelas mais significativas dos modelos. O parâmetro positivo para a componente relativa à Renda (CP3), por exemplo, faz sentido, pois, em escalara likert, a renda decresce de 1 para 8. Assim, o aumento numérico da componente Renda significa diminuição do poder aquisitivo, portanto maior probabilidade para escolhas do transporte coletivo. Já componente 11, que representa a quantidade de automóveis no domicílio, possui valor negativo no modelo calibrado, com relação inversa à escolha do transporte coletivo.
A componente que representa número total de viagens domiciliares (CP10) também possui valor negativo no modelo calibrado, com relação inversa à escolha do transporte coletivo. Sabe-se, portanto, que o automóvel viabiliza maior número de viagens domiciliares. Desta forma, pode-se afirmar que os resultados encontrados com os modelos logísticos binomiais com componentes principais corroboram com diversos resultados encontrados na literatura de escolha modal (Ahern e Tapley, 2008Ahern, A.; Tapley, N. (2008) The use of stated preference techniques to model modal choices on interurban trips in Ireland. Transportation Research Part A: Policy and Practice, vol 42, 1, 15-27. DOI: 10.1016/j.tra.2007.06.005.
https://doi.org/10.1016/j.tra.2007.06.00...
; Bhat, 1997Bhat, C. (1997). Work travel mode choice and number of non-work commute stops. Transpor-tation Research Part B: Methodological 31 (1), pp. 41–54. DOI: 10.1016/S0191-2615(96)00016-1; Brownstone et al., 2000Brownstone, D.; Bunch, D.; Train, K. (2000) Joint mixed logit models of stated and revealed preferences for alternative-fuel vehicles. Transportation Research Part B: Methodological, vol 34, 5, pp 315-338. DOI: 10.1016/S0191-2615(99)00031-4.
https://doi.org/10.1016/S0191-2615(99)00...
).
Finalmente, sabe-se que a multicolinearidade é um problema para o caso de modelos de regressão. Este trabalho apresenta uma forma de mitigar essa dificuldade, comum na modelagem de demanda por transportes, buscando melhorar estimativas relacionadas a escolhas discretas de transportes através do uso de componentes principais em modelos logísticos tradicionais.
Agradecimentos
Ao CNPq e CAPES pelo suporte financeiro fornecido à presente pesquisa.
Referências
- Aguilera, A.M.; Escabias, M.; Valderrama, M.J. (2006) Using principal components for estimating logistic regression with high-dimensional multicollinear data. Computational Statistics & Data Analysis 50, pp. 1905 – 1924. DOI:10.1016/j.csda.2005.03.011
» https://doi.org/10.1016/j.csda.2005.03.011 - Ahern, A.; Tapley, N. (2008) The use of stated preference techniques to model modal choices on interurban trips in Ireland. Transportation Research Part A: Policy and Practice, vol 42, 1, 15-27. DOI: 10.1016/j.tra.2007.06.005.
» https://doi.org/10.1016/j.tra.2007.06.005 - Bhat, C. (1997). Work travel mode choice and number of non-work commute stops. Transpor-tation Research Part B: Methodological 31 (1), pp. 41–54. DOI: 10.1016/S0191-2615(96)00016-1
- Brownstone, D.; Bunch, D.; Train, K. (2000) Joint mixed logit models of stated and revealed preferences for alternative-fuel vehicles. Transportation Research Part B: Methodological, vol 34, 5, pp 315-338. DOI: 10.1016/S0191-2615(99)00031-4.
» https://doi.org/10.1016/S0191-2615(99)00031-4 - Camminatiello, I.; Lucadamo, A. (2010) Estimating Multinomial Logit Model with Multicollinear Data. Asian Journal of Mathematics & Statistics, vol 3,2, 93-101. DOI: 10.3923/ajms.2010.93.101.
» https://doi.org/10.3923/ajms.2010.93.101 - Chen, Z. J., Cheng, L., Deng, H. N., Zhang, J. K (2010) Analyzing Residential Travel Mode Choice Based on Principal Component Analysis. Proceedings of the 10th International Conference of Chinese Transportation Professionals. 2010.
- Frank, I.E.; Friedman, J.H.; Wold, S.; Hastie, T.; Mallows, C. (1993) A statistical view of some chemometrics regression tools. Technometrics, c, pp 109-148.
- Hair Jr, J. F.; Black, W. C.; Babin, B. J.; Anderson, R. E. (2010) Multivariate Data Analysis. Prentice Hall. 7a ed. 785 p.
- Hoerl, A.; Kennard, R.W. (1970) Biased Estimation for Nonorthogonal Problems. Technometrics, Vol. 12, No. 1 (Feb., 1970), pp. 55-67.
- Jolliffe, I. T. (2002) Principal Component Analysis. 2a ed. Springer 2002. 518p.
- Neter, J.; Wasserman, W.; Kutner, M.H. (1989). Applied Linear regression Models. 2nd Edition. Irwin Homewood IL.
- Pitombo, C. S., Gomes, M. M. (2014). Study of Work-Travel Related Behavior Using Principal Component Analysis. Open Journal of Statistics,4(11), 889.
- Ortúzar, J. D.; Willumsen, L. G. (2011) Modelling Transport. Londres: Wiley. 4ª ed. 586p.
- Wold, H. (1985) Partial Least Squares, Encyclopedia of statistical sciences, vol 6. Pp 581-591.
Datas de Publicação
-
Publicação nesta coleção
Set 2016
Histórico
-
Recebido
09 Jul 2015 -
Aceito
12 Jul 2015