Revision da02fed50f656619ef5518cab97c0d00fb23d318 authored by Toni Giorgino on 08 January 2008, 00:00:00 UTC, committed by Gabor Csardi on 08 January 2008, 00:00:00 UTC
1 parent c8ecf75
dtw_test.Rout.save
``````
R version 2.6.1 (2007-11-26)
Copyright (C) 2007 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(dtw);

Attaching package: 'proxy'

The following object(s) are masked from package:stats :

as.dist,
dist

>
> ### Synthetic example: check indexes, distance
> ldist<-matrix(1,nrow=6,ncol=6);  # Matrix of ones
> ldist[2,]<-0; ldist[,5]<-0;      # Mark a clear path of zeroes
> ldist[2,5]<-.01;		 # Forcely cut the corner
>
> ds<-dtw(ldist);			 # DTW with user-supplied local cost matrix
> ds\$distance			 # 2
[1] 2
> ds\$index1			 # 1 2 2 2 2 3 4 5 6 6
[1] 1 2 2 2 2 3 4 5 6 6
> ds\$index2			 # 1 1 2 3 4 5 5 5 5 6
[1] 1 1 2 3 4 5 5 5 5 6
>
> da<-dtw(ldist,step="a");	 # Also compute the asymmetric
> da\$distance			 # 2
[1] 2
> da\$index1			 # 1 2 3 4 5 6
[1] 1 2 3 4 5 6
> da\$index2			 # 1 3 5 5 5 6
[1] 1 3 5 5 5 6
>
> da<-dtw(ldist,step=asymmetricP0);	 # The strange aP0 alignment
> da\$distance			 # 1
[1] 1
> da\$index1			 # 1 2 2 2 2 3 4 5 6 6
[1] 1 2 2 2 2 3 4 5 6 6
> da\$index2			 # 1 1 2 3 4 5 5 5 5 6
[1] 1 1 2 3 4 5 5 5 5 6
>
> ds<- globalCostMatrix(ldist)
> dsn<- globalCostMatrix(ldist,native=FALSE)
Warning message:
In globalCostMatrix(ldist, native = FALSE) :
Native dtw implementation not available: using (slow) interpreted fallback
> stopifnot(all.equal(ds,dsn))		 # TRUE
>
>
> ###  Sine/cosine example: verify native output
> ### there may be a random chance of failing due to rounding errors
> idx<-seq(0,6.28,len=100);
> query<-sin(idx)+runif(100)/10;
> template<-cos(idx)
> ldist<-outer(query,template,FUN=function(x,y){(x-y)^2})
> ds<- globalCostMatrix(ldist)
> dsn<- globalCostMatrix(ldist,native=FALSE)
Warning message:
In globalCostMatrix(ldist, native = FALSE) :
Native dtw implementation not available: using (slow) interpreted fallback
> all.equal(ds,dsn)		# TRUE
[1] TRUE
>
>
``````

Computing file changes ...