मैं इस समुदाय से चार प्रकार के तिरछे पर अधिक अंतर्दृष्टि की उम्मीद कर रहा हूं।
मेरे द्वारा निर्दिष्ट प्रकार http://www.inside-r.org/packages/cran/e1071/docs/skewness सहायता पृष्ठ में उल्लिखित हैं ।
मदद पृष्ठ में पुरानी पद्धति का उल्लेख नहीं किया गया था, लेकिन मैं इसे फिर भी शामिल करता हूं।
require(moments)
require(e1071)
x=rnorm(100)
n=length(x)
hist(x)
###############type=1
e1071::skewness(x,type=1)
sqrt(n) * sum((x-mean(x))^3)/(sum((x - mean(x))^2)^(3/2)) #from e1071::skewness source
m_r=function(x,r) {n=length(x); sum((x - mean(x))^r/n);} ##from e1071::skewness help
g_1=function(x) m_r(x,3)/m_r(x,2)^(3/2)
g_1(x) ##from e1071::skewness help
moments::skewness(x) ##from e1071::skewness help
(sum((x - mean(x))^3)/n)/(sum((x - mean(x))^2)/n)^(3/2) ##from moments::skewness code, exactly as skewness help page
###############type=2
e1071::skewness(x,type=2)
e1071::skewness(x,type=1) * sqrt(n * (n - 1))/(n - 2) #from e1071::skewness source
G_1=function(x) {n=length(x); g_1(x)*sqrt(n*(n-1))/(n-2);} #from e1071::help
G_1(x)
excel.skew=function(x) { n=length(x); return(n/((n-1)*(n-2))*sum(((x-mean(x))/sd(x))^3));}
excel.skew(x)
###############type=3
e1071::skewness(x,type=3)
e1071::skewness(x,type=1) * ((1 - 1/n))^(3/2) #from e1071::skewness source
b_1=function(x) {n=length(x); g_1(x)*((n-1)/n)^(3/2); } #from e1071::skewness help page
b_1(x);
prof.skew=function(x) sum((x-mean(x))^3)/(length(x)*sd(x)^3);
prof.skew(x)
###############very old method that fails in weird cases
(3*mean(x)-median(x))/sd(x)
#I found this to fail on certain data sets as well...
यहाँ वह कागज है जिसे e1071 के लेखक ने संदर्भित किया है: http://onlinelibrary.wiley.com/doi/10.1111/1467-9884.00122/pdf जॉयन्स और सीए गिल (1998), नमूना स्केचनेस और कुर्टोसिस के तुलनात्मक उपाय।
उस कागज के मेरे पढ़ने से, वे सुझाव देते हैं कि टाइप # 3 में सबसे कम त्रुटि है।
उपरोक्त कोड से तिरछेपन के उदाहरण इस प्रकार हैं:
e1071::skewness(x,type=1)
-0.1620332
e1071::skewness(x,type=2)
-0.1645113
e1071::skewness(x,type=3)
-0.1596088
#old type:
0.2694532
मैंने यह भी देखा कि ई 1071 के लेखक ने हेल्प पेज में नोटों से अलग तिरछा फंक्शन लिखा था। Sqrt को नोटिस करें:
sqrt(n) * sum((x-mean(x))^3)/(sum((x - mean(x))^2)^(3/2)) #from e1071::skewness source
(sum((x - mean(x))^3)/n)/(sum((x - mean(x))^2)/n)^(3/2) #from moments and e1071 help page
कोई विचार क्यों sqrt (n) पहले समीकरण में है? कौन सा समीकरण ओवरफ्लो / अंडरफ्लो को बेहतर तरीके से हैंडल करता है? किसी भी अन्य विचार क्यों वे अलग हैं (लेकिन समान परिणाम उत्पन्न करते हैं)?