A inferência estatística é um conjunto de técnicas que objetiva estudar a população através de evidências fornecidas por uma amostra. O verbo inferir tem o significado de deduzir através do raciocínio, na estatística o processo de dedução ocorre através da análise de dados presentes em uma amostra. Ao sortear uma amostra de uma população, a sua composição é por si só um fenômeno aleatório, ou seja, diferentes sorteios geram diferentes amostras. Por outro lado, diferentes amostras levam a diferentes resultados para as estatísticas de interesse tais como: média, proporção,variância,…Ao pensarmos no tradicional exemplo da amostragem de eleitores para estimação da proporção de votos de um candidato, diferentes institutos de pesquisa coletam diferentes amostras e, por conseqüência, geram diferentes estimativas. Então, cientes de que diferentes amostragens produzem diferentes valores, torna-se relevante questionar como é variabilidade das repostas obtidas em diferentes amostras.

O interesse em coletar uma amostra aleatória é fazer observações de um sequência n de váriaveis aleatórias que aqui será denotada por:

Para entender como as respostas em uma amostra aleatória podem variar, observe o exemplo em Magalhães e Lima (2004):

Exemplo 7.1: Uma empresa fabrica 100 equipamentos eletrônicos por semana e deseja verificar como se comporta a resistência deste tipo de equipamento em relação à alteração de voltagem. Um teste planejado pelo controle de qualidade da empresa, consiste em produzir sucessivas alterações padronizadas de voltagem e observar o efeito no aparelho. Como esses testes são demorados e demandam custos expressivos, apenas 5 desses aparelhos serão testados.

Neste exemplo, ao observar um aparelho, a resistência à alteração na voltagem pode ser considerada boa ou má. Isto corresponde a observar uma variável de Bernoulli, ou seja,

O procedimento do controle de qualidade consiste em observar 5 variáveis aleatórias de Bernoulli, o que pode ser representado pelo vetor:

Em cada semana a amostra de tamanho n =5 é coletada, ou seja, observa-se o vetor aleatório acima. Suponha que durante um mês este tenha sido o resultado das observações:

Caso nestas semanas o processo de fabricação esteja sob controle, e as peças tenham sido sorteadas de modo a representar bem os 100 equipamentos fabricados na semana, os valores acima representam 4 diferentes configurações para uma amostra aleatória.

Suponha que a estatística de interesse em cada amostra seja a quantidade de equipamentos com boa resistência, ou seja:

Esta estatística de interesse corresponde a uma variável aleatória com distribuição binomial. Assim, podemos calcular probabilidades para os seus possíveis valores.

(Inserir aqui Cap 10 Bussab)

Definição Parâmetro é uma quantidade de interesse que, em geral, é desconhecida na população.

Estas quantidades serão representadas, geralmente, por letras gregas e podem especificar uma distribuição de probabilidades. As técnicas de inferência estatística encarregam-se de deduzir o valor de um parâmetro por meio da avaliação de uma amostra.

Alguns parâmetros de interesse

Quando falarmos de forma geral sobre o parâmetro de interesse, vamos aqui utilizar a letra grega .

Definição Estimador é uma função dos valores de uma amostra construída para estimar um parâmetro de interesse. A notação apresentada a seguir traz uma forma padrão de representar um estimador, colocando um acento circunflexo no parâmetro a ser estimado.

Repare que o estimador é uma função de variáveis aleatórias e, portanto, também é uma variável aleatória.

Um exemplo clássico de estimador é a média amostral. Por meio desta função que soma n elementos de uma amostra e posteriormente divide o resultado por n, temos um estimador intuitivamente natural para a média populacional. Ao avaliar esta função em um conjunto de dados observados, o resultado é denominado estimativa. Com a aplicação de um estimador, podemos encontrar diferentes estimativas em diferentes amostras.

Veja o exemplo abaixo:

Exemplo 7.2 Para estudar o nível de colesterol em uma população de esportistas, coletamos uma amostra de 10 jovens atletas, obtendo os seguintes valores: 180 ; 196 ; 185 ; 165 ; 190 ; 180 ; 176 ; 165 ; 195.

O interesse nesta população é conhecer:

