forked from SVA-SE/SVASSS
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathVLT.r
62 lines (39 loc) · 1.92 KB
/
VLT.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
62
sp.label="VLT"
sp.acron="VLT"
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")
total <- CD.species
if(dim(total)[1]>0)(total$SYNDROMIC <- "Total")
CD.species <- rbind(CD.species,am,dd,total)
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"
# 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]])
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'"),")")))