r - 如何计算所有可能波段组合中的归一化比率指数，然后与 R 中的环境变量相关？

Question

我想计算所有可能波段组合中的归一化比率和简单比率指数，然后想与 R 中的环境变量相关联。然后我想确定具有最高相关性的最佳组合。该实现hsdar在 R 中的包中可用。但是对于大型数据集来说非常慢。我在这里附上一个小数据集

Environmental_variable  WV_400  WV_401  WV_402  WV_403  WV_404  WV_405
95.512  0.035   0.034   0.034   0.034   0.034   0.034
97.900  0.047   0.047   0.047   0.046   0.046   0.046
92.897  0.004   0.004   0.006   0.008   0.009   0.009
94.209  0.011   0.012   0.013   0.016   0.017   0.017
87.472  0.010   0.010   0.011   0.014   0.015   0.015
91.109  0.010   0.011   0.013   0.015   0.016   0.016
92.830  0.024   0.025   0.026   0.028   0.029   0.029

我从hsdar包中给出一个示例代码以供参考

library(hsdar)
data(spectral_data)
## Calculate normalised ratio indices in all possible combinations
nri_WV <- nri(spectral_data, recursive = TRUE)
## Build glm-models
glmnri <- glm.nri(nri_WV ~ chlorophyll, preddata = spectral_data)
## Return best 10 models
BM <- nri_best_performance(glmnri, n = 10, coefficient = "p.value")

任何以 R 代码形式作为hsdar包的快速替代方案的帮助都受到高度赞赏。

标签： rcombinations