Consideremos 3 estimadores para esta quantidade:

notando que X_(i) é o i-ésimo valor da amostra ordenada, o primeiro estimador emprega a média aritmética entre o mínimo e máximo da amostra como função para inferência.

O segundo estimador simplesmente utiliza o primeiro valor presente na amostra (não ordenada) como estimador para a média populacional.

O último dos estimadores propostos utiliza a média aritmética dos valores da amostra como estimador para a média populacional.

A aplicação dos três estimadores acima à amostra, vai gerar as estimativas :

Exemplo 7.3 Em uma cidade, os taxis estão numerados de 1 até , sendo é um parâmetro desconhecido que representa a quantidade de taxis na cidade. Supondo que os taxis circulam de modo uniforme por toda cidade, uma pessoal anotou a placa dos 5 primeiros taxis que passaram em uma determinada esquina. Estes números foram:

Esta amostra representa 5 observações da variável aleatória X, número do táxi que passa por aquela esquina. Para estimar o parâmetro, são propostos 3 estimadores:

Os três estimadores acima representam três propostas para estimar a quantidade total de taxis na cidade. As funções da amostra apresentadas acima são respectivamente: máximo, mediana e máximo+mínimo. Ao aplicarmos estes estimadores na amostra obtida teremos as seguintes estimativas:

Cada um dos exemplos acima propõe 3 estimadores, estes são utilizados em uma amostra observada da variável de interesse e são encontradas diferentes estimativas. A questão relevante neste momento é "Qual estimador é o mais apropriado ? ". A princípio esta questão parece não ter resposta, pois não conhecemos o valor do parâmetro de interesse.

Porém, o estimador é uma variável aleatória, logo podemos pensar em calcular probabilidades para seus possíveis valores e avaliar estatísticas como: valor esperado e variância. A partir deste fato são desenvolvidos princípios para qualificar e diferenciar os estimadores. Um estimador mais "preciso", por exemplo, é aquele que possui menor variabilidade de amostra para amostra. O valor esperado de um estimador deve ser o valor do parâmetro de interesse na população. Na sequência são apresentadas algumas propriedades desejáveis para um bom estimador.

A escolha entre vários estimadores estará norteada pelas propriedades apresentadas a seguir:

Um estimador é dito ser não viciado (não tendencioso) se o seu valor esperado é igual ao parâmetro de interesse. Matematicamente escrevemos esta propriedade como:

Caso o estimador seja viciado, a quantificação deste vício será representada por:

Exemplo 7.4: Uma variável aleatória X possui média e variância . Para uma amostra aleatória de tamanho n:

mostre que os estimadores :

e

são não viciados para o parâmetro

Estas demonstrações podem ser feitas utilizando somente as propriedades do operador esperança E(.). Para o primeiro estimador, é trivial demonstrar a não-tendenciosidade.

Para o segundo estimador, média amostral, basta lembrar que o valor esperado de uma soma de variáveis independentes é a soma dos valores esperados, logo,

O resultado acima é bem geral pois afirma que a média amostral é sempre um estimador não viciado da média populacional. Entretanto, se considerarmos para o mesmo exemplo o estimador :

para o parâmetro , constatamos que este é viciado !!!

O vício presente neste estimador é eliminado quando definimos um novo estimador:

Um estimador é dito ser consistente se o seu valor converge para o valor do parâmetro de interesse a medida que a amostra aumenta. Além disto, a sua variabilidade, medida pela variância, converge para 0. Matematicamente, as propriedades de um estimador consistente são expressas como:

Repare que a propriedade 1 é menos restritiva do que a propriedade de não-tendenciosidade. O que interessa é que o estimador seja não viciado para grandes amostras, ou seja, valores grandes de n. Eventualmente o estimador pode até ser viciado em pequenas amostras.

O vício do estimador

é dado por :

.

Entretanto, é direta a verificação de que:

e logo este estimador é consistente para a variância populacional.

Sejam dois estimadores , não viciados para . O estimador é dito ser mais eficiente do que caso tenha menor variabilidade, ou seja :

.

