Working with the Shootout 2016 data with R (Part 1)

wavelengths_A<-seq(from=730,to=1100,by=0.5)

wavelengths_B<-seq(from=570,to=1100,by=0.5)

wavelengths_C<-seq(from=850,to=1048,by=2)

######  WORKING WITH "MANUFACTURER A" DATA #########################

### A1 Instrument

CalSetA1_spec<-as.matrix(CalSetA1[,3:743])              #Defining the spectra range

    matplot(wavelengths_A,t(CalSetA1_spec),lty=1,type="l",

        pch=NULL,xlab="nm",ylab="abs",ylim=c(-0.5,6),col="red")

colsC<-seq(from=1,to=637,by=4)               #to convert from 0.5 to 2.0 nm

CalSetA1_2nm_spec<-CalSetA1_spec[,colsC]     # A1 with 2nm interval

CalSetA1_tC_spec<-CalSetA1_2nm_spec[,61:160] # A1 trimmed to C

    matplot(wavelengths_C,t(CalSetA1_tC_spec),lty=1,type="l",

        pch=NULL,xlab="nm",ylab="abs",ylim=c(-0.5,4),col="blue")

   

### A2 Instrument

CalSetA2_spec<-as.matrix(CalSetA2[,3:743])

    matplot(wavelengths_A,t(CalSetA2_spec),lty=1,type="l",

        pch=NULL,xlab="nm",ylab="abs",ylim=c(-0.5,6),col="red")

CalSetA2_2nm_spec<-CalSetA2_spec[,colsC]     # A2 with 2nm interval

CalSetA2_tC_spec<-CalSetA2_2nm_spec[,61:160] # A2 trimmed to C

    matplot(wavelengths_C,t(CalSetA2_tC_spec),lty=1,type="l",

        pch=NULL,xlab="nm",ylab="abs",ylim=c(-0.5,4),col="blue")

### A3 Instrument

CalSetA3_spec<-as.matrix(CalSetA3[,3:743])

    matplot(wavelengths_A,t(CalSetA3_spec),lty=1,type="l",

        pch=NULL,xlab="nm",ylab="abs",ylim=c(-0.5,6),col="red")

CalSetA3_2nm_spec<-CalSetA3_spec[,colsC]     # A3 with 2nm interval

CalSetA3_tC_spec<-CalSetA3_2nm_spec[,61:160] # A3 trimmed to C

    matplot(wavelengths_C,t(CalSetA3_tC_spec),lty=1,type="l",

        pch=NULL,xlab="nm",ylab="abs",ylim=c(-0.5,4),col="blue")

 This is the plot of the spectra in the calibration set for Manufacturer A scanned on the instruments A1, A2 and A3.

 This is the plot of the spectra in the calibration set for Manufacturer A scanned on the instruments A1, A2 and A3 timmed to the range of C (instrument with the lower range) .