Mercurial > repos > ecology > vigiechiro_bilanenrichirp
annotate BilanEnrichiRP.R @ 0:c55e09a8b4c8 draft
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
author | ecology |
---|---|
date | Wed, 13 Mar 2019 11:17:44 -0400 |
parents | |
children | be4e28da3919 |
rev | line source |
---|---|
0
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
1 suppressMessages(library(data.table)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
2 suppressMessages(library(DT)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
3 suppressMessages(library(htmlwidgets)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
4 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
5 args <- commandArgs(trailingOnly = TRUE) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
6 EchelleErreur=c("NA","POSSIBLE","PROBABLE","SUR") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
7 EchelleNumErreur=c(99,50,10,1) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
8 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
9 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
10 IdC2=fread(args[1]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
11 refRP=fread(args[2]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
12 GroupList=fread(args[3]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
13 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
14 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
15 if(substr(IdC2$`nom du fichier`[1],2,2)!="i") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
16 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
17 stop("Protocole non conforme, ce script doit etre lance pour un protocole Routier ou Pedestre",call.=FALSE) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
18 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
19 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
20 Routier=grepl("-",substr(IdC2$`nom du fichier`[1],4,7)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
21 #compute error risk by species (minimum error among files) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
22 #to be replaced by glm outputs if I'll have time |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
23 RisqueErreurT=aggregate(IdC2$IdProb,by=list(IdC2$IdExtrap),FUN=function(x) ceiling((1-max(x-0.0001))*100)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
24 #barplot(RisqueErreurT$x,names.arg=RisqueErreurT$Group.1,las=2) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
25 #compute error risk accoring to observer/validator (a little dirty because it relies on alphabetical order of confidence classes: POSSIBLE < PROBABLE < SUR) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
26 RisqueErreurOV0=match(IdC2$ConfV,EchelleErreur) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
27 RisqueErreurOV=aggregate(RisqueErreurOV0,by=list(IdC2$IdExtrap) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
28 ,FUN=max) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
29 RisqueErreurOV2=EchelleNumErreur[RisqueErreurOV$x] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
30 #compute minimum error risk between man and machine |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
31 RisqueErreur=pmin(RisqueErreurT$x,RisqueErreurOV2,na.rm=TRUE) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
32 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
33 #compute number of files validated per species |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
34 FichValid=aggregate(IdC2$IdV,by=list(IdC2$IdExtrap,IdC2$'nom du fichier') |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
35 ,FUN=function(x) sum(!is.na(x))) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
36 NbValid2=aggregate(FichValid$x,by=list(FichValid$Group.1),FUN=function(x) sum(x>0)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
37 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
38 DiffC50=vector() # to store the median of confidence difference between unvalidated records and validated ones |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
39 DiffT50=vector() # to store the median of time difference between unvalidated records and validated ones |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
40 for (j in 1:nlevels(as.factor(IdC2$IdExtrap))) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
41 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
42 IdSp=subset(IdC2 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
43 ,IdC2$IdExtrap==levels(as.factor(IdC2$IdExtrap))[j]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
44 IdSp$IdProb[is.na(IdSp$IdProb)]=0 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
45 IdSp=IdSp[order(IdSp$IdProb),] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
46 IdSpV=subset(IdSp,!is.na(IdSp$IdV)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
47 if(nrow(IdSpV)>0) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
48 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
49 cuts <- c(-Inf, IdSpV$IdProb[-1]-diff(IdSpV$IdProb)/2, Inf) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
50 CorrC=findInterval(IdSp$IdProb, cuts) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
51 CorrC2=IdSpV$IdProb[CorrC] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
52 DiffC=abs(IdSp$IdProb-CorrC2) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
53 DiffC50=c(DiffC50,median(DiffC)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
54 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
55 IdSp=IdSp[order(IdSp$TimeNum),] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
56 IdSpV=subset(IdSp,!is.na(IdSp$IdV)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
57 cuts <- c(-Inf, IdSpV$TimeNum[-1]-diff(IdSpV$TimeNum)/2, Inf) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
58 CorrT=findInterval(IdSp$TimeNum, cuts) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
59 CorrT2=IdSpV$TimeNum[CorrT] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
60 DiffT=abs(IdSp$TimeNum-CorrT2) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
61 DiffT50=c(DiffT50,median(DiffT)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
62 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
63 DiffC50=c(DiffC50,Inf) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
64 DiffT50=c(DiffT50,Inf) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
65 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
66 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
67 #compute an index of validation effort per species |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
68 EffortV=1/DiffC50/DiffT50 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
69 EffortClass=(EffortV>0.0005)+(EffortV>0.005)+RisqueErreurOV$x |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
70 #cbind(RisqueErreurOV,EffortV,DiffC50,DiffT50) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
71 #barplot(EffortClass-1,names.arg=NbValid2$Group.1,las=2) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
72 ClassEffortV=c("-","FAIBLE","SUFFISANT","SUFFISANT","FORT","FORT") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
73 EffortClassMot=ClassEffortV[EffortClass] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
74 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
75 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
76 #compare activity / reference frame |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
77 FileInfo=as.data.table(tstrsplit(IdC2$`nom du fichier`,"-")) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
78 IdC2$Tron=FileInfo$V4 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
79 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
80 MicTempsInfo=as.data.table(tstrsplit(as.data.frame(FileInfo)[,(ncol(FileInfo))],"_")) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
81 MicDroit=(as.data.frame(MicTempsInfo)[,(ncol(MicTempsInfo)-2)]=="1") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
82 IdC2$MicDroit=MicDroit |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
83 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
84 testTempsFin=aggregate(IdC2$temps_fin,by=list(MicDroit),FUN=max) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
85 testTempsFin$Direct=(testTempsFin$x>0.5) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
86 testTF2=sum((testTempsFin$x>0.5)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
87 if(testTF2>1){stop("Probleme stereo : les 2 canaux semblent etre en enregistrement direct")} |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
88 IdC2M=merge(IdC2,testTempsFin,by.x="MicDroit",by.y="Group.1") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
89 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
90 ActMoy=aggregate(IdC2$`nom du fichier` |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
91 ,by=list(IdC2M$IdExtrap,IdC2M$Direct),FUN=length) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
92 ListSpref=match(levels(as.factor(ActMoy$Group.1)),refRP$Espece) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
93 Subref=refRP[ListSpref] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
94 if(Routier) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
95 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
96 Subref=Subref[,c(1:17)] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
97 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
98 Subref=Subref[,c(1,18:33)] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
99 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
100 QualifActE=vector() |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
101 QualifActD=vector() |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
102 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
103 for (k in 1:nlevels(as.factor(ActMoy$Group.1))) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
104 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
105 Actsub=subset(ActMoy,ActMoy$Group.1==levels(as.factor(ActMoy$Group.1))[k]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
106 if(is.na(Subref[k,2])) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
107 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
108 QualifActE=c(QualifActE,NA) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
109 QualifActD=c(QualifActD,NA) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
110 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
111 ActE=subset(Actsub,Actsub$Group.2==F) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
112 if(nrow(ActE)==0) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
113 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
114 QualifActE=c(QualifActE,NA) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
115 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
116 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
117 cuts=cbind(-Inf,as.numeric(Subref[k,6]),as.numeric(Subref[k,7]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
118 ,as.numeric(Subref[k,8]),Inf) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
119 QualifActE=c(QualifActE,findInterval(ActE$x,cuts,left.open=T)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
120 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
121 ActD=subset(Actsub,Actsub$Group.2==T) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
122 if(nrow(ActD)==0) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
123 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
124 QualifActD=c(QualifActD,NA) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
125 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
126 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
127 cuts=cbind(-Inf,as.numeric(Subref[k,14]),as.numeric(Subref[k,15]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
128 ,as.numeric(Subref[k,16]),Inf) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
129 QualifActD=c(QualifActD,findInterval(ActD$x,cuts,left.open=T)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
130 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
131 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
132 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
133 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
134 ClassAct=c("FAIBLE","MODEREE","FORTE","TRES FORTE") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
135 QualifActMotE=ClassAct[QualifActE] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
136 QualifActMotD=ClassAct[QualifActD] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
137 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
138 #compute activity by nights (to be completed) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
139 #ActNuit=aggregate(IdC2M$`nom du fichier`,by=list(IdC2M$DateNuit,IdC2M$IdExtrap),FUN=length) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
140 ActED=dcast(data=ActMoy,formula=Group.1~Group.2,value=x) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
141 ActED[is.na(ActED)]=0 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
142 #organize the csv summary |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
143 SummPart0=cbind(Esp=levels(as.factor(IdC2M$IdExtrap)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
144 ,RisqueErreur,NbValid=NbValid2$x,EffortValid=EffortClassMot) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
145 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
146 test=match("FALSE",colnames(ActED)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
147 if(is.na(test)==F) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
148 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
149 SummPart0=cbind(SummPart0,Contacts_Expansion=ActED$'FALSE' |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
150 ,Niveau_Activite_Expansion=QualifActMotE) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
151 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
152 SummPart0=cbind(SummPart0,Contacts_Expansion="" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
153 ,Niveau_Activite_Expansion="") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
154 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
155 test=match("TRUE",colnames(ActED)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
156 if(is.na(test)==F) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
157 { |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
158 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
159 SummPart0=cbind(SummPart0,Contacts_Direct=ActED$'TRUE' |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
160 ,Niveau_Activite_Direct=QualifActMotD) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
161 }else{ |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
162 SummPart0=cbind(SummPart0,Contacts_Direct="" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
163 ,Niveau_Activite_Direct="") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
164 } |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
165 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
166 InfoSp=c("GroupFR","NomFR","Scientific name","Esp") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
167 GroupShort=GroupList[,..InfoSp] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
168 SummPart=merge(GroupShort,SummPart0,by="Esp") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
169 IndexGroupe=c("Autre","Sauterelle","Chauve-souris") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
170 SummPart$IndexSumm=match(SummPart$GroupFR,IndexGroupe) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
171 SummPart=SummPart[with(SummPart |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
172 ,order(IndexSumm,as.numeric(Contacts_Direct),as.numeric(Contacts_Expansion),decreasing=T)),] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
173 colnames(SummPart)=c("Code","Groupe","Nom francais","Nom scientifique" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
174 ,"Risque d'erreur (%)","Nb Validations" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
175 ,"Effort de validation","Nb de Contacts en expansion" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
176 ,"Niveau d'Activite en expansion" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
177 ,"Nb de Contacts en direct" |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
178 ,"Niveau d'Activite en direct","TriGroupe") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
179 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
180 #to do: extend colors to other columns to improve readability |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
181 SummHTML=datatable(SummPart, rownames = FALSE) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
182 formatStyle(columns = c("Code","Groupe","Nom francais","Nom scientifique","Risque d'erreur (%)"),valueColumns="Risque d'erreur (%)", |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
183 background = styleInterval(c(1, 10, 50), c("white", "khaki", "orange", "orangered"))) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
184 formatStyle(columns = "Effort de validation", |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
185 background = styleEqual(c("-","FAIBLE","SUFFISANT","FORT"), c("white", "cyan", "royalblue", "darkblue"))) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
186 formatStyle(columns = c("Nb de Contacts en expansion","Niveau d'Activite en expansion"),valueColumns="Niveau d'Activite en expansion", |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
187 background = styleEqual(c("FAIBLE","MODEREE","FORTE","TRES FORTE"), c("palegoldenrod", "greenyellow", "limegreen", "darkgreen"))) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
188 formatStyle(columns = c("Nb de Contacts en direct","Niveau d'Activite en direct"),valueColumns="Niveau d'Activite en direct", |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
189 background = styleEqual(c("FAIBLE","MODEREE","FORTE","TRES FORTE"), c("palegoldenrod", "greenyellow", "limegreen", "darkgreen"))) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
190 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
191 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
192 saveWidget(SummHTML,"output-summaryRP.html") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
193 write.table(SummPart,"output-summaryRP.tabular",row.names=F,sep="\t",quote=FALSE) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
194 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
195 #summary for each point/transect |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
196 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
197 #compute number of files validated per species |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
198 IdC2M$Canal=sapply(IdC2M$Direct,FUN=function(x) if(x){"Direct"}else{"Expansion"}) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
199 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
200 ActMoyTA=aggregate(IdC2M$`nom du fichier`,by=list(IdC2M$IdExtrap,IdC2M$Canal,IdC2M$Session),FUN=length) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
201 ActMoyT=dcast(data=IdC2M,formula=IdExtrap+Canal~Session |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
202 ,fun.aggregate=length) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
203 SummPartshort=cbind(SummPart[,c(1:5)],TriGroupe=SummPart[,TriGroupe]) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
204 SummPartTron=merge(SummPartshort,ActMoyT,by.x="Code",by.y="IdExtrap") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
205 SummPartTron=SummPartTron[order(TriGroupe,decreasing=T),] |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
206 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
207 ListSession=levels(as.factor(IdC2M$Session)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
208 brks <- quantile(ActMoyTA$x, probs = seq(.05, .95, .05), na.rm = TRUE)-1 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
209 clrs <- round(seq(255, 40, length.out = length(brks) + 1), 0) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
210 {paste0("rgb(255,", ., ",", ., ")")} |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
211 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
212 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
213 #to do: extend colors to other columns to improve readability |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
214 SummHTMLT=datatable(SummPartTron, rownames = FALSE) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
215 formatStyle(columns = c("Code","Groupe","Nom francais","Nom scientifique","Risque d'erreur (%)"),valueColumns="Risque d'erreur (%)", |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
216 background = styleInterval(c(1, 10, 50), c("white", "khaki", "orange", "orangered"))) %>% |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
217 formatStyle(columns=ListSession, backgroundColor = styleInterval(brks, clrs)) |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
218 |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
219 saveWidget(SummHTMLT,"output-detailRP.html") |
c55e09a8b4c8
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tools/vigiechiro commit d2de8e10c11bfa3b04729e59bba58e08d23b56aa
ecology
parents:
diff
changeset
|
220 write.table(SummPartTron,"output-detailRP.tabular",row.names=F,sep="\t",quote=FALSE) |