Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Effect of changing data collection parameters on statistical motor unit number estimates

The effect of number of samples and selection of data for analysis on the calculation of surface motor unit potential (SMUP) size in the statistical method of motor unit number estimates (MUNE) was determined in 10 normal subjects and 10 with amyotrophic lateral sclerosis (ALS). We recorded 500 sequential compound muscle action potentials (CMAPs) at three different stable stimulus intensities (10–50% of maximal CMAP). Estimated mean SMUP sizes were calculated using Poisson statistical assumptions from the variance of 500 sequential CMAP obtained at each stimulus intensity. The results with the 500 data points were compared with smaller subsets from the same data set. The results using a range of 50–80% of the 500 data points were compared with the full 500. The effect of restricting analysis to data between 5–20% of the CMAP and to standard deviation limits was also assessed. No differences in mean SMUP size were found with stimulus intensity or use of different ranges of data. Consistency was improved with a greater sample number. Data within 5% of CMAP size gave both increased consistency and reduced mean SMUP size in many subjects, but excluded valid responses present at that stimulus intensity. These changes were more prominent in ALS patients in whom the presence of isolated SMUP responses was a striking difference from normal subjects. Noise, spurious data, and large SMUP limited the Poisson assumptions. When these factors are considered, consistent statistical MUNE can be calculated from a continuous sequence of data points. A 2 to 2.5 SD or 10% window are reasonable methods of limiting data for analysis. Muscle Nerve 27: 320–331, 2003

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Usefulness of B-type natriuretic peptide in hypertensive patients with exertional dyspnea and normal left ventricular ejection fraction and correlation with new echocardiographic indexes of systolic and diastolic function

B-type natriuretic peptide (BNP) levels increase in systolic heart failure (HF). However, the value of BNP in hypertensive patients with suspected diastolic HF (symptoms suggestive of HF but normal ejection fraction) and its relation to myocardial function in these patients is unclear. We prospectively studied 72 ambulatory hypertensive subjects (40 women, mean age 58 +/- 8 years) with exertional dyspnea and ejection fraction greater than or equal to50%. Diastolic function was evaluated with transmitral and pulmonary venous Doppler, mitral annular velocities (pulsed-wave tissue Doppler), and flow propagation velocity (color M-mode). Systolic function was assessed with strain and strain rate derived from color tissue Doppler imaging. BNP was related to myocardial function and the presence or absence of global diastolic dysfunction. By conventional Doppler criteria, 34 patients had normal left ventricular diastolic function and 38 had isolated diastolic dysfunction. BNP values were higher in patients with diastolic dysfunction (46 +/- 48 vs 20 +/- 20 pg/ml, p = 0.004) and were related independently to blood pressure, systolic strain rate, left atrial function (p < 0.01 for all), and age (p = 0.015). Patients with diastolic dysfunction and pseudonormal filling had higher BNP levels compared with impaired relaxation (89 +/- 47 vs 35 +/- 42 pg/ml, p = 0.001). However, 79% of patients with diastolic dysfunction had BNP levels within the normal range. We conclude that in ambulatory hypertensive patients with symptoms suggestive of mild HF and normal ejection fraction, BNP is related to atrial and ventricular systolic parameters, blood pressure, and age. Although elevated in the presence of diastolic dysfunction, the BNP level mostly is in the normal range and, therefore, has limited diagnostic value in stable patients with suspected diastolic HF. (C) 2003 by Excerpta Medica, Inc.

Poor correlation of mid-femoral measurements by CT and hip measurements by DXA in the elderly

