Nothing
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, n), 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:")
print(total_t / t - genZ_t)
print(total_l / t)
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:")
print(total_t / t - genZ_t)
print(total_l / t)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.