Skip to content

Commit 5d5d794

Browse files
BF: forgot to rename lowlevel to base
1 parent fd0070a commit 5d5d794

File tree

1 file changed

+12
-12
lines changed

1 file changed

+12
-12
lines changed

selectiveInference/R/funs.lar.R

Lines changed: 12 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -381,18 +381,18 @@ larInf <- function(obj, sigma=NULL, alpha=0.1, k=NULL, type=c("active","all","ai
381381
vj = vj / mj # Standardize (divide by norm of vj)
382382

383383
limits.info = TG.limits(y, -Gj, -uj, vj, Sigma=diag(rep(sigma^2, n)))
384-
a = TG.pvalue.lowlevel(limits.info, bits=bits)
384+
a = TG.pvalue.base(limits.info, bits=bits)
385385
pv[j] = a$pv
386386
sxj = sx[vars[j]]
387387
vlo[j] = a$vlo * mj / sxj # Unstandardize (mult by norm of vj / sxj)
388388
vup[j] = a$vup * mj / sxj # Unstandardize (mult by norm of vj)
389389
vmat[j,] = vj * mj / sxj # Unstandardize (mult by norm of vj / sxj)
390390

391-
a = TG.interval.lowlevel(limits.info,
392-
alpha=alpha,
393-
gridrange=gridrange,
394-
flip=(sign[j]==-1),
395-
bits=bits)
391+
a = TG.interval.base(limits.info,
392+
alpha=alpha,
393+
gridrange=gridrange,
394+
flip=(sign[j]==-1),
395+
bits=bits)
396396
ci[j,] = a$int * mj / sxj # Unstandardize (mult by norm of vj / sxj)
397397
tailarea[j,] = a$tailarea
398398

@@ -439,19 +439,19 @@ larInf <- function(obj, sigma=NULL, alpha=0.1, k=NULL, type=c("active","all","ai
439439
uj = c(u,0)
440440

441441
limits.info = TG.limits(y, -Gj, -uj, vj, Sigma=diag(rep(sigma^2, n)))
442-
a = TG.pvalue.lowlevel(limits.info, bits=bits)
442+
a = TG.pvalue.base(limits.info, bits=bits)
443443

444444
pv[j] = a$pv
445445
sxj = sx[vars[j]]
446446
vlo[j] = a$vlo * mj / sxj # Unstandardize (mult by norm of vj / sxj)
447447
vup[j] = a$vup * mj / sxj # Unstandardize (mult by norm of vj / sxj)
448448
vmat[j,] = vj * mj / sxj # Unstandardize (mult by norm of vj / sxj)
449449

450-
a = TG.interval.lowlevel(limits.info,
451-
alpha=alpha,
452-
gridrange=gridrange,
453-
flip=(sign[j]==-1),
454-
bits=bits)
450+
a = TG.interval.base(limits.info,
451+
alpha=alpha,
452+
gridrange=gridrange,
453+
flip=(sign[j]==-1),
454+
bits=bits)
455455
ci[j,] = a$int * mj / sxj # Unstandardize (mult by norm of vj / sxj)
456456
tailarea[j,] = a$tailarea
457457
}

0 commit comments

Comments
 (0)