mxmehl/uni-ba-socialagenda
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Corrent status

Browse Source
This commit is contained in:
mxmehl
2015-03-04 11:36:31 +01:00
parent a388110551
commit 9bbf1b4f56
16 changed files with 2713 additions and 494 deletions
Download Patch File Download Diff File Expand all files Collapse all files
issuecomp-3-calc.R
+7 -2
View File
@@ -72,15 +72,20 @@ stats_entropy <- melt(stats_entropy, id="date")
g1 <- ggplot(data = stats_entropy, aes(x=date,y=value,colour=variable, group=variable)) +
geom_line() +
geom_smooth(size=1,formula = y ~ x, method="loess", se=FALSE, color=1)
g1
# g1
# VAR ---------------------------------------------------------------------
test <- VAR(issues[,2:32], p=3, type=c("const", "trend", "both", "none"), season=NULL, exogen = NULL, lag.max = NULL, ic = c("AIC", "HQ", "SC", "FPE"))
test <- VAR(issues[,2:32], p=1, type=c("const", "trend", "both", "none"), season=NULL, exogen = NULL, lag.max = NULL, ic = c("AIC", "HQ", "SC", "FPE"))
test <- VAR(issues_i[,2:22], p=1, type="none", exogen = issues_s[,2:3])
test <- VAR(issues_s[,2:11], p=1, type="none")
test <- VAR(issues[,2:32], p=1, type="none")
VAR(issues_s[,2:11], p=1, type=c("const", "trend", "both", "none"), season=NULL, exogen = issues_i[2:22])
plot(irf(test, impulse = names(issues_s[2:11]), response = names(issues_i[2:22])))
capture.output(print(summary(test), prmsd=TRUE, digits=1), file="out.txt")