Run Searchlight Analysis

Execute a searchlight analysis using multivariate pattern analysis.

Usage

run_searchlight(
  model_spec,
  radius,
  method = c("standard", "randomized"),
  niter = NULL,
  ...
)

Arguments

model_spec

A mvpa_model instance containing the model specifications

radius

The searchlight radius in millimeters

method

The type of searchlight, either 'randomized' or 'standard'

niter

The number of searchlight iterations (used only for 'randomized' method)

...

Extra arguments passed to specific searchlight methods

Value

A named list of NeuroVol objects containing performance metrics (e.g., AUC) at each voxel location

Examples

# \donttest{
  # Generate sample dataset with categorical response
  dataset <- gen_sample_dataset(
    D = c(8,8,8),           # 8x8x8 volume
    nobs = 100,             # 100 observations
    response_type = "categorical",
    data_mode = "image",
    blocks = 3,             # 3 blocks for cross-validation
    nlevels = 2             # binary classification
  )
  
  # Create cross-validation specification using blocks
  cval <- blocked_cross_validation(dataset$design$block_var)
  
  # Load the SDA classifier (Shrinkage Discriminant Analysis)
  model <- load_model("sda_notune")
  
  # Create MVPA model
  mspec <- mvpa_model(
    model = model,
    dataset = dataset$dataset,
    design = dataset$design,
    model_type = "classification",
    crossval = cval
  )
  
  # Run searchlight analysis
  results <- run_searchlight(
    mspec,
    radius = 8,            # 8mm radius
    method = "standard"    # Use standard searchlight
  )
#> function (model_spec, good_results, bad_results = NULL)  
#> [1] "combiner is "
#> INFO [2025-03-20 21:29:24] Running standard searchlight with radius = 8
#> INFO [2025-03-20 21:29:24] creating standard searchlight
#> INFO [2025-03-20 21:29:24] running standard searchlight iterator
#> INFO [2025-03-20 21:29:24] ⚡ Processing batch 1/11
#> INFO [2025-03-20 21:29:31] ⚡ Processing batch 2/11
#> INFO [2025-03-20 21:29:38] ⚡ Processing batch 3/11
#> INFO [2025-03-20 21:29:44] ⚡ Processing batch 4/11
#> INFO [2025-03-20 21:29:50] ⚡ Processing batch 5/11
#> INFO [2025-03-20 21:29:56] ⚡ Processing batch 6/11
#> INFO [2025-03-20 21:30:03] ⚡ Processing batch 7/11
#> INFO [2025-03-20 21:30:09] ⚡ Processing batch 8/11
#> INFO [2025-03-20 21:30:15] ⚡ Processing batch 9/11
#> INFO [2025-03-20 21:30:21] ⚡ Processing batch 10/11
#> INFO [2025-03-20 21:30:27] ⚡ Processing batch 11/11
#> INFO [2025-03-20 21:30:32] 
#> ✨ MVPA Iteration Complete
#> ├─ Total ROIs: 254
#> ├─ Processed: 254
#> └─ Skipped: 0
  
  # Results contain performance metrics
  # Access them with results$performance
# }