嵌套ifelse以输出3个响应R

Question

这是我的原始帖子中的一个相关问题：如何根据r

我从实验中有2个数据帧。第一个DF在40分钟内读取A（大致）连续信号。有5列，1：3是二进制的 - 说是否按下按钮。第四列是按下第2列或第3列的二进制。第五列是秒数的大约时间。 的示例：

启动	左右	下面	l或r	时间
0	0 0	1	1 2	2.8225
0	0	1 1	1	1.82375
0	0 0	1	1 1	2.82500
0	0	1	1	2.82625 1 2.82625
1	0 0	0 0	0 0 0 0	16.8200
1	0 0 0 0	0 0 0	0	16.8212

等

DF 信息的信息是每行是一个试验，通常取决于当天的100-150行。我有一列标记试验开始时间和另一列的列，该列以秒为单位的试验结束时间。我有另一列说明试验是否有干预。 Example from df below (I omitted several irrelevant columns):

trial	control	t start	t end
1	0	16.64709	35.49431
2	0	41.81843	57.74304
3	0	65.54510	71.16612
4	0	82.65743	87.30914
11	3	187.0787	193.5898
12	0	200.0486	203.1883
30	3	415.1710	418.0405

etc.

For the第一个数据框架，我想创建一个列，以指示在试用中是否按下按钮。如果确实在试验中按下了按钮，则我需要根据干预措施对其进行标签。这是基于第二个DF中的开始和结束时间以及控制信息。在此表中，0 =干预和3 =控制。

我希望它看起来像这样（iti = trial，wt_int =在试验中，wt_control = in trial＆amp; control）：

启动	左右	右	l或r	time time	trip trial_type
0	0	1	1	1.8225	ITI
0 0	0 0 0	1 1	1	2.82375	ITI
0	0	1	1 1	2.82500	0
16.82125	1	2.82625	1	0	WT_INT
1	0	0	0	0	0 0 0 0 0 0 16.82000 WT_INT
1	0	0	0	0	0
ITI	0	0	0	187.0800	WT_CONTROL

等

ITI 0 0 1 1 。我可以将其标记为所有试验，为“ ITI”或“ WT_INT”，以不同的失败尝试或第1037行的错误（当它从ITI更改为WT时）。从我最初的问题中，我现在在我的第一df中有一个“试用”列，我将用于以下代码。也许有一种更简单的方法结合了原始代码？

错误的一部分：

df %>% 
  rowwise() %>% 
  mutate(trial_type = ifelse(any(trial == "wt" & df2$control == 0,
                                 ifelse(trial == "wt" & df2$control == 3,
                                        "wt_omission", "iti"), "wt_odor")))

还尝试了此功能，它标记为WT_INT：

df$trial_type <- ifelse(df$trial == 'wt' && df2$control == 0,
                        ifelse(df$trial == 'wt' && df2$control == 3, 
                               "wt_control", "iti"), "wt_int")

谢谢！

Answer 1

您可以使用剪切来创建间隔并检查一个值是否落入其中：

library(dplyr)

df1 %>% 
  mutate(
    check_1 = cut(time, breaks = df2$t_start, labels = FALSE),
    check_2 = coalesce(cut(time, breaks = df2$t_end, labels = FALSE), 0),
    check_3 = df2$control[check_1],
    trial_type = case_when(
      check_1 - check_2 == 1 & check_3 == 0 ~ "wt_int",
      check_1 - check_2 == 1 & check_3 == 3 ~ "wt_control",
      TRUE ~ "iti"
      )
    ) %>% 
  select(-starts_with("check_"))

这将返回

# A tibble: 7 x 6
  initiate  left right l_or_r   time trial_type
     <dbl> <dbl> <dbl>  <dbl>  <dbl> <chr>     
1        0     0     1      1   2.82 iti       
2        0     0     1      1   2.82 iti       
3        0     0     1      1   2.82 iti       
4        0     0     1      1   2.83 iti       
5        1     0     0      0  16.8  wt_int    
6        1     0     0      0  16.8  wt_int    
7        1     0     0      0 187.   wt_control

数据

df1 <- structure(list(initiate = c(0, 0, 0, 0, 1, 1, 1), left = c(0, 
0, 0, 0, 0, 0, 0), right = c(1, 1, 1, 1, 0, 0, 0), l_or_r = c(1, 
1, 1, 1, 0, 0, 0), time = c(2.8225, 2.82375, 2.825, 2.82625, 
16.82, 16.8212, 187.08)), class = c("spec_tbl_df", "tbl_df", 
"tbl", "data.frame"), row.names = c(NA, -7L), spec = structure(list(
    cols = list(initiate = structure(list(), class = c("collector_double", 
    "collector")), left = structure(list(), class = c("collector_double", 
    "collector")), right = structure(list(), class = c("collector_double", 
    "collector")), l_or_r = structure(list(), class = c("collector_double", 
    "collector")), time = structure(list(), class = c("collector_double", 
    "collector"))), default = structure(list(), class = c("collector_guess", 
    "collector")), skip = 1L), class = "col_spec"))

df2 <- structure(list(trial = c(1, 2, 3, 4, 11, 12, 30), control = c(0, 
0, 0, 0, 3, 0, 3), t_start = c(16.64709, 41.81843, 65.5451, 82.65743, 
187.0787, 200.0486, 415.171), t_end = c(35.49431, 57.74304, 71.16612, 
87.30914, 193.5898, 203.1883, 418.0405)), class = c("spec_tbl_df", 
"tbl_df", "tbl", "data.frame"), row.names = c(NA, -7L), spec = structure(list(
    cols = list(trial = structure(list(), class = c("collector_double", 
    "collector")), control = structure(list(), class = c("collector_double", 
    "collector")), t_start = structure(list(), class = c("collector_double", 
    "collector")), t_end = structure(list(), class = c("collector_double", 
    "collector"))), default = structure(list(), class = c("collector_guess", 
    "collector")), skip = 1L), class = "col_spec"))