Análise AMMI com dados imputados em experimentos de interação genótipo x ambiente de algodão

O objetivo deste trabalho foi avaliar a conveniência de definir o número de componentes multiplicativos dos modelos de efeitos principais aditivos com interação multiplicativa (AMMI) em experimentos de interações genótipo x ambiente de algodão com dados imputados ou desbalanceados. Um estudo de simulação foi realizado com base em uma matriz de dados reais de produtividade de algodão em caroço, obtidos em ensaios de interação genótipo x ambiente, conduzidos com 15 cultivares em 27 locais no Brasil. A simulação foi feita com retiradas aleatórias de 10, 20 e 30% dos dados. O número ótimo de componentes multiplicativos para o modelo AMMI foi determinado usando o teste de Cornelius e o teste de razão de verossimilhança sobre as matrizes completadas por imputação. Para testar as hipóteses, quando a análise é feita a partir de médias e não são disponibilizadas as repetições, foi proposta uma correção com base nas observações ausentes no teste de Cornelius. Para a imputação de dados, foram considerados métodos usando submodelos robustos, mínimos quadrados alternados e imputação múltipla. Na análise de experimentos desbalanceados, é recomendável escolher o número de componentes multiplicativos do modelo AMMI somente a partir da informação observada e fazer a estimação clássica dos parâmetros com base nas matrizes completadas por imputação.

Curvas de crescimento em vacas de corte de diferentes tipos biológicos

O objetivo deste trabalho foi selecionar o modelo de curvas de crescimento mais adequado e avaliar a influência de efeitos de ambiente e de grupo genético sobre os parâmetros estimados do modelo. Cinco modelos não lineares, Brody, Gompertz, Logístico, Von Bertalanffy e Richards, foram ajustados a dados de peso-idade coletados de 316 vacas, de quatro grupos genéticos: G (Nelore, ½Canchim + ½Nelore, ½Angus + ½Nelore e ½Simental + ½Nelore), do nascimento até 100 meses de idade; em duas estações do ano: E (primavera e outono). As vacas foram submetidas a dois níveis de concentrado (S) durante quatro meses, pós-desmama. O ajuste dos modelos foi realizado por mínimos quadrados ordinários, usando os pesos ponderado e não ponderado pelo inverso da variância. Os modelos Brody e Von Bertalanffy convergiram para todos os grupos genéticos; porém, o Brody foi o mais adequado. As estimativas do peso assintótico (A) e da taxa de maturação (k) do modelo Brody ponderado pelo inverso da variância foram analisadas por modelo misto, que incluiu efeito médio global e efeitos principais de G, E e S, e suas interações. O parâmetro A foi influenciado pelo efeito de G e E, enquanto k foi influenciado por S, o que indica que melhorias no manejo alimentar resultam em menor variação na forma das curvas de crescimento e em altas taxas de maturação.

DESCRIÇÃO DA CURVA DE CRESCIMENTO DE FRUTOS DO CAJUEIRO POR MODELOS NÃO LINEARES

RESUMO O objetivo do trabalho foi avaliar o ajuste dos modelos Gompertz e Logístico, com estrutura de erros independentes e autoregressivos, no desenvolvimento de frutos de caju, com base em medidas de comprimento e largura do fruto, tomados ao longo do tempo. A estimação dos parâmetros foi feita por meio de rotinas no software R, utilizando-se o método dos mínimos quadrados e o processo iterativo de Gauss- Newton. Os modelos foram comparados pelos critérios: coeficiente de determinação ajustado (R2 aj), desvio padrão residual (DPR), critério de informação Akaike (AIC) ecritériobayesiano de Schwarz (BIC). Para os dois modelos, os dados apresentaram autocorrelação residual positiva, tanto para a variável comprimento como para a largura do fruto, descrita por processo autoregressivo de primeiraordem. Em geral, por todos critérios de avaliação da qualidade de ajuste, os dados se ajustaram ao modelo Logístico com uma estrutura autoregressiva da primeira ordem, havendo no entanto superestimação do tamanhodo fruto nas últimas idades, tanto no crescimento em comprimento (cm) e largura (cm).

