Regressão Linear

Ver Todos os Exercícios

Teoria

Introdução

Em química analítica, muitas vezes vamos estar lá no laboratório fazendo um experimento, e vamos precisar fazer uma análise em cima dos resultados que obtivemos. Uma das análises mais comuns nessa área é a regressão linear. Com ela, nós conseguiremos fazer com que os nossos dados seja aproximado por uma reta.

Imagine que você tenha que titular uma série de soluções padrões e anotar o volume de titulante usado. Os resultados você anotou na tabelinha abaixo:

Com esses valores, nós conseguimos colocar eles em um gráfico da seguinte forma:

O nosso objetivo, com a regressão linear, vai ser encontrar a reta que mais se aproxima de todos os pontos que encontramos. O gráfico abaixo mostra a melhor reta que representa os dados encontrados em nosso experimento:

Método dos mínimos quadrados

Para conseguir encontrar a equação dessa reta, podemos usar o método dos minímos quadrados. Esse método tem como principal objetivo encontrar a reta que minimiza o quadrado dos erros, isto é, a o quadrado da diferença entre os pontos experimentais e a reta aproximada. Sendo assim, precisaremos encontrar uma reta do tipo:

y = a x + b

Onde a é o coeficiente angular (inclinação da reta) e 1= b o coeficiente linear (interseção no eixo y ). Como eu disse, essa reta tem que ser aquela que minimiza os erros, ou resíduos. A soma dos quadrados dos resídiuos pode ser encontrada pela equação abaixo:

S S r e s = ∑ i = 1 N [ y i - a x i + b ] 2  

Precisamos então fazer com que S S r e s seja o menor possível, porque aí estaremos diminuindo os erros. Não vamos entrar no desenvolvimento matemático disso muito profundamente, mas vamos direto ao ponto.

A partir da equação acima, é possível chegar em três outras que vêm dela e que usaremos para calcular os coeficientes da nossa reta:

S x x = ∑ x i - x - 2

S y y = ∑ y i - y - 2

S x y = ∑ ( x 1 - x - ) ( y i - y - )

Sendo que x - e y - são os valores médios de x e y .

Com esses três valores, podemos calcular a e b da seguinte forma:

a = S x y S x x

b = y - - a x -

Vamos ver como ficaria para o nosso exemplo. Vamos fazer uma tabelinha para não nos perder:

Portanto:

S x x = 0,013

S y y = 169,040

S x y = 1,450

Vamos calcular os coeficientes da reta agora:

a = S x y S x x = 1,450 0,013 = 116

b = y - - a x - = 17,600 - 116 × 0,175 = - 2,7

Portanto, a reta que mais se aproxima dos dados obtidos no nosso experimento é:

y = 116 x - 2,7

Desvios padrões e intervalos de confiança

Com os valores de S x x , S x y e S y y também podemos calcular quatro tipos de desvio padrão da nossa reta:

  • Desvio padrão da regressão ( s r ) :
  • s r = S y y - a 2 S x x N - 2

  • Desvio padrão da inclinação ( s a ) e intervalo de confiança ( I C a ):
  • s a = s r 2 S x x

    Para encontrar o intervalo de confiança da inclinação basta aplicar a fórmula:

    I C a = a ± t ( n - 2 ) s a

  • Desvio padrão da interseção ( s b ) e intervalo de confiança ( I C b ):
  • s b = s r ∑ x i 2 N ∑ x i 2 - ∑ x i 2 = s r 1 N - ∑ x i 2 / ∑ x i 2

    Com intervalo de confiança:

    I C b = b ± t ( n - 2 ) s b

  • Desvio padrão dos resultados obtidos a partir da curva de calibração ( s x 0 ) e intervalo de confiança ( I C x 0 ):

s x 0 = s r a 1 M + 1 N + y - c - y - 2 a 2 S x x

Essa última é a expressão que indica o erro que obteremos caso utilizemos a equação de reta encontrada para encontrar um novo valor. Nessa equação y - c é o valor obtido na reta, M é o número de replicatas e N é o número de pontos na reta.

E um intervalo de confiança de um valor obtido através da equação da reta:

I C x 0 = x 0 ± t ( n - 2 ) s x 0

Coeficiente de correlação

