I would like to compare two absorption spectra ( oror interferograms) and conclude whether between these two there are statistically significant differences at particular wavelength intervals. At the moment, I have data of two experiments that look like this:
# A tibble: 6 x 5
t x1 y1 x2 y2
<dbl> <dbl> <dbl> <dbl> <dbl>
1 3999. 0.0124 0.0132 0.0122 0.0113
2 3998. 0.0125 0.0130 0.0122 0.0116
3 3997. 0.0122 0.0131 0.0122 0.0113
4 3996. 0.0121 0.0136 0.0122 0.0114
5 3995. 0.0124 0.0139 0.0122 0.0122
6 3994. 0.0125 0.0141 0.0122 0.0129
The first column represents the wavenumber, the xx
columns represent the absorbance of sample and the yy
columns represent the absorbance of irradiated sample ( beforebefore and after).
I was wondering whether I could compare these data ( xx
and yy
) as time series and if so, what could be the method to quantify the differences, if any, between the samples before and after irradiation. Maybe it's already been done and there is somewhere some information as to how to compare the spectra if the wavenumber is interpreted as time ( x$x$ axis).
I did the t test-test in R and in both experiments the null hypothesis could not be rejected, although for the second experiment (x2x2
, y2y2
) the p$p$ value was much lower than for the first.
If If I average the xx
and yy
, and then plot both data, I see that there are visible differences at certain wavelength intervals. But how can I check for sure the differences between spectrums?
Here
Here is a project with similar experiments: https://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-50532015001202571 , for by Zezell et al. [1]. For statistical analysis they use ANOVA and Tukey's test, but how do I do it for the vectored data?
Reference
- Zezell, D. M.; Benetti, C.; Veloso, M. N.; Castro, P. A. A.; Ana, P. A. FTIR Spectroscopy Revealing the Effects of Laser and Ionizing Radiation on Biological Hard Tissues. J. Braz. Chem. Soc 2015. DOI: 10.5935/0103-5053.20150246.