##############  ROC krive ##################

# Zadatak 1. Baza-sleep3ED. Da li obelezje hourwnit (broj sati spavanja tokom radnih dana) moze da bude marker
             # za pojavu problema sa spavanjem (probsleeprec)

install.packages("foreign") #instaliranje i ucitavanje paketa
library("foreign")

Sleep=read.spss("E:/baze/sleep3ED.sav", to.data.frame = TRUE) #ucitavanje spss baze
attach(Sleep)

#Primena ROC krivih

install.packages("cutpointr")
library(cutpointr)
ROCcurve=cutpointr(Sleep, hourwnit, probsleeprec, pos_class="yes", direction="<=", #naredba ne radi sa NA poljima
                   method=maximize_metric, metric=sum_sens_spec,na.rm=T) # na.rm=T je obavezan deo
summary(ROCcurve)
plot(ROCcurve)

#Zadatak 2. Baza-Baza; Da li KREATININpre moze da bude dobar marker za prezivljavanje pacijenata (ISHOD)

baza=read.spss("E:/baze/Baza.sav", to.data.frame = TRUE) #ucitavanje spss baze
attach(baza)

install.packages("cutpointr")
library(cutpointr)
ROCcurve=cutpointr(baza, KREATININpre, ISHOD, pos_class="mrtav", direction="<=", #naredba ne radi sa NA poljima
                   method=maximize_metric, metric=sum_sens_spec,na.rm=T) # na.rm=T je obavezan deo
summary(ROCcurve)
plot(ROCcurve)

rm(list=ls())                # Почисти променљиве.
cat("\014")