मैं सिर्फ निम्नलिखित पेपर में किए गए दावे को दोहराने की कोशिश कर रहा हूं , जो जीन एक्सप्रेशन डेटा से संबंधित सहसंबंधी बिकलस्टर्स ढूंढ रहा है , जो है:
प्रस्ताव 4. अगर । तो हमारे पास हैं:
मैं। यदि एडिटिव मॉडल के साथ एक सही बाइक्लेस्टर है, तो कॉलम पर सहसंबंध के साथ एक सही बाइक्लेस्टर है; ii। यदि एडिटिव मॉडल के साथ एक सही बाइक्लेस्टर है, तो पंक्तियों पर सहसंबंध के साथ एक सही बाइक्लेस्टर है; iii। यदि और दोनों एडिटिव मॉडल के साथ परफेक्ट हैं, तो एक परफेक्ट कॉरेलिटेड बाइक्लेस्टर है।
इन प्रस्तावों को आसानी से साबित किया जा सकता है ...
... लेकिन निश्चित रूप से, वे इसे साबित नहीं करते हैं।
मैं पेपर प्लस आधार + कस्टम आर कोड में कुछ सरल उदाहरणों का उपयोग कर रहा हूं यह देखने के लिए कि क्या मैं इस प्रस्ताव को प्रदर्शित कर सकता हूं।
corbic <- matrix(c(0,4,-4,2,2,-2,6,0,4,-8,16,-2,-2,10,-14,4), ncol=4)
(तालिका 1 एफ से)
कुछ कस्टम कोड मानक X = फॉर्म को में बदलने के लिए जैसा कि कागज में वर्णित है:
svdToRC <- function(x, ignoreRank = FALSE, r = length(x$d), zerothresh=1e-9) {
#convert standard SVD decomposed matrices UEV' to RC' form
#x -> output of svd(M)
#r -> rank of matrix (defaults to length of singular values vector)
# but really is the number of non-zero singular values
#ignoreRank -> return the full decomposition (ignore zero singular values)
#zerothresh -> how small is zero?
R <- with(x, t(t(u) * sqrt(d)))
C <- with(x, t(t(v) * sqrt(d)))
if (!ignoreRank) {
ind <- which(x$d >= zerothresh)
} else {
ind <- 1:r
}
return(list(R=as.matrix(R[,ind]), C=as.matrix(C[,ind])))
}
इस फ़ंक्शन को डेटासेट में लागू करें:
> svdToRC(svd(corbic))
$R
[,1] [,2]
[1,] 0.8727254 -0.9497284
[2,] -2.5789775 -1.1784221
[3,] 4.3244283 -0.7210346
[4,] -0.8531261 -1.0640752
$C
[,1] [,2]
[1,] -1.092343 -1.0037767
[2,] 1.223860 -0.9812343
[3,] 3.540063 -0.9586919
[4,] -3.408546 -1.0263191
जब तक मैं मतिभ्रम नहीं कर रहा हूं, तब तक यह मैट्रिस एडिटिव नहीं है, भले ही कॉर्बिक पंक्तियों और स्तंभों के बीच सही सहसंबंध प्रदर्शित करता है। यह अजीब लगता है कि वे जो उदाहरण प्रदान करते हैं वह उस संपत्ति को प्रदर्शित करता है जो उन्होंने कहा था कि ... जब तक मैं किसी प्रकार के पूर्व या बाद के परिवर्तन को नहीं भूल रहा हूँ - परिवर्तन परिवर्तन कदम?
4iii
कहता P(R), P(C), additivity => P(X)
? (मैं संक्षिप्त कर रहा हूं " Y
जैसा कि एक आदर्श बाइक्लस्टर है" P(Y)
)। ऐसा लगता है कि आप दूसरी दिशा में जा रहे हैं, उम्मीद है कि अन्य स्थितियों से सकारात्मकता बढ़ेगी। कृपया अधिक बताएं।
4iii
यह नहीं कहता है कि अगर X
पूरी तरह से सहसंबंधित बाइक्लेस्टर है R
और C
वह additive होगा। निहितार्थ दूसरी दिशा में जाता है। अब, मैं मानता हूं कि यह अजीब है कि वे जो उदाहरण देते हैं, वह उन प्रमेयों के साथ नहीं लगता है जो इसके बगल में हैं। शायद कुछ अन्य जानकारी है जो आप प्रदान कर सकते हैं? क्या कोई और प्रमेय है जो दूसरी दिशा में जाता है?