# \donttest{
#get simulation data
data(exampledata)
data=exampledata
# To get n.comp value, we can use NG_number function.
# use JB statistic as the measure of nongaussianity to run lngca with df=0
output_JB=singR(dX=exampledata$dX,dY=exampledata$dY,
df=0,rho_extent="small",distribution="JB",individual=TRUE)
est.Mj = aveM(output_JB$est.Mjx,output_JB$est.Mjy)
# }
Run the code above in your browser using DataLab