-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrechercheSimple.R
113 lines (100 loc) · 5.25 KB
/
rechercheSimple.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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
library(readr)
chiffres_cles_2018 <- read_delim("chiffres-cles-2018.csv",
";", escape_double = FALSE, trim_ws = TRUE)
View(chiffres_cles_2018)
chiffres_cles_2015 <- read_delim("chiffres-cles-2015.csv",
";", escape_double = FALSE, trim_ws = TRUE)
View(chiffres_cles_2015)
chiffre2018<-chiffres_cles_2018[,c(2,20,22:25)]
chiffre2017<-chiffres_cles_2018[,c(2,26,28:31)]
chiffre2016<-chiffres_cles_2018[,c(2,32,34:37)]
chiffre2015<-chiffres_cles_2015[,c(2,10,12:15)]
chiffre2014<-chiffres_cles_2015[,c(2,16,18:21)]
chiffre2013<-chiffres_cles_2015[,c(2,22,24:27)]
colnames(chiffre2018)<- c("Siren","Année","Durée","CA","Résultat","Effectif")
colnames(chiffre2017)<- c("Siren","Année","Durée","CA","Résultat","Effectif")
colnames(chiffre2016)<- c("Siren","Année","Durée","CA","Résultat","Effectif")
colnames(chiffre2015)<- c("Siren","Année","Durée","CA","Résultat","Effectif")
colnames(chiffre2014)<- c("Siren","Année","Durée","CA","Résultat","Effectif")
colnames(chiffre2013)<- c("Siren","Année","Durée","CA","Résultat","Effectif")
chiffre<-rbind(chiffre2013,chiffre2014,chiffre2015,chiffre2016,chiffre2017,chiffre2018)
chiffreClean <- chiffre[complete.cases(chiffre), ]
nrow(chiffreClean)/6
chiffre2018Clean<-chiffre2018[complete.cases(chiffre2018), ]
chiffre2017Clean<-chiffre2017[complete.cases(chiffre2017), ]
chiffre2016Clean<-chiffre2016[complete.cases(chiffre2016), ]
chiffre2015Clean<-chiffre2015[complete.cases(chiffre2015), ]
chiffre2014Clean<-chiffre2014[complete.cases(chiffre2014), ]
chiffre2013Clean<-chiffre2013[complete.cases(chiffre2013), ]
chiffretraj=data.frame();
#Prepare for traj
for(i in 1:nrow(chiffre2013Clean)){
siren<-chiffre2013Clean[i,]$Siren
siren2013<-chiffre2013Clean[i,]
siren2014<-unique(chiffre2014Clean[chiffre2014Clean$Siren==siren,])
siren2015<-unique(chiffre2015Clean[chiffre2015Clean$Siren==siren,])
siren2016<-unique(chiffre2016Clean[chiffre2016Clean$Siren==siren,])
siren2017<-unique(chiffre2017Clean[chiffre2017Clean$Siren==siren,])
siren2018<-unique(chiffre2018Clean[chiffre2018Clean$Siren==siren,])
if(nrow(siren2014)==1 && nrow(siren2015)==1 && nrow(siren2016)==1 && nrow(siren2017)==1 && nrow(siren2018)==1){
racc<-c(siren,
siren2013$CA/siren2013$Durée*12,
siren2014$CA/siren2014$Durée*12,
siren2015$CA/siren2015$Durée*12,
siren2016$CA/siren2016$Durée*12,
siren2017$CA/siren2017$Durée*12,
siren2018$CA/siren2018$Durée*12)
raccR<-c(
siren2013$Résultat/siren2013$Durée*12,
siren2014$Résultat/siren2014$Durée*12,
siren2015$Résultat/siren2015$Durée*12,
siren2016$Résultat/siren2016$Durée*12,
siren2017$Résultat/siren2017$Durée*12,
siren2018$Résultat/siren2018$Durée*12)
raccE<-c(
siren2013$Effectif/siren2013$Durée*12,
siren2014$Effectif/siren2014$Durée*12,
siren2015$Effectif/siren2015$Durée*12,
siren2016$Effectif/siren2016$Durée*12,
siren2017$Effectif/siren2017$Durée*12,
siren2018$Effectif/siren2018$Durée*12)
racc<-cbind(t(as.matrix(racc)),t(as.matrix(raccR)),t(as.matrix(raccE)))
chiffretraj<-rbind(chiffretraj,racc)
}
}
colnames(chiffretraj)<- c("Siren","CA2013","CA2014","CA2015","CA2016","CA2017","CA2018","Résultat2013","Résultat2014","Résultat2015","Résultat2016","Résultat2017","Résultat2018","Effectif2013","Effectif2014","Effectif2015","Effectif2016","Effectif2017","Effectif2018")
write_csv(chiffretraj,"chiffeAffaire.csv")
chiffretrajEnd<-chiffretraj
library(kml)
library(scales)
mtx <- matrix(t(1:6), ncol = 6)
mtx.rep <- (apply(t(mtx), 1, rep,nrow(chiffretraj) ))
time<-cbind(chiffretraj$Siren,mtx.rep)
data<-as.data.frame(cbind(chiffretraj$Siren,t(apply(chiffretraj[,2:7], 1, rescale))))
(ld2 <- clusterLongData(traj=data[,2:7],idAll=paste("I-",1:nrow(data),sep=""),time=(2:7)+0.5,varNames="R"))
kml(ld2, nbRedraw = 2, toPlot = "both")
chiffretrajEnd$clustersCA <- getClusters(ld2, 4)
data<-as.data.frame(cbind(chiffretraj$Siren,t(apply(chiffretraj[,8:13], 1, rescale))))
(ld2 <- clusterLongData(traj=data[,2:7],idAll=paste("I-",1:nrow(data),sep=""),time=(2:7)+0.5,varNames="R"))
kml(ld2, nbRedraw = 2, toPlot = "both")
chiffretrajEnd$clustersRES <- getClusters(ld2, 4)
data<-as.data.frame(cbind(chiffretraj$Siren,t(apply(chiffretraj[,14:19], 1, rescale))))
(ld2 <- clusterLongData(traj=data[,2:7],idAll=paste("I-",1:nrow(data),sep=""),time=(2:7)+0.5,varNames="R"))
kml(ld2, nbRedraw = 2, toPlot = "both")
chiffretrajEnd$clustersEFF <- getClusters(ld2, 4)
chiffretrajEnd$activite=""
chiffretrajEnd$nom=""
for(i in 1:nrow(chiffretrajEnd)){
siren<-chiffretrajEnd[i,1]
listAct<-chiffres_cles_2018[chiffres_cles_2018$Siren==siren,6]
caseAct<-unique(listAct[complete.cases(listAct),])
if(nrow(caseAct)==1){
chiffretrajEnd[i,23]=unique(listAct[complete.cases(listAct),])
}
listNom<-chiffres_cles_2018[chiffres_cles_2018$Siren==siren,1]
caseNom<-unique(listNom[complete.cases(listNom),])
if(nrow(caseNom)==1){
chiffretrajEnd[i,24]=unique(listNom[complete.cases(listNom),])
}
}
write_csv(chiffretrajEnd,"chiffeAffaireCluster.csv")