2017-10-13 54 views
1

我在R有4个不同日期的个人网站的状态数据。我希望为每个不在每个网站上看到的每个单独ID创建行,并附加一个存在/不存在列。根据在线状态数据添加缺席行

这是我存在的数据子集:

ID SP SITE DATE SMP 
975 1 1 41579 FJB 
997 1 1 41579 FAW 
114 1 1 41579 FJW 
926 2 1 41609 FJB 
101 2 1 41609 FJB 
108 2 1 41609 FAB 
129 3 2 41710 FUB 
131 3 2 41710 MAW 
132 3 2 41710 FAW 

这是我想创造什么:

ID SP SITE DATE SMP Present? 
975 1 1 41579 FJB Yes 
997 1 1 41579 FAW Yes 
114 1 1 41579 FJW Yes 
926 1 1 41579 FJB No 
101 1 1 41579 FJB No 
108 1 1 41579 FAB No 
129 1 1 41579 FUB No 
131 1 1 41579 MAW No 
132 1 1 41579 FAW No 
975 2 1 41609 FJB No 
997 2 1 41609 FAW No 
114 2 1 41609 FJW No 
926 2 1 41609 FJB Yes 
101 2 1 41609 FJB Yes 
108 2 1 41609 FAB Yes 
129 2 1 41609 FUB No 
131 2 1 41609 MAW No 
132 2 1 41609 FAW No 
975 3 2 41710 FJB No 
997 3 2 41710 FAW No 
114 3 2 41710 FJW No 
926 3 2 41710 FJB No 
101 3 2 41710 FJB No 
108 3 2 41710 FAB No 
129 3 2 41710 FUB Yes 
131 3 2 41710 MAW Yes 
132 3 2 41710 FAW Yes 

我希望有人也许能够帮助!

回答

0

首先,您需要准备一个包含ID,SP,SITE,DATE和SMP所有可能组合的数据集,我将其称为attendance_list

library(tidyverse) 
txt <- "ID SP SITE DATE SMP 
975 1 1 41579 FJB 
997 1 1 41579 FAW 
114 1 1 41579 FJW 
926 2 1 41609 FJB 
101 2 1 41609 FJB 
108 2 1 41609 FAB 
129 3 2 41710 FUB 
131 3 2 41710 MAW 
132 3 2 41710 FAW" 
presence_data <- txt %>% 
    gsub(" +", " ", x=.) %>% 
    read.delim(text=., sep=" ") 

# Attendance list 
## It looks like each ID has only one SMP 
id_smp <- presence_data %>% 
    select(ID, SMP) %>% 
    distinct() 
## Working site and date 
sp_site_date <- presence_data %>% 
    select(SP, SITE, DATE) %>% 
    distinct() 
attendance_list <- merge(id_smp, sp_site_date) 

然后,您可以留下一个“是”栏加入presence_dataattendance_list,并更换为“NA”(数据在attendance_list但不是在presence data)由“否”。

attendance_list <- attendance_list %>% 
    left_join(presence_data %>% mutate(PRESENT="Yes"), by=c("ID", "SMP", "SP", "SITE", "DATE")) %>% 
    mutate(PRESENT = ifelse(is.na(PRESENT), "No", PRESENT)) 

输出:

> attendance_list 
#  ID SMP SP SITE DATE PRESENT 
# 1 975 FJB 1 1 41579  Yes 
# 2 997 FAW 1 1 41579  Yes 
# 3 114 FJW 1 1 41579  Yes 
# 4 926 FJB 1 1 41579  No 
# 5 101 FJB 1 1 41579  No 
# 6 108 FAB 1 1 41579  No 
# 7 129 FUB 1 1 41579  No 
# 8 131 MAW 1 1 41579  No 
# 9 132 FAW 1 1 41579  No 
# 10 975 FJB 2 1 41609  No 
# 11 997 FAW 2 1 41609  No 
# 12 114 FJW 2 1 41609  No 
# 13 926 FJB 2 1 41609  Yes 
# 14 101 FJB 2 1 41609  Yes 
# 15 108 FAB 2 1 41609  Yes 
# 16 129 FUB 2 1 41609  No 
# 17 131 MAW 2 1 41609  No 
# 18 132 FAW 2 1 41609  No 
# 19 975 FJB 3 2 41710  No 
# 20 997 FAW 3 2 41710  No 
# 21 114 FJW 3 2 41710  No 
# 22 926 FJB 3 2 41710  No 
# 23 101 FJB 3 2 41710  No 
# 24 108 FAB 3 2 41710  No 
# 25 129 FUB 3 2 41710  Yes 
# 26 131 MAW 3 2 41710  Yes 
# 27 132 FAW 3 2 41710  Yes 
+0

哇! 这太不可思议了 - 非常感谢你! – RobP