Bom dia,
qual é o pacote mais adequado para teste de causalidade de granger?
o FIAR?
*
Type Package
Title Functional Integration Analysis in R
Version 0.4
Date 2011-04-01
Author Bjorn Roelstraete
Maintainer Bjorn Roelstraete
Pode ajudar.
http://en.wikipedia.org/wiki/Granger_causality
Here is an example of the function grangertest() in the lmtest library
of the R software:
Granger causality test
Model 1: fii ~ Lags(fii, 1:5) + Lags(rM, 1:5)
Model 2: fii ~ Lags(fii, 1:5)
Res.Df Df F Pr(F)
1629
2
Necessito gerar mapas a partir de X (LONGITUDE) Y (LATITUDE) Z(teor de umidade
no solo). O mapa deve ser 2D e a variável Z pretendo representar por cores
através de uma legenda, ex.: 12-16% azul ; 16,01 - 18% verde ; etc
Testei o GeoR e o adehabitatMA, porém nao tive os resultados pretendidos.
Dei uma googlada com maps with R e qualhou e resultados. Só na primeira
página tem vários que podem ser úteis para você, que incluem exemplos
reproduzíveis. Caso seu problema seja mais específico, envie um CMR.
À disposição.
Walmes.
Vitor,
você não mandou um CMR, mas assumindo que as colunas são comparadas
duas a duas (coluna1 x coluna2, coluna3 x coluna4, etc) o programa
abaixo deve resolver a situação.
Serão gerados 30 vetores na área de trabalho.
Faça os testes.
# Exemplo de dados
dados-trunc(rnorm(30*10,30,5))
Olá Benilton.
No manual encontrei instruções de como se criar uma função em um novo
pacote.
No meu caso, onde apenas modifiquei a função eigen.f para eigem.f, gostaria
de saber como fazer para carregá-la no momento da execução.
Quando chamo a função eigen:
eigen(cbind(c(1,-1),c(-1,1)))
Só alterando a função não quer dizer que elá estará disponível para
uso. Você precisa avisar o R de alguma maneira que existe uma
modificação em um pacote já compilado.
A primeira tentativa seria:
$ R CMD COMPILE eigem.f
e tente usar a função. Se não ser certo vc pode tentar
$ R CMD build
E, como ja citado aqui anteriormente, apos compilar, carregar a
biblioteca dinamica via dyn.load() e chamar a funcao propriamente dita
via .Fortran() ... adicionalmente, volto a repetir que nao eh uma boa
ideia criar funcoes com o mesmo nome de funcoes existentes no base ou
qq outro pacote
Olá Paulo,
muito obrigado.
Como eu estou apenas aprendendo programação, testei o seu código mas ainda
estou tentando compreender o resultado. Como assim gerar 30 vetores na área de
trabalho?
Pensei em outra coisa que também resolveria minha situação: substituir todos os
valores que aparecem