| 查看: 1913 | 回复: 0 | ||
[求助]
怎么做出文献中dssp的图
|
|
请教各位前辈 amber跑完轨迹，想看二级结构变，参照http://blog.sina.com.cn/s/blog_834872aa01019mnt.html这篇博文用vmd提取的信息 求问上面得到信息怎么才能做出文献中的图来？ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ #上一个博文中介绍用vmd提取dssp计算的二级结构，生成sec4traj.txt数据文件。然后就是需要做图，如同gmx的do_dssp生成的类似的图。 #导入数据 mysecdata <- read.table("sec4traj.txt", header=T, check.names=F, colClasses="character"  #定义为matrix mysecdata <- as.matrix(mysecdata) #拷贝一个进行数据处理 mysec <- mysecdata #现在的情况是mysec的data还是character类型，先转换一下，根据dssp的定义 #H = alpha helix #B = residue in isolated beta-bridge #E = extended strand, participates in beta ladder #G = 3-helix (3/10 helix) #I = 5 helix (pi helix) #T = hydrogen bonded turn #S = bend #我们的数据里的 X = loop or irregular #其实到这里为止就可以统计各种二级结构的百分含量了，例如我们要统计alpha helix的含量，就是计算每一行中字母“H”的含量； perH <- NULL for(i in 1:nrow(mysec)){perH=length(which(mysec[i,]=="H" )/ncol(mysec)}myhelix <- data.frame(rownames(mysec), perH) #对于myhelix就可以做图观察其含量变化了 #还是需要转换为数值型的 mysec[] <- match(mysecdata, c("X", "H", "B", "E", "G", "I", "T", "S" )mode(mysec) <- "numeric" mycol <- as.numeric(colnames(mysec)) myrow <- as.numeric(rownames(mysec)) #然后做图 #先定义一个颜色板吧，分别代表X H B E G I T S mycolor <- c("white", "blueviolet", "black", "yellow", "blue", "cyan", "red", "green",) image(min(myrow):max(myrow), min(mycol):max(mycol), mysec, col=mycolor) #这样应该差不多了,同样的做图适应于amber产生的数据 这是用的什么软件？ |
回复此楼» 猜你喜欢
什么是人一生最重要的？
已经有4人回复
-
版面费该交吗
已经有3人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复