Econometria Não-paramétrica no "R"

Sei que o que postarei é algo simples, e passível de conter algum erro, por isso peço desculpas aos mais avançados e motivação aos iniciantes "como eu".

O que segue abaixo é uma rotina simples no R, onde vou gerar duas variáveis e ajustar uma linha usando métodos não-paramétricos

em MAIUSCULO são feitos comentários que não fazem parte da rotina, apenas explicam o que está sendo realizado. Então, se você quer acompanhar o procedimento, copie o que está em "minusculo" e cole no programa "R"

Devemos primeiro carregar o pacote, em “pacotes”, “carregar pacote....”, “KernSmooth”

n = 200
NUMERO DE OBSERVAÇÕES

d = 0.1
DESVIO

e = rnorm(n,0,d)
ERRO, DISTRIBUIÇÃO RANDOMICA, MÉDIA ZERO E VARIÂNCIA DADA CONSTANTE = O.1

x = seq(1,10,length=n)
VARIÁVEL MONTADA ATRAVÉS DE UMA SEQUENCIA DE VALORES DE 1 A 10, DE TAMANHO N (200)

y = sin(x) + e
VARIÁVEL MONTADA COM UMA FUNÇÃO SENO DA ANTERIOR MAIS UM TERMO DE ERRO

plot(x,y)
FUNÇÃO PARA GERAR O GRÁFICO DAS DUAS VARIÁVEIS, PERCEBA O FORMATO “SENÓIDE” DESSE PLOT, DEVIDO AOS VALORES DE Y.

lines(locpoly(x,y,band=0.008),col="black")
LOCPOLY É A FUNÇÃO USADA PELO PACOTE KERNSMOOTH, BAND É O TAMANHO DA JANELA, COL = A COR DA LINHA. PERCEBA QUE A JANELA É MUITO PEQUENA, A LINHA QUE LIGA OS PONTOS NÃO APRESENTA UM FORMATO “SUAVE” ADAPTADO A FORMA QUE OS DADOS ESTÃO DISPOSTOS.

ABAIXO GRÁFICO, VALORES SIMULADOS SÃO OS PONTOS, ESTIMAÇÃO É A LINHA



lines(locpoly(x,y,band=0.5),col="blue")
COM O AUMENTO DA JANELA (DE 0.008 P/ O.5) A LINHA AJUSTA-SE DE MELHOR MANEIRA AOS DADOS.



lines(locpoly(x,y,band=0.3),col="yellow")
PERCEBA QUE COM A Diminuiçã DA JANELA A LINHA AJUSTA-SE AOS DADOS DE UMA FORMA "bOA" TAMB.



lines(locpoly(x,y,band=0.1),col="blue")
MAIS UMA DIMINUIÇÃO DA JANELA E A LINHA AJUSTA-SE DE UMA MANEIRA RUIM, POIS FICA "DEFORMADA". O QUE SUGERE QUE OS VALORES DE BOM AJUSTE ESTÃO ENTRE OS ANTERIORES ACIMA.




____________ESCOLHA DA JANELA - H____________

*QUANDO H VAI PARA INFINITO POLINÔMIO LOCAL IGUAL A AJUSTE LOCAL (EXEMPLO DA RETA)
EXEMPLO: PARA RETA, COM JANELA GRANDE
lines(locpoly(x,y,degree=1,band=10000),col="black")
CASO SEJA USADO UM POLINÔMIO DE GRAU DOIS E UMA JANELA GRANDE O GLOBAL SERIA UM POLINÔMIO E SEGUNDO GRAU, EXEMPLO:
lines(locpoly(x,y,degree=2,band=10000),col="red")
_______DEGREE É O GRAU DO POLINÔMIO QUE VAI SE AJUSTAR..

_____É POSSÍVEL INCLUIR OUTRO GRAU QUE NÃO O AUTOMÁTICO.....FICANDO, POR EXEMPLO, ASSIM:
lines(locpoly(x,y,degree=7,band=0.4),col="black")