library(XML) library(lattice) # # funkcja readHTMLTable wyciaga tabele z danymi z pliku HTML cancerdata <- readHTMLTable('http://www.guardian.co.uk/commentisfree/2011/oct/28/bad-science-diy-data-analysis', which=1) colnames(cancerdata) <- c('Area','Rate','Population','Number') for (i in 2:4) cancerdata[,i] <- as.numeric(as.character(cancerdata[,i])) summary(cancerdata) # # i wykres czestosci zachorowan od rozmiaru populacji png("zachorowalnoscRakJelita.png",1000,500) ktoreZNazwy <- with(cancerdata, Population^0.25 * Rate > 560 | Rate >= 27) xyplot(Rate~Population, data = cancerdata, pch=19, col="red", panel=function(...) { panel.xyplot(...) with(cancerdata, ltext(Population[ktoreZNazwy], Rate[ktoreZNazwy], Area[ktoreZNazwy], cex=0.8, adj=c(0,-0.5))) }) dev.off()