Newb нужна помощь в редактировании R скрипт-консультант написал

Question

Я редактирую скрипт, который написал мой начальник для анализа данных.Я смог запустить скрипт в R, но он добавил новый набор данных, который я пытаюсь заставить код работать правильно (используя два разных набора данных, не объединяя их).Я имею минимальное воздействие R, и я не уверен, как интерпретировать полученное сообщение об ошибке.Любая помощь приветствуется!

Я пробовал редактировать различные части скрипта, особенно ту, которая, как мне кажется, вызывает проблему.

# read data
downs <- read.csv(filename, header = TRUE)
dim(downs)
head(downs)
names(downs)

mri <- read.csv(mri_filename, header=TRUE)

######

lm_age_gender = lm(as.formula(paste0(interesting_rois[rois],"~","downs$Age_at_scan + downs$gender + downs$brain.average + downs$Site + mri$mri_rois[rois]")),na.action=na.omit)

Я ожидаю совсем немного от этого сценария, включая результаты линейной модели, которые дают средние значения, значения P, стандартные ошибки и т. Д., Но я получаю сообщение «traceback»:

7: eval(expr, envir, enclos)

6: eval(predvars, data, env)

5: model.frame.default(formula = as.formula(paste0(interesting_rois[rois], 
       "~", "downs$Age_at_scan + downs$gender + downs$brain.average + downs$Site + mri$mri_rois[rois]")), 
       na.action = na.omit, drop.unused.levels = TRUE)

4: stats::model.frame(formula = as.formula(paste0(interesting_rois[rois], 
       "~", "downs$Age_at_scan + downs$gender + downs$brain.average + downs$Site + mri$mri_rois[rois]")), 
       na.action = na.omit, drop.unused.levels = TRUE)

3: eval(expr, envir, enclos)

2: eval(mf, parent.frame())

1: lm(as.formula(paste0(interesting_rois[rois], "~", "downs$Age_at_scan + downs$gender + downs$brain.average + downs$Site + mri$mri_rois[rois]")), 
       na.action = na.omit) at #28

Строка 28 - это функция lm, которую я предоставил выше, и я попытался найти в строке 2 и 3 и прочитал ответы от людей на этом сайте, которые сталкивались с этими сообщениями об ошибках, но я не могу подключиться кчто означает сообщение об ошибке в этом скрипте.

ЛЮБАЯ проницательность была бы полезна, поскольку я чувствую, что перепробовал много вещей и продолжаю бить стену.