Background and aims: Hip fracture is a devastating event in terms of outcome in the elderly, and the best predictor of hip fracture risk is hip bone density, usually measured by dual X-ray absorptiometry (DXA). However, bone density can also be ascertained from computerized tomography (CT) scans, and mid-thigh scans are frequently employed to assess the muscle and fat composition of the lower limb. Therefore, we examined if it was possible to predict hip bone density using mid-femoral bone density. Methods: Subjects were 803 ambulatory white and black women and men, aged 70-79 years, participating in the Health, Aging and Body Composition (Health ABC) Study. Bone mineral content (BMC, g) and volumetric bone mineral density (vBMD, mg/cm(3)) of the mid-femur were obtained by CT, whereas BMC and areal bone mineral density (aBMD, g/cm(2)) of the hip (femoral neck and trochanter) were derived from DXA. Results: In regression analyses stratified by race and sex, the coefficient of determination was low with mid-femoral BMC, explaining 6-27% of the variance in hip BMC, with a standard error of estimate (SEE) ranging from 16 to 22% of the mean. For mid-femur vBMD, the variance explained in hip aBMD was 2-17% with a SEE ranging from 15 to 18%. Adjusting aBMD to approximate volumetric density did not improve the relationships. In addition, the utility of fracture prediction was examined. Forty-eight subjects had one or more fractures (various sites) during a mean follow-up of 4.07 years. In logistic regression analysis, there was no association between mid-femoral vBMD and fracture (all fractures), whereas a 1 SD increase in hip BMD was associated with reduced odds for fracture of similar to60%. Conclusions: These results do not support the use of CT-derived mid-femoral vBMD or BMC to predict DXA-measured hip bone mineral status, irrespective of race or sex in older adults. Further, in contrast to femoral neck and trochanter BMD, mid-femur vBMD was not able to predict fracture (all fractures). (C) 2003, Editrice Kurtis.

Análise Imuno-histoquímica de hipóxia renal em modelo de isquemia/reperfusão tratado com Fator Estimulador de Colônia de Granulócitos

O rim demonstra uma capacidade singular em reparar-se após danos locais, no entanto, depois de acometido, as chances de desenvolvimento de lesões renais elevam-se. A patofisiologia da isquemia/reperfusão (IR) é complexa porque há ocorrência simultânea de danos celulares e inflamação. O decréscimo na quantidade de oxigênio requer um sistema capaz de evitar seus efeitos prejudiciais e uma maquinaria molecular HIF (Hypoxia Inducible Factor), um complexo, atua como fator de transcrição de diversos genes desde os da regulação da proliferação celular e apoptose até a sinalização para angiogênese. O Fator Estimulador de Colônia de Granulócitos (G-CSF) é uma glicoproteína conhecida pela sua capacidade de promover a sobrevivência, proliferação e diferenciação de células estimulando a recuperação aos efeitos advindos da IR. Com o intuito de observar as influências dessas proteínas foi realizada uma análise semi-quantitativa de amostras renais submetidas ou não à IR, usando-se descrições microscópicas morfológicas e imunohistoquímicas, com os cálculos e gráficos estatísticos foram feitos no software GraphPad Prism®. Das análises morfológicas, constatou-se que as lesões características de IR foram observadas em espécimes não tratados: bolhas em epitélio tubular; vacuolização citoplasmática, distalização tubular e congestão luminal. De forma análoga, foi encontrada nos tratados, contudo em estágios menos avançados e em animais controle, não foi houve esta diferença tissular. As análises de microscopia eletrônica demonstraram alteração na barreira filtrante com concomitante perda de outras características glomerulares. Aos animais controle foi observada a arquitetura típica, ao passo que para os animais tratados notou-se conservação da barreira. A presença de HIF-1α nos rins contralaterais demonstrouse significante quando comparadas às amostras isquêmicas e tratadas (p<0,05). Já a ocorrência da mesma proteína em rins isquêmicos não apresentou qualquer diferença. Analisando-se a proteína VEGF foi comprovado que em rins contralaterais não há diferença estatística, contudo nos rins esquerdos há significância entre os três grupos (p<0,05). Já a correlação entre estas duas proteínas não se mostrou estatisticamente significante. Em relação às atividades de proliferação e morte celulares, todos os três grupos foram significantes entre si (p<0,05). Ao que concerne o tratamento, foi demonstrada a atividade protetora do medicamento e uma possível interação molecular com a HIF, enquanto que a ativação desta proteína corrobora sua rota metabólica já previamente descrita.

Controle estatístico de processo aplicado no monitoramento da absorção de água por carcaças de frango durante a etapa de pré-resfriamento

