Guest User

Untitled

a guest
Jun 21st, 2018
93
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.55 KB | None | 0 0
  1. import numpy as np
  2.  
  3. import rpy2.robjects.numpy2ri
  4. from rpy2.robjects.packages import importr
  5. rpy2.robjects.numpy2ri.activate()
  6.  
  7. # Set up our R namespaces
  8. R = rpy2.robjects.r
  9. DTW = importr('dtw')
  10.  
  11. # Generate our data
  12. idx = np.linspace(0, 2*np.pi, 100)
  13. template = np.cos(idx)
  14. query = np.sin(idx) + np.array(R.runif(100))/10
  15.  
  16. # Calculate the alignment vector and corresponding distance
  17. alignment = R.dtw(query, template, keep=True,window_type='sakoechiba',
  18. window_size=5)
  19.  
  20. >>> RRuntimeError: Error in window.function(row(wm), col(wm), query.size= n, reference.size = m, :
  21. argument "window.size" is missing, with no default
  22.  
  23. from rpy2.robjects.functions import SignatureTranslatedFunction
  24.  
  25. Rsession.dtw = SignatureTranslatedFunction(Rsession.dtw,
  26. init_prm_translate={'window_size': 'window.size'})
  27.  
  28. import numpy as np
  29.  
  30. import rpy2.robjects.numpy2ri
  31. from rpy2.robjects.packages import importr
  32. from rpy2.robjects.functions import SignatureTranslatedFunction
  33. rpy2.robjects.numpy2ri.activate()
  34.  
  35. # Set up our R namespaces
  36. R = rpy2.robjects.r
  37. DTW = importr('dtw')
  38. Rsession.dtw = SignatureTranslatedFunction(Rsession.dtw,init_prm_translate={'window_size': 'window.size'})
  39.  
  40. # Generate our data
  41. idx = np.linspace(0, 2*np.pi, 100)
  42. template = np.cos(idx)
  43. query = np.sin(idx) + np.array(R.runif(100))/10
  44.  
  45. # Calculate the alignment vector and corresponding distance
  46. alignment = R.dtw(query, template, keep=True,window_type='sakoechiba',
  47. window_size=10)
  48. dist = alignment.rx('distance')[0][0]
  49.  
  50. print(dist)
  51. >>> 117.348292359
Add Comment
Please, Sign In to add comment