anim.msd                     package: nenhum                   R Documentation

Animação dos Desvios Mínimos de uma Regressão Linear

Description:

A função apresenta e salva no diretório de trabalho uma animação em GIF demostrando a 
mudança no tamanho dos resíduos para diferentes retas de uma regressão linear, à medida 
que o coeficiente de inclinação muda. No final da animação é apresentada a reta com o 
melhor ajuste. Além disso, a função retorna um summary da regressão linear da variável 
y em função da variável x.

Usage:

anim.msd= (x, y, ...)

## Default
anim.msd(x, y, lim=1, int=0.05, cor= "tot", name="animation.gif", time=1, xlabel="X", 
ylabel="Y")

Arguments:

x	vetor numérico da variável preditora quantitativa.

y	vetor numérico da variável resposta quantitativa.

lim	valor positivo que será somado ou subtraído ao coeficiente de inclinação do 
	melhor modelo 	da regressão linear e que estabelecerá um intervalo com os 
	limites superior e inferior das inclinações mostradas pela animação. Para 
	melhor visualização 0 ≤ lim ≤ 20, ao usar valores próximo a 20 recomenda-
	se utilizar um int ≥ 0.5.  Default de lim=1.
  
int	valor positivo que determina a diferença entre as inclinações de cada reta 
	plotada dentro do intervalo estabelecido por lim. Para melhor visualização 
	0.01 ≤ int ≤ lim. Default de int=0.05.

time	tempo em segundos entre uma imagem (i.e. plot ou frame) e a próxima, 
	controlando a velocidade da animação. Para melhor visualização 
	0.01 ≤ time ≤ 5. O default é time=1.

cor	indica se as cores vermelha e verde serão computadas para o conjunto dos 
	tamanhos dos desvios (“tot”), ou se serão computadas para cada desvio 
	separadamente (“ind”). Default de cor="tot".

name	nome que será dado à animação em GIF. O default é name="animation.gif". Ao 
	mudar o nome sempre manter a extensão '.gif'.

xlabel	nome do eixo x, deve ser colocado entre aspas.

ylabel	nome do eixo y, deve ser colocado entre aspas.


Details:

O método de minimizar a soma dos desvios quadráticos, ou resíduos, entre os valores 
observados e esperados (i.e. aqueles providos pelo próprio modelo) de uma variável 
resposta é amplamente utilizado como técnica de ajuste dos modelos para um conjunto 
de dados. A regressão linear simples é uma análise bastante comum e a compreensão da 
sua lógica auxilia no entendimento de outros tipos de análises mais complexas.

Na animação quando o tamanho dos resíduos aumenta, ao mudar a inclinação da reta, os 
segmentos são plotados em vermelho e quando os resíduos diminuem, eles são plotados 
em verde. No primeiro plot os segmentos estão em preto, uma vez que não há um valor 
anterior para comparação. O ponto vermelho mostrado é ponto do fulcro, com os valores 
das médias para as variáveis x e y. A reta azul que aparece no final da animação é 
aquela com o melhor ajuste.

Essa função utiliza a função 'saveGIF' do pacote 'animation' para converter diferentes 
imagens em um único arquivo usando o software ImageMagick 
(disponível em http://www.imagemagick.org/script/index.php).

Value:

anim.msd retorna uma janela gráfica com uma animação em GIF e uma lista com o resultado 
da função summary() aplicada à função lm().

Warning:

A função é interrompida e mensagens de erro são retornadas quando:

1) o pacote 'animation' não é encontrado ao usar a função library() que procura os 
pacotes instalados.
2) os vetores x e y não tem o mesmo tamanho.
3) existem NAs nos vetores x e/ou y.

Author(s):

Karina Campos Tisovec Dufner
katisovec@gmail.com

References:

Gotelli, Nicholas J.; Ellison, Aaron M. (2004). Regression. In: A primer of ecological 
statistics. p. 240-287. Sinauer Associates. Sunderland, MA, U.S.A. 

Weisstein, Eric W. "Least Squares Fitting." From MathWorld--A Wolfram Web Resource.
http://mathworld.wolfram.com/LeastSquaresFitting.html

Least squares. From Wikipedia, the free encyclopedia.
https://en.wikipedia.org/wiki/Least_squares

Ordinary least squares. From Wikipedia, the free encyclopedia.
https://en.wikipedia.org/wiki/Ordinary_least_squares

See Also:

animation-package, saveGIF(), summary(), lm().

Examples:

## criando os objetos x e y
x= rnorm(100, 10, 1)
y= rnorm(100, 100, 2)

## rodando a função
anim.msd(x, y)

## mudando o limite superior e inferior da gama de inclinações apresentadas
anim.msd(x, y, lim=0.5)
anim.msd(x, y, lim=5)
anim.msd(x, y, lim=20, int=0.5)

## mudando a magnitude da diferença entre uma inclinação e a seguinte
anim.msd(x, y, int=0.01)
anim.msd(x, y, lim=10, int=10)
anim.msd(x, y, lim=10, int=1)

## mudando a velocidade da animação
anim.msd(x, y, time=0.01)
anim.msd(x, y, time=2)

## visualizando a diferença dos resíduos para cada um deles
anim.msd(x, y, cor="ind")