Description
df <- cbind(cs, dip)
head(df)
## region dep libdep wave ids idt pop
## 1 82 01 Ain 1968 82_01 82_01_1968 89180
## 2 22 02 Aisne 1968 22_02 22_02_1968 115968
## 3 83 03 Allier 1968 83_03 83_03_1968 90240
## 4 93 04 Alpes-de-Haute-Provence 1968 93_04 93_04_1968 24600
## 5 93 05 Hautes-Alpes 1968 93_05 93_05_1968 22808
## 6 93 06 Alpes-Maritimes 1968 93_06 93_06_1968 178236
## p_cs1 p_cs2 p_cs3 p_cs4 p_cs5 p_cs6 region dep
## 1 19.533530 11.751510 3.758690 11.28056 16.48800 37.18771 82 01
## 2 9.078366 9.716474 3.959713 13.10706 17.44274 46.69564 22 02
## 3 15.806740 11.728720 4.312943 12.07890 18.86082 37.21188 83 03
## 4 13.528460 14.000000 5.756098 13.86992 17.07317 35.77236 93 04
## 5 19.729920 13.170820 5.559453 12.92529 19.69484 28.91968 93 05
## 6 3.456092 16.214460 7.199443 12.74041 26.68148 33.70812 93 06
## libdep wave ids idt pop p_dip0 p_dip1
## 1 Ain 1968 82_01 82_01_1968 89180 36.33102 37.26844
## 2 Aisne 1968 22_02 22_02_1968 115968 40.22834 33.73344
## 3 Allier 1968 83_03 83_03_1968 90240 33.03635 37.79255
## 4 Alpes-de-Haute-Provence 1968 93_04 93_04_1968 24600 36.86179 32.16260
## 5 Hautes-Alpes 1968 93_05 93_05_1968 22808 32.88320 37.82883
## 6 Alpes-Maritimes 1968 93_06 93_06_1968 178236 35.44738 29.71566
## p_dip2 p_dip3 p_dip4 p_dip7
## 1 3.740749 12.93115 6.310832 3.417807
## 2 3.338852 12.43791 6.857064 3.404387
## 3 4.485816 12.55762 8.280142 3.847518
## 4 5.739837 11.02439 7.804878 6.406504
## 5 5.243774 11.60996 7.821817 4.612417
## 6 6.712449 12.28932 9.461613 6.373572
p <- ggplot(df, aes(p_dip7, p_cs4))
p + geom_path(aes( colour=factor(dep)))
p + geom_path(aes( colour=wave) ) + facet_wrap(~dep)
p <- ggplot(df, aes(p_cs3, p_dip7))
p + geom_path(aes( colour=wave) ) + facet_wrap(~dep)
dx <- cbind(df, ptamult$Tli)
head(dx)
## region dep libdep wave ids idt pop
## 1.Ana1 82 01 Ain 1968 82_01 82_01_1968 89180
## 2.Ana1 22 02 Aisne 1968 22_02 22_02_1968 115968
## 3.Ana1 83 03 Allier 1968 83_03 83_03_1968 90240
## 4.Ana1 93 04 Alpes-de-Haute-Provence 1968 93_04 93_04_1968 24600
## 5.Ana1 93 05 Hautes-Alpes 1968 93_05 93_05_1968 22808
## 6.Ana1 93 06 Alpes-Maritimes 1968 93_06 93_06_1968 178236
## p_cs1 p_cs2 p_cs3 p_cs4 p_cs5 p_cs6 region dep
## 1.Ana1 19.533530 11.751510 3.758690 11.28056 16.48800 37.18771 82 01
## 2.Ana1 9.078366 9.716474 3.959713 13.10706 17.44274 46.69564 22 02
## 3.Ana1 15.806740 11.728720 4.312943 12.07890 18.86082 37.21188 83 03
## 4.Ana1 13.528460 14.000000 5.756098 13.86992 17.07317 35.77236 93 04
## 5.Ana1 19.729920 13.170820 5.559453 12.92529 19.69484 28.91968 93 05
## 6.Ana1 3.456092 16.214460 7.199443 12.74041 26.68148 33.70812 93 06
## libdep wave ids idt pop p_dip0
## 1.Ana1 Ain 1968 82_01 82_01_1968 89180 36.33102
## 2.Ana1 Aisne 1968 22_02 22_02_1968 115968 40.22834
## 3.Ana1 Allier 1968 83_03 83_03_1968 90240 33.03635
## 4.Ana1 Alpes-de-Haute-Provence 1968 93_04 93_04_1968 24600 36.86179
## 5.Ana1 Hautes-Alpes 1968 93_05 93_05_1968 22808 32.88320
## 6.Ana1 Alpes-Maritimes 1968 93_06 93_06_1968 178236 35.44738
## p_dip1 p_dip2 p_dip3 p_dip4 p_dip7 CS1 CS2
## 1.Ana1 37.26844 3.740749 12.93115 6.310832 3.417807 -0.4915543 -1.9701348
## 2.Ana1 33.73344 3.338852 12.43791 6.857064 3.404387 -1.1271310 -2.2096100
## 3.Ana1 37.79255 4.485816 12.55762 8.280142 3.847518 1.8217249 0.8739154
## 4.Ana1 32.16260 5.739837 11.02439 7.804878 6.406504 0.8474948 2.4853602
## 5.Ana1 37.82883 5.243774 11.60996 7.821817 4.612417 1.7468425 1.1674771
## 6.Ana1 29.71566 6.712449 12.28932 9.461613 6.373572 3.3472343 5.6006847
p <- ggplot(dx, aes(CS1, CS2))
p + geom_path(aes( colour=dep) )
p + geom_path(aes( colour=dep) ) + facet_wrap(~dep)
dx1 <- cbind(df, ptawdip$Tli)
p <- ggplot(dx1, aes(CS1, CS2))
p + geom_path(aes( colour=dep) )
p + geom_path(aes( colour=dep) ) + facet_wrap(~dep)
dx2 <- cbind(df, ptawcs$Tli)
p <- ggplot(dx1, aes(CS1, CS2))
p + geom_path(aes( colour=dep) )
p + geom_path(aes( colour=dep) ) + facet_wrap(~dep)
a <- ptawcs$Tli
b <- ptawdip$Tli
a$CSA <- a$CS1
a$CSB <- a$CS2
head(a)
## CS1 CS2 CSA CSB
## 1.Ana1 -0.40098469 -0.02385029 -0.40098469 -0.02385029
## 2.Ana1 0.01648166 0.72266447 0.01648166 0.72266447
## 3.Ana1 -0.16129948 -0.08463197 -0.16129948 -0.08463197
## 4.Ana1 -0.01651174 -0.43113204 -0.01651174 -0.43113204
## 5.Ana1 -0.03646098 -0.75798454 -0.03646098 -0.75798454
## 6.Ana1 0.30454036 -1.08101069 0.30454036 -1.08101069
a <- select(a, CSA, CSB)
head(a)
## CSA CSB
## 1.Ana1 -0.40098469 -0.02385029
## 2.Ana1 0.01648166 0.72266447
## 3.Ana1 -0.16129948 -0.08463197
## 4.Ana1 -0.01651174 -0.43113204
## 5.Ana1 -0.03646098 -0.75798454
## 6.Ana1 0.30454036 -1.08101069
## CSA CSB CS1 CS2
## 1.Ana1 -0.40098469 -0.02385029 -0.49736101 0.19337406
## 2.Ana1 0.01648166 0.72266447 -0.33018361 0.07554173
## 3.Ana1 -0.16129948 -0.08463197 -0.08635054 0.08141458
## 4.Ana1 -0.01651174 -0.43113204 0.52178265 -0.18401691
## 5.Ana1 -0.03646098 -0.75798454 0.04208958 -0.03415939
## 6.Ana1 0.30454036 -1.08101069 0.89155132 -0.37386548
pca <- dudi.pca(c, scannf=FALSE, nf=2)
x <- ggplot(pca$li, aes(x=Axis1, y=Axis2))
x + geom_point()
dfx <- cbind(pca$li, df )
colnames(dfx)
## [1] "Axis1" "Axis2" "region" "dep" "libdep" "wave" "ids"
## [8] "idt" "pop" "p_cs1" "p_cs2" "p_cs3" "p_cs4" "p_cs5"
## [15] "p_cs6" "region" "dep" "libdep" "wave" "ids" "idt"
## [22] "pop" "p_dip0" "p_dip1" "p_dip2" "p_dip3" "p_dip4" "p_dip7"
dfx$wave <- as.factor(dfx$wave)
s.traject(pca$li, dfx$wave)
s.traject(pca$li, dfx$dep)
s.traject(pca$li, dfx$dep, col=rainbow(100))
xx <- ggplot(dfx, aes(x=Axis1, y=Axis2))
xx + geom_path(aes(colour=dep))