data(atlas)
op <- par(no.readonly = TRUE)
par(mfrow = c(2,2))
area.plot(atlas$area, cpoin = 1.5)
area.plot(atlas$area, lab = atlas$names.district, clab = 1)
x <- atlas$meteo$mini.jan
names(x) <- row.names(atlas$meteo)
area.plot(atlas$area, val = x, ncl = 12, sub = "Temp Mini January",
csub = 2, cleg = 1)
s.corcircle((dudi.pca(atlas$meteo, scann = FALSE)$co),
clab = 1)
area.plot(atlas$area, val = dudi.pca(atlas$meteo,scann=FALSE)$li[,1],
ncl = 12, sub = "Principal Component Analysis analysis", csub = 1.5,
cleg = 1)
birds.coa <- dudi.coa(atlas$birds, sca = FALSE, nf = 1)
x <- birds.coa$li$Axis1
area.plot(atlas$area, val = x, ncl = 12,
sub = "Correspondence analysis", csub = 1.5, cleg = 1)
s.value(atlas$xy, x, contour = atlas$contour, csi = 2,
incl = FALSE, addax = FALSE)
triangle.plot(atlas$alti)
par(op)
par(mfrow=c(1,1))
Run the code above in your browser using DataLab