plot(DATA[patient$WALKING,]$AnkleLeft.y,
     type = "l",
     ylab = "Hoogte",
     col = ifelse(DATA$AnkleLeft.y > -1.1, "green", "red"),
     ylim = c(-1.15, -0.9)
)
par(new=TRUE)
plot(DATA[patient$WALKING,]$AnkleRight.y,
     type = "l",
     ylab = "Hoogte",
     col = ifelse(DATA$AnkleRight.y > -1.1, "blue", "purple"),
     ylim = c(-1.15, -0.9)
)

#Calculates the mean height of AnkleLeft.y
meanLine <- mean(DATA[patient$WALKING,]$AnkleLeft.y)

##-LEFT-##
#Creates vector which contains the values above the mean height (1e quarter)
Q1 <- NULL
for(i in patient$WALKING){
  if (DATA[i,]$AnkleLeft.y > meanLine){
    Q1 <- c(Q1, DATA[i,]$AnkleLeft.y)
  }
}

#Gives the mean value of the first quarter
meanQ1 <- mean(Q1)

#Creates vector which containts the values under the mean height value (3e quarter)
Q3 <- NULL
for(i in patient$WALKING){
  if(DATA[i,]$AnkleLeft.y < meanLine){
    Q3 <- c(Q3, DATA[i,]$AnkleLeft.y)
  }
}

#Gives the mean value of the third quarter
meanQ3 <- mean(Q3)

DiffAnkleLeft <- abs(meanQ3 - meanQ1)


##-RIGHT-##
#Creates vector which contains the values above the mean height (1e quarter)
Q1 <- NULL
for(i in patient$WALKING){
  if(DATA[i,]$AnkleRight.y > meanLine){
    Q1 <- c(Q1, DATA[i,]$AnkleRight.y)
  }
}

#Gives the mean value of the first quarter
meanQ1 <- mean(Q1)

#Creates vector which contains the values above the mean height (3e quarter)
Q3 <- NULL
for(i in patient$WALKING){
  if(DATA[i,]$AnkleRight.y < meanLine){
    Q3 <- c(Q3, DATA[i,]$AnkleRight.y)
  }
}

#Gives the mean value of the third quarter, And the difference
meanQ3 <- mean(Q3)

DiffAnkleRight <- abs(meanQ3 - meanQ1)

DiffAnkleLeft
DiffAnkleRight

if(DiffAnkleLeft < 0.05 || DiffAnkleRight < 0.05){
  print("The patient barely lift his foot up, there is a potentional falling risk")
}else {
  print("The patient walks just fine according his ankles")
}

#Creates Smoothline of the left ankle
SmoothAnkleLeft <- smooth.spline(DATA[patient$WALKING,]$AnkleLeft.y, spar=0.35)
plot(SmoothAnkleLeft, col = "green",
     type = "l",
     ylim = c(-1.17, -0.9),
     ylab = "Hoogte",
     xlab = "Index")


par(new=TRUE)
SmoothAnkleRight <- smooth.spline(DATA[patient$WALKING,]$AnkleRight.y, spar=0.35)
plot(SmoothAnkleRight, col = "blue",
     type = "l",
     ylim = c(-1.17, -0.9),
     ylab = "Hoogte",
     xlab = "Index")