Momento R do Dia – Ver a página do “Family Guy” na Wikipedia te leva a ver a do “The Simpsons”, mas não o contrário

Momento R do Dia – Crossover?

Seguindo a brincadeira das últimas horas (originalmente a coisa toda começou aqui), resolvi fazer um experimento com dois verbetes favoritos do pessoal de língua inglesa: “The_Simpsons” e “Family_Guy”. Adianto que o meu favorito é o último. A inspeção visual inicial nos mostra uma possível correlação entre as séries.

guysimpsons1

Como sabemos, a correlação não quer dizer causalidade (e nem a correlação parece ser algo lá muito importante, né?). Assim, verifiquei a estacionaridade das séries e, para a possível surpresa de alguns, descobri que ambas possuem uma raiz unitária. Logo, tive que diferenciá-las (o gráfico das mesmas está abaixo).

diffs

A brincadeira foir fazer uns testes de causalidade de Granger entre as duas séries diárias de pageviews. Arbitrariamente, usei 4, 8 e 12 defasagens. Adivinhem só o resultado?

A hipótese nula de que Simpsons não-Granger causam Family Guy não é rejeitada sob qualquer grau de significância padrão. Já a hipótese nula de que Family Guy não-Granger causa Simpsons é rejeitada (a 5% para a especificação com quatro defasagens e a 1% para as com oito e doze defasagens).

O que isto significa, meu Deus???

Grosso modo, quem visualiza o verbete do Family Guy vai lá ver mais sobre os Simpsons, mas não o contrário. Por que será que isto ocorre? Não sei. Palpites?

Agora, o que fazer? Meu mundo acabou? O que faltou?

O que eu não fiz? Não fiz várias coisas. Dentre outras, por exemplo, um teste de cointegração (pensando bem, teria sido divertido…). Poderia ter feito mas, enfim, não era a minha intenção. O código? Veja aí.

library(wikipediatrend)

# exemplo 1

page_views <- 
  wp_trend( 
    page = c( "Family_Guy", "The_Simpsons") ,
    from = "2008-01-01",
    to   = prev_month_end(),
    file = "C:/Users/cdshi_000/Documents/Meus Documentos/Meus Documentos/guysimpsons.csv"
  )

guys<- read.zoo("C:/Users/cdshi_000/Documents/Meus Documentos/Meus Documentos/guysimpsons.csv",header=TRUE,sep=",",format = "%m/%d/%Y")
head(guys)
summary(guys)

library(forecast)
ndiffs(guys$family_guy)
ndiffs(guys$the_simpsons)

plot(guys)
plot(diff(guys))

# Simple Scatterplot (não apresentado)

plot(guys$family_guy, guys$the_simpsons, main="Scatterplot Example", 
     xlab="Family Guy ", ylab="The Simpsons", pch=19)

abline(lm(guys$the_simpsons~guys$family_guy), col="red") # regression line (y~x) 
lines(lowess(guys$family_guy,guys$the_simpsons), col="blue") # lowess line (x,y)

# a correlação apresentada

library(car) 
scatterplot(the_simpsons~family_guy, data=guys, 
            main="Enhanced Scatter Plot")

# Os testes de causalidade...

library(lmtest)
grangertest(diff(family_guy) ~diff(the_simpsons), order=4, data=guys)
grangertest(diff(the_simpsons) ~ diff(family_guy) , order=4, data=guys)

grangertest(diff(family_guy) ~diff(the_simpsons), order=8, data=guys)
grangertest(diff(the_simpsons) ~ diff(family_guy) , order=8, data=guys)

grangertest(diff(family_guy) ~diff(the_simpsons), order=12, data=guys)
grangertest(diff(the_simpsons) ~ diff(family_guy) , order=12, data=guys)

Algum comentário? Talvez Peter tenha um.

Anúncios

Deixe um comentário

Preencha os seus dados abaixo ou clique em um ícone para log in:

Logotipo do WordPress.com

Você está comentando utilizando sua conta WordPress.com. Sair / Alterar )

Imagem do Twitter

Você está comentando utilizando sua conta Twitter. Sair / Alterar )

Foto do Facebook

Você está comentando utilizando sua conta Facebook. Sair / Alterar )

Foto do Google+

Você está comentando utilizando sua conta Google+. Sair / Alterar )

Conectando a %s