我使用这些变量对运行回归:
regression_pairs=c("A~B", "C~D", "E~F", "G~H","I~J","K~L","M~N","O~P","Q~R")
regression_pairs
并将残基存储在列表residuals.new
中:
residuals.new <- list()
a=c(1,2,3,4,5)
b=c(6,7,8,9,10)
c=c(11,12,13,14,15)
d=c(16, 17, 18, 19, 20)
e=c(21,22,23,24,25)
f=c(26,27,28,29,30)
g=c(31,32,33,34,35)
h=c(36,37,38,39,40)
i=c(41,42,43,44,45)
residuals.new[[1]]=a
residuals.new[[2]]=b
residuals.new[[3]]=c
residuals.new[[4]]=d
residuals.new[[5]]=e
residuals.new[[6]]=f
residuals.new[[7]]=g
residuals.new[[8]]=h
residuals.new[[9]]=i
residuals.new
因此,回归A~B
的残差为1,2,3,4,5
,回归C〜D的残差为6,7,8,9,10
,依此类推。
经过一些标准,我得到了最好的配对:
bestresults<-list()
best_pairs=c("A~B", "G~H", "Q~R")
现在,我需要找到一种创建列表best.residuals
的方式,该列表将最佳对与自己的残差相关联。
我可以使用什么功能?有任何建议吗?
所以我的输出列表将是:
the.bestresiduals<-list()
the.bestresiduals[[1]]=c(1,2,3,4,5)
the.bestresiduals[[2]]=c(16, 17, 18, 19, 20)
the.bestresiduals[[3]]=c(41,42,43,44,45)
the.bestresiduals
[[1]]
[1] 1 2 3 4 5
[[2]]
[1] 16 17 18 19 20
[[3]]
[1] 41 42 43 44 45
我该怎么做?