Guest User

Untitled

a guest
Oct 18th, 2017
107
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.92 KB | None | 0 0
  1. a=c(0,15,30,45,60,90,120)
  2. b=c(4.3,6.3,7.9,5.3,4.1,4.6,4.9)
  3.  
  4. incrementalAUC <- vector()
  5.  
  6. for (i in 1:(length(a)-1)){
  7.  
  8. if((b[i+1]-b[1] >= 0) && (b[i]-b[1] >= 0))
  9. {incrementalAUC[i] <-((b[i]-b[1]+b[i+1]-b[1])/2)*(a[i+1]-a[i])}
  10.  
  11. else if((b[i+1]-b[1] < 0) && (b[i]-b[1] >= 0))
  12. {incrementalAUC[i] <-(b[i]-b[1])*((b[i]-b[1])/(b[i]-b[i+1])*(a[i+1]-a[i])/2)}
  13.  
  14. else if((b[i+1]-b[1] >= 0) && (b[i]-b[1] < 0))
  15. {incrementalAUC[i] <-(b[i+1]-b[1])*((b[i+1]-b[1])/(b[i+1]-b[i])*(a[i+1]-a[i])/2)}
  16.  
  17. else if((b[i]-b[1] < 0) && (b[i+1]-b[1] < 0))
  18. {incrementalAUC[i] <- 0}
  19. }
  20. incrementalAUC
  21. sum(incrementalAUC)
  22.  
  23. ID Group calculated AUC
  24. ID=1 Group=1
  25. ID=1 Group=2
  26. ID=2 Group=1
  27.  
  28. a <-rep(c(0,15,30,45,60,90,120), 4)
  29. Id<-rep(1:2,each=14)
  30. group<-rep(1:2,each=7, times = 2)
  31. b=c(4.3,6.3,7.9,5.3,4.1,4.6,4.9,4.3,6.3,7.9,5.6,4.4,4.5,4.8,4.3,6.3,8.9,9.3,4.1,4.6,5.9,4.3,7.3,8.9,7.3,4.1,5.6,4.9)
  32. data<-data.frame(Id,group,a,b)
Add Comment
Please, Sign In to add comment