Uma das principais formar de determinar se uma reta é uma boa representação do conjunto de pontos obtidos é através do coeficiente de correlação ( R 2 ). Quanto mais perto o valor de R 2 for de 1 melhor terá sido a regressão linear. O valor de R 2 pode ser encontrado da seguinte forma:

R 2 = 1 - S S r e s i d S S t o t

Onde S S r e s i d é a soma dos quadrados dos resíduos que definimos no início:

S S r e s = ∑ i = 1 N [ y i - a x i + b ] 2  

E S S t o t é a soma dos quadrados dos erros em y :

S S t o t = S y y

Vamos adicionar uma linha na nossa tabela para calcular S S r e s :

Sendo assim, temos:

S S r e s = 0,84

S S t o t = 169,040

Portanto, podemos calcular R 2 :

R 2 = 1 - 0,84 169,04 = 0,995

O valor encontrado é muito próximo de 1 e bastante satisfatório também. Geralmente, valores de R 2 maiores que 0,95 já são ótimos.

Desvios padrões e intervalos de confiança

Coeficiente de correlação

Exercícios Resolvidos

Exercício Resolvido #1

Os dados na tabela que segue foram obtidos durante uma determinação colorimétrica de glicose em soro sangüíneo.

Faça a regressão linear, encontre a equação da reta pelo método dos mínimos quadrado e determine o coeficiente de correlação.

Passo 1

Para começar, vamos primeiro definir que a concentração de glicose será nosso x e os valores de absorbância serão o nosso y :

Para fazer a análise dos mínimos quadrados, também precisamos das médias de x e y , então vamos calculá-las:

x - = 0,0 + 2,0 + 4,0 + 6,0 + 8,0 + 10,0 6 = 5,0

y - = 0,002 + 0,150 + 0,294 + 0,434 + 0,570 + 0,704 6 = 0,359

Passo 2

Além das médias, também é importante calcular os quadrados e os quadrados das diferenças dos valores, assim como suas somas. Sendo assim, vamos expandir nossa tabela de apoio:

Com o auxílio dessa tabela podemos determinar:

S x x = ∑ x i - x - 2 = 70,00

S y y = ∑ y i - y - 2 = 0,34

S x y = ∑ x 1 - x - y i - y - = 4,91

E com isso, já conseguiresmo determinar os coeficientes da reta:

a = S x y S x x = 4,91 70,0 = 0,07

b = y - - a x - = 0,352 - 0,007 × 5 = 0,008

Como nossa reta é do tipo:

y = a x + b

Ela será então:

y = 0,07 x + 0,008

Passo 3

Para calcular R 2 precisaremos de mais uma coluna na nossa tabela, a coluna de S S r e s i d .

Como:

S S r e s i d = ∑ y i - a x i + b 2

Logo:

S S r e s i d = 0,0001

O coeficiente de correlação R 2 pode ser calculado através da seguinte fórmula:

R 2 = 1 - S S r e s i d S S t o t

Onde:

S S t o t = S y y = 0,34

Portanto:

R 2 = 1 - 0,0001 0,34 = 0,9997

Resposta

y = 0,07 x + 0,008

R 2 = 0,9997

Exercício Resolvido #2

Os dados na tabela a seguir representam o potencial de eletro E verus a concentração c .

  1. Transforme os dados em valores de E versus - log ⁡ c
  2. Acha a equação da reta da regressãp linear de E versus - log ⁡ c
  3. Encontre o coeficiente de correlação.

Passo 1

Primeiramente, devemos transformar nossos dados na coluna da concentração. Para obter uma relação linear entre os dois dados, precisaremos aplicar - log ⁡ c em todos os valores de concentração. Sendo assim, nossa nova tabela, será:

Passo 2

Agora, já poderemos começar a aplicar o método dos mínimos quadrados. Vamos chamar o potencial do eletro de y e - log ⁡ c de x .

Podemos também calcular os valores de x - e y - :

x - = ∑ i = 0 N x i N = 1,94974

y - = ∑ i = 0 N y i N = 166

Passo 3

Além das médias, também é importante calcular os quadrados e os quadrados das diferenças dos valores, assim como suas somas. Sendo assim, vamos expandir nossa tabela de apoio:

Com o auxílio dessa tabela podemos determinar:

S x x = ∑ x i - x - 2 = 5,97

S y y = ∑ y i - y - 2 = 18490,00

S x y = ∑ x 1 - x - y i - y - = 330,60

