Διαγράμματα με την R
22/11/2023
Απλό διάγραμμα διασποράς σημείων
Ας υποθέσουμε πως θέλουμε να παρουσιάσουμε σε διάγραμμα τιμές μιας μεταβλητής. Εδώ έχουμε ένα παράδειγμα με 5 τιμές τις οποίες τοποθετούμε σε διάνυσμα με όνομα y.
Η εντολή plot μπορεί να χρησιμοποιηθεί για τη γραφική αναπαράσταση:
Οριζόντιος άξονας
Στο οριζόντιο άξονα υπάρχει απλά ένας σειριακός δείκτης, 1,2,3,4,5 που αντιστοιχεί στο 1:length(y). Αν επιθυμούμε μπορούμε να το τοποθετήσουμε απευθείας στην εντολή:
Τύπος σημείων γραμμών
Αν δεν ορίσουμε τον τύπο γραφήματος, τότε θα τοποθετηθούν σημεία.
Αυτό υπονοεί την επιλογή type = “p”.
Μπορούμε να τα αλλάξουμε σε γραμμές που ενώνουν τα σημεία.
Παράδειγμα
26 διαφορετικά σύμβολα σημείου
Code
k <- 0:25
plot(k, k, pch = k, cex = 2, axes = FALSE, xlab = "", ylab = "", ylim = c(-1, 25))
text(k, k, labels = k, pos = 1, cex = 1.5)
6 τύποι γραμμών
Code
plot(0:1, c(1, 1), type = "l", lwd = 2, lty = 1, ylim = c(1, 6.5), axes = FALSE, xlab = "", ylab = "")
for (i in 2:6) {
lines(0:1, c(i, i), lwd = 2, lty = i)
}
text(0, 1:6, labels = 1:6, pos = 2, cex = 2)
Γράφημα με δύο σειρές δεδομένων
Code
y1 <- c(12, 34, 21, 15, 20)
y2 <- c(10, 30, 18, 19, 22)
d <- c("ΔΕΥ", "ΤΡΙ", "ΤΕΤ", "ΠΕΜ", "ΠΑΡ")
plot(y1, col = "green", ylab = "Ποσότητα", type = "b",
ylim = c(0, 40), xaxt = "n", xlab = "Ημέρα")
lines(y2, col = "red", type = "b", pch = 2)
axis(1, labels = d, at = 1:length(d))
legend("topright", c("Amstel", "Heineken"), col = c("red", "green"), pch = c(1, 2), lty = 1)
Τροποποίηση
Code
y1 <- c(12, 34, 21, 15, 20)
y2 <- c(10, 30, 18, 19, 22)
d <- c("ΔΕΥ", "ΤΡΙ", "ΤΕΤ", "ΠΕΜ", "ΠΑΡ")
plot(y1, type = "n", xaxt = "n", ylim = c(0,40),
xlab = "Ημέρα", ylab = "Ποσότητα")
lines(y1, col = "red", type = "b", pch = 1, lwd = 1.5, lty = 2)
lines(y2, col = "green", type = "b", pch = 2, lwd = 1.5, lty = 1)
axis(1, labels = d, at = 1:length(d))
legend(3.5, 40, c("Amstel", "Heineken"), col = c("red","green"),
pch = c(1, 2), lwd = 2, lty = c(2, 1))
Απλό ραβδόγραμμα
Τροποποίηση, διαφορά από τη μέση τιμή
Ραβδόγραμμα με δύο σειρές δεδομένων
Υπολογισμοί και οπτικοποίηση σε θηκόγραμμα
summary(mtcars$mpg) Min. 1st Qu. Median Mean 3rd Qu. Max.
10.40 15.43 19.20 20.09 22.80 33.90 IQR(mtcars$mpg)[1] 7.375median(mtcars$mpg) - 1.5*IQR(mtcars$mpg)[1] 8.1375median(mtcars$mpg) + 1.5*IQR(mtcars$mpg)[1] 30.2625boxplot(mtcars$mpg)
Ένα απλό ιστόγραμμα
Το ιστόγραμμα γίνεται με την εντολή hist.
Για παράδειγμα, το ιστόγραμμα μιας μεταβλητής που έχει 1000 τυχαίους αριθμούς της τυπικής κανονικής κατανομής είναι:
y <- rnorm(1000)
hist(y)
Παράμετροι ιστογράμματος
Όπως και σε άλλες περιπτώσεις, υπάρχουν διάφορες επιλογές που μπορεί ο χρήστης να “επέμβει”. Πχ, να αλλάξει το πλήθος των των τάξεων (nclass), να αλλάξει τα όρια των αξόνων (xlim, ylim) ή να αλλάξει την απεικόνιση από συχνότητα σε πυκνότητα (freq = FALSE).
y <- rnorm(1000)
hist(y, nclass = 15, xlim = c(-4, 4), ylim = c(0, 0.5), freq = FALSE)
Πολύγωνο συχνοτήτων
