आर: एकीकृत करें: उपखंडों की अधिकतम संख्या तक पहुंच गई, राउंडऑफ त्रुटि

Question

मैं एक दिलचस्प लेकिन बदसूरत समस्या में आया।

मैं एक ऐसे फ़ंक्शन को एकीकृत करने की कोशिश कर रहा हूं जिसे डेटासेट से गणना की गई है। डेटा यहां पाया जा सकता है: Link to sample.txt।

मैं अपने डेटा पर एक लाइन फिट करके शुरू करता हूं। इसे approxfun या splinefun के साथ गैर-रैखिक के साथ रैखिक किया जा सकता है। नीचे दिए गए मेरे उदाहरण में मैं बाद वाले का उपयोग करता हूं। अब, जब मैं फिट समारोह को एकीकृत करने की कोशिश मैं त्रुटि

• maximum number of subdivisions reached

में चलाने लेकिन जब मैं उपखंड वृद्धि मैं

• roundoff error

मिल मेरे नमूना कोड में मानों से, आप इस विशिष्ट डेटा के लिए देख सकते हैं सीमा निर्धारित करें 754-> 755 है।

मेरे सहयोगी को मैटलैब में इस डेटासेट को एकीकृत करने में कोई समस्या नहीं है। एकीकृत करने के लिए मेरे डेटा में हेरफेर करने का कोई तरीका है? आर में संख्यात्मक एकीकरण के लिए एक और तरीका है?

data<-read.table('sample.txt',sep=',') 
colnames(data)<-c('wave','trans') 
plot(data$wave,data$trans,type='l') 

trans<- -1 * log(data$trans) 
plot(data$wave,trans,type='l') 

fx.spline<-splinefun(data$wave,trans) 

#Try either 
Fx.spline<-integrate(fx.spline,min(data$wave),max(data$wave)) 
#Above: Number of subdivision reached 
Fx.spline<-integrate(fx.spline,min(data$wave),max(data$wave),subdivisions=754) 
#Above: Number of subdivision reached 
Fx.spline<-integrate(fx.spline,min(data$wave),max(data$wave),subdivisions=755) 
#Above: Roundoff error 

Answer 1

आर में कई एकीकरण दिनचर्या रहे हैं, और आप 'RSiteSearch'ing या का उपयोग करके' द्वारा एसओएस 'पैकेज उनमें से कुछ पा सकते हैं।

उदाहरण के लिए, पैकेज pracma, कई कार्यान्वयन है उदाहरण के लिए

quad(fx.spline,min(data$wave),max(data$wave)) # adaptive Simpson 
# [1] 2.170449         # 2.5 sec 
quadgk(fx.spline,min(data$wave),max(data$wave)) # adaptive Gauss-Kronrod 
# [1] 2.170449         # 0.9 sec 
quadl(fx.spline,min(data$wave),max(data$wave)) # adaptive Lobatto 
# [1] 2.170449         # 0.8 sec 

कृपया नहीं है कि इन शुद्ध आर स्क्रिप्ट और इसलिए की तुलना में धीमी है, जैसे, इस तरह के एक दोलन समारोह के साथ संकलित integrate दिनचर्या है।