Enunciado

As falhas decorrentes de fadiga de alto ciclo em motores de avião de turbina a gás são um problema que tende a se propagar. O artigo "Effect of Crystal Orientation on Fatigue Failure of Single Crystal Nickel Base Turbine Blade Superalloys (J. of Engineering for Gas Turbines and Power, 2002, p. 161–176) apresentou os dados a seguir e ajustou um modelo de regressão não-linear a fim de prever a amplitude do esforço dos ciclos até a fratura. Ajuste um modelo apropriado, investigue a qualidade desse ajuste e preveja a amplitude quando o número de ciclos até a fratura = 5.000.

Passo 1

Um modelo apropriado, após verificar os gráficos de dispersão (o que você deve sempre fazer), é

y = α + β ln ⁡ x

A transformação para linearizar é

x ' = ln ⁡ x

e obteria uma forma linear, que pode ser estimada pelo método dos mínimos quadrados

y = α + β x '

Um gráfico de dispersão de dados transformados sugere um padrão linear; assim, um modelo de regressão linear simples é um bom ajuste.

Passo 2

A estimativa de mínimos quadrados da inclinação no modelo de regressão linear simples é

A estimativa de mínimos quadrados da interceptação é

As seguintes somas são necessárias para calcular as estimativas de mínimos quadrados e são fornecidas acima (estes valores são arredondados):

Usando os resultados anteriores

Encontrando y ¯   e   x ¯ usando as somas dadas acima, a estimativa do intercepto é

A linha de regressão linear é

y ^ = 0.0197 - 0.0013 ∙ ln ⁡ x

Para calcular as estimativas de α e β um modelo deve ser transformado; assim, as estimativas são

α ^ = 0.0197

β ^ = - 0.0013

Para investigar a qualidade do ajuste, devem ser calculados gráficos residuais, um teste t ou coeficiente de determinação. A maneira mais fácil de investigar a qualidade é testando as hipóteses apropriadas.

Ao testar H 0   : β 1 = β 10 contra uma das alternativas, o valor da estatística de teste é

t = β 1 ^ - β 10 s β 1 ^

Dependendo da hipótese alternativa, o valor P pode ser determinado como área correspondente sob t n - 2 . O teste de utilidade do modelo:

H 0   : β 1 = 0   v s   H a   : β 1 ≠ 0

O teste de uso é o teste de utilidade do modelo.

A variância e o desvio padrão do estimador   β ^ 1 são

Ao estimar o desvio padrão σ com s você obtém estimativa de σ β 1 ^

A soma dos quadrados dos erros, SSE, é

A estimativa do desvio padrão σ é

σ ^ = s = S S E n - 2

A soma dos quadrados dos erros, SSE, é

A fórmula computacional para SSE é

Pela fórmula computacional, o valor de SSE é

Passo 3

Portanto, a estimativa de σ é

s = S S E n - 2 = 0.0001 17 = 0.0027

A estimativa do desvio padrão do estimador   β ^ 1 é

Portanto, o valor da estatística de teste é

t = - 0.0013 - 0 0.0003 = - 4.0956

O valor PP é a área correspondente sob o t 17 curva

P = P T < t = 0.0008

que foi calculado por meio de um software. Finalmente, desde

P = 0.0008 < 0.01 = α

Rejeitar a hipótese nula em determinado nível de significância. O modelo é adequado.

Um intervalo de previsão de 95% para y quando x * =   5000 deve ser calculado.

Intervalo de previsão

é um intervalo de previsão de 100 ( 1 - α ) % para uma observação futura Y, quando x = x * .

O valor de   y ^ quando x * =   5000 , usando o modelo de regressão linear simples, é

A variância de uma variável aleatória   Y ^ é

O desvio padrão é a raiz quadrada da expressão. Ao substituir (estimar) σ por s você obtém a estimativa da variância e/ou desvio padrão da variável aleatória   Y ^ .

A estimativa do desvio padrão de   Y ^ é

Os valores necessários para calcular um intervalo de previsão de 95% podem ser encontrados acima; assim, um intervalo de previsão de 95% quando ln ⁡ x * = 8.517   é

Este é um intervalo de previsão de y quando x * = 5000 .

Há muitas outras maneiras de ver a validade do ajuste; no entanto, eles não serão apresentados aqui.

Resposta

Transformar para linearizar; Use teste t e gráficos; Calcular intervalo de previsão: (0.0030,0.0146).

Exercícios de Livros Relacionados

Refira-se ao Exercício 5. Calcule o IC de 95% t para
Considere os dados a seguir sobre o crescimento de plantas após a aplicação de diferentes tipos de hormônio de crescimento. a. Execute um teste F ao nível 0,05. b. O que ocorre quando se aplica o proc
Analise novamente os dados de rigidez axial fornecidos no Exercício 8. Seguem os resultados ANOVA obtidos no MINITAB:a. É razoável afirmar que as variâncias das cinco distribuições do índice de rigide
Considere os dados resumidos a seguir sobre o módulo de elasticidade (106 libras por polegada quadrada ou psi) para madeiras de três diferentes categorias (bem próximos dos valores usados no artigo "B
Repita o Exercício 12, supondo x2 502,8 e x3 427,5