Functions wcc
and wac
are meant for profile data
(intensities measured at equidistant time points), whereas
wcc.st
and wac.st
are meant for peak lists. In general,
wcc values calculated for profiles will be higher since they will also
include the large similarity in the empty spaces, i.e., parts of the
profiles where no peaks are present (and that will appear to be
perfectly aligned), whereas the peak-based version concentrates only
on the peaks.
wcc(pattern1, pattern2, trwdth, wghts = NULL, acors1 = NULL, acors2 = NULL)
wac(pattern1, trwdth, wghts = NULL)
wcc
and wac
are defined such that the
triangle width stands for
the number of points on one side of and including the current
point. Thus, a trwdth
of 0 signifies a non-existent triangle and
results in an error; a trwdth
equal to 1 only includes the current
point with weight 1 and no neighbouring points. For the stick-based
equivalents, the units of the time axis are used for the triangle width.
data(gaschrom)
wcc(gaschrom[1,], gaschrom[2,], 20)
wcc.st(gaschrom.st[[1]], gaschrom.st[[2]], 20)
Run the code above in your browser using DataLab