R/Psetspecific.R

In mhajihos/GPAQ: Global Physical Activity Questionnaire

 #psetspecific (unweighted mean & median values)

Psetspecific<-function(data)
{
data%>%

  mutate(p1t3 = ifelse(p1t3cln==1, p2*p3, NA)) %>% 
  mutate(p4t6 = ifelse(p4t6cln==1, p5*p6, NA)) %>% 
  mutate(p7t9 = ifelse(p7t9cln==1, p8*p9, NA)) %>% 
  mutate(p10t12 = ifelse(p10t12cln==1, p11*p12, NA)) %>% 
  mutate(p13t15 = ifelse(p13t15cln==1, p14*p15, NA)) %>% 

  mutate(pworkday = (p1t3+p4t6)/7) %>% 
  mutate(ptravelday = p7t9/7) %>% 
  mutate(precday = (p10t12+p13t15)/7)%>% 
  mutate(cln_specific = if_else(p1t3cln==1 & p4t6cln==1 & p7t9cln==1 & p10t12cln==1 & p13t15cln==1, 1, 2, missing = 2)) %>% 
  mutate(cln_specific = replace(cln_specific, is.na(p1) & is.na(p4) & is.na(p7) & is.na(p10) & is.na(p13), 2))%>%
  
  return(data)

}