# NOT RUN {
library(rTensor)
indices = c(2,3,4)
noise = rand_tensor(indices)@data
Theta = array(runif(prod(indices),min=-3,max = 3),indices)
# The signal plus noise model
Y = Theta + noise
# Estimate Theta from CP low rank tensor method
hatTheta = fit_continuous(Y,3)
print(hatTheta$est)
# }
Run the code above in your browser using DataLab