[R] alternative for multiple if_else statements
S Ellison
S.Ellison at LGCGroup.com
Mon Feb 26 12:56:39 CET 2018
That many ifelse statements is obviously rather a pain.
Would you not have got what you want with
... paste("survey", year, sep="_")
?
If that is not what you're looking for (eg because 'year' is the observation year and not the study start year), perhaps something that picks the minimum year for a subject or other relevant group might work? For example
paste("survey", ave(year, studyno, FUN=min), sep="_")
S Ellison
> -----Original Message-----
> From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Kevin
> Wamae
> Sent: 21 February 2018 20:34
> To: R-help at r-project.org
> Subject: [R] alternative for multiple if_else statements
>
> Hi, I am having trouble trying to figure out why if_else is behaving the way it is,
> it may be my code or the way the data is structured.
>
> Below is a snapshot of a database am working on and it represents a
> longitudinal survey of study participants in a trial with weekly follow up.
>
> The variable "survey_start" represents the start of the study-defined one year
> follow up (which we called "survey_year").
>
> I am trying to populate all subsequent entries for each participant, per survey
> year, with the entry "survey" followed by an underscore and the respective
> year, eg. survey_2014.
>
> There are missing entries such as the participant represented here, wasn't
> available at the start of the 2015 survey. Also, some participants don’t have
> complete one-year follow ups but I still need to include them.
>
> I have written two codes, first one fails while the second works, the only
> difference being I have reversed the order in which the entries are populated in
> the second code (from 2007-2016 to 2016-2007) and removed the if_else
> statement for 2015. Also noticed, that for the second code, which spans the
> years 2007-2016 (less 2015), if a participants entries start from 2010-2016, the
> code fails.
>
> Kindly assist in figuring this out...or better yet, an alternative.
>
> trialData <- structure(list(study = c("site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1", "site_1",
> "site_1"), studyno = c("child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1", "child_1", "child_1", "child_1", "child_1", "child_1", "child_1",
> "child_1"), date = structure(c(16078, 16085, 16092, 16098, 16104, 16115,
> 16121, 16129, 16135, 16140, 16146, 16156, 16162, 16168, 16177, 16185,
> 16191, 16195, 16203, 16210, 16217, 16225, 16234, 16237, 16246, 16253,
> 16262, 16269, 16278, 16283, 16288, 16297, 16304, 16311, 16319, 16326,
> 16332, 16337, 16346, 16353, 16360, 16366, 16370, 16381, 16384, 16395,
> 16399, 16407, 16415, 16422, 16444, 16452, 16454, 16467, 16474, 16477,
> 16484, 16490, 16501, 16508, 16514, 16520, 16529, 16533, 16539, 16550,
> 16556, 16564, 16566, 16578, 16582, 16593, 16599, 16604, 16613, 16620,
> 16623, 16635, 16636, 16654, 16660, 16666, 16673, 16681, 16688, 16693,
> 16702, 16706, 16714, 16721, 16728, 16734, 16745, 16749, 16757, 16764,
> 16769, 16778, 16785, 16792, 16805, 16812, 16819, 16830, 16832, 16839,
> 16846, 16856, 16862, 16867, 16877, 16884, 16890, 16898, 16904, 16912,
> 16917, 16923, 16936, 16938, 16953, 16960, 16966, 16973, 16980), class =
> "Date"), year = c(2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L,
> 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L,
> 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L,
> 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L,
> 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2014L, 2015L, 2015L,
> 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L,
> 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L,
> 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L,
> 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L, 2015L,
> 2015L, 2015L, 2015L, 2015L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L,
> 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L,
> 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L), month = c(1L, 1L, 1L, 1L, 2L,
> 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L,
> 7L, 7L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 11L, 11L, 11L,
> 11L, 12L, 12L, 12L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L,
> 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 10L,
> 10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
> 2L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L), survey_start =
> c("", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "Y", "", "", "", "", "", "", "",
> "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
> "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
> "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "", "",
> "", "", "", "", "", "", "", "Y", "", "", "", "", "", "", "", "", "", "", "", "", "", "")), class =
> "data.frame", row.names = c(NA, -125L), .Names = c("study", "studyno", "date",
> "year", "month",
> "survey_start"))
>
>
> code 1 fails:
>
> trialData <- trialData %>% arrange(studyno, date) %>% group_by(studyno) %>%
> mutate(survey_year = if_else(date >= date[survey_start == "Y" & year == 2007
> & study == "site_1"][1] & date < date[month == 5 & year == 2008 & study ==
> "site_1"][1], "survey_2007",
> if_else(date >= date[survey_start == "Y" & year == 2008 & study
> == "site_1"][1] & date < date[month == 4 & year == 2009 & study ==
> "site_1"][1], "survey_2008",
> if_else(date >= date[survey_start == "Y" & year == 2009 & study
> == "site_1"][1] & date < date[month == 5 & year == 2010 & study ==
> "site_1"][1], "survey_2009",
> if_else(date >= date[survey_start == "Y" & year == 2010 & study
> == "site_1"][1] & date < date[month == 5 & year == 2011 & study ==
> "site_1"][1], "survey_2010",
> if_else(date >= date[survey_start == "Y" & year == 2011 & study
> == "site_1"][1] & date < date[month == 4 & year == 2012 & study ==
> "site_1"][1], "survey_2011",
> if_else(date >= date[survey_start == "Y" & year == 2012 & study
> == "site_1"][1] & date < date[month == 4 & year == 2013 & study ==
> "site_1"][1], "survey_2012",
> if_else(date >= date[survey_start == "Y" & year == 2013 & study
> == "site_1"][1] & date < date[month == 4 & year == 2014 & study ==
> "site_1"][1], "survey_2013",
> if_else(date >= date[survey_start == "Y" & year == 2014 & study
> == "site_1"][1] & date < date[month == 4 & year == 2015 & study ==
> "site_1"][1], "survey_2014",
> if_else(date >= date[survey_start == "Y" & year == 2015 & study
> == "site_1"][1] & date < date[month == 3 & year == 2016 & study ==
> "site_1"][1], "survey_2015",
> if_else(date >= date[survey_start == "Y" & year == 2016 & study
> == "site_1"][1], "survey_2016","")))))))))))
>
> code 2 works:
>
> trialData <- trialData %>% arrange(studyno, date) %>% group_by(studyno)
> %>%
> mutate(survey_year = if_else(date >= date[survey_start == "Y" & year == 2016
> & study == "site_1"][1] , "survey_2016",
> if_else(date >= date[survey_start == "Y" & year == 2014 &
> study == "site_1"][1] & date < date[month == 4 & year == 2015 & study ==
> "site_1"][1], "survey_2014",
> if_else(date >= date[survey_start == "Y" & year == 2013 &
> study == "site_1"][1] & date < date[month == 4 & year == 2014 & study ==
> "site_1"][1], "survey_2013",
> if_else(date >= date[survey_start == "Y" & year == 2012 &
> study == "site_1"][1] & date < date[month == 4 & year == 2013 & study ==
> "site_1"][1], "survey_2012",
> if_else(date >= date[survey_start == "Y" & year == 2011 &
> study == "site_1"][1] & date < date[month == 4 & year == 2012 & study ==
> "site_1"][1], "survey_2011",
> if_else(date >= date[survey_start == "Y" & year == 2010 &
> study == "site_1"][1] & date < date[month == 5 & year == 2011 & study ==
> "site_1"][1], "survey_2010",
> if_else(date >= date[survey_start == "Y" & year == 2009 &
> study == "site_1"][1] & date < date[month == 5 & year == 2010 & study ==
> "site_1"][1], "survey_2009",
> if_else(date >= date[survey_start == "Y" & year == 2008 &
> study == "site_1"][1] & date < date[month == 4 & year == 2009 & study ==
> "site_1"][1], "survey_2008",
> if_else(date >= date[survey_start == "Y" & year == 2007 &
> study == "site_1"][1] & date < date[month == 5 & year == 2008 & study ==
> "site_1"][1], "survey_2007",""))))))))))
>
> ________________________________________________________________
> ______
>
> This e-mail contains information which is confidential. It is intended only for the
> use of the named recipient. If you have received this e-mail in error, please let
> us know by replying to the sender, and immediately delete it from your system.
> Please note, that in these circumstances, the use, disclosure, distribution or
> copying of this information is strictly prohibited. KEMRI-Wellcome Trust
> Programme cannot accept any responsibility for the accuracy or completeness
> of this message as it has been transmitted over a public network. Although the
> Programme has taken reasonable precautions to ensure no viruses are present
> in emails, it cannot accept responsibility for any loss or damage arising from
> the use of the email or attachments. Any views expressed in this message are
> those of the individual sender, except where the sender specifically states them
> to be the views of KEMRI-Wellcome Trust Programme.
> ________________________________________________________________
> ______
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
*******************************************************************
This email and any attachments are confidential. Any use, copying or
disclosure other than by the intended recipient is unauthorised. If
you have received this message in error, please notify the sender
immediately via +44(0)20 8943 7000 or notify postmaster at lgcgroup.com
and delete this message and any copies from your computer and network.
LGC Limited. Registered in England 2991879.
Registered office: Queens Road, Teddington, Middlesex, TW11 0LY, UK
More information about the R-help
mailing list