Movatterモバイル変換

Title:	Impact Measurement Toolkit
Version:	1.0.0
Description:	A toolkit for causal inference in experimental and observational studies. Implements various simple Bayesian models including linear, negative binomial, and logistic regression for impact estimation. Provides functionality for randomization and checking baseline equivalence in experimental designs. The package aims to simplify the process of impact measurement for researchers and analysts across different fields. Examples and detailed usage instructions are available athttps://book.martinez.fyi.
License:	Apache License (≥ 2.0)
Encoding:	UTF-8
RoxygenNote:	7.3.2
Biarch:	true
Depends:	R (≥ 3.4.0)
Imports:	R6, bayesplot, caret, dplyr, ggplot2, glue, magrittr, methods,Rcpp, RcppParallel, rlang, rstan, rstantools, tibble, tidyr,vizdraws, purrr, scales, tidyselect
LinkingTo:	BH (≥ 1.66.0), Rcpp (≥ 0.12.0), RcppEigen (≥ 0.3.3.3.0),RcppParallel (≥ 5.0.1), rstan (≥ 2.18.1), StanHeaders (≥2.18.0)
Suggests:	knitr, rmarkdown, roxygen2, testthat
SystemRequirements:	GNU make
URL:	https://github.com/google/imt
BugReports:	https://github.com/google/imt/issues
NeedsCompilation:	yes
Packaged:	2024-08-28 13:19:01 UTC; ignacio
Author:	Ignacio Martinez [aut, cre]
Maintainer:	Ignacio Martinez <martinezig@google.com>
Repository:	CRAN
Date/Publication:	2024-09-02 21:40:10 UTC

The 'imt' package

Description

A DESCRIPTION OF THE PACKAGE

Author(s)

Maintainer: Ignacio Martinezmartinezig@google.com (ORCID)

References

Stan Development Team (NA). RStan: the R interface to Stan. R package version 2.32.6. https://mc-stan.org

See Also

Useful links:

• https://github.com/google/imt

• Report bugs athttps://github.com/google/imt/issues

Pipe operator

Description

Seemagrittr::%>% for details.

Usage

lhs %>% rhs

Arguments

Value

The result of callingrhs(lhs).

Combine and Unite Columns

Description

This function takes a data frame, identifies columns that are not specifiedin the exclusion list, and combines them into a new column called 'group'.The original columns used for the combination are then removed.Finally, it returns a data frame with only the 'group','variables', 'std_diff', and 'balance' columns.

Usage

.combineColumns(df)

Arguments

Value

A data frame with the combined 'group' column and specified columns.

Add a Random Treatment Indicator Column to a Data Frame

Description

This function takes a data frame and adds a new column named "treated" withrandomly assigned TRUE/FALSE values. Randomization can be done either onthe entire data frame or stratified by specified columns. The probability ofbeing assigned to the treatment group can be specified, with adefault of 0.5.

Usage

.randomize_internal(data, group_by = NULL, seed = NULL, pr_treated = 0.5)

Arguments

Value

A new data frame with the added "treated" column.

Create a Baseline Balance Plot.

Description

Create a Baseline Balance Plot.

Usage

balancePlot(data)

Arguments

Value

ggplot2 baseline balance plot.

Examples

library(imt)set.seed(123)  # for reproducibilityN <- 1000fake_data <- tibble::tibble(x1 = rnorm(N), x2 = rnorm(N), t = rbinom(N, 1, 0.5))baseline <- checkBaseline(data = fake_data, variables = c("x1", "x2"), treatment = "t")balancePlot(data = baseline)

Bayesian Linear Model Factory

Description

Bayesian Linear Model Factory

Bayesian Linear Model Factory

Active bindings

Methods

Public methods

• blm$new()

• blm$ppcDensOverlay()

• blm$tracePlot()

• blm$posteriorProb()

• blm$pointEstimate()

• blm$credibleInterval()

• blm$priorProb()

• blm$priorInterval()

• blm$vizdraws()

• blm$clone()

Methodnew()

Get the package version

Get the posterior draws for eta

Get the MCMC diagnostics

Get the credible interval

Create a new Bayesian Linear Model object.

Usage

blm$new(  data,  y,  x,  treatment,  eta_mean,  eta_sd,  generate_fake_data = 0,  seed = 1982,  ...)

Arguments

Returns

invisible

MethodppcDensOverlay()