E com isso, já conseguiresmo determinar os coeficientes da reta:

a = S x y S x x = 330,60 5,97 = 55,37

b = y - - a x - = 166 - 55,37 × 1,94974 = 58,04

Como nossa reta é do tipo:

y = a x + b

Ela será então:

y = 55,37 x + 58,04

Passo 4

Para calcular R 2 precisaremos de mais uma coluna na nossa tabela, a coluna de S S r e s i d :

Como:

S S r e s i d = ∑ y i - a x i + b 2

Logo:

S S r e s i d = 183,9996

O coeficiente de correlação R 2 pode ser calculado através da seguinte fórmula:

R 2 = 1 - S S r e s i d S S t o t

Onde:

S S t o t = S y y = 18490,00

Portanto:

R 2 = 1 - 183,9996 18490,00 = 0,990

Resposta

y = 55,37 x + 58,04

R 2 = 0,990

Exercício Resolvido #3

A concentração do íon sulfato em águas naturais pode ser determinada pela medida da turbidez que resulta quando um excesso de BaCl2 é adicionado a uma quantidade medida da amostra. Um turbidímetro, instrumento usado para essa análise, foi calibrado com uma série de padrões de soluções padrão de Na2SO4. Os seguintes dados foram obtidos na calibração:

  1. Encontre a equação da reta através do método dos mínimos quadrados.
  2. Encontre o coeficiente de correlação, R 2
  3. Para uma leitura no turbidímetro de 3,50 , qual seria a concentração de sulfato?

Passo 1

Precisamos então aplicar o método dos mínimos quadros nos dados que o enunciado nos deu para conseguir calcular nossa regressão linear. Antes de começar a fazer qualquer coisa, vamos definir a concentração de sulfato como nossa variável independente ( x ) e a leitura do turbidímetro como a variável dependente ( y ):

Com isso, vamos calcular nossos valores médios, x - e y - :

x - = 0,00 + 5,00 + 10,00 + 15,0 + 20,0 5 = 10,00

y - = 0,06 + 1,48 + 2,28 + 3,98 + 4,61 5 = 2,482

Passo 2

Além das médias, também é importante calcular os quadrados e os quadrados das diferenças dos valores, assim como suas somas. Sendo assim, vamos expandir nossa tabela de apoio:

Com o auxílio dessa tabela podemos determinar:

S x x = ∑ x i - x - 2 = 250,0000

S y y = ∑ y i - y - 2 = 13,6833

S x y = ∑ x 1 - x - y i - y - = 58,0000

E com isso, já conseguiresmo determinar os coeficientes da reta:

a = S x y S x x = 58,00 250,0 = 0,232

b = y - - a x - = 2,482 - 0,232 × 10,00 = 0,162

Como nossa reta é do tipo:

y = a x + b

Ela será então:

y = 0,232 x + 0,162

Passo 3

Para calcular R 2 precisaremos de mais uma coluna na nossa tabela, a coluna de S S r e s i d .

Como:

S S r e s i d = ∑ y i - a x i + b 2

Logo:

S S r e s i d = 0,2273

O coeficiente de correlação R 2 pode ser calculado através da seguinte fórmula:

R 2 = 1 - S S r e s i d S S t o t

Onde:

S S t o t = S y y = 13,6833

Portanto:

R 2 = 1 - 0,2273 13,6833 = 0,983

Passo 4

Por fim, precisamos estimar qual seria a concentração da amostra se a leitura no turbidímetro fosse de 3,50 . Para isso, vamos usar a equação da reta que encontramos:

y = 0,232 x + 0,162

Como a leitura é 3,50 , temos que:

y = 3,50

Jpgando na equação da reta:

3,50 = 0,232 x + 0,162

x = 3,50 - 0,162 0,232 = 14,39   m g   L - 1

Resposta

  1. y = 0,232 x + 0,162
  2. R 2 = 0,983
  3. x = 14,39   m g   L - 1

Exercício Resolvido #4

Na determinação de um fármaco por fluorescência em meio ácido, foram obtidos os seguintes resultados:

  1. Calcule o valor de R 2 e a equação da reta
  2. Considerando que a intensidade de uma amostra foi de 16,1 , calcule o valor de quinina na amostra
  3. Calcule os desvio padrões e os intervalos de confiança a 95 % da inclinação e da interseção da reta.

