test_function {DetR} | R Documentation |
Functions to test the cpp
codes in the package.
test_function()
This is a series of R
functions that, together, implement the c++
codes used
in the package and which can be used to test those.
Vakili Kaveh.
n<-100 p<-5 #set.seed(123) #for repoducibility. Z<-matrix(rnorm(n*(p+1)),nc=p+1) x<-Z[,1:p] y<-Z[,p+1] datao<-cbind(x,y) alpha<-0.6; test_R_0<-DetR:::test_fxOGK(x0=x,y0=y,cent_est='scaleTau2_test',scal_est='scaleTau2_test', alpha=alpha) h<-DetR:::quanf(alpha,n=n,p=p+1) #intercept=1 test_cpp<-DetR:::fxOGK(Data=datao,scale_est="scaleTau2",intercept=1,h=h,doCsteps=1) ####should be the same sort(test_cpp$bestRaw) sort(as.numeric(test_R_0$bestRaw)) ############# test_R_1<-DetR:::test_Cstep(x=x,y=y,h=h,z0=test_R_0$bestRaw) ####should be the same sort(test_R_1$bestCStep) sort(test_cpp$bestCStep[1:h]) #################################### n<-100 p<-5 set.seed(123) #for repoducibility. Z<-matrix(rnorm(n*(p+1)),nc=p+1) x<-Z[,1:p] y<-Z[,p+1] datao<-cbind(x,y) alpha<-0.6; test_R_0<-DetR:::test_fxOGK(x0=x,y0=y,cent_est='median',scal_est='qn', alpha=alpha) h<-DetR:::quanf(alpha,n=n,p=p+1) #intercept=1 test_cpp<-DetR:::fxOGK(Data=datao,scale_est="qn",intercept=1,h=h,doCsteps=1) ####should be the same sort(test_cpp$bestRaw) sort(as.numeric(test_R_0$bestRaw)) ############# test_R_1<-DetR:::test_Cstep(x=x,y=y,h=h,z0=test_R_0$bestRaw) ####should be the same sort(test_R_1$bestCStep) sort(test_cpp$bestCStep[1:h])