-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathDOG.r
61 lines (41 loc) · 2.11 KB
/
DOG.r
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
50
51
52
53
54
55
56
57
58
59
60
61
sp.label="Dog"
sp.acron="DOG"
sp.position=which(species.acronyms==sp.acron)
# SVA data: ----
# refining syndromic classification for species specific syndromes ----
CD.species <-current.SVA.data[current.SVA.data$SPECIES==sp.label,]
am <- CD.species[CD.species$General=="AntimResistance",]
if(dim(am)[1]>0)(am$SYNDROMIC <- "AntimResistance")
dd <- CD.species[CD.species$Doddes==1,]
if(dim(dd)[1]>0)(dd$SYNDROMIC <- "Doddes")
CD.species <- rbind(CD.species,am,dd)
CD.species$SYNDROMIC[CD.species$Specific=="Sarcoptes"] <- "Sarcoptes"
CD.species$SYNDROMIC[CD.species$SYNDROMIC=="GIT"&
CD.species$General=="Macroparasites"] <- "Macropar"
CD.species$SYNDROMIC[CD.species$SYNDROMIC=="GIT"&
CD.species$General=="Microparasites"] <- "Micropar"
CD.species$SYNDROMIC[CD.species$SYNDROMIC=="GIT"&
CD.species$General=="Bacterial"] <- "Bacterial"
CD.species$SYNDROMIC[CD.species$SYNDROMIC=="GIT"&
CD.species$General=="Viral"] <- "Viral"
classified.species.data <- bind_rows(classified.species.data,CD.species)
# running vetsyn-based-engine ----
#(running the codes that are the same for every syndromic object, namely:
# 1. update ths yndromic object;
# 2. run detection algorithms;
# 3. report alarms)
source("vetsyn-based-engine.r",local=TRUE,encoding="native.enc")
# svaga data ----
source("svaga-data-filtering.r",local=TRUE,encoding="native.enc")
# save data ----
assign(paste0(sp.acron,".svaga"),svaga.object[[sp.position]])
#assign(paste0(sp.acron,".non.svaga"),non.svaga.object[[sp.position]])
assign(paste0(sp.acron,".weekly"),weekly.object[[sp.position]])
assign(paste0(sp.acron,".daily"),daily.object[[sp.position]])
eval(parse(text=paste0("save(",
paste0(sp.acron,'.daily,'),
paste0(sp.acron,'.weekly,'),
paste0(sp.acron,'.svaga,'),
#paste0(sp.acron,'.non.svaga,'),
"file='",
paste0(wd.history,sp.acron,".RData'"),")")))