我正在R中使用基线包,只是在阅读文档。他们的示例代码产生以下图像。
data(milk)
bc.irls <- baseline(milk$spectra[1,, drop=FALSE])
plot(bc.irls)
在我的情况下,x刻度不应以0开头,而应以5000开头(因此x轴的注释偏离了5000)。是否可以在基线构造函数中更改范围?还是可以以某种方式在plot方法中调整x标签?
根据文档,baseline
是带有访问器方法getCorrected
的S4类,该方法返回一个矩阵。
在您的示例中,此矩阵只有一行,共有21,452列。因此,您可以将其转换为向量,并删除前4999个元素,以使光谱从5000开始。您还需要创建5,000到21,452的序列作为x轴值。然后,您可以像绘制任何其他x,y向量一样进行绘图:
new_spec <- as.vector(getCorrected(bc.irls))
new_spec <- new_spec[5000:length(new_spec)]
new_x <- seq_along(new_spec) + 4999
plot(new_x, new_spec, type = "l", main = "Baseline corrected spectrum")