यहां कुछ आर कोड / उदाहरण दिए गए हैं जो आपको एक फिट और एक तंदुरुस्त फिट के लिए फिट की तुलना करने देंगे:
library(TeachingDemos)
library(splines)
tmpfun <- function(x,y,span=.75,df=3) {
plot(x,y)
fit1 <- lm(y ~ ns(x,df))
xx <- seq( min(x), max(x), length.out=250 )
yy <- predict(fit1, data.frame(x=xx))
lines(xx,yy, col='blue')
fit2 <- loess(y~x, span=span)
yy <- predict(fit2, data.frame(x=xx))
lines(xx,yy, col='green')
invisible(NULL)
}
tmplst <- list(
span=list('slider', from=0.1, to=1.5, resolution=0.05, init=0.75),
df=list('slider', from=3, to=25, resolution=1, init=3))
tkexamp( tmpfun(ethanol$E, ethanol$NOx), tmplst )
आप इसे अपने डेटा के साथ आज़मा सकते हैं और अन्य प्रकारों या विकल्पों को आज़माने के लिए कोड बदल सकते हैं। आप यह भी देख सकते हैं कि loess.demo
टीचिंगडेमोस पैकेज में फंक्शन को बेहतर तरीके से समझने के लिए लूप अल्गोरिथम क्या करता है। ध्यान दें कि आप जिस चीज को देखते हैं, वह अक्सर दूसरे इंटरपोलेशन स्मूदी (कभी-कभी खुद की एक पट्टी) के साथ लॉयस का एक संयोजन होता है, loess.demo
फ़ंक्शन वास्तव में स्मूथ और कच्ची दोनों को फिट दिखाता है।
सैद्धांतिक रूप से आप हमेशा एक ऐसी पट्टी पा सकते हैं जो किसी अन्य निरंतर फ़ंक्शन को जितना आप चाहें उतना करीब से पता लगा सकती है, लेकिन यह संभावना नहीं है कि समुद्री मील का एक सरल विकल्प होगा जो मज़बूती से किसी भी डेटा सेट के लिए एक करीब फिट के लिए एक अनुमानित सन्निकटन देगा।