C_R <- log10(-NMF_spectra)修改这个错误，将吸光度数据转换为反射率数据

C_R <- log10(-NMF_spectra)

这个表达式看起来是在尝试使用对数变换（log10），但它直接应用于负的NMF_spectra值，这在吸光度到反射率的转换中是不正确的。吸光度通常是以正数表示的，其值越高代表吸收越强；而反射率则是从0到1的范围内，通常表示为负的吸光度。

正确的步骤应该是首先将吸光度转换为浓度（例如，对于比色法，浓度可以是负的吸光度除以特定常数），然后再取对数。然后，为了得到反射率，你需要将浓度转换成反射率比例，而不是取对数。这里是一个简单的例子：

1. 如果\( A_{sample} \)是样品的吸光度，\( A_{blank} \)是空白溶液的吸光度（作为参考），那么浓度\( C \)通常是:
\( C = \frac{-A_{sample}}{A_{blank}} \)

2. 然后将浓度转换为反射率:
\( R = 100 \times (1 - e^{-C}) \) 或 \( R = 100 \times (1 + C) \)（取决于原始吸光度是否已减去空白）

所以，正确的转换步骤可能是这样：

# 假设A_sample是你的NMF_spectra
Concentration <- -NMF_spectra / A_blank
Reflectivity <- 100 * (1 - exp(-Concentration)) # 或者 100 * (1 + Concentration)

请注意，具体的转换公式可能依赖于实验设置，比如使用的单位和校准方式。