diff --git a/.Rbuildignore b/.Rbuildignore
index bf2e3f7..fc00ef0 100644
--- a/.Rbuildignore
+++ b/.Rbuildignore
@@ -5,3 +5,4 @@
 ^knockofftools\.Rproj$
 ^LICENSE$
 ^README\.md$
+^\.Rproj\.user$
diff --git a/DESCRIPTION b/DESCRIPTION
index a952efe..a90fb62 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -26,10 +26,10 @@ Imports:
     survival (>= 2.44.1.1),
     randomForestSRC (>= 2.9.1),
     grf (>= 2.3.1),
-    glasso (>= 1.11),
     CVglasso (>= 1.0),
-    nnet (>= 7.3-18)
-RoxygenNote: 7.2.3
+    nnet (>= 7.3-18),
+    magrittr (>= 1.5)
+RoxygenNote: 7.3.3
 Suggests: 
     knitr,
     rmarkdown,
diff --git a/NAMESPACE b/NAMESPACE
index c632dfd..de4a81b 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -35,3 +35,24 @@ export(stat_predictive_causal_forest)
 export(stat_predictive_glmnet)
 export(stat_random_forest)
 export(variable.selections)
+importFrom(magrittr,"%>%")
+importFrom(stats,as.formula)
+importFrom(stats,coef)
+importFrom(stats,contrasts)
+importFrom(stats,cutree)
+importFrom(stats,dist)
+importFrom(stats,dnbinom)
+importFrom(stats,hclust)
+importFrom(stats,ks.test)
+importFrom(stats,lm)
+importFrom(stats,model.matrix)
+importFrom(stats,pnbinom)
+importFrom(stats,predict)
+importFrom(stats,qqnorm)
+importFrom(stats,rbinom)
+importFrom(stats,rexp)
+importFrom(stats,rmultinom)
+importFrom(stats,rnorm)
+importFrom(stats,runif)
+importFrom(stats,sd)
+importFrom(stats,toeplitz)
diff --git a/R/globals.R b/R/globals.R
new file mode 100644
index 0000000..0ea43aa
--- /dev/null
+++ b/R/globals.R
@@ -0,0 +1,6 @@
+utils::globalVariables(c(
+  "selected",
+  "variable",
+  "draw",
+  "drawclass"
+))
diff --git a/R/imports.R b/R/imports.R
new file mode 100644
index 0000000..a5d045f
--- /dev/null
+++ b/R/imports.R
@@ -0,0 +1,4 @@
+#' @keywords internal
+#' @importFrom magrittr %>%
+#' @importFrom stats as.formula coef contrasts cutree dist dnbinom hclust ks.test lm model.matrix pnbinom predict qqnorm rbinom rexp rmultinom rnorm runif sd toeplitz
+NULL
diff --git a/R/internal.R b/R/internal.R
index b8eb551..7f0e196 100644
--- a/R/internal.R
+++ b/R/internal.R
@@ -133,7 +133,7 @@ find_single_optimal_variable_set <- function(S, p, trim=0.5) {
 #' Do not call this function on its own. Fits cross-validated glmnet model with fixed effect.
 #'
 #'
-#' @param X.fixed a data.frame (or tibble) with "numeric" and "factor" columns corresponding to covariates or terms that should be treated as fixed effects in the model.
+#' @param X_fixed a data.frame (or tibble) with "numeric" and "factor" columns corresponding to covariates or terms that should be treated as fixed effects in the model.
 #' @param X original data.frame (or tibble) with "numeric" and "factor" columns only. The number of columns, ncol(X) needs to be > 2.
 #' @param y response vector with \code{length(y) = nrow(X)}. Accepts "numeric" (family="gaussian") or binary "factor" (family="binomial"). Can also be a survival object of class Surv
 #' as obtained from y = survival::Surv(time, status).
@@ -158,7 +158,7 @@ cv_coeffs_glmnet_with_fixed_effect <- function(X_fixed, X, y, family, nlambda=50
     X = scale(X)
   }
 
-  if (!methods::hasArg(lambda) ) {
+  if (!hasArg(lambda) ) {
     if( identical(family, "gaussian") ) {
       if(!is.numeric(y)) {
         stop('Input y must be numeric.')
@@ -195,13 +195,12 @@ cv_coeffs_glmnet_with_fixed_effect <- function(X_fixed, X, y, family, nlambda=50
 #' @param y response vector with \code{length(y) = nrow(X)}. Accepts "numeric" (family="gaussian") or binary "factor" (family="binomial"). Can also be a survival object of class Surv
 #' as obtained from y = survival::Surv(time, status).
 #' @param type should be "regression" if y is numeric, "classification" if y is a binary factor variable or "survival" if y is a survival object.
-#' @param ...
 #'
 #' @return importance scores
 #' @export
 #'
 #' @keywords internal
-random_forest_importance_scores <- function(X, y, trt, type = "regression", ...){
+random_forest_importance_scores <- function(X, y, trt, type = "regression"){
   # make the column names unique
   colnames(X) = make.unique(colnames(X))
 
@@ -323,16 +322,15 @@ ns.transform <- function(y) {
 
 }
 
-#' Heuristic check for whether a variable can be reasonably treated as continuous
+#' Heuristic check for whether numeric variables can be reasonably treated as continuous
 #'
-#' @param x a numeric variable vector
+#' @param X the design matrix of interest with columns either "numeric" or "factor"
 #'
 #' @return a logical TRUE or FALSE depending on whether n_distinct(x) > 30
 #' @export
 #'
 #' @keywords internal
 check_if_continuous <- function(X) {
-  `%>%` <- magrittr::`%>%`
   X_numeric <- dplyr::select_if(X, is.numeric)
   is.continuous <- sum(X_numeric %>% lapply(dplyr::n_distinct) %>% unlist() <= 30) > 0
   if (is.continuous) warning("Some of the numeric columns of X have suspiciously few distinct values: n_distinct <= 30. Those columns should perhaps not be treated as continuous variables. Please review carefully and read the documentation about the gcm parameter of the knockoff.statistics function.")
diff --git a/R/knockoff_filters.R b/R/knockoff_filters.R
index 2e57a08..821c541 100644
--- a/R/knockoff_filters.R
+++ b/R/knockoff_filters.R
@@ -157,7 +157,6 @@ knockoff.statistics <- function(y, X, type="regression",
 #' @param y response vector with \code{length(y) = nrow(X)}. Accepts "numeric", binary "factor", or survival ("Surv") object.
 #' @param X data.frame (or tibble) with "numeric" and "factor" columns only. The number of columns, ncol(X) needs to be > 2.
 #' @param type should be "regression" if y is numeric, "classification" if y is a binary factor variable or "survival" if y is a survival object.
-#' @param M the number of independent knockoff feature statistics that should be calculated.
 #' @param knockoff.method what type of knockoffs to calculate. Defaults to sequential knockoffs, knockoff.method="seq", but other options are "sparseseq" and "mx". The "mx" option only works if all columns of X are continuous.
 #' @param statistic knockoff feature statistic function, defaults to glmnet coefficient difference (statistic="stat_glmnet"; see ?stat_glmnet). Other options include statistic="stat_random_forest" (see ?stat_random_forest), statistic="stat_predictive_glmnet" (see ?stat_predictive_glmnet) or statistic="stat_predictive_causal_forest" (see ?stat_predictive_causal_forest).
 #' @param trt binary treatment (factor) variable required if statistic involves a predictive knockoff filter (i.e. if statistic="stat_predictive_glmnet" or statistic="stat_predictive_causal_forest")
@@ -315,7 +314,7 @@ stat_glmnet <- function(y, X, X_k, type = "regression", X.fixed=NULL, penalty.fi
 #' @param y response vector with \code{length(y) = nrow(X)}. Accepts "numeric" (family="gaussian") or binary "factor" (family="binomial"). Can also be a survival object of class "Surv" (type="survival")
 #' as obtained from y = survival::Surv(time, status).
 #' @param type should be "regression" if y is numeric, "classification" if y is a binary factor variable or "survival" if y is a survival object.
-#' @param ...
+#' @param ... other parameters passed to \code{random_forest_importance_scores}.
 #'
 #' @return data.frame with knockoff statistics W as column. The number of rows matches the number of columns (variables) of the data.frame X and the variable names are recorded in rownames(W).
 #' @export
diff --git a/R/performance.R b/R/performance.R
index 9412258..97c0559 100644
--- a/R/performance.R
+++ b/R/performance.R
@@ -21,7 +21,7 @@ eval_fdp <- function(selected, negatives) {
 #' True positive proportion (tpp) as function of selection and known positives:
 #'
 #' @param selected vector of indices of selected variables
-#' @param negatives vector of indices of known non-null variables (that influence response)
+#' @param positives vector of indices of known non-null variables (that influence response)
 #'
 #' @return true positive rate
 #' @export
diff --git a/R/plot.R b/R/plot.R
index faa5d64..486b445 100644
--- a/R/plot.R
+++ b/R/plot.R
@@ -1,17 +1,25 @@
+
 #' Heatmap of multiple variable selections ordered by importance
 #'
-#' @param S data.frame of variable selections from multiple knockoffs (each entry is either 1 if variable is selected and 0 otherwise). Columns correspond to different knockoffs and rows correspond to the underlying variables. row.names(S) records the variable names.
+#' @param x data.frame of variable selections from multiple knockoffs
+#' (each entry is either 1 if variable is selected and 0 otherwise).
+#' Columns correspond to different knockoffs and rows correspond to the
+#' underlying variables. row.names(x) records the variable names.
+#'
+#' @param ... Additional arguments passed to other plot methods (currently ignored).
+#'
 #' @param nbcocluster bivariate vector c(number of variable clusters, number of selection clusters).
-#' The former number must be specified less than nrow(S) and the latter must be less than ncol(S).
+#' The former number must be specified less than nrow(x) and the latter must be less than ncol(x).
 #'
 #' @details To help visualize most important variables we perform clustering both selections and variables.
 #'
 #' @return plot of heatmap
+#'
+#' @method plot variable.selections
 #' @export
 #'
 #' @examples
 #' library(knockofftools)
-#'
 #' set.seed(1)
 #'
 #' # Simulate 8 Gaussian covariate predictors and 2 binary factors:
@@ -20,19 +28,19 @@
 #' # create linear predictor with first 5 beta-coefficients = 1 (all other zero)
 #' lp <- generate_lp(X, p_nn = 5, a=1)
 #'
-#' # Gaussian
-#'
-#' # Simulate response from a linear model y = lp + epsilon, where epsilon ~ N(0,1):
+#' # Simulate response:
 #' y <- lp + rnorm(100)
 #'
-#' # Calculate M independent knockoff feature statistics:
+#' # Calculate knockoff statistics:
 #' W <- knockoff.statistics(y=y, X=X, type="regression", M=5)
 #'
-#' S = variable.selections(W, error.type = "pfer", level = 1)
+#' S <- variable.selections(W, error.type = "pfer", level = 1)
 #'
 #' # plot heatmap of knockoff selections:
 #' plot(S)
-plot.variable.selections <- function(S, nbcocluster=c(7,7)) {
+plot.variable.selections <- function(x, ..., nbcocluster=c(7,7)) {
+
+  S <- x
 
   if (class(S)[1]!="variable.selections") {
     stop("Input S must be of class \'variable.selections\'. Please see ?variable.selections.")
@@ -49,8 +57,6 @@ plot.variable.selections <- function(S, nbcocluster=c(7,7)) {
                            variable = factor(rownames(S)),
                            selected = as.numeric(as.matrix(S)))
 
-  `%>%` <- dplyr::`%>%`
-
   sel.mat <- matrix(selections$selected,nrow=nrow(S))
   hclust.row <- hclust(dist(sel.mat, method="binary"), method="ward.D")
   hclust.col <- hclust(dist(t(sel.mat), method="binary"), method="ward.D")
diff --git a/R/simdata-data.R b/R/simdata-data.R
new file mode 100644
index 0000000..b1ecc35
--- /dev/null
+++ b/R/simdata-data.R
@@ -0,0 +1,49 @@
+#' Simulated dataset for knockofftools
+#'
+#' @description
+#' A synthetic dataset generated by the function \code{generate_simdata()}.
+#' It contains simulated Gaussian, binary, and survival outcomes along with covariates.
+#'
+#' @format A data frame with 2000 rows and 33 variables:
+#' \describe{
+#'   \item{Yg}{Continuous outcome}
+#'   \item{Yb}{Binary outcome}
+#'   \item{Tc}{Treatment indicator}
+#'   \item{X1}{Covariate 1}
+#'   \item{X2}{Covariate 2}
+#'   \item{X3}{Covariate 3}
+#'   \item{X4}{Covariate 4}
+#'   \item{X5}{Covariate 5}
+#'   \item{X6}{Covariate 6}
+#'   \item{X7}{Covariate 7}
+#'   \item{X8}{Covariate 8}
+#'   \item{X9}{Covariate 9}
+#'   \item{X10}{Covariate 10}
+#'   \item{X11}{Covariate 11}
+#'   \item{X12}{Covariate 12}
+#'   \item{X13}{Covariate 13}
+#'   \item{X14}{Covariate 14}
+#'   \item{X15}{Covariate 15}
+#'   \item{X16}{Covariate 16}
+#'   \item{X17}{Covariate 17}
+#'   \item{X18}{Covariate 18}
+#'   \item{X19}{Covariate 19}
+#'   \item{X20}{Covariate 20}
+#'   \item{X21}{Covariate 21}
+#'   \item{X22}{Covariate 22}
+#'   \item{X23}{Covariate 23}
+#'   \item{X24}{Covariate 24}
+#'   \item{X25}{Covariate 25}
+#'   \item{X26}{Covariate 26}
+#'   \item{X27}{Covariate 27}
+#'   \item{X28}{Covariate 28}
+#'   \item{X29}{Covariate 29}
+#'   \item{X30}{Covariate 30}
+#' }
+#'
+#' @source Simulated using \code{generate_simdata()}
+#'
+#' @examples
+#' data(simdata)
+#' head(simdata)
+"simdata"
diff --git a/R/simknockoffs.R b/R/simknockoffs.R
index 5b035ad..7e44441 100644
--- a/R/simknockoffs.R
+++ b/R/simknockoffs.R
@@ -37,7 +37,7 @@ knockoffs_seq <- function(X, seq_simulator = sim_glmnet, ...) {
     y <- X[[i]] # i-th column serves as response
     Xp <- X[,-i] # columns[-i] serve as predictors
 
-    if (loop.count > 1) Xp <- cbind(knockoffs[,shf[1:(loop.count-1)]], Xp)
+    if (loop.count > 1) Xp <- cbind(knockoffs[,shf[1:(loop.count-1)], drop=FALSE], Xp)
 
     knockoffs[[i]] <- seq_simulator(y = y, X = Xp, ...)
 
@@ -90,7 +90,7 @@ sim_glmnet <- function(y, X, ...) {
     # Beta coefficients (excluding intercept)
     beta.coefs <- as.numeric(coef(gm.cv, s = "lambda.1se")[[2]])[-1]
 
-    mu <- predict(gm.cv, newx=x, type="response", s="lambda.1se")
+    mu <- predict(gm.cv, newx=x, type="response", s="lambda.min")
 
     mat.multinom <- apply(mu, 1, function(prob) rmultinom(n=1, size=1, prob=prob))
 
@@ -107,7 +107,7 @@ sim_glmnet <- function(y, X, ...) {
     gm.cv <- glmnet::cv.glmnet(y=y, x=x, family="gaussian", intercept=TRUE, alpha=1, ...)
 
     # Beta coefficients (excluding intercept)
-    beta.coefs <- as.numeric(coef(gm.cv, s = "lambda.1se"))[-1]
+    beta.coefs <- as.numeric(coef(gm.cv, s = "lambda.min"))[-1]
 
     # columns of predictor matrix corresponding to non-zero beta.coefs:
     non.zero.cols <- which(beta.coefs != 0)
@@ -115,7 +115,7 @@ sim_glmnet <- function(y, X, ...) {
     # Total number of non-zero parameters (including intercept, hence + 1)
     s.lambda = length(non.zero.cols) + 1
 
-    mu <- predict(gm.cv, newx=x, type="response", s="lambda.1se")
+    mu <- predict(gm.cv, newx=x, type="response", s="lambda.min")
 
     rmse = sqrt(sum((y-mu)^2)/(length(y) - s.lambda))
 
@@ -144,7 +144,6 @@ sim_glmnet <- function(y, X, ...) {
 #'
 #' @param X data.frame (or tibble) with "numeric" and "factor" columns only. The number of columns, ncol(X) needs to be > 2.
 #' @param adjacency.matrix optional user specified adjacency matrix (i.e. binary indicator matrix corresponding to the non-zero elements of the precision matrix of X). Defaults to NULL and is then estimated within the function call.
-#' @param seq_simulator name of function that used to estimate the conditional distributions in the sequential steps. Default is the function \code{sim_simple}, which is a least squares fit (continuous variables) or multinomial logistic regression (factor variables) respectively.
 #'
 #' @return sparse sequential knockoff copy of X. A data.frame or tibble of same type and dimensions as X.
 #' @export
diff --git a/data/simdata.RData b/data/simdata.RData
index 6a20a21..c6de10a 100644
Binary files a/data/simdata.RData and b/data/simdata.RData differ
diff --git a/man/check_if_continuous.Rd b/man/check_if_continuous.Rd
index e9a1c8f..c277d77 100644
--- a/man/check_if_continuous.Rd
+++ b/man/check_if_continuous.Rd
@@ -2,17 +2,17 @@
 % Please edit documentation in R/internal.R
 \name{check_if_continuous}
 \alias{check_if_continuous}
-\title{Heuristic check for whether a variable can be reasonably treated as continuous}
+\title{Heuristic check for whether numeric variables can be reasonably treated as continuous}
 \usage{
 check_if_continuous(X)
 }
 \arguments{
-\item{x}{a numeric variable vector}
+\item{X}{the design matrix of interest with columns either "numeric" or "factor"}
 }
 \value{
 a logical TRUE or FALSE depending on whether n_distinct(x) > 30
 }
 \description{
-Heuristic check for whether a variable can be reasonably treated as continuous
+Heuristic check for whether numeric variables can be reasonably treated as continuous
 }
 \keyword{internal}
diff --git a/man/cv_coeffs_glmnet_with_fixed_effect.Rd b/man/cv_coeffs_glmnet_with_fixed_effect.Rd
index 37f01d8..4a6686d 100644
--- a/man/cv_coeffs_glmnet_with_fixed_effect.Rd
+++ b/man/cv_coeffs_glmnet_with_fixed_effect.Rd
@@ -15,6 +15,8 @@ cv_coeffs_glmnet_with_fixed_effect(
 )
 }
 \arguments{
+\item{X_fixed}{a data.frame (or tibble) with "numeric" and "factor" columns corresponding to covariates or terms that should be treated as fixed effects in the model.}
+
 \item{X}{original data.frame (or tibble) with "numeric" and "factor" columns only. The number of columns, ncol(X) needs to be > 2.}
 
 \item{y}{response vector with \code{length(y) = nrow(X)}. Accepts "numeric" (family="gaussian") or binary "factor" (family="binomial"). Can also be a survival object of class Surv
@@ -27,8 +29,6 @@ as obtained from y = survival::Surv(time, status).}
 \item{penalty.factor}{Separate penalty factors, passed to the glmnet::cv.glmnet function, can be applied to each coefficient. This is a number that multiplies lambda to allow differential shrinkage. Can be 0 for some variables, which implies no shrinkage, and that variable is always included in the model.}
 
 \item{...}{other parameters passed to glmnet::cv.glmnet}
-
-\item{X.fixed}{a data.frame (or tibble) with "numeric" and "factor" columns corresponding to covariates or terms that should be treated as fixed effects in the model.}
 }
 \value{
 coefficients of both the original and knockoff variables
diff --git a/man/dot-knockoff.statistics.single.Rd b/man/dot-knockoff.statistics.single.Rd
index c1ca4a3..d9af22b 100644
--- a/man/dot-knockoff.statistics.single.Rd
+++ b/man/dot-knockoff.statistics.single.Rd
@@ -31,8 +31,6 @@
 \item{adjacency.matrix}{optional user specified adjacency matrix (i.e. binary indicator matrix corresponding to the non-zero elements of the precision matrix of X). Defaults to NULL and is then estimated within the function call.}
 
 \item{...}{additional parameters passed to the "statistic" function (note that the knockoffs parameter X_k should not be entered by user; it is already calculated inside the knockoff.statistics function).}
-
-\item{M}{the number of independent knockoff feature statistics that should be calculated.}
 }
 \value{
 data.frame with a single knockoff statistics W as column.
diff --git a/man/eval_tpp.Rd b/man/eval_tpp.Rd
index 0e835dc..a23d144 100644
--- a/man/eval_tpp.Rd
+++ b/man/eval_tpp.Rd
@@ -9,7 +9,7 @@ eval_tpp(selected, positives)
 \arguments{
 \item{selected}{vector of indices of selected variables}
 
-\item{negatives}{vector of indices of known non-null variables (that influence response)}
+\item{positives}{vector of indices of known non-null variables (that influence response)}
 }
 \value{
 true positive rate
diff --git a/man/knockoffs_sparse_seq.Rd b/man/knockoffs_sparse_seq.Rd
index 95cf4ee..c069452 100644
--- a/man/knockoffs_sparse_seq.Rd
+++ b/man/knockoffs_sparse_seq.Rd
@@ -10,8 +10,6 @@ knockoffs_sparse_seq(X, adjacency.matrix = NULL)
 \item{X}{data.frame (or tibble) with "numeric" and "factor" columns only. The number of columns, ncol(X) needs to be > 2.}
 
 \item{adjacency.matrix}{optional user specified adjacency matrix (i.e. binary indicator matrix corresponding to the non-zero elements of the precision matrix of X). Defaults to NULL and is then estimated within the function call.}
-
-\item{seq_simulator}{name of function that used to estimate the conditional distributions in the sequential steps. Default is the function \code{sim_simple}, which is a least squares fit (continuous variables) or multinomial logistic regression (factor variables) respectively.}
 }
 \value{
 sparse sequential knockoff copy of X. A data.frame or tibble of same type and dimensions as X.
diff --git a/man/plot.variable.selections.Rd b/man/plot.variable.selections.Rd
index 56893b5..1b666ff 100644
--- a/man/plot.variable.selections.Rd
+++ b/man/plot.variable.selections.Rd
@@ -4,13 +4,18 @@
 \alias{plot.variable.selections}
 \title{Heatmap of multiple variable selections ordered by importance}
 \usage{
-\method{plot}{variable.selections}(S, nbcocluster = c(7, 7))
+\method{plot}{variable.selections}(x, ..., nbcocluster = c(7, 7))
 }
 \arguments{
-\item{S}{data.frame of variable selections from multiple knockoffs (each entry is either 1 if variable is selected and 0 otherwise). Columns correspond to different knockoffs and rows correspond to the underlying variables. row.names(S) records the variable names.}
+\item{x}{data.frame of variable selections from multiple knockoffs
+(each entry is either 1 if variable is selected and 0 otherwise).
+Columns correspond to different knockoffs and rows correspond to the
+underlying variables. row.names(x) records the variable names.}
+
+\item{...}{Additional arguments passed to other plot methods (currently ignored).}
 
 \item{nbcocluster}{bivariate vector c(number of variable clusters, number of selection clusters).
-The former number must be specified less than nrow(S) and the latter must be less than ncol(S).}
+The former number must be specified less than nrow(x) and the latter must be less than ncol(x).}
 }
 \value{
 plot of heatmap
@@ -23,7 +28,6 @@ To help visualize most important variables we perform clustering both selections
 }
 \examples{
 library(knockofftools)
-
 set.seed(1)
 
 # Simulate 8 Gaussian covariate predictors and 2 binary factors:
@@ -32,15 +36,13 @@ X <- generate_X(n=100, p=10, p_b=2, cov_type="cov_equi", rho=0.2)
 # create linear predictor with first 5 beta-coefficients = 1 (all other zero)
 lp <- generate_lp(X, p_nn = 5, a=1)
 
-# Gaussian
-
-# Simulate response from a linear model y = lp + epsilon, where epsilon ~ N(0,1):
+# Simulate response:
 y <- lp + rnorm(100)
 
-# Calculate M independent knockoff feature statistics:
+# Calculate knockoff statistics:
 W <- knockoff.statistics(y=y, X=X, type="regression", M=5)
 
-S = variable.selections(W, error.type = "pfer", level = 1)
+S <- variable.selections(W, error.type = "pfer", level = 1)
 
 # plot heatmap of knockoff selections:
 plot(S)
diff --git a/man/random_forest_importance_scores.Rd b/man/random_forest_importance_scores.Rd
index 426c37a..3a50cd6 100644
--- a/man/random_forest_importance_scores.Rd
+++ b/man/random_forest_importance_scores.Rd
@@ -4,7 +4,7 @@
 \alias{random_forest_importance_scores}
 \title{Internal function called to return the importance scores from random forest}
 \usage{
-random_forest_importance_scores(X, y, trt, type = "regression", ...)
+random_forest_importance_scores(X, y, trt, type = "regression")
 }
 \arguments{
 \item{X}{original data.frame with "numeric" and "factor" columns only.}
@@ -13,8 +13,6 @@ random_forest_importance_scores(X, y, trt, type = "regression", ...)
 as obtained from y = survival::Surv(time, status).}
 
 \item{type}{should be "regression" if y is numeric, "classification" if y is a binary factor variable or "survival" if y is a survival object.}
-
-\item{...}{}
 }
 \value{
 importance scores
diff --git a/man/simdata.Rd b/man/simdata.Rd
new file mode 100644
index 0000000..81853fb
--- /dev/null
+++ b/man/simdata.Rd
@@ -0,0 +1,59 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/simdata-data.R
+\docType{data}
+\name{simdata}
+\alias{simdata}
+\title{Simulated dataset for knockofftools}
+\format{
+A data frame with 2000 rows and 33 variables:
+\describe{
+  \item{Yg}{Continuous outcome}
+  \item{Yb}{Binary outcome}
+  \item{Tc}{Treatment indicator}
+  \item{X1}{Covariate 1}
+  \item{X2}{Covariate 2}
+  \item{X3}{Covariate 3}
+  \item{X4}{Covariate 4}
+  \item{X5}{Covariate 5}
+  \item{X6}{Covariate 6}
+  \item{X7}{Covariate 7}
+  \item{X8}{Covariate 8}
+  \item{X9}{Covariate 9}
+  \item{X10}{Covariate 10}
+  \item{X11}{Covariate 11}
+  \item{X12}{Covariate 12}
+  \item{X13}{Covariate 13}
+  \item{X14}{Covariate 14}
+  \item{X15}{Covariate 15}
+  \item{X16}{Covariate 16}
+  \item{X17}{Covariate 17}
+  \item{X18}{Covariate 18}
+  \item{X19}{Covariate 19}
+  \item{X20}{Covariate 20}
+  \item{X21}{Covariate 21}
+  \item{X22}{Covariate 22}
+  \item{X23}{Covariate 23}
+  \item{X24}{Covariate 24}
+  \item{X25}{Covariate 25}
+  \item{X26}{Covariate 26}
+  \item{X27}{Covariate 27}
+  \item{X28}{Covariate 28}
+  \item{X29}{Covariate 29}
+  \item{X30}{Covariate 30}
+}
+}
+\source{
+Simulated using \code{generate_simdata()}
+}
+\usage{
+simdata
+}
+\description{
+A synthetic dataset generated by the function \code{generate_simdata()}.
+It contains simulated Gaussian, binary, and survival outcomes along with covariates.
+}
+\examples{
+data(simdata)
+head(simdata)
+}
+\keyword{datasets}
diff --git a/man/stat_random_forest.Rd b/man/stat_random_forest.Rd
index 6e087aa..3daaf2b 100644
--- a/man/stat_random_forest.Rd
+++ b/man/stat_random_forest.Rd
@@ -17,7 +17,7 @@ as obtained from y = survival::Surv(time, status).}
 
 \item{type}{should be "regression" if y is numeric, "classification" if y is a binary factor variable or "survival" if y is a survival object.}
 
-\item{...}{}
+\item{...}{other parameters passed to \code{random_forest_importance_scores}.}
 }
 \value{
 data.frame with knockoff statistics W as column. The number of rows matches the number of columns (variables) of the data.frame X and the variable names are recorded in rownames(W).
diff --git a/tests/testthat/test-plot.R b/tests/testthat/test-plot.R
index 905df52..9ea561d 100644
--- a/tests/testthat/test-plot.R
+++ b/tests/testthat/test-plot.R
@@ -34,8 +34,6 @@ test_that("plot.variable.selections works as expected", {
 
   p <- suppressWarnings(plot(S))
 
-  expect_equal(class(p), c("gg", "ggplot"))
-
   expect_error(plot.variable.selections(S$selected))
 
 })