如何在R中通过循环实现cNORM包predictNorm函数逐行处理列变量批量计算T值？

Absolutely! You can easily compute T scores for all 999 subjects by iterating over your data row-wise. Here are two straightforward approaches in R, depending on your preference:

Approach 1: For Loop (Intuitive for Beginners) #

If you prefer an explicit, easy-to-follow method, a for loop works perfectly. Assuming your data is stored in a data frame named df with columns Raw (your raw scores) and Age (your age values), here's how to do it:

# Create a new column to store the computed T scores
df$T_score <- NA

# Loop through each row in the data frame
for (i in 1:nrow(df)) {
  # Extract the raw score and age for the current subject
  current_raw <- df$Raw[i]
  current_age <- df$Age[i]
  
  # Calculate the T score using predictNorm and store it
  df$T_score[i] <- predictNorm(Raw = current_raw, Age = current_age, model = model1, minNorm = 0, maxNorm = 100)
}

# View the first few rows of your updated data frame
head(df)

Approach 2: Functional Programming with purrr (Cleaner & More Efficient) #

For a more idiomatic R approach, use the pmap_dbl function from the purrr package. This avoids explicit loops and is faster for large datasets:

First, install and load the purrr package if you haven't already:

install.packages("purrr")
library(purrr)

Then compute the T scores row-wise:

# Apply predictNorm to each pair of Raw and Age values
df$T_score <- pmap_dbl(df[, c("Raw", "Age")], function(Raw, Age) {
  predictNorm(Raw = Raw, Age = Age, model = model1, minNorm = 0, maxNorm = 100)
})

# Or use a shorter formula syntax if you prefer
df$T_score <- pmap_dbl(df[, c("Raw", "Age")], ~ predictNorm(Raw = ..1, Age = ..2, model = model1, minNorm = 0, maxNorm = 100))

Quick Check: Does predictNorm Accept Vectors? #

Before diving into iteration, it’s worth testing if predictNorm can handle vector inputs directly (sometimes documentation or initial testing can be misleading). If this works, it’s the simplest method:

df$T_score <- predictNorm(Raw = df$Raw, Age = df$Age, model = model1, minNorm = 0, maxNorm = 100)

If this throws an error, stick with the loop or purrr approaches above.

Notes #

• Ensure your model1 is properly fitted using cNORM::cnorm() before running these calculations.
• If you have missing values in Raw or Age, you might want to add checks (e.g., if (!is.na(current_raw) && !is.na(current_age))) to avoid errors during iteration.

内容的提问来源于stack exchange，提问作者I'm not a robot