--- /dev/null
+expertsArray = c(
+ "CLM", #1
+ "GAM", #2
+ "CLM1", #3
+ "CLM2", #4
+ "S_AIRPARIF", #5
+ "S_INERIS", #6
+ "D_ESMERALDA", #7
+ "D_PREVAIR", #8
+ "D_PREVAIR2", #9
+ "PERSIST", #10
+ #new additions only for PQV_2014
+ #TODO: default behavior on station != PQV_2014 ?
+ "GAM_sepMar", #11
+ "GAM_aprAug", #12
+ "GAM_highPollution", #13
+ "GAM_lowPollution", #14
+ "GAM_hotTemperature", #15
+ "GAM_coldTemperature", #16
+ "GAM_eastWind", #17
+ "GAM_westWind", #18
+ "GAM_noRain", #19
+ "GAM_rain" #20
+)
+
+stationsArray = c(
+ "AIL", #1
+ "ALE", #2
+ "CAE", #3
+ "CHD", #4
+ "EVT", #5
+ "HRI", #6
+ "IFS", #7
+ "JUS", #8
+ "LIS", #9
+ "MAS", #10
+ "MRA", #11
+ "NEI", #12
+ "POS", #13
+ "PQV", #14
+ "SLO", #15
+ "HRI_2014", #16
+ "LIS_2014", #17
+ "PQV_2014", #18
+ "PQV2" #19
+)
+
+#' @title Get forecasts + observations
+#'
+#' @description Get forecasts of all specified experts for all specified stations, also with (ordered) dates and (unordered) stations indices.
+#'
+#' @param experts Names of the experts. Default: all
+#' @param station Names of the stations. Default: all
+#'
+#' @export
+getData = function(experts=expertsArray, stations=stationsArray)
+{
+ #no need because of "LazyData: true" in DESCRIPTION
+ #data(list=stations, package="aggexp")
+ data = as.data.frame(matrix(nrow=0, ncol=1 + length(experts) + 2))
+ names(data) = c("Date", experts, "Measure", "Station")
+ for (i in 1:length(stations))
+ {
+ stationInfo = get(stations[i])
+ #date index is sufficient; also add station index
+ stationInfo = cbind(Date = 1:nrow(stationInfo), stationInfo[,names(stationInfo) %in% experts], Measure = stationInfo[,"Measure"], Station = i)
+ data = rbind(data, stationInfo)
+ }
+
+ #extra step: order by date (would be a DB request)
+ data = data[order(data[,"Date"]),]
+
+ return (data)
+}