Code: Alles auswählen
data <- structure(list(Jahr = 1989:2014, X0 = c(0.003023015, 0.002927452,
0.002297123, 0.002815327, 0.002763109, 0.002838143, 0.003664018,
0.003178668, 0.004350883, 0.003850793, 0.003509275, 0.003250701,
0.003345808, 0.00308408, 0.002670496, 0.003213876, 0.002821363,
0.002765209, 0.002883229, 0.003257329, 0.002789505, 0.002373275,
0.002696899, 0.002846016, 0.002207287, 0.002817528), X1 = c(0.000652742,
0.000653723, 0.000544034, 0.000927038, 0.001038739, 0.000935872,
0.000867124, 0.000804148, 0.001057385, 0.000658921, 0.0007864,
0.000617871, 0.000990958, 0.000884782, 0.001689317, 0.001002813,
0.001014877, 0.000989908, 0.00128063, 0.000865994, 0.000754945,
0.001154167, 0.000881523, 0.00112217, 0.000850712, 0.001031472
), X2 = c(0.00042304, 0.000844688, 0.000600304, 0.000363535,
0.000561097, 0.000485496, 0.000690552, 0.000679026, 0.000592129,
0.000638541, 0.000591872, 0.000782869, 0.000995662, 0.000910635,
0.000902703, 0.001273949, 0.000895537, 0.000809995, 0.001178394,
0.000833681, 0.000960142, 0.001103172, 0.000783296, 0.000926111,
0.000895559, 0.001037169), X3 = c(0.000461802, 0.000814171, 0.001182525,
0.000652673, 0.000547784, 0.000557483, 0.000564596, 0.000627238,
0.000555007, 0.000486309, 0.000527241, 0.00065475, 0.001217803,
0.001132075, 0.001027046, 0.000919007, 0.000989144, 0.001008864,
0.000856164, 0.000791595, 0.000806511, 0.001251564, 0.001025872,
0.000949991, 0.001043764, 0.00062494), X4 = c(0.000630476, 0.000501122,
0.000738104, 0.000818031, 0.000602898, 0.00054448, 0.000639426,
0.000661973, 0.000545907, 0.000514097, 0.000654533, 0.00074408,
0.000805854, 0.000868949, 0.001474064, 0.001091888, 0.001366087,
0.001446397, 0.001219804, 0.001290574, 0.00100356, 0.001121158,
0.000881144, 0.000996314, 0.000896383, 0.000815721), X5 = c(0.000603513,
0.000691055, 0.00056014, 0.000512908, 0.000766006, 0.000684505,
0.000647236, 0.000597064, 0.000721674, 0.00068741, 0.000677632,
0.000778816, 0.000852757, 0.00108279, 0.001248553, 0.001091982,
0.000793574, 0.00138092, 0.001270741, 0.001145639, 0.001505768,
0.001264403, 0.001169591, 0.001195151, 0.001215398, 0.000987833
), X6 = c(0.000560969, 0.000558024, 0.000985095, 0.000510345,
0.000635042, 0.000784813, 0.000581287, 0.000604393, 0.000411286,
0.000621841, 0.000605938, 0.000778051, 0.000925439, 0.001219618,
0.000688987, 0.001015893, 0.001106858, 0.001194172, 0.00144759,
0.001393064, 0.001214633, 0.001018577, 0.001120381, 0.001093534,
0.001021525, 0.000987484)), .Names = c("Jahr", "X0", "X1", "X2",
"X3", "X4", "X5", "X6"), class = "data.frame", row.names = c(NA,
-26L))
###
limJahr <- 2000
pred.Jahr <- 2001:2014
mypred <- function(x, predJ) {
m <- lm(x ~ Jahr, data=data.frame(x=x, Jahr=limData$Jahr))
predict(m, newdata=data.frame(Jahr=predJ))
}
limData <- data[data$Jahr<=limJahr,]
d.pred <- sapply(limData[-1], mypred, predJ=pred.Jahr)
cbind(data[data$Jahr %in% pred.Jahr,], d.pred) # a)
rbind(data[data$Jahr %in% pred.Jahr,], cbind(Jahr=pred.Jahr, d.pred)) # b)