Ajuste de parâmetros de interação em modelos termodinâmicos

Para misturas nas quais existem interações entre moléculas diferentes aparece um parâmetro energético de interação ajustável na relação das variáveis. As grandezas termodinâmicas de excesso dão informação qualitativa dos possíveis tipos de interação. Duas grandezas importantes neste sentido são o volume de excesso e a entalpia de excesso. Neste artigo se estudam diferentes métodos, regressão com mínimos quadrados e por minimização do valor absoluto e relativo dos resíduos, para determinar a possibilidade de aplicação dos mesmos no cálculo de parâmetros de interação molecular em modelos termodinâmicos de soluções. Usam-se dados experimentais de volume molar de excesso para todo o intervalo de frações molares de 10 sistemas binários formados por acetonitrila + um álcool primário (do metanol ao n-decanol). Comparam-se os valores calculados e os dados experimentais das grandezas, e são feitas análises dos resultados.

Comparação do ajuste dos modelos de morgan-mercer-flodin e de johnson-mehl-avrami na reação de precipitação na liga cu-3%al-5%ag

Em geral, os modelos de crescimento empregados na descrição de processos cinéticos do estado sólido utilizam as funções de Morgan-Mercer-Flodin (MMF) e Johnson-Mehl-Avrami (JMA). Neste trabalho comparou-se o comportamento das estimativas dos parâmetros dessas funções, em duas parametrizações, quando ajustadas aos dados experimentais de variação isotérmica da microdureza da liga Cu-3%Al-5%Ag com o tempo de envelhecimento, em cinco temperaturas diferentes. Na estimativa dos parâmetros aplicou-se o método dos mínimos quadrados à função linearizada do modelo e para o refinamento da solução, os procedimentos da regressão não linear. Foram obtidas as distribuições de freqüências das estimativas dos parâmetros, por simulação, e determinados seus erros relativos. Observou-se que, em uma das parametrizações das funções, a estimativa do parâmetro cinético apresentou maior estabilidade. Apesar dessas funções serem anteriormente consideradas distintas, verificou-se que a função de MMF é uma aproximação da função de JMA, sendo que esta, com a parametrização adequada, é a mais indicada para descrever o processo cinético considerado.

Variabilidade espacial de Colletotrichum truncatum em campo de soja sob três níveis de sanidade de sementes

Objetivou-se caracterizar a variabilidade espacial da incidência de Colletotrichum truncatum (Schwein) em sementes colhidas de soja e determinar o melhor modelo e método de semivariograma que represente a incidência de C. truncatum dessas sementes. O experimento foi realizado em condições de campo na safra 2009/10, em 3 parcelas de 9,9 x 10 m, com sementes inoculadas com C.truncatum. O inóculo correspondeu a 0,8; 1,6 e 2,4% do total semeado. Foram demarcadas 3 malhas com receptores GNSS, totalizando 112 pontos em cada parcela distanciados a 1,3 m na linha. No final do ciclo da soja, realizaram-se a colheita, a secagem e a análise sanitária das sementes pelo método 'blotter test', referente aos 336 pontos demarcados. Quatro modelos de semivariogramas foram ajustados aos dados coletados utilizando os métodos mínimos quadrados ordinários (OLS), mínimos quadrados ponderados (WLS), máxima verossimilhança (ML) e máxima verossimilhança restrita (REML). A validação cruzada foi empregada para escolha final do modelo e método do semivariograma. Em seguida, efetuaram-se a krigagem e o desvio padrão da krigagem. Os mapas de krigagem ilustraram a transmissão da semente para a semente e a sua variância. Verificou-se estrutura de dependência espacial da transmissão de C. truncatum via semente. O melhor modelo de semivariograma foi o esférico e o melhor ajuste foi o REML. Houve alcance de 0,95, 4,03 e 7,05 m para as parcelas com 0,8, 1,6 e 2,4% de sementes inoculadas respectivamente. Quanto maior o inóculo primário da parcela, maior foi a transmissão para as sementes próximas à fonte de inóculo.

Uso de análise exploratória de dados e de regressão robusta na avaliação do crescimento de espécies comerciais de terra firme da Amazônia