This method compares the empirical distribution of the data 'y'to the distributions of simulated/replicated data 'yrep' from theposterior predictive distribution. This is done by creating a densityoverlay plot using thebayesplot::ppc_dens_overlay function.

Usage

blm$ppcDensOverlay(n = 50)

Arguments

Returns

ggplot2 visualization

MethodtracePlot()

Plot MCMC trace for the eta and sigma parameters.

Usage

blm$tracePlot(...)

Arguments

Returns

A ggplot object.

MethodposteriorProb()

Calculate posterior probability of effect exceeding a threshold

This function calculates the posterior probability of the effectbeing larger or smaller than a specified threshold.

Usage

blm$posteriorProb(threshold = 0, gt = TRUE)

Arguments

Details

This function uses the private$..eta_draws internal variableto obtain draws from the posterior distribution of the effect size.Based on the specified arguments, the function calculates theproportion of draws exceeding/falling below the threshold andreturns a formatted statement describing the estimated probability.

Calculate point estimate of the effect

This R6 method calculates the point estimate of the effect sizebased on the posterior draws of the eta parameter.

Returns

A character string summarizing the estimated probability.

MethodpointEstimate()

Usage

blm$pointEstimate(median = TRUE)

Arguments

Details

This method uses the private$..eta_draws internal variablewhich contains MCMC draws of the eta parameter representing theeffect size. Based on the specified median argument, the methodcalculates and returns either the median or the mean of the draws.Calculates credible interval for the effect of the intervention

This R6 method calculates and returns a formatted statement summarizingthe credible interval of a specified width for the effect of the intervention.

Returns

A numeric value representing the point estimate.

MethodcredibleInterval()

Usage

blm$credibleInterval(width = 0.75, round = 0)

Arguments

Details

This method uses the private$..eta_draws internal variablecontaining MCMC draws of the eta parameter representing the effectsize. It calculates the credible interval, stores it internally, andreturns a formatted statement summarizing the findings.

Calculate and format probability statement based on prior

This method calculates the probability that the effect is greater than orless than a given threshold based on the prior distribution of the effect.The results are formatted into a statement suitable for reporting.

Returns

A character string with the following information:

• The probability associated with the specified width

• The lower and upper bounds of the credible interval, rounded to thespecified number of decimal places

MethodpriorProb()

Usage

blm$priorProb(threshold = 0, gt = TRUE)

Arguments

Returns

A character string containing the formatted probability statement.Calculates the prior for the effedt of the intervention based onthe hyperpriors.

This method computes and formats a statement about the probability intervalof the effect based on the prior distribution.

MethodpriorInterval()

Usage

blm$priorInterval(width = 0.75, round = 0)

Arguments

Details

This method calculates the lower and upper bounds of the interval basedon the specified probability width and the prior distribution of the effect.It then formats the results into a clear and informative statement.

Note that the probability is checked to be within valid range (0-1) withconsideration of machine precision using .Machine$double.eps.

Returns

A character string containing the formatted probability intervalstatement.

Methodvizdraws()

Plots impact's prior and posterior distributions.

For more details seevizdraws::vizdraws().

Usage

blm$vizdraws(...)

Arguments

Returns

An interactive plot of the prior and posterior distributions.

Methodclone()

The objects of this class are cloneable with this method.

Usage

blm$clone(deep = FALSE)

Arguments

Calculate Probability of Posterior Draws Falling Within a Range

Description

This function estimates the probability that a vector of posterior draws,represented by the parametereta, falls within a specified range.It provides flexibility to use either prior distributions or posterior draws,and to specify one-sided or two-sided probability calculations.

Usage

calcProb(x, a = NULL, b = NULL, prior = FALSE, group_name = "group average")

Arguments

Details

This function checks the following cases:

• If botha andb areNULL, it returns an empty string.

• Ifb is less than or equal toa, it throws an error.

The calculated probability and range are presented in a human-readable stringusing theglue package for formatting.

Value

A formatted string stating the calculated probability and thespecified range. The probability is the proportion of samples(either prior or posterior) that fall within the defined range.

Calculate Effect Sizes for Treatment vs. Control

Description

This function calculates effect sizes for each variable in a data frame, comparing treatment and control groups. It handles continuous, binary, and categorical variables using appropriate effect size measures.

Usage

calculateEffectSizes(data, treatment_column, to_check = NULL)

Arguments

