У меня есть одна непрерывная переменная с 76 наблюдениями и еще 29877 непрерывными переменными, против которых я хотел бы провести линейную регрессию.Код для одного может выглядеть примерно так:
lm(data=spread_gene_exp, spread_gene_exp[,i] ~ spread_gene_exp[,2])
, где i - одна из 29877 переменных Y.
Затем я хотел бы извлечь оценку коэффициента и P-значения из результатов и добавить их к некоторому виду данных, где я могу легко определить, какие переменные были наиболее значимыми.Для этого, вероятно, потребуется какой-то цикл, но я не уверен, с чего начать.
Я бы с удовольствием запустил это на своем рабочем столе, но это, скорее всего, займет довольно много времени, и поэтому, вероятно, придетсяЗапустите его на университетском сервере, и любая оценка того, сколько времени потребуется на среднечастотный / полужесткий рабочий стол, будет полезна.