Nothing
#*********************************************************************************
#** function: prodVA product array to matrix **
#** compute: 'array * vector', which is (p,p) matrix. **
#** variables: ary: is (n*p*p) array. Diffrent from before, then change. **
#** vector: is a vector with dimension (n) **
#*********************************************************************************
prodVA <- function(ary, vector)
{
if(class(ary)=="numeric") return(sum(ary*vector))
p <- length(ary[ 1, 1,])
n <- length(ary[,1, 1 ])
prd <- matrix(rep(0, times=p * p),nrow=p)
for (k1 in 1:p)
for(k2 in 1:p)
{
t1 <- ary[,k1 ,k2 ]
t2 <- sum(vector * t1)
prd[k1,k2] <- t2
}
return(prd)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.