This function calculates the Distribution descriptor of the CTD descriptors (dim: 105).
Note
For this descriptor type, users need to intelligently evaluate the underlying details of the descriptors provided, instead of using this function with their data blindly. It would be wise to use some negative and positive control comparisons where relevant to help guide interpretation of the results.
References
Inna Dubchak, Ilya Muchink, Stephen R. Holbrook and Sung-Hou Kim. Prediction of protein folding class using global description of amino acid sequence. Proceedings of the National Academy of Sciences. USA, 1995, 92, 8700-8704.
Inna Dubchak, Ilya Muchink, Christopher Mayor, Igor Dralyuk and Sung-Hou Kim. Recognition of a Protein Fold in the Context of the SCOP classification. Proteins: Structure, Function and Genetics, 1999, 35, 401-407.
See also
See extractCTDC and extractCTDT
for Composition and Transition of the CTD descriptors.
Author
Nan Xiao <https://nanx.me>
Examples
x <- readFASTA(system.file("protseq/P00750.fasta", package = "protr"))[[1]]
extractCTDD(x)
#> prop1.G1.residue0 prop1.G1.residue25 prop1.G1.residue50 prop1.G1.residue75
#> 0.3558719 23.1316726 50.1779359 73.8434164
#> prop1.G1.residue100 prop1.G2.residue0 prop1.G2.residue25 prop1.G2.residue50
#> 99.8220641 0.5338078 27.4021352 47.3309609
#> prop1.G2.residue75 prop1.G2.residue100 prop1.G3.residue0 prop1.G3.residue25
#> 75.2669039 100.0000000 0.1779359 19.5729537
#> prop1.G3.residue50 prop1.G3.residue75 prop1.G3.residue100 prop2.G1.residue0
#> 51.7793594 75.6227758 99.6441281 0.3558719
#> prop2.G1.residue25 prop2.G1.residue50 prop2.G1.residue75 prop2.G1.residue100
#> 25.6227758 48.0427046 75.4448399 100.0000000
#> prop2.G2.residue0 prop2.G2.residue25 prop2.G2.residue50 prop2.G2.residue75
#> 1.4234875 23.3096085 54.4483986 76.3345196
#> prop2.G2.residue100 prop2.G3.residue0 prop2.G3.residue25 prop2.G3.residue50
#> 99.4661922 0.1779359 22.7758007 48.9323843
#> prop2.G3.residue75 prop2.G3.residue100 prop3.G1.residue0 prop3.G1.residue25
#> 69.5729537 99.8220641 0.1779359 20.9964413
#> prop3.G1.residue50 prop3.G1.residue75 prop3.G1.residue100 prop3.G2.residue0
#> 50.8896797 74.5551601 99.6441281 0.5338078
#> prop3.G2.residue25 prop3.G2.residue50 prop3.G2.residue75 prop3.G2.residue100
#> 26.5124555 46.2633452 75.4448399 100.0000000
#> prop3.G3.residue0 prop3.G3.residue25 prop3.G3.residue50 prop3.G3.residue75
#> 0.3558719 24.1992883 50.5338078 73.8434164
#> prop3.G3.residue100 prop4.G1.residue0 prop4.G1.residue25 prop4.G1.residue50
#> 99.8220641 0.3558719 26.5124555 48.3985765
#> prop4.G1.residue75 prop4.G1.residue100 prop4.G2.residue0 prop4.G2.residue25
#> 76.1565836 99.2882562 1.4234875 21.5302491
#> prop4.G2.residue50 prop4.G2.residue75 prop4.G2.residue100 prop4.G3.residue0
#> 51.4234875 75.8007117 100.0000000 0.1779359
#> prop4.G3.residue25 prop4.G3.residue50 prop4.G3.residue75 prop4.G3.residue100
#> 22.7758007 48.9323843 69.5729537 99.8220641
#> prop5.G1.residue0 prop5.G1.residue25 prop5.G1.residue50 prop5.G1.residue75
#> 0.8896797 20.8185053 48.9323843 69.5729537
#> prop5.G1.residue100 prop5.G2.residue0 prop5.G2.residue25 prop5.G2.residue50
#> 99.8220641 0.1779359 24.9110320 49.1103203
#> prop5.G2.residue75 prop5.G2.residue100 prop5.G3.residue0 prop5.G3.residue25
#> 75.2669039 100.0000000 0.3558719 26.1565836
#> prop5.G3.residue50 prop5.G3.residue75 prop5.G3.residue100 prop6.G1.residue0
#> 64.2348754 77.4021352 99.2882562 0.1779359
#> prop6.G1.residue25 prop6.G1.residue50 prop6.G1.residue75 prop6.G1.residue100
#> 22.9537367 50.8896797 74.3772242 99.8220641
#> prop6.G2.residue0 prop6.G2.residue25 prop6.G2.residue50 prop6.G2.residue75
#> 1.6014235 21.5302491 49.2882562 70.8185053
#> prop6.G2.residue100 prop6.G3.residue0 prop6.G3.residue25 prop6.G3.residue50
#> 98.9323843 0.3558719 29.0035587 48.2206406
#> prop6.G3.residue75 prop6.G3.residue100 prop7.G1.residue0 prop7.G1.residue25
#> 77.4021352 100.0000000 0.5338078 23.4875445
#> prop7.G1.residue50 prop7.G1.residue75 prop7.G1.residue100 prop7.G2.residue0
#> 50.0000000 74.5551601 98.9323843 0.3558719
#> prop7.G2.residue25 prop7.G2.residue50 prop7.G2.residue75 prop7.G2.residue100
#> 23.1316726 50.1779359 73.8434164 99.8220641
#> prop7.G3.residue0 prop7.G3.residue25 prop7.G3.residue50 prop7.G3.residue75
#> 0.1779359 27.2241993 48.0427046 75.4448399
#> prop7.G3.residue100
#> 100.0000000