-1
我正在處理一個簡單的奇異值分解10x10的矩陣,並且所有條目爲零。我想添加一個硬幣翻轉到這個矩陣並運行SVD。該代碼是下面:越界下標錯誤
# Simple Example
set.seed(1234)
par(mar=rep(0.2,4))
dataMatrix = matrix(rnorm(100), nrow=10)
dataMatrix
getMatrix = dataMatrix*0
getMatrix
svd1 = svd(getMatrix)
svd1
par(mfrow=c(1,5))
image(t(getMatrix)[,nrow(getMatrix):1], col='blue')
plot(svd1$d,xlab='Column',ylab='Singular Values', pch=19, col='red')
## Add a pattern to the data set , flipping a coin
set.seed(78910)
for (i in 1:40){
coinFlip1 = rbinom(1,size=1, prob=0.5)
coinFlip2 = rbinom(1,size=1,prob=0.5)
if (coinFlip1){
getMatrix[i,]=getMatrix[i,]+rep(c(0,5),each=5)
}
if (coinFlip2){
getMatrix[i,]=getMatrix[i,]+rep(c(0,5),each=5)
}
}
hh = hclust(dist(getMatrix))
getMatrixOrdered= getMatrix[hh$order,]
當我運行R中的程序時,它給我一個錯誤下面:
錯誤getMatrix [I,]:下標出界
的我知道界限是無序的,所以我該如何解決?
謝謝!