Question

我长期遇到“矩阵上的下标数不正确”的问题，请帮帮我。我试图找到错误但无法找到，我的数据集很大，但我使用较小的初始数据进行程序运行。我尽我所能，但一切都是徒劳的。错误出现在最后一行代码中。我是编码的新手，所以请原谅我在缩进中的错误。

########### inputs
#Y <- as.matrix(read.table("data2.dat"));
#X <- as.matrix(read.table("data1.dat"));
###########
# smaller X and Y for initial code run.
X <- matrix(c(71,22,53,14,75,68,74,93,72),9,1)
Y <- matrix(c(7,1,4,1,7,6,6,8,5),9,1)



########### Constant
Xlen <- length(X[,1])  #### rows of X
Yhat <- as.matrix(mat.or.vec(Xlen,1))
error <- as.matrix(mat.or.vec(Xlen,1))
Rmax <- 30
Z <- as.matrix(cbind(X,Y)) 
Zlen <- length(Z[1,]) #### columns of Z
P <- as.matrix(mat.or.vec(1,Rmax))
U <- as.matrix(mat.or.vec(1,Rmax)) # mu calculation
lambda <- as.matrix(mat.or.vec(1,Rmax)) # lambda calculation
C <- as.matrix(mat.or.vec(Zlen,Rmax))
R <- 0
sigma <- 0
beta <- as.matrix(mat.or.vec(1,Zlen))
Zp <- mat.or.vec(1,Zlen)
V1 <- 0
V2 <- 0
V3 <- 0
Pz <- 0
Dmin <-0
Dind <-0
Pmax <-0
Pind <-0
r <- 0.4
alpha <- 4/(r^2)
P_mat <- 10^4 * diag(2)
gamma <- as.matrix(mat.or.vec(2,Rmax))
w <- as.matrix(mat.or.vec(2,Rmax))
Xe <- as.matrix(cbind(1,X))
##### 


for (i in 1:Xlen)
{   
if (i==1) 
{R <- R+1 ; C[,R] <- Z[i,] ; P[R] <- 1; Yhat <- Y[i,]}
else 

{V1 <- sum(Z[i,]^2)
sigma <- sigma + sum(Z[i-1,]^2)
beta <- beta + Z[i-1,]
V2 <- sum(Z[i,]*t(beta))
Pz <- (i-1)/((i-1)*(V1+1)+(sigma-(2*V2)))
for (j in 1:R)
{V3 <- sum((Z[i,]-Z[i-1,])^2);P[j] <- ((i-1)*P[j])/((i-2)+(P[j])+(P[j]*V3))}
dist <- as.matrix(mat.or.vec(1,R))

for (k in 1:R)
{dist[k] <- abs(sqrt(sum((Z[i,]-C[,k])^2)))}


Dmin <- min(dist)
Dind <- which.min(dist)


Pmax <- max(P)
Pind <- which.max(P)


if (r > Dmin/(1-(Pz/Pmax)) && Pz > Pmax)
{C[,Dind] <- Z[i,]; P[Dind] <- Pz}

if (r < (Dmin/(1-(Pz/Pmax))) && Pz > Pmax && R< Rmax)
{R <- R+1; C[,R] <- Z[i,]; P[R] <- Pz}

for(n in 1:R)
{U[,n] <- exp(-alpha*(X[i]-C[1,n])^2)}

for(n in 1:R)
{lambda[,n] <- U[,n]/sum(U)}


A1 <- as.matrix(P_mat%*%Xe[i,]%*%lambda)
B1 <- as.matrix(t(lambda)%*%Xe[i,])
C1 <- as.matrix(diag(Rmax))
D1 <- as.matrix(P_mat)
E1 <- as.matrix(B1%*%D1%*%t(B1))
F1 <- C1+E1
gamma <- A1%*%solve(F1)
P_mat <- (diag(2)%*%gamma%*%t(lambda)%*%Xe[i,])%*%P_mat
G1 <- as.matrix(lambda%*%t(Xe[i,]%*%w))
H1 <- as.matrix(Y[i,])
I1 <- as.vector(H1-G1)
J1 <- I1*gamma
w <- w + J1; #as.matrix(I1%*%gamma)
K1 <- Xe[i,]%*%w
L1 <- lambda%*%t(K1)

这是Yhat的错误

Yhat[i,] <- L1  
#error[i,] <- Y[i,]-Yhat[i,]
#b <- b+1
}
}
#Yhat <- as.vector(Yhat)
#plot(Y,type="l",col="red")
#lines(Yhat,col="green")
#plot(Yhat)
#plot(error)

Answer 1

如果不通过所有代码，让我们看看Yhat会发生什么：

首先，你把它变成一个矩阵：

Yhat <- as.matrix(mat.or.vec(Xlen,1))
> Yhat
      [,1]
 [1,]    0
 [2,]    0
 [3,]    0
 [4,]    0
 [5,]    0
 [6,]    0
 [7,]    0
 [8,]    0
 [9,]    0

在循环for (i in 1:Xlen)的第一次迭代中，您将Y[i,]分配给Yhat：

if (i==1)
{R <- R+1 ; C[,R] <- Z[i,] ; P[R] <- 1; Yhat <- Y[i,]}

此时，Y[i,]为7。 Yhat不再是矩阵，现在它的类型为numeric。因此，Yhat[i,] <- L1会抛出错误。

我收到“矩阵上的下标数不正确”的问题？

1 个答案: