Intermediate dot indexing produced a comma-separated list with 34 values, but it must produce a single value when followed by subsequent indexing operations
10 次查看(过去 30 天)
显示 更早的评论
% Starting from pdb file
Proteinname='1UXD';
uxd = getpdb(Proteinname);
No_of_residues = uxd.Sequence.NumOfResidues;
uxd.Model.Atom.resName
isTYR = ({uxd.Model.Atom.resName} == "TYR")
uxdTYR = uxd.Model.Atom(isTYR)
In the above code, I am trying to extract the details of only the tyrosine residues from the protein '1UXD'. I am guessing the error is due to the presence of many models of the same protein in the pdb file which might have caused the output to be a comma separated list. I am in need of only the first model's tyrosine (TYR) details.. There is a possibility to handle this by turning the output into a cell array and extracting the first output which corresponds to the first model but I am uncertain of how to approach the code for this. Can anyone help :) Thankyou!
0 个评论
采纳的回答
Stephen23
2023-3-18
编辑:Stephen23
2023-3-18
The documentation
states the the MODEL field can be a structure array. Your code does not take that into account.
"I am in need of only the first model's tyrosine (TYR) details"
If you only want the first model, then just use indexing into that array:
uxd.Model(1).Atom.resName
For more information:
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!