मुझे लगता है कि उत्तर हां होना चाहिए, लेकिन मुझे अभी भी लगता है कि कुछ सही नहीं है। साहित्य में कुछ सामान्य परिणाम होने चाहिए, क्या कोई मेरी मदद कर सकता है?
मुझे लगता है कि उत्तर हां होना चाहिए, लेकिन मुझे अभी भी लगता है कि कुछ सही नहीं है। साहित्य में कुछ सामान्य परिणाम होने चाहिए, क्या कोई मेरी मदद कर सकता है?
जवाबों:
नहीं।
तीन चर, , Y और Z = X + Y पर विचार करें । उनके सहप्रसरण मैट्रिक्स, एम , वहाँ एक वेक्टर है, क्योंकि, सकारात्मक निश्चित नहीं है z ( = ( 1 , 1 , - 1 ) ' जिसके लिए) जेड ' एम जेड सकारात्मक नहीं है।
जनसंख्या सहसंयोजक matrices सकारात्मक अर्द्ध निश्चित हैं।
( यहां संपत्ति 2 देखें )
आम तौर पर पूर्ण नमूनों (कोई लापता मूल्यों) के सहसंयोजक मैट्रिक्स पर लागू नहीं होना चाहिए, क्योंकि उन्हें असतत जनसंख्या सहसंयोजक के रूप में भी देखा जा सकता है।
हालाँकि, अस्थायी बिंदु संख्यात्मक संगणना की अक्षमता के कारण, यहां तक कि बीजगणितीय रूप से सकारात्मक निश्चित मामलों को कभी-कभी सकारात्मक अर्ध-निश्चित नहीं होने के लिए गणना की जा सकती है; एल्गोरिदम का अच्छा विकल्प इसमें मदद कर सकता है।
आम तौर पर, नमूना सहसंयोजक परिपक्वता - इस बात पर निर्भर करता है कि वे कुछ चर में लापता मूल्यों से कैसे निपटते हैं - सिद्धांत में भी सकारात्मक अर्ध-निश्चित या नहीं हो सकता है। उदाहरण के लिए, यदि जोड़ीदार विलोपन का उपयोग किया जाता है, तो सकारात्मक अर्ध-निश्चितता की कोई गारंटी नहीं है। इसके अलावा, संचित संख्यात्मक त्रुटि के कारण नमूना कोविरियस मैट्रिस हो सकते हैं जो होने के लिए विफल होने के लिए सकारात्मक रूप से सकारात्मक अर्ध-निश्चित होना चाहिए।
इस तरह:
x <- rnorm(30)
y <- rnorm(30) - x/10 # it doesn't matter for this if x and y are correlated or not
z <- x+y
M <- cov(data.frame(x=x,y=y,z=z))
z <- rbind(1,1,-1)
t(z)%*%M%*%z
[,1]
[1,] -1.110223e-16
यह मेरे द्वारा किए गए पहले उदाहरण पर हुआ (मुझे संभवतः एक बीज की आपूर्ति करनी चाहिए लेकिन यह इतना दुर्लभ नहीं है कि आपको एक प्राप्त करने से पहले बहुत सारे उदाहरणों का प्रयास करना चाहिए)।
परिणाम नकारात्मक निकला , भले ही यह बीजगणितीय रूप से शून्य होना चाहिए। संख्याओं का एक अलग सेट एक सकारात्मक संख्या या "सटीक" शून्य प्राप्त कर सकता है।
-
उदारवादी लापता होने का उदाहरण, युग्मक विलोपन के माध्यम से सकारात्मक अर्धचालकता के नुकसान की ओर जाता है:
z <- x + y + rnorm(30)/50 # same x and y as before.
xyz1 <- data.frame(x=x,y=y,z=z) # high correlation but definitely of full rank
xyz1$x[sample(1:30,5)] <- NA # make 5 x's missing
xyz1$y[sample(1:30,5)] <- NA # make 5 y's missing
xyz1$z[sample(1:30,5)] <- NA # make 5 z's missing
cov(xyz1,use="pairwise") # the individual pairwise covars are fine ...
x y z
x 1.2107760 -0.2552947 1.255868
y -0.2552947 1.2728156 1.037446
z 1.2558683 1.0374456 2.367978
chol(cov(xyz1,use="pairwise")) # ... but leave the matrix not positive semi-definite
Error in chol.default(cov(xyz1, use = "pairwise")) :
the leading minor of order 3 is not positive definite
chol(cov(xyz1,use="complete")) # but deleting even more rows leaves it PSD
x y z
x 0.8760209 -0.2253484 0.64303448
y 0.0000000 1.1088741 1.11270078
z 0.0000000 0.0000000 0.01345364