#############################################################################
###
###  labII.txt - Mauro Gasparini, Ottobre 2002
###              secondo laboratorio di R per Ingegneri
###      
###              Riferimento: Devore, J.L. (2000)
###                           Probability and Statistics
###                           for Engineering and the Sciences
###                           Duxbury, 5th edition
###
#############################################################################

# A differenza del primo laboratorio, tutti i commenti sono in questa forma.
# L'intestazione e il formato del testo costituiscono infatti  
# un semplice prototipo di codice di un generico progetto, 
# leggibile dentro R con il comando source().


library(Devore5)  # chiama la libreria con i dati del libro  di Devore

# in alternativa, si provi a scaricare su un file di testo i dati 
# per i singoli esempi o esercizi dal sito web 
# http://calvino.polito.it/~gasparin/teaching/devore/
# facendo quanto illustrato nel primo laboratorio

# carica i dati  di xmp01.01, cioè l'esempio 1.1 
data(xmp01.01)

# fai una breve descrizione di cosa è xmp01.01
str(xmp01.01)

# fai la stessa cosa con l'esercizio 1.23 (svolto durante le esercitazioni)
data(ex01.23)
str(ex01.23)

# questa non funziona, perchè la variabile cycles non è riconosciuta
hist(cycles, prob=T) 

# porta in prima fila (nel path interno) i dati
attach(ex01.23)

# stavolta l'istogramma di cycles funziona
# prob=T è una opzione che  normalizza l'area totale a 1 come desiderato
hist(cycles, prob=T) 

# leggi la miriade di opzioni per l'istogramma
help(hist)
# e applicalo agli altri esercizi svolte in classe.

# assegna cycles a x (vedi primo laboratorio)
x <- cycles

# Altri strumenti di descrittiva:

# un grafico per punti
rug(x)

# un grafico a barre
barplot(x) 

# media e mediana
mean(x) median(x) 

# il riassunto delle cinque statistiche
summary(x) 

# un grafico a rami e foglie
stem(x) 

# un boxplot
boxplot(x) 

# misure di variabilità
var(x); sd(x) 

# una trasformazione di x
log(x) 

# altre trasformazioni di x
log(x, 10)
sqrt(x)
x^(1/2)

# ora si rifacciano gli esercizi del primo capitolo in R.
# Si usino le opzioni ricavate dagli help().
# Per esempio, come si cambiano le classi in un istogramma?
# Come si fanno boxplot appaiati?


# si può anche trarre ispirazione da questo demo
demo(graphics)

# per salvare i grafici, usare il menu della finestra grafica.
# Spesso è meglio salvare i comandi che producono il grafico
# invece del grafico stesso. Questo file ne è un esempio.