Блог

Я пытаюсь построить кривую ROC с точками, но отображаются только 2 из 5 точек, и я не уверен, почему? Что я упускаю? Я неправильно рассчитал объем.значение 3 переменных?( roc_max_inf_vol , roc_med_adult_vol , и roc_med_inf_vol не появляются). Я строю график ( 1 - prop.adults ) против ( 1 - prop.infants ). Каждая точка должна соответствовать volume.value одной из следующих переменных ( roc_max_diff и roc_a1_i только 2 точки отображаются правильно).:

Эти 2 отображаются правильно:

 roc_max_diff <- volume.value[which.max(smooth.difference)]
roc_a1_i <- volume.value[volume.value > max(mydata[mydata$CLASS == "A1" amp; mydata$TYPE == "I", "VOLUME"])]
 

Эти 3 не отображаются:

 roc_max_inf_vol <- volume.value[max(mydata$VOLUME[mydata$TYPE == "I"])] 
roc_med_adult_vol <- volume.value[median(mydata$VOLUME[mydata$TYPE=="ADULT"])]
roc_med_inf_vol <- volume.value[median(mydata$VOLUME[mydata$TYPE=="I"])]
 

Вот что у меня есть:

 cuts <- which(volume.value %in% c(roc_max_inf_vol, roc_med_adult_vol, max_diff, a1_i, roc_med_inf_vol))

ggplot(mapping = aes((1 - prop.infants), (1 - prop.adults)))   
    geom_line(color = 'darkblue', size = 1)   

geom_abline(intercept = 0, slope = 1, color = 'darkred', linetype = 2)  
  geom_point(mapping = aes((1 - prop.infants[cuts]), (1 - prop.adults[cuts])),
             shape = 21, size = 3)  
  annotate("text", label = paste('Median infant nvol =', round(med_inf_vol, 1)),
           x = 0.58, y = 0.90, size = 3)  
  annotate("text", label = paste('Zero A1 infants nvol =', round(a1_i, 1)), 
           x = 0.25, y = 0.90, size = 3)  
  annotate("text", label = paste('Max difference nvol =', round(max_diff, 1)), 
           x = 0.25, y = 0.68, size = 3)  
  annotate("text", label = paste('Median adult nvol =', round(med_adult_vol, 1)), 
           x = 0.12, y = 0.46, size = 3)  
  annotate("text", label = paste('Protect all infants nvol =', round(max_inf_vol, 1)),
           x = 0.09, y = 0.20, size = 3)  
  labs(title = "ROC curve of adult and infant harvest proportions",
       x = "Infant harvest proportion", y = "Adult harvest proportion")   
  theme(axis.text.y = element_text(angle = 90, vjust = 0.5, hjust = 0.5))
 

Для справки, вот исходные данные, которые я использую: https://raw.githubusercontent.com/pandahappy204/Sample/main/mydata.csv

И вот некоторые фоновые переменные, которые, я думаю, необходимы, чтобы понять вышесказанное (например, проп.младенцы и проп.взрослые). Мы очень ценим любую помощь.

 mydata$TYPE <- factor(ifelse(mydata$SEX == "I", "I", "ADULT"))
mydata$TYPE[mydata$CLASS == "A4" | mydata$CLASS == "A5"] <- "ADULT"

idxi <- mydata$TYPE == "I"
idxa <- mydata$TYPE == "ADULT"

max.v <- max(mydata$VOLUME)
min.v <- min(mydata$VOLUME)
delta <- (max.v - min.v)/10000
prop.infants <- numeric(10000)
prop.adults <- numeric(10000)
volume.value <- numeric(10000)

total.infants <- sum(idxi)  
total.adults <- sum(idxa)

for (k in 1:10000) { 
    value <- min.v   k*delta
    volume.value[k] <- value
    prop.infants[k] <- sum(mydata$VOLUME[idxi] <= value)/total.infants
    prop.adults[k] <-  sum(mydata$VOLUME[idxa] <= value)/total.adults
}

y.loess.a <- loess(1 - prop.adults ~ volume.value, span = 0.25,
    family = c("symmetric"))
y.loess.i <- loess(1 - prop.infants ~ volume.value, span = 0.25,
    family = c("symmetric"))
smooth.difference <- predict(y.loess.a) - predict(y.loess.i)