Economia
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")
- Tristão Da Silva
Aquela janela virada para o mar... Cem anos que eu viva não posso esquecer-me Daquele navio que eu vi naufragar Na boca da barra tentando perder-me E aquela janela virada para o mar Crédito: Música Anos 60 ...
- Uma Dúvida
A indicação de hiperligações pode ser feita abrindo a nova janela em substituição da que se está a consultar (exemplo), ou abrindo uma nova janela (exemplo). Tenho encontrado os dois sistemas em blogues e gostaria de saber a opinião dos leitores...
- Flutuações Macroeconômicas No Brasil
O gráfico apresenta o desvio percentual do PIB corrente brasileiro em relação ao PIB potencial medido segundo o método de Hodrik-Prescott. Os dados cobrem o período 1947-2006. A linha horizontal corresponde ao pleno emprego. Observações abaixo...
- Material Escolar E O Custo TributÁrio Na Volta Ás Aulas
A volta as aulas é motivo de preocupação para os pais. Afinal de contas, o custo do material escolar pesa no bolso dos brasileiros. Como se não bastasse a quantidade consumida nas escolas, é imortante que os responsáveis tenham conhecimento de...
- Flutuações Macroeconômicas No Brasil
O gráfico apresenta o desvio percentual do PIB corrente brasileiro em relação ao PIB potencial medido segundo o método de Hodrik-Prescott. Os dados cobrem o período 1947-2006. A linha horizontal corresponde ao pleno emprego. Observações abaixo...
Economia
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")
loading...
- Tristão Da Silva
Aquela janela virada para o mar... Cem anos que eu viva não posso esquecer-me Daquele navio que eu vi naufragar Na boca da barra tentando perder-me E aquela janela virada para o mar Crédito: Música Anos 60 ...
- Uma Dúvida
A indicação de hiperligações pode ser feita abrindo a nova janela em substituição da que se está a consultar (exemplo), ou abrindo uma nova janela (exemplo). Tenho encontrado os dois sistemas em blogues e gostaria de saber a opinião dos leitores...
- Flutuações Macroeconômicas No Brasil
O gráfico apresenta o desvio percentual do PIB corrente brasileiro em relação ao PIB potencial medido segundo o método de Hodrik-Prescott. Os dados cobrem o período 1947-2006. A linha horizontal corresponde ao pleno emprego. Observações abaixo...
- Material Escolar E O Custo TributÁrio Na Volta Ás Aulas
A volta as aulas é motivo de preocupação para os pais. Afinal de contas, o custo do material escolar pesa no bolso dos brasileiros. Como se não bastasse a quantidade consumida nas escolas, é imortante que os responsáveis tenham conhecimento de...
- Flutuações Macroeconômicas No Brasil
O gráfico apresenta o desvio percentual do PIB corrente brasileiro em relação ao PIB potencial medido segundo o método de Hodrik-Prescott. Os dados cobrem o período 1947-2006. A linha horizontal corresponde ao pleno emprego. Observações abaixo...