R Under development (unstable) (2026-02-18 r89435 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2026 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> set.seed(19970224)
> nt = 100
> n = 100
> d = 200
> p = 3
> m = d * p
> t = 1
> 
> library(splines)
> 
> X = 0.5*matrix(runif(n*d),n,d) + matrix(rep(0.5*runif(n),d),n,d)
> y = sign(((X[,1]-0.5)^2 + (X[,2]-0.5)^2)-0.06)
> ## flipping about 5 percent of y
> y = y*sign(runif(n)-0.05) 
> y = sign(y==1)
> 
> genZ_t = proc.time()
> Z = matrix(0,n,m)
> for(j in 1:d){
+   tmp = (j-1)*p + c(1:p)
+   tmp0 = ns(X[,j],df=p)
+   Z[,tmp] = tmp0
+ }
> 
> Z <-cbind(rep(1, n), Z)
> colnames(Z) <- c("Intercept", paste("X", 1:m, sep = ""))
> 
> genZ_t = proc.time() - genZ_t
> 
> Xt = 0.5*matrix(runif(nt*d),nt,d) + matrix(rep(0.5*runif(nt),d),nt,d)
> yt = sign(((Xt[,1]-0.5)^2 + (Xt[,2]-0.5)^2)-0.06)
> ## flipping about 5 percent of y
> yt = yt*sign(runif(nt)-0.05) 
> yt = sign(yt==1)
> 
> Zt = matrix(0,nt,m)
> for(j in 1:d){
+   tmp = (j-1)*p + c(1:p)
+   tmp0 = ns(Xt[,j],df=p)
+   Zt[,tmp] = tmp0
+ }
> 
> Zt <-cbind(rep(1, nt), Zt)
> colnames(Zt) <- c("Intercept", paste("X", 1:m, sep = ""))
> 
> index <- NA
> for (i in 1:d) {
+   index <- c(index, rep(i, p))
+ }
> 
> # SAM
> 
> library(SAM)
> 
> total_t = 0
> total_l = 0
> nlamb = 20
> for (i in 1:t) {
+   t0 = proc.time()
+   out.trn = samLL(X, y, nlambda=nlamb)
+   total_t = total_t + proc.time() - t0
+   out.tst = predict(out.trn, Xt)
+   total_l = total_l + mean(out.tst$labels[,nlamb]==yt)
+ }
> print("sam log-reg:")
[1] "sam log-reg:"
> print(total_t / t - genZ_t)
   user  system elapsed 
   0.44   -0.03    0.41 
> print(total_l / t)
[1] 0.8
> 
> total_t = 0
> total_l = 0
> nlamb = 20
> for (i in 1:t) {
+   t0 = proc.time()
+   out.trn = samLL(X, y, nlambda=nlamb, regfunc="MCP")
+   total_t = total_t + proc.time() - t0
+   out.tst = predict(out.trn, Xt)
+   total_l = total_l + mean(out.tst$labels[,nlamb]==yt)
+ }
> print("sam log-reg with MCP:")
[1] "sam log-reg with MCP:"
> print(total_t / t - genZ_t)
   user  system elapsed 
   1.33   -0.02    1.31 
> print(total_l / t)
[1] 0.66
> 
> 
> proc.time()
   user  system elapsed 
   2.76    0.10    2.86