Passo 1

Para começar, vamos primeiro definir que a concentração de quinina será nosso x e os valores de fluorescência serão o nosso y :

Para fazer a análise dos Mínimos quadrados, também precisamos das médias de x e y , então vamos calculá-las:

x - = 0,00 + 0,10 + 0,20 + 0,30 + 0,40 5 = 0,20

y - = 0,00 + 5,20 + 9,90 + 15,30 + 19,10 5 = 9,90

Passo 2

Além das médias, também é importante calcular os quadrados e os quadrados das diferenças dos valores, assim como suas somas. Sendo assim, vamos expandir nossa tabela de apoio:

Com o auxílio dessa tabela podemos determinar:

S x x = ∑ x i - x - 2 = 0,10

S y y = ∑ y i - y - 2 = 233,90

S x y = ∑ x 1 - x - y i - y - = 4,83

E com isso, já conseguiresmo determinar os coeficientes da reta:

a = S x y S x x = 4,83 0,10 = 48,3

b = y - - a x - = 9,90 - 48,3 × 0,20 = 0,24

Como nossa reta é do tipo:

y = a x + b

Ela será então:

y = 48,3 x + 0,24

Passo 3

Para calcular R 2 precisaremos de mais uma coluna na nossa tabela, a coluna de S S r e s i d :

Como:

S S r e s i d = ∑ y i - a x i + b 2

Logo:

S S r e s i d = 0,61

O coeficiente de correlação R 2 pode ser calculado através da seguinte fórmula:

R 2 = 1 - S S r e s i d S S t o t

Onde:

S S t o t = S y y = 233,90

Portanto:

R 2 = 1 - 0,61 233,90 = 0,997

Passo 4

Caso a intensidade de fluorescência em uma amostra fosse de 16,1 , ou seja, nosso valor de y , para achar a concentração da amostra, precisamos apenas aplicar esse valor na equação da reta encontrada, que foi:

y = 48,3 x + 0,24

16,1 = 48,3 x + 0,24

x = 16,1 - 0,24 16,1 = 0,328   m g   L - 1

Passo 5

Agora devemos calcular o desvio padrão e intervalo de confiança para a inclinação e a interseção. Para calular isso, primeiro precisamos calcular o desvio padrão da regressão ( s r ) conforme a fórmula abaixo:

s r = S y y - a 2 S x x N - 2

Com N = 5 , pois é o número de pontos. Sendo assim, o desvio padrão da regressão será:

s r = 233,90 - 48,3 2 × 0,10 3 = 0,451

Para a inclinação, o desvio padrão será:

s a = s r 2 S x x = 0 , 451 2 0,1 = 1,427

E o intervalo de confiança:

I C a = a ± t ( n - 2 ) s a

Onde t n - 2 é o coeficiente t com grau de liberdade 2 graus menor que o número de pontos e nível de confiança de 95 % . Consultando uma tabela, teremos:

t n - 2 = 3,182

Portanto:

I C a = 48,3 ± 3,182 × 1,427 = 48,3 ± 4,54

Para a interseção, a fórmula do desvio padrão é:

s b = s r 1 N - ∑ x i 2 / ∑ x i 2

Usando os valores da tabela:

s b = 0,451 1 5 - 1 2 / 0,30 = 0,349

E o intervalo de confiança é calculado de maneira análoga:

I C b = b ± t ( n - 2 ) s b

I C b = 0,24 ± 3,182 × 0,349 = 0,24 ± 1,11

Resposta

  1. y = 48,3 x + 0,24
  2. R 2 = 0,99 7

  3. x = 0,38   m g   L - 1
  4. s a = 1,427

I C a = 48,3 ± 4,54

s b = 0,349

I C b = 0,24 ± 1,11

Exercício Resolvido #5

Desenvolva uma análise de mínimos quadrados dos dados fornecidos pelas duas primeiras colunas da tabela abaixo que representa os dados obtidos na calibração para a determinação cromotográfica do isoctano em uma mistura de hidrocarbonetos

Passo 1

Bom pessoal, quando a questão pede para a gente fazer a análise do método dos mínimos quadrado, isso significa que a gente precisa fazer TUDO. Ou seja, calular a reta, o coeficiente de correlação, as variâncias e os intervalos de confiança.

Vamos começar então pelos médias:

