setwd("c:/") library(vioplot) library(car) library(MASS) library(R2HTML) library(Hmisc) # # wczytujemy dane # wyniki = read.table("mimLista14.csv", header=T, sep=";") poziomy = c(0.5,0.75,0.9,0.95) kwantyle = quantile(wyniki[,9],poziomy,na.rm=T) # # przygotujmy rysunki z podsumowaniami # png("gr1.png",700,500) plot(ecdf(wyniki[,9]), main="Dystrybuanta empiryczna") for (i in seq_along(poziomy)) { lines(rep(kwantyle[i],2), c(0,poziomy[i]), col="red") lines(c(-5,kwantyle[i]), rep(poziomy[i],2), col="red") text(kwantyle[i]+1,0.2,paste("centyl",poziomy[i]*100,"=",kwantyle[i]),srt=270, adj=c(0.5,0)) } dev.off() png("gr2.png",600,400) boxplot(wyniki[,5:8],notch=T,varwidth=T) dev.off() png("gr3.png",600,400) boxplot(wyniki[,9]~wyniki[,4],notch=T,varwidth=T) dev.off() png("gr4.png",600,400) plot(wyniki[,3], wyniki[,9], xlab="nr. indeksu", ylab="punkty", main = "Punkty a rok studiow") abline(rlm(wyniki[,9]~ wyniki[,3])) dev.off() plec = factor(sapply(as.character(wyniki[,1]), function(x) substr(x, nchar(x), nchar(x))=="a"),labels = c("M","K")) png("gr5.png",600,400) out <- histbackback(split(wyniki[,9], plec), brks=seq(0,40,2), xlim=c(-10,10), main="Punkty a plec") barplot(-out$left, col="red" , horiz=TRUE, space=0, add=TRUE, axes=FALSE) barplot(out$right, col="blue", horiz=TRUE, space=0, add=TRUE, axes=FALSE) dev.off() tabelka = wyniki[order(wyniki[,3]),c(3,5,6,7,8,9)] tabelka = tabelka[!is.na(tabelka[,2]),] # # wyswietlamy wyniki w formacie HTML # cat("

Wyniki kolokwium ze Statystyka I MIMUW 2007/2008

",file="index.html") HTML(tabelka,file="index.html",row.names=F,innerBorder=1,Border=2) sink(file="index.html",append=T) cat("

Podsumowania

") cat("

Skrypt wyswietlajacy podsumowania") sink()