Details

- For continuous variables, Hedges' g effect size is calculated.  - For binary variables, Cox's Proportional Hazards Index (Cox's C) is calculated.  - For categorical variables, the variable is converted into multiple indicator    (dummy) variables, and the average Cox's C across these indicators is reported.

Value

A data frame with two columns:  *Variable: The name of each variable in the original data frame.  *EffectSize: The calculated effect size for each variable.

Check Baseline Equivalency.

Description

Check Baseline Equivalency.

Usage

checkBaseline(data, variables, treatment)

Arguments

Value

tibble with the standardized difference of the pre-interventionvariables. The tibble includesvariables: the variable name,std_diff: the standardized difference for that variable as a numberbalance: the standardized difference for that variable as a factorvariable.For more details about this methodology checkhttps://ies.ed.gov/ncee/wwc/Docs/OnlineTraining/wwc_training_m3.pdf.

Cleans and prepares data for analysis

Description

This function performs a series of data cleaning and preprocessing stepsto ensure the data is suitable for analysis. This includes:

• Missing data handling

• Variable type checks

• Collinearity and zero-variance feature removal

Usage

cleanData(data, y, treatment, x = NULL, binary = FALSE)

Arguments

Value

A list containing the cleaned dataset and relevant metadata:

• N: The number of observations after cleaning.

• K The number of covariates after cleaning.

• X The cleaned covariate matrix.

• treat_vec: Treatment vector as integers (1 for TRUE, 0 for FALSE).

• Y: The dependent variable vector.

Count missing values (NA) in a dataframe

Description

This function takes a dataframe as input and returns a tibblesummarizing the number of missing values (NA) in each column andthe number of rows with at least one missing value.

Usage

countMissing(df)

Arguments

Value

A tibble with the following columns:

• NA_<column_name>: Number of NAs in each original column

• missing_rows: Number of rows with at least one NA across all columns

Cox's Proportional Hazards Index (Cox's C)

Description

Calculates Cox's C, a standardized effect size measure for comparing hazard ratesbetween two groups in survival analysis.

Usage

coxsIndex(p_t, p_c, n_t, n_c)

Arguments

Details

Cox's C is a useful effect size for survival analysis when hazard ratios are not constant over time.It's calculated based on the log odds ratio of events and includes a small sample size correction.The value 1.65 is used to approximate a conversion to a Cohen's d-like scale.

Value

The calculated Cox's C effect size.

References

Cox, D. R. (1972). Regression models and life-tables. Journal of the Royal Statistical Society: Series B (Methodological), 34(2), 187-202.

Converts a dataframe into a named list to provide data to a Stan model

Description

Converts a dataframe into a named list to provide data to a Stan model

Usage

createData(  data,  y,  treatment,  x = NULL,  eta_mean = 0,  eta_sd = 1,  run_estimation = 1)

Arguments

Value

Returns
stan_data a named list providing the data for the stan model

Calculate credible interval from MCMC draws

Description

This function calculates a credible interval of the specifiedwidth from a vector of MCMC draws.

Usage

credibleInterval(draws, width)

Arguments

Details

The function calculates the lower and upper bounds of thecredible interval using the quantile function based on thespecified width.

Value

A named list containing three elements:

• width: The specified width of the credible interval.

• lower_bound: The lower bound of the credible interval.

• upper_bound: The upper bound of the credible interval.

Examples

# Generate example drawsdraws <- rnorm(1000)# Calculate 95% credible intervalcredibleInterval(draws, width = 0.95)

Fits Stan model.

Description

Fits Stan model.

Usage

fit(stan_data, model = "blm", ...)

Arguments

Value

The complete StanFit object from the fitted Stan model.

Extracts parameter from Stan model.

Description

Extracts parameter from Stan model.

Usage

getStanParameter(fit, par)

Arguments

Value

A long tibble with draws for the pamete.

Hedges' g Effect Size with Pooled Standard Deviation

Description

Calculates Hedges' g, a standardized effect size for comparing means.This version includes a small-sample correction factor (omega)and uses the pooled standard deviation.

Usage

hedgesG(n_t, n_c, y_t, y_c, s_t, s_c)

Arguments

Details

Hedges' g is a variation of Cohen's d that adjusts for small-sample bias.It is calculated as the difference in means divided by the pooled standarddeviation, then multiplied by a correction factor.

Value

The calculated Hedges' g effect size.