Esta propriedade estabelece termos de comparação entre dois estimadores. Casos ambos sejam não viciados, a escolha deve recair sobre aquele que for mais eficiente. Como em geral a variabilidade de um estimador está em função de n, tamanho de amostra, um estimador mais eficiente requer uma amostra menor para gerar a mesma precisão do estimador "menos" eficiente.

Exemplo 7.5 : Vamos ilustrar a propriedade de eficiência considerando uma amostra de tamanho n da variável aleatória do Exemplo 7.4 e dois estimadores não viciados para a média populacional:

e

Utilizando as propriedades do operador variância, temos que :

e

de acordo com a expressão das duas variâncias, o estimador é mais eficiente nas situações em que n>1.

Na inferência estatística, procura-se avaliar como as estimativas produzidas por um estimador variam de uma amostra para outra. A situação ideal é aquela em que há conhecimento da distribuição de probabilidade do estimador. Por intermédio dela, é possível avaliar o valor esperado do estimador, assim como sua variância. Nesta seção, a distribuição de probabilidade do estimador é chamada de distribuição amostral.

Veja o exemplo em Magalhães e Lima (2004) que será discutido posteriormente

Exemplo 7.12: Um jogo consiste em lançar uma moeda honesta 3 vezes. Para cada lançamento, se cair cara você ganha 1 ponto, caso saia coroa você perde 1 ponto.

A variável aleatória em questão pode ser expressa da seguinte maneira :

X: resultado, em pontos, após o lançamento de uma moeda.

Sendo assim,

Uma amostra de tamanho n=3 corresponde a observação do seguinte vetor de variáveis aleatórias:

onde cada uma delas tem a distribuição de probabilidades :

Nesta situação, a variável pode assumir somente dois valores e podemos representar todas as possíveis amostras de tamanho 3 , conforme mostrado abaixo:

O valor da média amostral varia conforme a amostra e, neste caso, temos conhecimento das probabilidades das possíveis médias amostrais. Com estas probabilidades, é direto verificar que:

Dessa forma, uma vez que

então o estimador é não viciado para o parâmetro estimado.

Veja agora um resultado geral para a distribuição de probabilidade da média amostral: o Teorema Central do Limite

Exemplo 7.15: Uma variável aleatória X assume os valores 3, 6 e 8 com probabilidades 0,4; 0,3 e 0,3. Uma amostra com 40 observações é sorteada. A variável X não tem distribuição Normal e obtemos µ=5,4 e σ²=4,44. Apesar de não ser simétrica, consideramos que 40 observações é uma amostra grande o suficiente para usar o Teorema Central do Limite.

Para calcular a probabilidade da média amostral superar o valor 5, temos:

Uma aplicação importante do Teorema Central do Limite relaciona-se com a distribuição da proporção amostral:

Se construirmos para o i-ésimo indivíduo uma variável aleatória Y_i tal que Y_i=1 se o indivíduo apresenta a característica, e Y_i=0 caso contrário, podemos reescrever a proporção amostral como:

Logo, a proporção amostral nada mais é do que a média de variáveis aleatórias convenientemente definidas.

Assim, da mesma forma que um conjunto de médias amostrais são distribuídas nas proximidades da média populacional, as proporções amostrais são distribuídas ao redor da verdadeira proporção populacional p. Devido ao Teorema Central do Limite, para n grande e p não muito próximo de 0 ou 1, a distribuição de será aproximadamente normalmente distribuída com média p e variância dada por , ou seja, .

Os estimadores discutidos até aqui são estimadores pontuais, pois fornecem como estimativa um único valor numérico para o parâmetro de interesse. Por serem variáveis aleatórias, os estimadores possuem uma distribuição de probabilidade e, levando este fato em consideração, podemos apresentar uma estimativa mais informativa para o parâmetro de interesse que inclua uma medida de precisão do valor obtido.

Esse método de estimação, denominado intervalo de confiança, incorpora à estimativa pontual do parâmetro, informações a respeito de sua variabilidade.

Consideremos, inicialmente, o intervalo de confiança para a média µ de uma certa população Normal, com variância conhecida σ². Supondo uma amostra de tamanho n dada por (X₁, X₂,…, X_n), vimos que a média amostral tem distribuição Normal com a mesma média µ e variância σ²/n.

