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值，标准误差等的线性模型结果，但是我收到一条“回溯”消息：

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函数，并且我尝试使用Google搜索2和3行，并阅读了该网站上遇到这些错误消息的人员的回复，但是我无法以某种方式连接含义该脚本的错误消息。

任何见解都会有所帮助，因为我觉得我已经尝试了很多事情并不断碰壁。

Newb需要帮助编辑R脚本顾问写的

0 个答案: