as.numeric(covidNZ$Case_date) cumsum(covidNZ) hcases["x"][hcases["selected"]] Female <- covidNZ$Sex == "Female" nFemale <- sum(Female) cFemale <- nFemale:1 Male <- covidNZ$Sex == "Male" nMale <- sum(Male) cMale <- nMale:1 length(hcases["x"]) cumulativeCases <- numeric(length = nMale + nFemale) cumulativeCases[Female] <- cFemale cumulativeCases[Male] <- cMale p <- l_plot(covidNZ$Case_date, cumulativeCases, linkingGroup = "CovidNZ") colFemale <- hcases["color"][Female][1] l_layer_line(p, x = covidNZ$Case_date[Female], y = cumulativeCases[Female], color = colFemale, linewidth = 5) colMale <- hcases["color"][Male][1] l_layer_line(p, x = covidNZ$Case_date[Male], y = cumulativeCases[Male], color = colMale, linewidth = 5, dash = c(10, 10)) # To see what states are "linked" l_getLinkedStates(p) l_getLinkedStates(hcases) l_getLinkedStates(hsex)