Esta pesquisa originou-se de uma base de dados de dois inventários a 100%, realizados em 1984 e 2000 sobre uma mesma área de 576 ha de floresta tropical primária, localizada na Floresta Nacional do Tapajós, Belterra, Pará. O objetivo da pesquisa foi utilizar a análise exploratória de dados e a regressão robusta para modelar o crescimento em diâmetro e área basal. Nos dois inventários as circunferências à altura do peito (CAP) foram medidas com fitamétrica, enquanto as alturas comerciais em 1984 foram medidas com o hipsômetro de Weiss e em 2000, estimadas com o uso de varas; o DAP mínimo em 1984 foi de 55 cm para todas as espécies e em 2000 foi adotado o DAP mínimo de 35 cm. As análises estatísticas foram precedidas de análise exploratória de dados (AED), em que foram utilizados o box plot (caixa-de-bigodes) na detecção de outliers (observações discrepantes) e o gráfico stem-and-leaf (tronco-e-folhas) para filtrar as observações extremas. Utilizou-se a regressão robusta para ajustar os modelos na presença de outliers. A pesquisa mostrou que, apesar das variações intra e interespecíficas, as seis espécies, em conjunto, apresentaram taxas de crescimento medianas observadas e esperadas em diâmetro e área basal de 4,56 e 4,28 mm/ano e 13,00 e 13,09 cm²/ha/ano, respectivamente. A pesquisa também mostrou que o uso da análise exploratória de dados e da regressão robusta viabilizou a análise e a determinação dos incrementos periódicos em diâmetro e área basal em bases consistentes. A metodologia empregada no caso específico dos dados disponíveis e tipos de variáveis mostrou ser mais eficiente com o uso de regressão robusta pelo método Least Trimmed Square do que pelo método dos mínimos quadrados ordinários.

A função de produção da indústria brasileira de celulose, em 2004

O segmento de celulose e papel tem-se destacado como o mais importante do setor florestal no Brasil, e o conhecimento das relações de produção e das características da tecnologia pode gerar informações úteis para o setor. Este trabalho teve como objetivos especificar e ajustar a função de produção da indústria de celulose no Brasil, no ano de 2004, utilizando-se como referência os fundamentos da teoria da produção. Foi empregado um modelo Cobb-Douglas, em que as variáveis explicativas foram terra, trabalho e capital. As estimativas confirmaram que essas variáveis afetam, significativamente, a produção de celulose. Os resultados, usando o método dos Mínimos Quadrados Ordinários (MQO), apontaram que os fatores de produção analisados estão sendo empregados racionalmente na indústria de celulose no Brasil. Essa indústria apresenta retorno constantes à escala, terra e trabalho constituem nos fatores mais importantes na formação da quantidade produzida de celulose e as empresas mostram-se intensivas em capital.

Análise do mercado brasileiro de celulose, 1969-2005

O segmento de celulose e papel contribui de forma significativa para a economia brasileira. O conhecimento da estrutura do mercado de celulose é fundamental para fins de planejamento da produção, comercialização e previsão para o setor. Este trabalho teve como objetivo principal especificar e estimar a oferta e demanda brasileira de celulose, bem como a oferta e a demanda de exportação, no período de 1969 a 2005. Os resultados, usando o método dos Mínimos Quadrados de Dois Estágios (MQ2E), indicaram que a oferta e a demanda interna de celulose, bem como a oferta e a demanda de exportação, são preço-inelástico. A demanda interna é inelástica com relação à renda. Já a demanda de exportação é elástica com relação à renda. A oferta interna é mais sensível a variações no preço das máquinas e equipamentos. A oferta de exportação é sensível a variações na produção e no consumo interno e a demanda de exportação, a variações na renda dos países importadores.

Influência do modelo de análise estatística e da forma das parcelas experimentais na seleção de clones de Eucalyptus spp.

