R使用内置函数避免循环

Question

我在R语言工作。

我有两个数组，分别表示为E和y。 E被指数为E [状态]，y被指数为y [时间，样本]。因此，E是1维，y是2维。

我必须计算以下3维数组，表示为P.

P（状态，时间，样本）= f（E（状态），y [时间，样本]），其中f是某个函数。我想在不使用循环的情况下这样做。

示例代码如下。

# EM algorithm for HMM  
#y is array containing data
  T #length of time series
  nsamples #number of sequences observed
  nstates # no. of states in HMM
  P_forward<-array(0,dim=c(nstates,T,nsamples)) #forward probabilities in 
                                         #EM algo. for HMM
P_forward[state,time,sample] = 1/abs(E[state]-y[time,sample])