nanxstats / protr

🧬 Toolkit for generating various numerical features of protein sequences
https://nanx.me/protr/
Other
52 stars 13 forks source link

extractCTDC fails when no amino acids map to any one group #4

Closed koefoed closed 8 years ago

koefoed commented 9 years ago
> packageVersion('protr')
[1] ‘0.5.1’
> extractCTDC("NAKADQASSDAQTANAKADQASNDANAARSDAQAAKDDAARANQRADNAA")
Error in names(CTDC) = paste("prop", rep(1:7, each = 3), ".", names(CTDC),  : 
  'names' attribute [21] must be the same length as the vector [19]
> 

In the above example no amino acids map to property 1, group 3 or property 3 group 1. Thus the entire group is left out in lines 32-34 of the source code (which returns nothing).