Foram avaliadas três formas de parcelas experimentais (retangular, uma linha - linear e parcela de uma árvore - STP) em testes clonais de Eucalyptus spp, utilizando-se três experimentos, cada um com 18 clones. Foram usados três modelos de análise (mínimos quadrados ordinários - ANOVA tradicional, modelos mistos com fator clone fixo ou com fator clone aleatório - REML/BLUP). Os dois primeiros modelos apresentaram resultados similares. Com REML/BLUP houve estreitamento das predições em relação às amplitudes obtidas com as médias, e essa redução foi proporcionalmente maior com parcelas retangulares e STP. O ordenamento dos clones também foi similar com esses dois tipos de parcelas. É provável que com parcelas STP haja um balanço compensatório das alocompetições, pois se pode trabalhar com maior número de repetições e menor custo. Portanto, com parcelas STP haverá economia de recursos e sem prejuízos para o Programa de Melhoramento Florestal.

Estimativa da relação hipsométrica em clones de Eucalyptus sp. com o modelo de curtis ajustado por métodos bayesianos empíricos

Neste trabalho foi considerado o modelo de Curtis para a relação hipsométrica em clones de Eucalyptus sp. com os parâmetros sujeitos a restrições. Para fazer a inferência dos parâmetros do modelo com restrições, utilizou-se uma abordagem bayesiana com densidade a priori construída empiricamente. As estimativas bayesianas são calculadas com a técnica de simulação de Monte Carlo em Cadeia de Markov (MCMC). O método proposto foi aplicado a diferentes conjuntos de dados reais, dos quais foram selecionados cinco para exemplificar os resultados. Estes foram comparados com os resultados obtidos pelo método de mínimos quadrados, destacando-se a superioridade da abordagem bayesiana proposta.

Variabilidade espacial e temporal da precipitação pluvial na bacia hidrográfica do Ribeirão João Leite-GO

A precipitação pluvial constitui-se na principal variável hidrológica de entrada do ciclo hidrológico e para conhecimento da variabilidade espacial e temporal como informação básica para estudos hidrológicos, manejo de bacias hidrográficas, gestão de recursos hídricos, dentre outros. Os recursos hídricos superficiais produzidos nessa bacia são essenciais ao desenvolvimento do eixo Goiânia-Anápolis-Brasília e, recentemente, com o aumento da demanda hídrica, tem sido relatados conflitos pelo uso da água. Dessa maneira, tem-se que o desenvolvimento de estudos hidrológicos visando ao melhor entendimento e aproveitamento dos recursos hídricos é estratégico para o desenvolvimento da região. Nesse sentido, objetivou-se com este estudo mapear a precipitação média mensal e anual na bacia hidrográfica do Ribeirão João Leite, com o auxílio de técnicas geoestatísticas. Foram avaliados os modelos de semivariograma esférico, exponencial e gaussiano, ajustados pelo método dos mínimos quadrados ponderados, sendo que o modelo utilizado no mapeamento por krigagem foi o que apresentou o menor erro médio indicado pela validação cruzada. Foi constatado bom desempenho das técnicas geoestatísticas no mapeamento da precipitação média mensal e anual, indicado pelos pequenos erros encontrados, podendo-se destacar o modelo de semivariograma exponencial, que se sobressaiu na maioria dos eventos estudados.

Emprego de ferramentas numéricas na avaliação do módulo de elasticidade em vigas roliças de madeira

Vigas são elementos estruturais encontrados na maioria das construções civis. Dentre os materiais de engenharia, destaca-se a madeira, por ter resistência mecânica satisfatória aliada a baixa densidade. A madeira roliça apresenta-se como boa solução na confecção de vigas, uma vez que não precisa ser processada, como é o caso da madeira serrada. O projeto de elementos estruturais de madeira requer o conhecimento de suas propriedades físicas e mecânicas, obtidas segundo as premissas de documentos normativos. Em se tratando da madeira roliça, os documentos normativos nacionais que tratam da determinação das propriedades de resistência e rigidez estão vigentes há mais de vinte anos sem revisão técnica. De forma geral, tanto as normas nacionais como as internacionais idealizam geometria troncocônica para as peças roliças de madeira, implicando equações simplificadas incapazes de prever a influência das irregularidades da forma na determinação do módulo de elasticidade longitudinal. Este trabalho objetiva avaliar a influência das irregularidades da geometria em peças roliças de madeira Corymbia citriodora e Pinus caribaea no cálculo do módulo de elasticidade longitudinal. Para tanto, utilizou-se do ensaio de flexão estática a três pontos, considerando também um modelo matemático simplificado, assumindo seção circular constante para a forma do elemento. As irregularidades das peças são consideradas nos modelos numéricos, constituídos de elementos finitos de barra e tridimensionais. Os resultados encontrados revelam equivalência estatística entre os módulos de elasticidade para ambas as formas de cálculo, indicando ser plausível a consideração de seção circular constante para as peças de madeira aqui avaliadas.

