# Creiamo una sequenza di valori da -5 a 5 con step 0.01
x <- seq(-5, 5, 0.01)
# Calcoliamo la distribuzione normale con media=0 e sd=1 a quei valori
y <- dnorm(x, 0, 1)

# Plot senza assi
# xaxt/yaxt sta per "x (o y) axis type", n sta per "none"
plot(x, y, t="l", xaxt="n", yaxt="n", xlab="y", ylab="f(y)")
# Creiamo le linee verticali tratteggiate (line type=2) a -1, 0, 1
pos <- c(-1, 0, 1)
segments(pos, -1, pos, dnorm(pos, 0, 1), lty=2)
# Alternativamente potevamo usare abline, 
# ma sarebbero arrivate fino in alto
# Ad es. abline(v=c(-1, 0, 1), lty=2)

# Creiamo l'asse x usando axis
# Il primo parametro è l'asse (1 in basso, 2 a sx, 3 in alto, 4 a dx)
# at specifica dove mettere i valori e labels le etichette
# expression ci permette di scrivere (tra le altre cose) lettere greche!
axis(1, at=pos, labels=c(expression(sigma - mu), 
                         expression(mu), expression(sigma + mu)))

Citazione:

---

Messaggio inserito da chick80

# Creiamo l'asse x usando axis
# Il primo parametro è l'asse (1 in basso, 2 a sx, 3 in alto, 4 a dx)
# at specifica dove mettere i valori e labels le etichette
# expression ci permette di scrivere (tra le altre cose) lettere greche!
axis(1, at=pos, labels=c(expression(sigma - mu), 
                         expression(mu), expression(sigma + mu)))

---

Citazione:

---

ascolta che c'è un modo per aumentare la dimensione delle didascalie tipo m-sigma ecc.?
mi serve per altri grafici

---

axis(1, at=pos, labels=c(expression(sigma - mu), 
                         expression(mu), expression(sigma + mu)), cex.axis=1.5)

Citazione:

---

Un po' di tempo fa mi serviva zoomare un plot su un range di valori ben definito

---

Citazione:

---

Inoltre volevo aggiungere i minor tics sugli assi

---

Citazione:

---

Inoltre volevo aggiungere i minor tics sugli assi
O usi minor.ticks del pacchetto hMisc, oppure chiami due volte axis...

In generale io esporto il plot in PDF o in svg e poi me lo spapocchio in Inkscape quindi non mi metto troppo ad abbellirlo in R.
Oppure probablimente ggplot2 ha delle funzioni apposite (ma non lo amo molto)

---

Citazione:

---

tipo con rect io faccio rect(0,0,1,dnorm(x,0,1))
ma mi esce una schifezza :D

---

x.fill <- seq(0, 1, 0.01)
y.fill <- dnorm(x.fill, 0, 1)
# Mettiamo gli estremi a 0 per riempire fino alla base
y.fill[1] <- 0
y.fill[length(y.fill)] <- 0
polygon(x.fill, y.fill, col=rgb(0, 0.2, 1, .5))