References

Hedges, L. V. (1981). Distribution theory for Glass's estimator of effect size and related estimators.Journal of Educational Statistics, 6(2), 107-128.

Bayesian Logit Model Factory

Description

A class for creating and managing Bayesian Logit Models

Active bindings

Methods

Public methods

• logit$new()

• logit$tracePlot()

• logit$calcProb()

• logit$pointEstimate()

• logit$credibleInterval()

• logit$vizdraws()

• logit$lollipop()

• logit$plotPrior()

• logit$predict()

• logit$getPred()

• logit$predSummary()

• logit$predCompare()

• logit$clone()

Methodnew()

Get the package version

Get the posterior draws for tau

Get the MCMC diagnostics

Get the credible interval

Get the prior for eta

Get the prior for tau

Get the prior for mean y

Get the posterior draws for eta

Get the list of predictions

Create a new Bayesian Logit Model object.

Usage

logit$new(  data,  y,  x,  treatment,  mean_alpha,  sd_alpha,  mean_beta,  sd_beta,  tau_mean,  tau_sd,  seed = 1982,  fit = TRUE,  ...)

Arguments

Returns

invisible

MethodtracePlot()

Plot MCMC trace for the eta and sigma parameters.

Usage

logit$tracePlot(...)

Arguments

Returns

A ggplot object.

MethodcalcProb()

Calculates the posterior of an effect being greater than, less than,or within a range defined by thresholds.

Usage

logit$calcProb(a = 0, b = NULL, prior = FALSE)

Arguments

Returns

A character string summarizing the estimated probability

Calculate point estimate of the effect

This R6 method calculates the point estimate of the effect sizebased on the posterior draws of the eta parameter.

MethodpointEstimate()

Usage

logit$pointEstimate(median = TRUE)

Arguments

Details

This method uses the private$..eta_draws internal variablewhich contains MCMC draws of the eta parameter representing theeffect size. Based on the specified median argument, the methodcalculates and returns either the median or the mean of the draws.Calculates credible interval for the effect of the intervention

This R6 method calculates and returns a formatted statement summarizingthe credible interval of a specified width for the effect of the intervention.

Returns

A numeric value representing the point estimate.

MethodcredibleInterval()

Usage

logit$credibleInterval(width = 0.75, round = 0)

Arguments

Details

This method uses the private$..eta_draws internal variablecontaining MCMC draws of the eta parameter representing the effectsize. It calculates the credible interval, stores it internally, andreturns a formatted statement summarizing the findings.

Returns

A character string with the following information:

• The probability associated with the specified width

• The lower and upper bounds of the credible interval, rounded to thespecified number of decimal places

Methodvizdraws()

Plots impact's prior and posterior distributions.

Usage

logit$vizdraws(tau = FALSE, ...)

Arguments

Returns

An interactive plot of the prior and posterior distributions.

Methodlollipop()

Plots lollipop chart for the prior and posterior of the impact beinggreater or less than a threshold.

For more details seevizdraws::lollipops().

Usage

logit$lollipop(threshold = 0, ...)

Arguments

Returns

A lollipop chart with the prior and posterior probability ofthe impact being above or below a threshold.

MethodplotPrior()

Plots draws from the prior distribution of the outcome, tau, and impactin percentage points.

Usage

logit$plotPrior()

Methodpredict()

Predict new data

Usage

logit$predict(new_data, name = NULL, M = NULL, ...)

Arguments

Returns

invisible(self)

MethodgetPred()

Get posterior predictive draws

Usage

logit$getPred(name = NULL, ...)

Arguments

Returns

Matrix of posterior predictive draws

MethodpredSummary()

Get point estimate, credible interval and prob summary of predictive draws

Usage

logit$predSummary(  name = NULL,  subgroup = NULL,  median = TRUE,  width = 0.75,  round = 0,  a = NULL,  b = NULL,  ...)

Arguments

Returns

A character string with summary information

MethodpredCompare()

Compare the average of the posterior draws of two groups

Usage

logit$predCompare(  name1,  name2,  subgroup1 = NULL,  subgroup2 = NULL,  median = TRUE,  width = 0.75,  round = 0,  a = NULL,  b = NULL,  ...)

Arguments

Returns

A character string with comparison summary

Methodclone()

The objects of this class are cloneable with this method.

Usage

logit$clone(deep = FALSE)

Arguments