A absorção de água por carcaças de frango na etapa de pré-resfriamento da linha abate representa uma característica de qualidade importante relacionada ao rendimento do produto final. Uma forma de manter o padrão de qualidade de um produto é garantir que as etapas do processo sejam estáveis e replicáveis. Ao empregar o Controle Estatístico de Processo (CEP) é possível obter estabilidade e melhorias nos processos, por meio da redução da variabilidade. Neste contexto, o objetivo deste trabalho foi a aplicação de gráficos de controle, análise de correlação, estatística descritiva, testes de hipóteses e regressão linear múltipla na linha de abate de um abatedouro-frigorífico de aves para monitorar a variabilidade da absorção de água pelas carcaças de frango após a etapa de pré-resfriamento. Como resultado, verificou-se que o teor de absorção de água das carcaças de frango apresentou elevada variabilidade, sendo que 10% (8/80) das carcaças apresentaram absorção de água superior ao limite de 8% definido pela legislação brasileira. Do total de 16 variáveis de entrada analisadas, as mais impactantes no teor de absorção de água foram o “tempo de retenção da carcaça no pré-chiller” e o “tempo de espera da carcaça após a etapa de gotejamento”. Entretanto, o modelo de regressão obtido apresentou baixa correlação (R²=0,16) que foi associada à elevada variabilidade da variável-resposta. Os resultados da estatística descritiva demonstraram que as variáveis de entrada também apresentaram elevada variabilidade, com coeficiente de variação entre 7,95 e 63,5%. Verificou-se, pela análise dos gráficos de controle de medida individual e da amplitude móvel, que 15 das 16 variáveis de entrada se apresentaram fora de controle estatístico assim como a variável-resposta. Baseado no fluxograma e na descrição das etapas da linha de abate, previamente realizados, atribuiu-se à falta de padronização na condução das etapas e de procedimentos para o controle de qualidade das operações na linha de abate como fatores relevantes que poderiam estar associados à presença de causas especiais no processo. Concluiu-se que para reduzir a elevada variabilidade das variáveis e eliminar as causas especiais presentes são necessários ajustes operacionais para, dessa forma, obter um processo mais estável e mais uniforme garantindo o padrão de qualidade das carcaças de frango em relação ao teor de absorção de água.

Variabilidade genética em famílias de meios-irmãos de cafeeiro robusta

O cultivo do café é uma das atividades do agronegócio de maior importância socioeconômica dentre as diferentes atividades ligadas ao comércio agrícola mundial. Uma das maiores contribuições da genética quantitativa para o melhoramento genético é a possibilidade de prever ganhos genéticos. Quando diferentes critérios de seleção são considerados, a predição de ganhos referentes a cada critério tem grande importância, pois indica os melhoristas sobre como utilizar o material genético disponível, visando obter o máximo de ganhos possível para as características de interesse. O presente trabalho foi instalado em julho de 2004, na Fazenda Experimental de Bananal do Norte, conduzida pelo Incaper, no distrito de Pacotuba, município de Cachoeiro de Itapemirim, região Sul do Estado, com o objetivo de selecionar as melhores plantas entre e dentro de progênies de meios- irmãos de Coffea canephora, por meio de diferentes critérios de seleção. Foram realizadas análises de variância individuais e conjuntas para 26 progênies de meios- irmãos Coffea canephora. O delineamento experimental utilizado foi em blocos ao acaso com quatro testemunhas adicionais com quatro repetições e parcela composta por cinco plantas, com o espaçamento de 3,0 m x 1,2 m. Neste trabalho, considerou-se os dados das últimas cinco colheitas. As características mensuradas foram: florescimento, maturação, tamanho do grão, peso, porte, vigor, ferrugem, mancha cercóspora, seca de ponteiros, escala geral, porcentagem de frutos boia e bicho mineiro. Todas as análises estatísticas foram realizadas com o aplicativo computacional em genética e estatística (GENES). Foram estimados os ganhos de seleção em função da porcentagem de seleção de 20% entre e dentro, sendo as mesmas mantidas para todas as características. Todas as características foram submetidas a seleção no sentido positivo, exceto para florescimento, porte, ferrugem, mancha cercóspora, seca de ponteiros, porcentagem de frutos boia e bicho mineiro, para obter decréscimo em suas médias originais. Os critérios de seleção estudados foram: seleção convencional entre e dentro das famílias, índice de seleção combinada, seleção massal e seleção massal estratificada. Esta dissertação é composta por dois capítulos, em que foram realizadas análises biométricas, como a obtenção de estimativas de parâmetros genéticos. Na maioria das características estudadas, verificaram-se diferenças significativas (P<0,05) para genótipos que, associados aos coeficientes de variação genotípicos e também ao coeficiente de determinação genotípico e à relação CVg/CVe, indicam a existência de variabilidade genética nos materiais genéticos para a maioria das características e condições favoráveis para obtenção de ganhos genéticos pela seleção. Essas características também foram correlacionadas. Os dados foram submetidos às análises de variância e multivariada, aplicando-se a técnica de agrupamento e UPGMA, teste de médias e estudo de correlações. Na técnica de agrupamento, foi utilizada a distância generalizada de Mahalanobis como medida de dissimilaridade, e na delimitação dos grupos, o método de Tocher. Foi encontrada diversidade genética para as características associadas à qualidade fisiológica, mobilização de reserva das sementes, dimensões e biomassa das plântulas. Quatro grupos de genótipos puderam ser formados. Peso de massa seca de sementes, redução de reserva de sementes e peso de massa seca de plântulas estão positivamente correlacionados entre si, enquanto a redução de reserva das sementes e a eficiência na conversão dessas reservas em plântulas estão negativamente correlacionadas. De acordo com os resultados obtidos, verificou-se que todas as características apresentaram níveis diferenciados de variabilidade genética e os critérios de seleção utilizados mostraram-se eficientes para o melhoramento, no qual o índice de seleção combinada é o critério de seleção que apresentou os melhores resultados em termos de ganhos, sendo indicado como critério mais apropriado para o melhoramento genético da população estudada. Nos estudos de correlações, em 70% dos casos, a correlação fenotípica foi superior à genotípica, mostrando maior influência dos fatores ambientais em relação aos genotípicos e condições propícias ao melhoramento dos diferentes caracteres. No estudo de divergência genética, observou-se que pelo agrupamento de genótipos, pela técnica de Tocher, indicou que os genótipos foram distribuídos em três grupos.

