Istruzione in R
Ciao a tutti. Sono alle prese con il pacchetto R, poichè devo presentare un progetto all'esame di statistica. E non ho frequentato i laboratori.
So tratta delle seguenti operazioni, che mi consentono con un grafico di controllare il grado di "normalità° dei miei dati.
I dati sono delle temperature
# importo il data set
#verifica dell' assunzione di normalità
Ottengo un grafico con una serie di punti disposti approssimativamente lungo una linea con una certa pendenza.
Ora voglio inserire una linea rossa teorica e lo faccio in questo modo:
quest' ultima istruzione mi fornisce l'errore:
Errore in Im(z_j ~ temp.ord) : argomento non numerico alla funzione
A questo punto avrei dovuto usare l'istruzione:
Ma ovviamente non va in porto poichè gli mancano i valori delle ordinate.
Ripeto, l'istruzione
Ovviamente ho gia fatto uso dell'help, ma non mi ha aiutato.
Non sono assolutamente esperto di R, forse dovrei esserelo. Ma mi servirà solamente per presentere un progetto. Qualcuno puo aiutarmi?
Grazie
So tratta delle seguenti operazioni, che mi consentono con un grafico di controllare il grado di "normalità° dei miei dati.
I dati sono delle temperature
# importo il data set
dati<-read.table('d:/polimi/Statistica/R/Dataset/temperatura.txt', header=T) n<-dim(dati) attach(dati)
#verifica dell' assunzione di normalità
qqnorm(Temperatura,datax=T, main="Norman Probability Plot", xlab="Quantili empirici", ylab="Quantili teorici")
Ottengo un grafico con una serie di punti disposti approssimativamente lungo una linea con una certa pendenza.
Ora voglio inserire una linea rossa teorica e lo faccio in questo modo:
temp.ord<-sort(Temperatura) ranghi<- 1:n F.emp<-(ranghi -0.5)/n z_j<-qnorm(F.emp)) y_j<-Im(z_j ~ temp.ord)$fitted.values
quest' ultima istruzione mi fornisce l'errore:
Errore in Im(z_j ~ temp.ord) : argomento non numerico alla funzione
A questo punto avrei dovuto usare l'istruzione:
lines(temp.ord,y_j, col='red',lwd=2)
Ma ovviamente non va in porto poichè gli mancano i valori delle ordinate.
Ripeto, l'istruzione
y_j<-Im(z_j ~ temp.ord)$fitted.valuesche cosa ha di sbagliato?
Ovviamente ho gia fatto uso dell'help, ma non mi ha aiutato.
Non sono assolutamente esperto di R, forse dovrei esserelo. Ma mi servirà solamente per presentere un progetto. Qualcuno puo aiutarmi?
Grazie
Risposte
Il testo che sto seguendo mi sembra che complichi parecchio le cose.
Cosi è molto piu intuitivo. Grazie.
Cosi è molto piu intuitivo. Grazie.