我刚接触R,需要使用单核苷酸多态性(SNP)列表(也是离散值0或1)对表型列表进行逻辑回归,每个表型的离散值为0或1。我已经用这两个列表做了一个矩阵,现在我想做一个单变量分析,涵盖两组中每个实体的每一个组合
phenotypes <- colnames(regressioninput)[c(2:32, 47:56)]lapply(glm.formula, function(f) glm(as.formula(f), family = binomial(link='logit'), data = d