# 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 nonparametic completion method via sign series
hatTheta = SignT(Y,truer = 3,H = 3,Lmin = -3,Lmax = 3, option =1)
print(hatTheta$est)
# }
Run the code above in your browser using DataLab