我认为这会产生您正在寻找的东西。 param_df 将是一个数据框,其中包含每对分布的 N、均值和 sd。
# Generate some fake parameters
param_df <- data.frame(n_1 = sample(50:100, 10, replace = T),
mean_1 = sample(5:10, 10, replace = T),
sd_1 = sample(1:9, 10, replace = T),
n_2 = sample(50:100, 10, replace = T),
mean_2 = sample(5:10, 10, replace = T),
sd_2 = sample(1:9, 10, replace = T),
variable_name = paste0("variable_", 1:10),
stringsAsFactors = F)
# Function that will generate the distributions and calculate statistics
my_function <- function(n_1, mean_1, sd_1, n_2, mean_2, sd_2, var_name){
samp_A <- rnorm(as.numeric(n_1), as.numeric(mean_1), as.numeric(sd_1))
samp_B <- rnorm(as.numeric(n_2), as.numeric(mean_2), as.numeric(sd_2))
var_score <- (mean(samp_A) + mean(samp_B)) / 2
var_range <- sd(samp_A)
return(c(Variable = var_name, Score = var_score, Range = var_range))
}
# Apply the function to the param_df by row, passing each column as an argument to the function
df <- apply(param_df, 1, function(params) my_function(params[1],
params[2],
params[3],
params[4],
params[5],
params[6],
params[7]))
# Convert to a dataframe
df <- as.data.frame(t(df))
请注意,此处的代码计算“分数”的方式与您上面的代码不同。我假设您想除以平均值的总和,然后将 A 的平均值和 B 的平均值相加。