plotline<-function(x,y) { #Description # Da un grafico sintetico della variabilità delle medie di una serie di dati suddivisa in fattori #Usage # plotline(x,y) #Arguments # x Vettore dei dati # y Vettore numerico indicante i fattori in cui sono suddivisi i dati #Value # Grafico rappresentante per ogni livello del fattore y una linea compresa fra media del livello # e +/- 2 volte lo scarto quadratico medio sempre dello stesso livello yy<-unique(y) l<-length(yy) media<-tapply(x,y,mean) sqm<-sqrt(tapply(x,y,var)) z<-sort(media) z<-sort(sqm) plot(yy,media,xlab="Fattori",ylab="Valori",ylim=c(media[1]-2*sqm[l],media[l]+2*sqm[l]),main="variabilita' delle medie") rug(y) for (i in 1:l) { segments(yy[i],media[i]-sqm[i],yy[i],media[i]+sqm[i]) } }