Calculate logit link and sample from binomial distribution

Description

This function takes prior / posterior draws of the Bayesian logit model,applies theinverse logit (logistic) transformation to obtain probabilities,and then generates random samples from binomial distributions.

Usage

logitRng(alpha, tau, beta, treat, X, N)

Arguments

Value

A vector of N samples of preditive draws

MCMC Checks

Description

Checks convergence, mixing, effective sample size, and divergent transitions

Methods

$new(fit, pars)

Runs diagnostics on the suppliedstanfitobject, restricted to parameters identified by the character vectorpars.
Tests include:
Share of specified parameters with an Rhat less than 1.1. If anyhave an Rhat > 1.1,everything_looks_fine is set toFALSE.
Share of specified parameters with an n_eff at least 0.1% of thetotal number of posterior draws. If any have n_eff < 0.001 * N,everything_looks_fine is set toFALSE.
Share of specified parameters with an n_eff of at least 100. If anyhave n_eff < 100,everything_looks_fine is set toFALSE.
Number of divergent transitions during posterior sampling. If thereare any whatsoever,everything_looks_fine is set toFALSE.
Share of posterior iterations where the sampler reached themaximum treedepth. If more than 25\everything_looks_fine is set toFALSE.

Active bindings

Methods

Public methods

• mcmcChecks$new()

• mcmcChecks$clone()

Methodnew()

Initialize a new mcmcChecks object and run diagnostics

Usage

mcmcChecks$new(fit, pars)

Arguments

Methodclone()

The objects of this class are cloneable with this method.

Usage

mcmcChecks$clone(deep = FALSE)

Arguments

Create a Meta-Analysis Object Using Data From Previous Studies

Description

Create a Meta-Analysis Object Using Data From Previous Studies

Create a Meta-Analysis Object Using Data From Previous Studies

Details

A meta analysis has raw data and draws from the lift's posteriordistribution. This is represented by an R6 Class.

Active bindings

Methods

Public methods

• metaAnalysis$new()

• metaAnalysis$PlotRawData()

• metaAnalysis$PlotLift()

• metaAnalysis$UpdateCI()

• metaAnalysis$probability()

• metaAnalysis$findings()

• metaAnalysis$clone()

Methodnew()

Create a new meta analysis object.

Usage

metaAnalysis$new(  data,  point_estimates,  standard_errors,  id,  mean_mu = 0,  sd_mu = 0.05,  ci_width = 0.75,  X = NULL,  run_estimation = 1,  ...)

Arguments

Returns

A newmeta_analysis object.

MethodPlotRawData()

Plots the raw data.

Usage

metaAnalysis$PlotRawData()

Returns

A plot with point estimates and 95% confidence intervals.

MethodPlotLift()

Plots lift's prior and posterior distributions.

For more details seevizdraws::vizdraws().

Usage

metaAnalysis$PlotLift(...)

Arguments

Returns

An interactive plot of the prior and posterior distributions.

MethodUpdateCI()

Update the width of the credible interval.

Usage

metaAnalysis$UpdateCI(ci_width)

Arguments

Methodprobability()

Calculates that probability that lift is between a and b.

Usage

metaAnalysis$probability(a = -Inf, b = Inf, percent = TRUE)

Arguments

Returns

A string with the probability.

Methodfindings()

Calculates the point estimate a credible interval for the meta analysis.

Usage

metaAnalysis$findings(percent = TRUE)

Arguments

Returns

A string with the findings

Methodclone()

The objects of this class are cloneable with this method.

Usage

metaAnalysis$clone(deep = FALSE)

Arguments

Bayesian Negative Binomial Model Factory

Description

Bayesian Negative Binomial Model Factory

Bayesian Negative Binomial Model Factory

Active bindings

Methods

Public methods

• negativeBinomial$new()

• negativeBinomial$tracePlot()

• negativeBinomial$posteriorProb()

• negativeBinomial$pointEstimate()

• negativeBinomial$credibleInterval()

• negativeBinomial$vizdraws()

• negativeBinomial$clone()

Methodnew()

Create a new Bayesian Negative Binomial Model object.

Usage

negativeBinomial$new(  data,  y,  x,  treatment,  tau_mean,  tau_sd,  run_estimation = 1,  seed = 1982,  ...)

Arguments

Returns

invisible

MethodtracePlot()

Plot MCMC trace for the eta and sigma parameters.