Cenários de mudanças climáticas para a América do Sul : impacto no zoneamento agroclimático da cana-de-açúcar

O Painel Intergovernamental sobre Mudanças Climáticas (IPCC) através do seu Quarto Relatório de Avaliação das Mudanças Climáticas Globais (IPCC-AR4), publicado em 2007, atribui as emissões de gases de efeito estufa como a principal causa do aumento médio das temperaturas e alerta para uma elevação entre 1,8 ºC e 6,4 ºC até 2100, podendo modificar assim a aptidão climática para as culturas agrícolas em diversas regiões do planeta. Diante disso, existe a necessidade de substituição dos combustíveis fósseis por fontes renováveis e limpas de energia, como o etanol. A cana-de-açúcar apresenta-se, portanto, como uma cultura estratégica na produção do etanol. O presente trabalho teve como objetivos: 1) avaliar o desempenho dos Modelos Climáticos Globais (MCGs) do IPCC-AR4 na simulação de dados climáticos de temperatura do ar e precipitação pluviométrica para o período anual e mensal; 2) elaborar o zoneamento agroclimático da cana-de-açúcar para a América do Sul considerando o clima referência e o futuro para as décadas de 2020, 2050 e 2080 em função do cenário de emissão A1B considerado pessimista e que usa um equilíbrio entre todas as fontes de energia. Para a avaliação do desempenho dos MCGs, foram utilizados dados climáticos médios mensais observados de precipitação e temperatura do ar provenientes do Climatic Research Unit (CRU) e dados simulados oriundos dos 22 MCGs do IPCC (cenário 20c3m) compreendidos entre o período de 1961-1990, além do Multimodel (ensemble) – MM que é a média da combinação dos dados de todos os modelos. O desempenho dos MCGs foi avaliado pelos índices estatísticos: desvio padrão, correlação, raiz quadrada da média do quadrado das diferenças centralizadas e o “bias” dos dados simulados com os observados, que foram representados no diagrama de Taylor. Para a etapa da elaboração do zoneamento agroclimático procedeu-se o cálculo dos balanços hídricos (referência e futuros) da cultura, pelo método de Thornthwaite & Mather (1955). Para o cenário referência, utilizaram-se dados das médias mensais da precipitação e temperatura provenientes do CRU, enquanto que para as projeções futuras, dados provenientes das anomalias do Multimodel (ensemble) – MM para as décadas de 2020, 2050 e 2080, que foram ajustados, obtendo-se assim as projeções futuras para cada período analisado. Baseado nos mapas temáticos reclassificados de deficiência hídrica anual, temperatura média anual, excedente hídrico anual e no índice de satisfação das necessidades de água (ISNA), realizou-se uma sobreposição dessas informações obtendo assim, os mapas finais do zoneamento agroclimático da cana-de-açúcar. Posteriormente ao zoneamento, realizou-se a análise das transições (ganhos, perdas e persistências) entre as classes de aptidão climática da cultura. Os resultados mostram que o Multimodel (ensemble) – MM para o período mensal apresenta o melhor desempenho entre os modelos analisados. As áreas inaptas correspondem a maior parte da América do Sul e uma expressiva transição entre as classes de aptidão climática da cultura.

