1

I have some data on concentrations observed for chemicals by tasks performed. There are three people, each performing two tasks, and repeating each task twice. concentrations are measured simultaneously for the three different chemicals, at five time points. Person A has only completed two repetitions, and some other concentrations were missing. The data looks like this:

    test_dt <- data.table(person = rep(LETTERS[1:3],each = 20),
           task = rep(LETTERS[24:25], each = 10),
           reps = rep(c(1,2),each = 5),
           time = 1:5, chem1 = rnorm(60,1,0.2),
           chem2 = rnorm(60,4,1.1),chem3 = rnorm(60,2,0.75))
   test_dt[person == "A" & reps == 2,`:=`(chem3 = NA_real_)]
   test_dt[person == "B" & task == "X" & reps == 1 &time %in% 3:5,chem1 := NA_real_]
   test_dt[person == "C" & task == "Y" & reps == 2 &time %in% 3:4,chem2 := NA_real_]
   

I want to get the time when the data is first NA and the time NAs end, for each person, by task and repetition. I tried doing this:

lapply(c("chem1","chem2","chem3"),function(var){
  start_var = paste0("na_start_",var)
  end_var = paste0("na_end_",var)
  test_dt[is.na(get(var)), 
          .(deparse(substitute(start_var)) = min(time),
            deparse(substitute(end_var)) = max(time)),
          .(person,task,reps)]
})

But ended up with this error:

"  test_dt[is.na(get(var)), 
          .(deparse(substitute(start_var)) ="
>             deparse(substitute(end_var)) = max(time)),
Error: unexpected ')' in "            deparse(substitute(end_var)) = max(time))"
>           .(person,task,reps)]
Error: unexpected ']' in "          .(person,task,reps)]"
> })
Error: unexpected '}' in "}"

How can I do this in data.table in R?

1 Answer 1

2

Instead of deparse/substitute, use setNames or setnames. Note that = wouldn't allow evaluation on the lhs

lapply(c("chem1","chem2","chem3"),function(var){
  start_var = paste0("na_start_",var)
  end_var = paste0("na_end_",var)
 test_dt[is.na(get(var)), setNames(.(min(time), max(time)),
     c(start_var, end_var)), .(person, task, reps)]
})

-output

[[1]]
   person task reps na_start_chem1 na_end_chem1
1:      B    X    1              3            5

[[2]]
   person task reps na_start_chem2 na_end_chem2
1:      C    Y    2              3            4

[[3]]
   person task reps na_start_chem3 na_end_chem3
1:      A    X    2              1            5
2:      A    Y    2              1            5
Sign up to request clarification or add additional context in comments.

Comments

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Start asking to get answers

Find the answer to your question by asking.

Ask question

Explore related questions

See similar questions with these tags.