kryzysyDem1 <- function(ramka,stopienGrDolna = -Inf, stopienGrGorna = Inf) { czas <- ramka[,1] zdarzenia <- ramka[,2] koniecWektoraIter <- length(zdarzenia)-10 jednostkiCzasuBez10pierwszych <- c() for(i in 1 : koniecWektoraIter) { jednostkiCzasuBez10pierwszych[i] <- czas[i+10] } jednostkiCzasuBez10pierwszych daneRok11iNastepneLata <- c() for(i in 1 : koniecWektoraIter) { daneRok11iNastepneLata[i] <- zdarzenia[i+10] } daneRok11iNastepneLata wynikiSredniejAryt <- c() for(i in 1 : koniecWektoraIter) { wynikiSredniejAryt[i] <- mean(zdarzenia[i:(i+9)]) } wynikiSredniejAryt wynikiOdchyleniaStand <- c() for(i in 1: koniecWektoraIter) { wynikiOdchyleniaStand[i] <- sd(zdarzenia[i:(i+9)]) } wynikiOdchyleniaStand wynikiDane <- (daneRok11iNastepneLata - wynikiSredniejAryt) / wynikiOdchyleniaStand wynikiWybor <- ifelse(wynikiDane >= stopienGrDolna & wynikiDane < stopienGrGorna, round(wynikiDane, 3),"") data.frame(jednostkiCzasuBez10pierwszych, wynikiWybor) }