sFName<-c('0007165','0000902','0007010','0007005','0007033','0006364','0007047','0016568','0006260','0007067','0007126','0000074','0006605','0006811','0045045','0016192','0000747','0006066','0009309','0006520','0008610','0006412','0005975','0006511','0015980','0006281','0006397','0008380','0006357','0016310','0006464','0009628','0009607','0006950',);
for(j in 0:4)
{
for(i in 1:34)
{
filename<-paste("/home/hyunjul/public_html/klr//ppipfamCVBPYeast/LRFeatures/TrainFeature",sFName[i],"_",j,".txt",sep="")
features<-scan(filename,list(protein="",label=0,feature1=0,feature2=0,feature3=0,feature4=0),sep='\t');
myfit<-glm(formula=label ~ feature1+feature2+feature3+feature4,family=binomial(link=logit),data = features, na.action = na.exclude);
para<-matrix(myfit$coefficients);
filename<-paste("/home/hyunjul/public_html/klr//ppipfamCVBPYeast/Para/para",sFName[i],"_",j,".txt",sep="")
write.table(para,file=filename,row.names=FALSE,col.names=FALSE);
}
}