x - = 0,352 + 0,803 + 1,08 + 1,38 + 1,75 5 = 1,07

y - = 1,09 + 1,78 + 2,60 + 3,03 + 4,01 5 = 2,50

Passo 2

Com os valores das médias, conseguiremos calcular tudo que precisamos na nossa tabela para então começar a fazer os nosso cálculos:

Com o auxílio dessa tabela podemos determinar:

S x x = ∑ x i - x - 2 = 1,1454

S y y = ∑ y i - y - 2 = 5,0775

S x y = ∑ x 1 - x - y i - y - = 2,3697

E com isso, já conseguiresmo determinar os coeficientes da reta:

a = S x y S x x = 2,3697 1,1454 = 2,09

b = y - - a x - = 2,50 - 2,09254 × 1,07 = 0,26

Como nossa reta é do tipo:

y = a x + b

Ela será então:

y = 2,09 x + 0,26

Passo 3

Vamos partir agora para o cálculo nosso coeficiente de correlação, R 2 . Para isso, precisaremos de mais uma coluna na nossa tabela, a coluna de S S r e s i d .

Como:

S S r e s i d = ∑ y i - a x i + b 2

Logo:

S S r e s i d = 0 , 0625

O coeficiente de correlação R 2 pode ser calculado através da seguinte fórmula:

R 2 = 1 - S S r e s i d S S t o t

Onde:

S S t o t = S y y = 5,0775

Portanto:

R 2 = 1 - 0 , 0625 5,0775 = 0,98 8

Passo 4

Agora precisaremos calcular as variâncias e os intervalos de confiança do coeficiente angular e do coeficiente linear.

Porém, antes disso, precisaremos da variância da regressão, ela é dada por:

s r = S y y - a 2 S x x N - 2

Substituindo os valores, teremos:

s r = 5,0775 - 2,09 2 × 1,1454 3 = 0,144

Passo 5

Vamos calular agora a variância do coeficiente angular ( a ) :

s a = s r 2 S x x

Substituindo tudo:

s a = 0,144 2 1,1454 = 0,135

Com a variância de a , podemos calcular seu intervalo de confiança, pela fórmula abaixo:

I C a = a ± t ( n - 2 ) s a

Onde t é tem grau de liberdade duas unidades menor que oo número de pontos, ou seja, três. Como o enunciado não determinou o nível de confiança, vamos fazer para 95 % que é o nível de confiança mais usual. O valor de t para grau de liberdade 3 e nível de confiança 95 % é de 3,182 . Portanto, nosso intervalo de confiança será:

I C a = 2,09 ± 3,182 × 0,135 = 2,09 ± 0,43

Passo 6

Vamos fazer o mesmo agora para o coeficiente angular b :

s b = s r ∑ x i 2 N ∑ x i 2 - ∑ x i 2 = s r 1 N - ∑ x i 2 / ∑ x i 2

Substituindo na fórmula:

s b = 0,144 1 5 - 5,37 2 / 6,9020 = 0,158

O intervalo de confiança será portanto:

I C b = b ± t ( n - 2 ) s b

I C b = 0,26 ± 3,182 × 0,158 = 0,26 ± 0,50

Resposta

y = 2,09 x + 0,26

R 2 = 0,988

s r = 0,144

s a = 0,135

s b = 0,158

I C a = 2,09 ± 0,43

I C b = 0,26 ± 0,50

Exercícios de Livros Relacionados

Problema Desafiador. Zwanziger e Sârbu conduziram um estudo para validar os métodos analíticos e instrumentos. Os seguintes dados são os resultados obtidos na determinação de mercúrio em resíduos sóli
Os dados na tabela que segue foram obtidos durante uma determinação colorimétrica de glicose em soro sanguíneo.Quais os desvios padrão para a inclinação e o intercepto? Qual o erro padrão para a estim
Os dados que seguem representam áreas relativas de picos obtidas para cromatogramas de soluções padrão de metilvinilcetona (MVC).Uma amostra contendo MVC gerou uma área relativa de pico de 10,3. Calcu
Os dados na tabela a seguir representam o potencial de eletro E versus a concentração c .Transforme os dados em valores de E versus – l o gc .
Os dados que seguem representam áreas relativas de picos obtidas para cromatogramas de soluções padrão de metilvinilcetona (MVC).Faça uma tabela ANOVA.