#### Statistik 1 : Vorlesung am 21.12.2015 ## Paket Epi sollte geladen werden ### library(Epi) ROC(c(1,2,3,4,5),c(0,0,1,0,1)) ROC(c(1,2,3,4,5),c(0,0,0,1,1)) ### extremer Fall daten zu bei 1 alle größer als bei x=0 ROC(c(1,3,3,4,5),c(0,0,1,0,1)) ROC(c(1,1,2,2,3,3),c(0,1,0,1,0,1)) ### Anderes Extrem. Kein Zusammenhang! ### Staubdaten einlesen #### staub <- read.table("dust.asc", header=TRUE, sep="", na.strings="NA", dec=".", strip.white=TRUE) attach(staub) head(staub) ## die ersten werte ## ROC(dust,cbr) ROC(smoking,cbr) ldust<- log(1+dust) ROC(ldust,cbr)