some small changes

issuecomp-3-calc.R

@@ -81,13 +81,19 @@ g1
 # 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:44], p=1, type="none")
 # VAR(issues_s[,2:23], p=1, type=c("const", "trend", "both", "none"), season=NULL, exogen = issues_i[2:22])
 
+issues_ts <- as.ts(issues)
+vIssues <- VAR(issues_ts[,2:44], p=1, type="both")
+
 plot(irf(test, impulse = names(issues_s[2:23]), response = names(issues_i[2:22])))
 
 capture.output(print(summary(test), prmsd=TRUE, digits=1), file="out.txt")
 
+# Tests
+issues_ts <- as.ts(issues)
+VARselect(issues[2:44], lag.max = 8, type = "both")
+summary(ur.df(issues_ts[, 2], type ="none", lags=1))
 
 # SOME TESTS --------------------------------------------------------------
 
@@ -123,6 +129,16 @@ pie(acc_parties$twitter, col=c("black", "red", "purple", "green"), labels = c("C
 rm(acc_parties, p)
 
 
+# Count all tags
+num <- 0
+for(i in 1:length(issuelist)) {
+  j <- length(issuelist[[i]])
+  num <- num + j
+  rm(j)
+}
+num
+
+
 # VISUALS -----------------------------------------------------------------