R语言中列表上的向量迭代：RNA数据蒙特卡洛结果分组问题

Hey there! Sounds like you're deep into RNA feature extraction with Monte Carlo simulations—let's get this grouping sorted out for you. The core idea here is to create a group identifier that maps each value in your result vector back to its original sequence, then use that to split your data.

Let's break this down with concrete examples (using R, since your code snippet uses n <- 10) #

First, let's define our variables to mirror your scenario:

• n_sequences: Number of RNA sequences (2 in your case)
• n_iterations: Monte Carlo iterations per sequence (10 in your case)
• result_vector: Your combined 20-value vector from all iterations

Method 1: Base R (Quick & Simple) #

This is perfect if you want a list where each element holds the iterations for one sequence.

1. Create grouping IDs: Repeat each sequence's identifier (e.g., 1, 2) exactly n_iterations times
2. Split the vector using these IDs

# Simulate your result vector (replace this with your actual data)
result_vector <- rnorm(20)  # 2 sequences × 10 iterations = 20 values

# Define your parameters
n_sequences <- 2
n_iterations <- 10

# Generate group IDs: repeat each sequence ID n_iterations times
group_ids <- rep(1:n_sequences, each = n_iterations)

# Split the result vector into groups
split_results <- split(result_vector, group_ids)

# Access results for each sequence
split_results[[1]]  # Iterations for sequence 1
split_results[[2]]  # Iterations for sequence 2

If your sequences have custom names (like "RNA_seq_001" instead of 1/2), just swap the sequence IDs with your names:

seq_names <- c("RNA_seq_001", "RNA_seq_002")
group_ids <- rep(seq_names, each = n_iterations)
split_results <- split(result_vector, group_ids)

# Now you can access by sequence name directly
split_results$RNA_seq_001

Method 2: Tidyverse (For Data Frame Workflows) #

If you prefer working with data frames (great for downstream analysis like calculating per-sequence stats), use dplyr and tidyr:

library(tidyverse)

# Simulate your result vector
result_vector <- rnorm(20)
n_sequences <- 2
n_iterations <- 10

# Build a structured data frame
result_df <- tibble(
  iteration_value = result_vector,
  sequence_id = rep(1:n_sequences, each = n_iterations),
  iteration_num = rep(1:n_iterations, times = n_sequences)  # Track iteration number if needed
)

# Group by sequence and analyze (example: calculate mean/sd per sequence)
sequence_stats <- result_df %>%
  group_by(sequence_id) %>%
  summarize(
    mean_iter = mean(iteration_value),
    sd_iter = sd(iteration_value),
    all_iterations = list(iteration_value)  # Store all iterations as a list column
  )

# View the result
print(sequence_stats)

This approach keeps your data organized and makes it easy to run further statistical analyses on each sequence's Monte Carlo results.

内容的提问来源于stack exchange，提问作者Emilio Mármol Sánchez