Analisi statistica multivariata - Modulo R

Guida alle funzioni e comandi di R

Funzione di aiuto

La funzione ?funzione serve per vedere che cosa fanno le varie funzioni, ad esempio ?sqrt.

Funzioni matematiche

Ecco alcune delle principali funzioni matematiche utilizzabili in R:

• exp()
• sin(), cos(), tan()
• asin(), acos(), atan()

Operatori logici

• AND: &
• OR: |
• NOT: !

Operazioni matematiche

Alcune operazioni utili in matematica:

• Valore assoluto: abs(x)
• Radice quadrata: sqrt(x)
• Logaritmo: log(x, base), ad esempio log(5, 10) è il logaritmo di 5 in base 10

Contare gli elementi uguali a qualcosa: sum(x==6). Non fa la somma dei 6 (6+6+6...) ma conta quanti elementi sono uguali a sei. Ad esempio, in un vettore (1,2,6,7,6,8), sum(x==6) restituisce 2.

Arrotondamenti

Per arrotondare:

• Round(n) arrotonda il numero n senza decimali.
• Round(n, digits=3) arrotonda a tre decimali.

Gestione delle variabili

• rm(x) rimuove la variabile x dalla memoria di R.
• objects() serve per vedere le variabili in uso.

Vettori

Creare e gestire vettori:

• c(elementi del vettore), ad esempio x=c(1,2,3)
• length(x) lunghezza del vettore
• mode(x) tipo del vettore
• names(x) ad esempio, con x=c(altezza=171, peso=60) restituisce "altezza" e "peso"
• Per creare un vettore vuoto uso il comando vector, ad esempio x=vector("numeric", len=3) => vettore numerico di lunghezza 3.

Liste

Gestione delle liste:

• x=list()
• Con x[[1]] trovo il primo elemento della lista.
• Con [[1]][2] il secondo elemento del primo elemento della lista.
• Con NA lascio un elemento vuoto.
• Con is.na() vedo quali elementi sono vuoti.

Uso di which

Condizione di tipo logico:

Se avessi x=c(2,5,8,1,2) e mettessi indici = which(x < 5) mi restituisce gli indici degli elementi del vettore <5, ovvero 1 4 5.

Factor

Variabili categoriali non ordinate:

• as.factor(x) rende la variabile factor anche se all’inizio R non la riconosce come tale.
• levels: mi dice quali sono i livelli della variabile. Per riconoscere i livelli la variabile deve essere factor (per vedere se è factor uso: is.factor(x)).

Ripetizioni con REP

Ripete qualcosa un determinato numero di volte:

• Rep(1,2) ripete l’1 per due volte.
• Rep(c(2,3), c(6,7)) ripete il 2 per sei volte e il 3 per sette volte.

Ordinamento

• is.ordered(x): per vedere se il vettore x è ordinato.
• as.ordered(x): mi fa vedere i livelli ordinati.

Frequenze

• Assoluta (fi): numero di volte che si presenta un certo dato. Funzione table()
• Relativa (pi): rapporto tra la assoluta e il numero totale di dati. pi=fi/sum(fi).
• Cumulata (Fi): Fi = cumsum(pi)

Data frame

Crea una matrice in cui posso mettere variabili di natura diversa:

• Per crearla uso: data.frame("nome colonna" = nomi righe)
• Subset(): per selezionare alcune righe del data frame. Ad esempio maschi.alti=subset(classe, subset= genere=="M" & cm>190)
• Subset(classe, select=genere) per selezionare una colonna, posso anche usare semplicemente "classe$genere"

Ordinamento di un vettore

• sort(x): ordina in senso crescente
• sort(x, decreasing = TRUE): ordina in senso decrescente
• order(x) mi restituisce gli indici ordinati, ad esempio con il vettore v=c(2,5,8,1,2), order(v) restituisce 4 1 5 2 3

Grafici

Per creare un data frame e calcolare le frequenze assolute fi con table():

• Per le relative devo dividere per il totale: pi=fi/sum(fi).
• Grafico a torta: comando pie(pi) dove pi è il vettore delle frequenze relative.
• Per aggiungere delle etichette: calcolo per esempio le percentuali con perc=pi*100
• paste() per mettere accanto un nome alla sua percentuale, ad esempio paste(categorie, perc).
• pie(pi, labels = categorie) dove labels serve per dire quali etichette voglio aggiungere.