Assim,

Fixado um valor γ tal que 0<γ<1, podemos encontrar um valor z_γ/2 tal que

O índice de z_γ/2 apresenta o valor de γ dividido por 2 uma vez que a "massa" γ deve ser distribuída igualmente em torno de 0.

O valor z_γ/2 pode ser obtido da tabela da Normal padrão, localizando o valor de γ/2 no corpo da tabela e obtendo o valor z_γ/2 nas margens correspondentes. Feito isso, temos o intervalo

que pode reescrito como

Assim, o intervalo de confiança para μ, com coeficiente de confiança γ, é dado por

A interpretação do intervalo de confiança deve ser feita com cuidado. A expressão IC(μ;γ) envolve a quantidade que é uma variável aleatória e, portanto, o intervalo obtido também é aleatório. A probabilidade que ele contenha o verdadeiro valor da média populacional μ é dada por γ.

No entanto, note que ao coletar a amostra, torna-se e, como conhecemos σ, n e z_γ/2, o intervalo passa a ser numérico.

Desta forma, uma interpretação conveniente é a seguinte:

se obtivermos várias amostras de mesmo tamanho e, para cada uma delas, calculamos os correspondentes intervalos de confiança com coeficientes de confiança γ, esperamos que a proporção de intervalos que contenham o valor de μ seja igual a γ.

O exemplo de Magalhães e Lima (2004) a seguir ilustra os conceitos discutidos.

Exemplo 7.18: Suponha que os comprimentos de jacarés adultos de uma certa raça siga o modelo Normal com média μ desconhecida e variância σ²=0,01 m². Uma amostra de dez animais foi sorteada e forneceu média 1,69 m.

Desejamos uma estimativa para o parâmetro desconhecido μ.

Uma vez que podemos proceder de forma análoga ao desenvolvimento apresentado acima e obter uma estimativa por intervalo para μ.

Estabelecendo γ=95% obtemos da tabela da Normal z_γ/2=z_0,475=1,96.

Segue então que

A aplicação do Teorema Central do Limite (TCL) permite a obtenção de intervalos de confiança para μ e p. Neste caso, o intervalo construído terá um coeficiente de confiança aproximadamente igual a γ, sendo que esta aproximação melhora à medida que aumenta o tamanho da amostra.

Segundo o TCL, por exemplo, um intervalo de confiança de aproximadamente 95% para p é portanto

em que

Note que não sabemos o verdadeiro valor de p, e portanto usamos na fórmula acima para estimar SE. Chamamos SE de erro padrão da proporção amostral.

Exemplo: Um ensaio clínico foi realizado para determinar a preferência entre dois analgésicos, A e B, contra dor de cabeça. Cem pacientes que sofrem de dor de cabeça crônica receberam em dois tempos diferentes o analgésico A e o analgésico B. A ordem na qual os pacientes receberam os analgésicos foi determinada ao acaso. Os pacientes desconheciam esta ordem. Ao final do estudo foi perguntado a cada paciente qual analgésico lhe proporcionou maior alívio: o primeiro ou o segundo. Dos 100 pacientes, 45 preferiram A e 55 preferiram B. Com base nestas informações podemos dizer que há prefência por algum dos analgésicos?

Dizemos que não há preferência por um dos analgésicos quando a proporção dos que preferem A (p_A), é igual a proporção dos que preferem B (p_B). Como temos dois resultados possíveis, p_A e p_B são iguais quando p_A=p_B=0,5.

Um intervalo de aproximadamente 95% de confiança para a verdadeira proporção de pacientes que preferem o analgésico A é:

Então com aproximadamente 95% de confiança, a verdadeira proporção de pacientes que preferem o analgésico A está entre 0,35 e 0,55. Observe que este intervalo contem o valor 0,5 então concluímos que não existem evidências amostrais de preferência por um dos analgésicos.

A amplitude do intervalo de confiança para μ é dada por ou seja, depende da confiança γ, do desvio padrão σ e do tamanho da amostra n.

É usual referir-se à semi-amplitude, como o erro envolvido na estimação. Podemos utilizar esta medida para fazer o dimensionamento de amostras.