एमटीएस ऑब्जेक्ट्स पर फ़ंक्शंस के परिवार को लागू करने का उपयोग करना

Question

किसी mts ऑब्जेक्ट पर लागू (या sapply) का उपयोग करना फ़ंक्शन पर भेजते समय इसकी टाइम श्रृंखला गुणों को हटा देता है। मुझे एमटीएस ऑब्जेक्ट में प्रत्येक बार श्रृंखला पर एक ही फ़ंक्शन (टीएस इनपुट और टीएस आउटपुट के साथ) कैसे लागू करना चाहिए और इसे वापस करना (अधिमानतः एमटीएस के रूप में) [मेरा मतलब है लूप के लिए उपयोग करने के अलावा]?

उदाहरण के लिए मैं एक समारोह है कि एक नमूना मीटर

z <- ts(matrix(rnorm(90), 30, 3), start=c(1961, 1), frequency=4) 
class(z) 

केवल समय श्रृंखला कार्यों में से एक के लिए भेजा जा रहा है अब एक समय श्रृंखला की प्रवृत्ति (एसटीएल का उपयोग कर)

myfunc <- function(x) { 
     return(stl(x,"per")$time.series[,2]) 
} 

रिटर्न लिखने मान लें कि सही:

myfunc(z[,1]) # works correctly, returns the trend of first series 

मेरा फ़ंक्शन एकाधिक समय श्रृंखला के लिए डिज़ाइन नहीं किया गया है:

myfunc(z) # will not work returning the error below 

Error in stl(x, "per") : only univariate series are allowed 

मीटर पर लागू होते हैं का उपयोग करते हुए आपत्ति वेक्टर के रूप में समय श्रृंखला के प्रत्येक भेजने के लिए, अपने समय श्रृंखला गुण (चम्मच) संरक्षण नहीं:

apply(z,2,myfunc) # will not work returning the error below 

Error in stl(x, "per") : 
series is not periodic or has less than two periods 

Answer 1

इस के चारों ओर एक सरल तरीका है, एक साफ apply के बजाय सूचकांक के साथ काम करने के लिए है:

sapply(seq_len(ncol(z)),function(i) myfunc(z[,i])) 

apply समारोह के अंदर साफ वैक्टर डालता है, क्योंकि यह पहले एक मैट्रिक्स को कोई चीज बदल देता है। टाइम श्रृंखला ऑब्जेक्ट्स के लिए परिभाषित [ फ़ंक्शन का उपयोग करके, आप सुनिश्चित हैं कि आप प्रत्येक बार एक वैध समय श्रृंखला निकालते हैं।

Answer 2

मैं अगर यह एक है की जाँच करने के myfunc बदल पैरामीटर एक्स के रूप में ts वस्तु।

यदि एक्स टीएस नहीं है, तो इसे 0 ऑब्जेक्ट में बदल दिया गया है क्योंकि stl इस पैरामीटर प्रकार की आवश्यकता है।

myfunc <- function(x,...){ 
     y <- x 
     if(class(x) != 'ts') { 
     dots <- c(...) 
     y <- ts(x,start=c(dots[1], dots[2]), frequency=dots[3]) 
     } 
     return(stl(y,"per")$time.series[,2]) 
    } 
    ## no need to conversion (already ts object) 
    myfunc(z[,1]) 


    ## mts object (here we give parameter necessary for conversion) 
    apply(z,2,myfunc,1961,1,4)