r ऑनल वक्र के तहत क्षेत्र के 95% को कवर करने वाले दो नंबर ढूंढने के लिए खोजें



मोबाइल नंबर लोकेशन ऐप्स (1)

पृष्ठभूमि:

मेरे पास दो कॉलम नंबर हैं (स्तंभ " dd ", और स्तंभ " yy ", इस डेटा तक पहुंचने के लिए नीचे देखें )। जब मैं " yy " ( वाई-अक्ष ) के विरुद्ध "डीडी" ( एक्स-अक्ष ) का प्रयोग करता हूं, तो इसका उपयोग करते हुए:

plot(dd, yy, ty="l", xlab = "dd", ylab ="yy", lwd=2, xaxs="i")

मैंने इसे प्राप्त किया:

सवाल कोडिंग:

"-इन्फ़" से "0" में आगे बढ़कर जब तक अन्य बिंदु को इस वक्र के नीचे क्षेत्र का 95% तक शामिल नहीं किया जाता है? दूसरे शब्दों में, यदि हम इस वक्र के नीचे तीर का 95% हिस्सा ढकेल रहे हैं , तो हम दो तीरों के दो संख्याओं को देखना चाहते हैं, जो कि तीर के शीर्षकों से संबंधित हैं ?

आपके Rstudio में आसानी से इस छोटे डेटा को आयात करने के लिए आर कोड है :

id <- "0B5V8AyEFBTmXM1VIYUYxSG5tSjQ"

Data <- read.csv(paste0("https://docs.google.com/uc?id=",id,"&export=download")) 

संपादित करें: मेरे पिछले सन्निकटन बहुत सटीक नहीं थे। यह पता चला है कि pracma पैकेज में एक फ़ंक्शन, cumtrapz , जो cumtrapz साथ समता का अनुमान लगाता है और बाएं से शुरू होने वाले प्रत्येक बिंदु पर योग को आउटपुट करता है। उस फ़ंक्शन का उपयोग करना एक अधिक सटीक उत्तर देता है।

library(pracma)
points <- read.csv(paste0("https://docs.google.com/uc?id=",id,"&export=download")) 
area   <- cumtrapz(points$dd, points$yy)
ind    <- rev(which(area[length(area)] - area > .95))[1]

point  <- points$dd[ind]
point

# [1] -0.06406406




confidence-interval