批量计算多肽段的理化性质----R
主程序
library("Peptides")
ww<-"data_1.csv"
data<-read.csv(ww,1)
cc<-dim(data)
rowname<-(data$Sequence)
colname<-c("mw","hydrophobicity","Boman index","Net charge","GRAVY")
numcol<-5
jj<-matrix(nrow=cc[1],ncol=numcol,dimnames=list(rowname,colname))
for (i in 1:cc[1]){
jj[i,1]=mw(seq=data[i,1])
jj[i,2]=hydrophobicity(seq =data[i,1] , scale = "Eisenberg")
jj[i,3]=boman(seq =data[i,1])
jj[i,4]=charge(seq =data[i,1], pH = 7, pKscale = "EMBOSS")
jj[i,5]=hydrophobicity(seq =data[i,1] , scale = "Eisenberg")
}
write.csv(jj,"protein-peptides_5.csv")
输入文件数据
![在这里插入图片描述](https://img-blog.csdnimg.cn/bd759800cedd4ff298d91f3af86d911b.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA54Ot5Yqb5a2m6LSf5LiA5a6a55CG,size_10,color_FFFFFF,t_70,g_se,x_16)
输出结果示例
![在这里插入图片描述](https://img-blog.csdnimg.cn/bd0df338b86f4c48afc9812c8ce89b2f.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA54Ot5Yqb5a2m6LSf5LiA5a6a55CG,size_15,color_FFFFFF,t_70,g_se,x_16)