Usage

negativeBinomial$tracePlot(...)

Arguments

Returns

A ggplot object.

MethodposteriorProb()

Calculate posterior probability of effect exceeding a threshold

This function calculates the posterior probability of the effectbeing larger or smaller than a specified threshold.

Usage

negativeBinomial$posteriorProb(threshold = 0, gt = TRUE)

Arguments

Details

This function uses the private$..tau_draws internal variableto obtain draws from the posterior distribution of the effect size.Based on the specified arguments, the function calculates theproportion of draws exceeding/falling below the threshold andreturns a formatted statement describing the estimated probability.

Calculate point estimate of the effect

This R6 method calculates the point estimate of the effect sizebased on the posterior draws of the tau parameter.

Returns

A character string summarizing the estimated probability.

MethodpointEstimate()

Usage

negativeBinomial$pointEstimate(median = TRUE)

Arguments

Details

This method uses the private$..tau_draws internal variablewhich contains MCMC draws of the tau parameter representing theeffect size. Based on the specified median argument, the methodcalculates and returns either the median or the mean of the draws.Calculates credible interval for the effect of the intervention

This R6 method calculates and returns a formatted statement summarizingthe credible interval of a specified width for the effect of the intervention.

Returns

A numeric value representing the point estimate.

MethodcredibleInterval()

Usage

negativeBinomial$credibleInterval(width = 0.75, round = 0)

Arguments

Details

This method uses the private$..tau_draws internal variablecontaining MCMC draws of the tau parameter representing the effectsize. It calculates the credible interval, stores it internally, andreturns a formatted statement summarizing the findings.

Returns

A character string with the following information:

• The probability associated with the specified width

• The lower and upper bounds of the credible interval, rounded to thespecified number of decimal places

Methodvizdraws()

Plots impact's prior and posterior distributions.

For more details seevizdraws::vizdraws().

Usage

negativeBinomial$vizdraws(...)

Arguments

Returns

An interactive plot of the prior and posterior distributions.

Methodclone()

The objects of this class are cloneable with this method.

Usage

negativeBinomial$clone(deep = FALSE)

Arguments

Calculate Point Estimate (Median or Mean) as Percentage

Description

This function computes a point estimate from a numeric vector,returning either the median or the mean as a percentage.

Usage

pointEstimate(x, median = TRUE)

Arguments

Details

This function provides a simple way to obtain either the median or meanof a numeric vector as a percentage. The choice between these two measuresof central tendency can be controlled by themedian argument.

Value

A numeric value representing the chosen point estimate (median or mean)of the input vectorx, multiplied by 100 to express it as a percentage.

Randomly Assign Treatment While Controlling for Baseline Equivalency

Description

This function repeatedly randomizes treatmentassignment (using.randomizer) untilbaseline equivalency is achieved across specifiedvariables, as measured by thecheckBaseline function fromtheim package. It can optionally stratify the randomization byspecified groups.

Usage

randomize(  data,  variables,  standard = "Not Concerned",  seed = NULL,  max_attempts = 100,  pr_treated = 0.5,  group_by = NULL)

Arguments

Value

A new data frame with the added "treated"column, if baseline equivalency is achieved within the specified number ofattempts. Otherwise, an error is thrown.

Randomization Class for Treatment Assignment

Description

This class provides methods to randomly assign treatmentsto a dataset while ensuring baseline covariate balance. It canhandle both simple and stratified randomization.

Active bindings

Methods

Public methods

• randomizer$new()

• randomizer$clone()

Methodnew()

Initialize a new Randomizer object.

Usage

randomizer$new(  data,  variables,  standard = "Not Concerned",  seed = NULL,  max_attempts = 100,  group_by = NULL)

Arguments

Returns

A newrandomizer object.

Methodclone()

The objects of this class are cloneable with this method.

Usage

randomizer$clone(deep = FALSE)

Arguments

Validate a Logical Subgroup Vector

Description

This function checks if a given vector is a logical vector (TRUE/FALSE)and whether its length matches the number of rows in a specified matrix.It is designed to validate subgroup vectors used for subsettin data.

Usage

validate_logical_vector(subgroup, N, name = NULL)

Arguments

Details

This function performs two key validations:

1. Checks if thesubgroup vector is logical.

2. Checks if the length of thesubgroup vector matches the N.

Value

The originalsubgroup vector if it passes all validation checks.