Estudo sobre acetilação da isoniazida em pacientes com tuberculose pulmonar e da sua implicação na redução ou eliminação da carga bacilar no escarro

A N-acetiltransferase 2 é a principal enzima responsável pelo metabolismo e inativação da isoniazida no organismo humano. Mutações no gene NAT2 levam a 3 perfis genotípicos de acetilação que alteram os níveis séricos do fármaco: acetiladores lentos, intermediários e rápidos, o que pode alterar o desfecho terapêutico. O objetivo do estudo foi investigar se os diferentes perfis podem influenciar no tempo de negativação da cultura de escarro, e se existe correlação entre carga bacilar e gravidade da doença com tempo de conversão da cultura. A população de estudo foi composta por 62 pacientes, que tiveram seus DNAs sequenciados para identificação de mutações no gene NAT2 e seus perfis de acetilação determinados. A análise genotípica detectou 10 SNPs, sendo as mutações 341 T>C (39,65%) e 481 C>T (38,71%) as mais frequentes. A determinação das variantes alélicas identificou NAT2*5B (29,03%), NAT2*6A (23,39%) e NAT2*4 (24,19%) como os alelos mais frequentes e NAT2*5B/*5B como o genótipo mais frequente (20,4%). Dentre os 62 pacientes, foi possível correlacionar tempo de negativação da cultura e perfil de acetilação entre 43 deles, os quais 58,3% e 55,6% tiveram o genótipo lento com maior frequência no mês 1 e mês 3, respectivamente. Por meio de dados microbiológicos, a carga bacilar e a gravidade da doença também foram comparadas com o tempo de negativação, indicando que os pacientes com doença moderada ou avançada (76,7%) e aqueles com carga bacilar alta (60,4%), não tiveram associação estatística com o tempo de conversão da cultura. Por último, curvas de crescimento de isolados de M. tuberculosis de pacientes foram construídas para verificar possíveis diferenças na duração da fase lag entre os isolados, porém não foi observada diferença estatística entre elas. Com base nos resultados encontrados, verifica-se que não existe associação entre o perfil de acetilação do paciente, a carga bacilar, a gravidade da doença e o tempo de negativação da cultura de escarro

Correlation between grain nutritional content and pasting properties of pre-gelatinized red rice flour

As with any variety of rice, red rice characteristics are subject to varietal differences, growing conditions, types of processing, and nutritional and rheological properties. This study determined the nutritional characteristics (centesimal composition and minerals) and paste viscosity properties of raw grains of four red rice genotypes (Tradicional MNAPB0405, MNACE0501 and MNACH0501) and the paste viscosity properties of pre-gelatinized flours obtained at different cooking times (20, 30 and 40 min). The main nutritional properties were correlated with the pasting properties of the pre-gelatinized flours. The samples showed differences in nutritional properties and paste viscosity. MNAPB0405 and MNACE0501 showed higher levels of fiber and fat and provided higher caloric energy than Tradicional and MNACH0501, which, in turn, showed higher levels of amylose. MNACH0501 showed higher peak viscosity (2402 cP), higher breakdown viscosity (696 cP) and a greater tendency to retrogradation (1510 cP), while Tradicional, MNAPB0405 and MNACE0501 had pasting profiles with peak viscosities varying between 855 and 1093 cP, breaking viscosity below 85 cP and retrogradation tendency between 376 and 1206 cP. The factors genotype and cooking time influenced the rheological behavior of pre-gelatinized flours, decreasing their pasting properties. The protein and amylose levels are correlated with the pasting properties and can be used as indicators of these properties in different genotypes of red rice, whether raw or processed into pre-gelatinized flours.

Comparative analysis between statistical and artificial intelligence models in business failure prediction

A growing number of predicting corporate failure models has emerged since 60s. Economic and social consequences of business failure can be dramatic, thus it is not surprise that the issue has been of growing interest in academic research as well as in business context. The main purpose of this study is to compare the predictive ability of five developed models based on three statistical techniques (Discriminant Analysis, Logit and Probit) and two models based on Artificial Intelligence (Neural Networks and Rough Sets). The five models were employed to a dataset of 420 non-bankrupt firms and 125 bankrupt firms belonging to the textile and clothing industry, over the period 2003–09. Results show that all the models performed well, with an overall correct classification level higher than 90%, and a type II error always less than 2%. The type I error increases as we move away from the year prior to failure. Our models contribute to the discussion of corporate financial distress causes. Moreover it can be used to assist decisions of creditors, investors and auditors. Additionally, this research can be of great contribution to devisers of national economic policies that aim to reduce industrial unemployment.

Comparative Analysis between Statistical and Artificial Intelligence Models in Business Failure Prediction

A growing number of predicting corporate failure models has emerged since 60s. Economic and social consequences of business failure can be dramatic, thus it is not surprise that the issue has been of growing interest in academic research as well as in business context. The main purpose of this study is to compare the predictive ability of five developed models based on three statistical techniques (Discriminant Analysis, Logit and Probit) and two models based on Artificial Intelligence (Neural Networks and Rough Sets). The five models were employed to a dataset of 420 non-bankrupt firms and 125 bankrupt firms belonging to the textile and clothing industry, over the period 2003–09. Results show that all the models performed well, with an overall correct classification level higher than 90%, and a type II error always less than 2%. The type I error increases as we move away from the year prior to failure. Our models contribute to the discussion of corporate financial distress causes. Moreover it can be used to assist decisions of creditors, investors and auditors. Additionally, this research can be of great contribution to devisers of national economic policies that aim to reduce industrial unemployment.

A relevância do relatório de auditoria para os investidores institucionais, analistas financeiros e analistas de crédito

In the 70s, a new line of research focused on the study of the influence of the audit report on the decision process of investors, financial analysts and credit analysts. Notwithstanding the numerous studies that have been carried out, results have not been consistent. Given the above, and considering the lack, in Portugal, of a research of this nature, it seems urgent to carry out a study that allows the analysis of the use of the audit report, as well as its influence on the decision making process of Portuguese stakeholders. For that purpose, in the light of the positivist research paradigm, a questionnaire was designed, which was administered by mail and on the Survey Monkey platform to a sample of institutional investors, financial analysts and credit analysts. The statistical analysis of the data obtained was undertaken with resource to the Statistical Package for the Social Sciences and SmartPLS 2.0. Corroborating the literature review and the assumptions of the Agency Theory and the Stakeholder Theory, used in the theoretical framework of analysis, empirical evidence has shown that the audit report influences the decision of institutional investors, financial analysts and credit analysts, and that the opinion expressed in that document is the most determinant factor of this influence. In addition to this factor, it was found that the degree of utilization of the audit report, as well as the value ascribed to this document, determine its influence in the decision process of research groups studied. Only in the case of institutional investors, the results did not reveal a correlation between the utility ascribed to the audit report and the influence of this document in their decision making process. In turn, the statistical inference of the model explaining the degree of use of the audit report revealed that it is conditioned by the perceived quality of the information enclosed in the audit report, the utility assigned to the audit report on the decision process, as well as the relevance of the other sources of information used by stakeholders. Therefore, this study allowed proving the importance of the audit report to its users. As a result, we believe to have filled a gap in national literature and to have contributed to the enhancement of international literature. The importance that this document has for the development of any country is, therefore, shown, and it is urgent to maintain rigor in the selection of its staff, in the development of its standards, and especially in the development of audits. Moreover, we also consider that this research may contribute to the improvement of the audit report, insofar as it will help professional bodies to understand the information needs and perceptions of stakeholders.