Question

我有一个for循环，它适用于（相对）少量重复次数= 10,100。但是对于“时间”的较大值，我通过填充矩阵得到错误：下标超出界限...（参见下面的代码和代码说明）

错误： Error in M_zp_var[j, (1:n)] : subscript out of bounds

要填充矩阵，我在第一个中使用第二个for循环

M_zp<-matrix(numeric(1),B,N)
for(j in 1:dim(M_zp)[1]){ M_zp[j,]<-Z_pi(x,y) }

我甚至尝试过

M_zp<-t(replicate(B, Z_pi(x,y)))

相反，但它也不起作用。

正如我所说，如果我的外循环很小，代码就可以工作。我不明白为什么如果我选择变量“倍数”（1000,5000），它的工作方式会有所不同。

我知道“下标越界”的含义

我希望你能帮助我，我会非常棒！：（

代码：

T_p<- matrix(numeric(1),times,Bvar)
  H0 <- numeric(times)
  for(i in 1:times){
    # Each loop starts with a new random vector
    x<-rnorm(n)   #x<- rdist(n,...); #
y<-rnorm(m)   #y<- rdist(m,...); #

# Order statstic with T_pi's T(Z(pi))
n<-length(x)
m<-length(y)
N<-n+m

#Permutate (x,y) B times
**M_zp<-matrix(numeric(1),B,N)
for(j in 1:dim(M_zp)[1]){ M_zp[j,]<-Z_pi(x,y) }** #see below for Z_pi() funct.
#M_zp<-t(replicate(B, Z_pi(x,y)))
M_zp_var<-unique(M_zp)
Bf<-dim(unique(M_zp))[1]

#Test statistic computation for each one of the Permutations in M_zp_var
T_pvec<-numeric(Bvar)
for(j in 1:Bvar){
  xp<-M_zp_var[j,(1:n)]
  yp<-M_zp_var[j,((n+1):(n+m))]
  m_xp<- mstern(xp) 
  m_yp<- mstern(yp) 
  T_pvec[j]<-(sd(xp)-sd(yp))/sqrt(m_xp/n+m_yp/m)
}
T_p[i,]<-sort(T_pvec)
}


Z_pi<-function(x,y){
  n<-length(x)
  m<-length(y)
  x<-sort(x)
  y<-sort(y)
  wicy<-function(s){ wi<-which(y==s); return(wi)}
  wicx<-function(s){ wi<-which(x==s); return(wi)}
  r<-sample(1:min(n,m),1)   #interger beween 1 and min(m+n) = length of entries to interchange between x and y
  zwy<-sample(y,r)          #vector of length r out of entries in y (randomly) =intries to interchange with x
  wy<-unlist(lapply(zwy,wicy))
  zwx<-sample(x,r)
  wx<-unlist(lapply(zwx,wicx))

  yp<-y
  xp<-x
  yp[wy]<-zwx
  xp[wx]<-zwy
  zp<-c(sort(xp), sort(yp))
  return(zp)
}

说明：

循环模拟“时间” - 排列测试（蒙特卡罗）。一个循环生成2个随机样本，然后随机置换它们B = 16000次（用函数Z_pi()进行置换），然后它只需要不同的permuatations（unique()）并且每个置换它计算一个测试统计信息......代码停止在标记的行（粗体）中工作，该行涉及函数Z_pi()。

对于大型for循环，下标超出界限R.

0 个答案: