Modelagem de atributos físico-hídricos do solo numa bacia hidrográfica da região do Alto Rio Grande, MG

O presente trabalho objetivou desenvolver modelos para estimativa dos atributos físico-hídricos do solo porosidade drenável (PD) e capacidade total de retenção de água (CTA), com base em atributos físicos de fácil e rápida obtenção (textura, argila dispersa em água, volume total de poros, densidade do solo e densidade de partícula) e a determinação de matéria orgânica. Os atributos foram coletados na bacia hidrográfica do Ribeirão Marcela, na camada de 0 a 15 cm, obedecendo aos grids de 240 x 240 e 60 x 60 m, totalizando 165 pontos amostrados. Trabalhou-se com regressão múltipla linear, constituindo-se variáveis por meio da combinação dos diferentes atributos entre si e estimando os respectivos coeficientes pelo método dos mínimos quadrados, utilizando-se o programa SAS for Windows. Os modelos gerados para CTA e PD apresentaram boa qualidade estatística, com elevados coeficientes de determinação, baixos erros de estimativa e variáveis significativas, o que permitiu sua aplicação na estimativa desses atributos na bacia hidrográfica do Ribeirão Marcela, representativa do domínio dos Latossolos na região do Alto Rio Grande, fornecendo assim subsídios para implantação de projetos agrícolas e ambientais, especialmente no auxílio à parametrização de modelos de simulação hidrossedimentológica, os quais estão em desenvolvimento para aplicação na região.

Seleção de modelos de variabilidade espacial para elaboração de mapas temáticos de atributos físicos do solo e produtividade da soja

Pesquisas sobre a variabilidade espacial dos atributos do solo que influenciam a produtividade são de uma grande importância para o desenvolvimento de novas técnicas que beneficiam a agricultura. A variabilidade desses atributos pode ser avaliada por técnicas de geoestatística e auxiliar no mapeamento e manejo do solo. Este trabalho teve por objetivo avaliar a qualidade do ajuste dos modelos teóricos espaciais segundo o Critério de Informação de Akaike, de Filliben, de Validação Cruzada e o valor máximo do logaritmo da função verossimilhança, de dados da umidade do solo, da densidade do solo e da resistência do solo à penetração, nas camadas de 0 a 0,1, 0,1 a 0,2 e 0,2 a 0,3 m, e de produtividade da soja do ano agrícola 2004-2005. Os parâmetros dos modelos de variabilidade espacial foram estimados por meio dos métodos de mínimos quadrados ordinários, mínimos quadrados ponderados e máxima verossimilhança. A pesquisa foi desenvolvida em uma área de 57 ha de um Latossolo Vermelho distroférrico, utilizando-se uma malha de 75 x 75 m georreferenciada. Concluiu-se que, dos métodos de avaliação de ajustes estudados, o da Validação Cruzada foi o mais adequado para escolha do melhor ajuste do modelo de variabilidade espacial; conseqüentemente têm-se mapas temáticos mais acurados.

A Lei de Lotka na bibliometria brasileira

Usando os dados reportados em artigos publicados em revistas brasileiras e trabalhos apresentados em congressos nacionais, replicaram-se as aplicações da Lei de Lotka à literatura brasileira em 10 campos diferentes. Utilizou-se o modelo do poder inverso pelos métodos do mínimo quadrado e probabilidade máxima. Das 10 literaturas nacionais analisadas, somente a literatura de medicina, siderurgia, jaca e biblioteconomia ajustaram-se ao modelo do poder inverso generalizado pelo método dos mínimos quadrados. No entanto, só duas literaturas (veterinária e cartas do Arquivo Privado de Getúlio Vargas) não se ajustaram ao modelo quando se usou o método da máxima probabilidade. Para ambas literaturas, tentaram-se modelos diferentes. A literatura de veterinária ajustou-se à distribuição binomial negativa, e as cartas do Arquivo Privado de Getúlio Vargas ajustaram-se melhor à distribuição Gauss-Poisson Inversa Generalizada.

Heterose sobre os pesos de bovinos Canchim e Aberdeen Angus e de seus cruzamentos recíprocos

O trabalho foi conduzido para estimar a heterose sobre os pesos ao nascimento (PNT), à desmama (P210) e ao ano (P365) e sobre os ganhos de pesos médios diários do nascimento à desmama (G210) e da desmama ao ano (G365) nas quatro primeiras gerações do sistema de cruzamentos alternados entre as raças Canchim (C) e Aberdeen Angus (A). Os dados de 1.147 bezerros nascidos de 1981 a 1998 foram analisados pelo método dos mínimos quadrados, ajustando-se um modelo linear que incluiu os efeitos linear e quadrático da idade da mãe do bezerro e os efeitos fixos de sexo, grupo genético, mês e ano de nascimento do bezerro. Estimativas de heterose e de outras diferenças genéticas foram estimadas por contrastes entre médias e testadas pelo teste t. O contraste "CA" foi positivo e significativo (P<0,001) para as cinco características. O contraste F1CAF1AC teve sinal negativo e foi altamente significativo (P<0,001) para P210 e G210 e significativo (P<0,05) para P365. A geração F1 exibiu heterose de 4,8% para P210 e de 4,9% para G210. A heterose materna foi de 3,7%, 5,8%, 6,3% e 20,4%, respectivamente, para P210, G210, P365 e G365. A heterose média das terceira e quarta gerações do cruzamento alternado entre C e A foi de 4,6% para P210, 5,3% para G210 e de 3,5% para P365.

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).

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.

Determinação simultânea de NADH e ácido ascórbico usando voltametria de onda quadrada com eletrodo de carbono vítreo e calibração multivariada

O método dos mínimos quadrados parciais (PLS) foi aplicado aos dados obtidos por voltametria de onda quadrada para a determinação simultânea de ácido ascórbico (AA) e do co-fator b-nicotinamida adenina dinucleotídeo (NADH) em misturas sintéticas. As curvas voltamétricas foram obtidas em tampão fosfato 0,2 mol L-1 (pH 8,0). Foi possível verificar que o método PLS permite a determinação destes compostos simultaneamente nas condições escolhidas. Os resultados mostraram um erro relativo máximo de 1,7% para o NADH e 2,1 % para o AA. A metodologia proposta é simples e desenvolvimentos posteriores podem torná-la útil para análises in vivo.

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.

Determinação de proteína em café cru por espectroscopia NIR e regressão PLS

A combinação da espectroscopia no infravermelho próximo (NIR) e calibração multivariada (método dos mínimos quadrados parciais - PLS) para a determinação do teor de proteína total em amostras de café cru, foi investigada. Os teores de proteína total foram inicialmente determinados usando-se como método de referência o de Kjeldhal, e, posteriormente foram construídos modelos de regressão a partir dos espectros na região do infravermelho próximo das amostras de café cru. Foram coletados 159 espectros das amostras de café cru utilizando um acessório de reflectância difusa, na faixa espectral de 4500 a 10000cm-1. Os espectros originais no NIR sofreram diferentes transformações e pré-tratamento matemático, como a transformação Kubelka-Munk; correção multiplicativa de sinal (MSC); alisamento (SPLINE); derivada primeira; média móvel e o pré-tratamento dos dados escalados pela variância. O método analítico proposto possibilitou a determinação direta, sem destruição da amostra, com obtenção de resultados rápidos e sem o consumo de reagentes químicos de forma a preservar o meio ambiente. O método proposto forneceu resultados com boa capacidade de previsão do teor de proteína total, sendo que os erros médios foram inferiores a 6,7%.

Determinação de umidade em café cru usando espectroscopia NIR e regressão multivariada

A espectroscopia na região do infravermelho próximo (NIR) foi usada para determinar o teor de umidade em amostras de café cru. Foram construídos modelos de regressão usando o método dos mínimos quadrados parciais (PLS) com diferentes pré-tratamentos de dados e 157 espectros NIR coletados de amostras de café usando um acessório de reflectância difusa, na região entre 4500 e 10000 cm-1. Os espectros originais passaram por diferentes transformações e pré-tratamentos matemáticos, como a transformação Kubelka-Munk; a correção multiplicativa de sinal (MSC); o alisamento com SPLINE e a média móvel, e os dados foram escalados pela variância. O modelo de regressão permitiu determinar o teor de umidade nas amostras de café cru com erro quadrático médio de calibração (SEC) de 0,569 g.100 g -1; erro quadrático médio de validação de 0,298 g.100 g -1; coeficiente de correlação (r) 0,712 e 0,818 para calibração e validação, respectivamente; e erro relativo médio de 4,1% para amostras de validação.