1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49
|
BR139 <- str_c("BR", 1:139, sep ="")
AE139 <- str_c("AE", 1:139, sep ="")
df2 <- data.frame(CC = rep(c("CC1","CC2", "CC3"), each = 139*139*2),
SI = rep(c("SI1","SI2"), each = 139*139),
BR = rep(BR139 , 139),
AE = rep(AE139 ,each = 139),
TA = rep(c("TA1","TA2"), each = 139*139*2*3),
VAL = sample(c(NA, 1:10), 139*139*2*3, replace = T, prob = c(1/3, rep(1/15,10)))
)
df1 <- filter(df2,
CC == "CC1",
SI %in% c("SI1", "SI2"),
BR %in% BR139,
AE %in% AE139,
TA == "TA1")
df2 <- drop_na(df2, VAL)
df1$VAL2 <- 0
df1$VAL2 <- ifelse (str_sub(df1$AE, start = 3, end = 5) == str_sub(df1$BR, start = 3, end = 5),
rowSums(cbind(df2[df2$AE %in% df1$AE &
df2$BR %in% df1$BR &
df2$SI %in% df1$SI&
df2$CC == "CC1"&
df2$TA == "TA1", "VAL"],
df2[df2$AE %in% df1$AE &
df2$BR == "BR87" &
df2$SI %in% df1$SI&
df2$CC == "CC1"&
df2$TA == "TA1", "VAL"],
df2[df2$AE %in% df1$AE &
df2$BR == "BR88" &
df2$SI %in% df1$SI&
df2$CC == "CC1"&
df2$TA == "TA1", "VAL"],
df2[df2$AE %in% df1$AE &
df2$BR == "BR94" &
df2$SI %in% df1$SI&
df2$CC == "CC1"&
df2$TA == "TA1", "VAL"], na.rm = T)),
df2[df2$AE %in% df1$AE &
df2$BR %in% df1$BR &
df2$SI %in% df1$SI&
df2$CC == "CC1"&
df2$TA == "TA1", "VAL"]) |
Partager