所有成績單的線性回歸

Question

我有一個名為 timecourse1 的 dataframe ，其中包含時間的名稱和不同成績單的行，如下所示。 timecourse1我可以知道如何一次對每個成績單進行線性回歸嗎？ 我已經嘗試了下面的代碼，但沒有成功。 真的很感激！

  data = timecourse1,
  grouping.vars = c(transcript),           
  formula = expression_level ~ time
)```

Answer 1

unique_transcripts <- unique(timecourse1$transcript)

lm_list <- list()

for(u in unique_transcripts){
  
  lm_list[u] <- tryCatch({lm(expression_level~time,data=timecourse1[timecourse1$transcript == u,])},
                         error = function(e){NA})
  
}