data(clementines)
op <- par(no.readonly = TRUE)
par(mfrow = c(5,4)) ; par(mar = c(2,2,1,1))
for (i in 1:20) {
w0 <- 1:15 ; plot(w0, clementines[,i], type = "b")
abline(lm(clementines[,i] ~ w0))
}
par(op)
pca1 <- dudi.pca(clementines, scan = FALSE)
s.corcircle(pca1$co, clab = 0.75)
barplot(pca1$li[,1])
op <- par(no.readonly = TRUE)
par(mfrow = c(5,4)) ; par(mar = c(2,2,1,1))
clem0 <- pca1$tab
croi <- 1:15
alter <- c(rep(c(1,-1),7),1)
for (i in 1:20) {
y <- clem0[,i]
plot(w0, y, type = "b", ylim = c(-2,2))
z <- predict(lm(clem0[,i] ~ croi*alter))
points(w0, z, pch = 20, cex = 2)
for(j in 1:15) segments (j,y[j],j,z[j])
}
par(op)
Run the code above in your browser using DataLab