आर

2011-10-18 11 views
5

में एयूसी अनुपात की गणना करना मैं प्रजातियों के एक सेट के लिए पारिस्थितिकीय विशिष्ट मॉडल उत्पन्न कर रहा हूं और मैं एयूसी को पारिस्थितिकीय विशिष्ट गुणवत्ता के लिए मीट्रिक के रूप में उपयोग करना चाहता हूं। मैक्सेंट विकसित करने वाले स्टीवन फिलिप्स ने आर में एयूसी की गणना के लिए अपने मैक्सेंट मैनुअल में कोड प्रदान किया है। हालांकि, मैं उन कागजात पढ़ रहा हूं जो आंशिक एयूसी अनुपात को अधिक मजबूत और अवधारणात्मक ध्वनि मीट्रिक के रूप में रिपोर्ट करते हैं। मुझे लगता है कि मैं समझता हूं कि आरओसीआर आर पैकेज का उपयोग करके आंशिक एयूसी की गणना कैसे करें, लेकिन एयूसी अनुपात की गणना कैसे करता है?आर

यहाँ फिलिप्स से ट्यूटोरियल स्क्रिप्ट है:

presence<-read.csv("bradypus_variegatus_samplePredictions.csv") 
background<-read.csv("bradypus_variegatus_backgroundPredictions.csv") 
pp<-presence$Logistic.prediction 
testpp<-pp[presence$Test.or.train=="test"] 
trainpp<-pp[presence$Test.or.train=="train"] 
bb<-background$logistic 

combined<-c(testpp,bb) 
label<-c(rep(1,length(testpp)),rep(0,length(bb))) 
pred<-prediction(combined,label) 
perf<-performance(pred,"tpr","fpr") 
plot(perf,colorize=TRUE) 
performance(pred,"auc")@y.values[[1]] #RETURNS AUC 

AUC<-function(p,ind){ 
    pres<-p[ind] 
    combined<-c(pres,bb) 
    label<-c(rep(1,length(pres)),rep(0,length(bb))) 
    predic<-prediction(combined,label) 
    return(performance(predic,'auc')@y.values[[1]]) 
} 

b1<-boot(testpp,AUC,100) #RETURNS AUC WITH STANDARD ERROR 
b1 

कोई सलाह या सुझाव बहुत सराहना की जाएगी! धन्यवाद।

उत्तर

2
आपके डेटासेट और आवेदन की बारीकियों को जानने के बिना

,

  • आंशिक एयूसी: वक्र के केवल एक हिस्से के तहत क्षेत्र। (आमतौर पर चुना जाता है क्योंकि यह अधिक मजबूत या अन्यथा वांछनीय है, जैसा कि आपने कहा था)
  • एयूसी अनुपात: एक एयूसी का अनुपात दूसरे में। (आमतौर पर किसी प्रकार का संदर्भ)

सू ...

  • आंशिक एयूसी अनुपात: एक से दूसरे आंशिक एयूसी के अनुपात।
+0

तो तुम कह रहे हैं कि एयूसी अनुपात वास्तव में एक आला मॉडल की गुणवत्ता का एक उपाय नहीं है, लेकिन एक मीट्रिक के और अधिक जब कुछ कुछ पैरामीटर बदलती के साथ कई मॉडलिंग आउटपुट की तुलना? – Pascal

+0

सही, जैसा कि मैंने देखा है कि यह आमतौर पर उपयोग किया जाता है। यह अभी भी संभव है कि आपके क्षेत्र में कुछ और विशिष्ट अर्थ है। मुझे पारिस्थितिक आला मॉडल के बारे में कुछ भी पता नहीं है, इसलिए दुर्भाग्य से मैं वहां कोई मदद नहीं कर रहा हूं। वैसे भी, यहां एक उदाहरण दिया गया है कि मैं इसे आम तौर पर कैसे देखता हूं, अगर यह मदद करता है: http://www.drug-interactions.eu/eng/AUC.htm –

+0

ठीक है, तो आप शायद सही हैं। मुझे लगता है कि मैं जो पेपर पढ़ रहा था उसमें एयूसी अनुपात के साथ दिखाने की कोशिश कर रहा था, मुझे गलतफहमी हो रही थी। स्पष्टीकरण के लिए धन्यवाद! – Pascal

1

पैकेज आरओसीआर fpr.stop= पैरामीटर का उपयोग करके आंशिक एयूसी मानों की गणना कर सकता है। जैसा कि जॉन ने कहा था कि अनुपात सिर्फ यह मान आपके संदर्भ मॉडल के लिए समान गणना से विभाजित है।