Variabilidade espacial da força de desprendimento de frutos do cafeeiro

O objetivo deste trabalho foi utilizar técnicas de georreferenciamento e de geoestatística para analisar a variabilidade espacial da força de desprendimento de frutos de cafeeiros por meio de semivariogramas e pela interpolação por krigagem. O trabalho foi conduzido no município de Três Pontas - MG, Brasil. A força de desprendimento dos frutos verdes e de cerejas dos cafeeiros foi obtida por meio de um protótipo de dinamômetro em pontos georreferenciados. A dependência espacial dos dados foi analisada por meio de ajustes de semivariogramas, clássico e robusto, para o método dos mínimos quadrados ordinários e ponderados, e apenas o ajuste clássico, para os métodos da máxima verossimilhança e máxima verossimilhança restrita. Testaram-se, para cada um dos métodos, os modelos esférico, exponencial e gaussiano. Os mapas de isolinhas obtidos por krigagem foram produzidos, baseados no melhor método e modelo de ajuste da função semivariograma, que foram obtidos pelas estatísticas de validação. As variáveis em estudo apresentaram estruturas de dependência espacial, as quais foram modeladas pelos semivariogramas, o que possibilitou a confecção dos mapas de isolinhas de distribuição espacial, obtidos por krigagem. Foi possível identificar os locais mais propícios para se iniciar a colheita seletiva e mecanizada dos frutos do cafeeiro.

Concentração Sérica Materna da Proteína C Reativa em Gestações Complicadas pela Pré-eclâmpsia

Objetivo: investigar a associação entre a concentração sérica da proteína C reativa e a ocorrência de pré-eclâmpsia, bem como sua relação com a gravidade da doença. Métodos: foram avaliadas, em estudo caso-controle transversal, 27 gestantes portadoras de pré-eclâmpsia e outras 27 sem nenhuma intercorrência clínica, no terceiro trimestre gestacional. As pacientes tiveram a dosagem sérica da proteína C reativa realizada no período antenatal, além de exames clínico e laboratoriais para diagnóstico da doença. Foram investigadas a associação entre a concentração sérica da proteína C reativa e a presença da pré-eclâmpsia e a correlação entre os valores desta proteína plasmática com os níveis da pressão arterial e excreção urinária de proteína. Empregou-se o teste de significância (chi²) e análise de regressão pela técnica dos mínimos quadrados, considerando-se significância estatística quando p<0,05. Resultados: as gestantes portadoras de pré-eclâmpsia apresentaram níveis da pressão arterial média maiores do que seus controles (129,9±12,1 e 87,2±6,5 mmHg, respectivamente) e valores médios da proteína C reativa significativamente superiores aos das normotensas (18,9±4,9 e 1,5±0,8 mg/L, respectivamente). Houve associação significativa entre a elevação da concentração de proteína C reativa e a ocorrência da pré-eclâmpsia (p<0,0001, "odds ratio": 20,1). Verificou-se também que a medida da pressão arterial média e a proteinúria apresentam correlação direta com a concentração da proteína C reativa circulante no sangue materno (p=0,001 e p=0,018, respectivamente). Conclusão: a proteína C reativa mostrou-se um marcador efetivo da ocorrência da pré-eclâmpsia e tem significativa correlação com a gravidade da doença. O uso deste exame para diagnóstico diferencial entre os diversos quadros hipertensivos da gestante e sua utilização como marcador de prognóstico da pré-eclâmpsia merecem novos estudos.