Hi Friends,
I am getting this error when I am running my R scripts. Can u please
tell me how this error will happened. My script is given below.
Errot: "Error in do.call("expand.grid", dimnames(x)) :"
Script:
************************************************************************************************************************
#!/usr/bin/R
s<-c("GG","GA","GV","GL","GI","GM","GF","GW","GP","GS","GT","GC","GY","GN","GQ","GD","GE","GK","GR","GH","AG","s[l]","AV","AL","AI","AM","AF","AW","AP","AS","AT","AC","AY","AN","AQ","AD","AE","AK","AR","AH","VG","VA","VV","VL","VI","VM","VF","VW","VP","VS","VT","VC","VY","VN","VQ","VD","VE","VK","VR","VH","LG","LA","LV","LL","LI","LM","LF","LW","LP","LS","LT","LC","LY","LN","LQ","LD","LE","LK","LR","LH","IG","IA","IV","IL","II","IM","IF","IW","IP","IS","IT","IC","IY","IN","IQ","ID","IE","IK","IR","IH","MG","MA","MV","ML","MI","MM","MF","MW","MP","MS","MT","MC","MY","MN","MQ","MD","ME","MK","MR","MH","FG","FA","FV","FL","FI","FM","FF","FW","FP","FS","FT","FC","FY","FN","FQ","FD","FE","FK","FR","FH","WG","WA","WV","WL","WI","WM","WF","WW","WP","WS","WT","WC","WY","WN","WQ","WD","WE","WK","WR","WH","PG","PA","PV","PL","PI","PM","PF","PW","PP","PS","PT","PC","PY","PN","PQ","PD","PE","PK","PR","PH","SG","SA","SV","SL","SI","SM","SF","SW","SP","SS","ST","SC","SY","SN","SQ","SD","SE","SK","SR","SH","TG","TA","TV","TL","TI","TM","TF","TW","TP","TS","TT","TC","TY","TN","TQ","TD","TE","TK","TR","TH","CG","CA","CV","CL","CI","CM","CF","CW","CP","CS","CT","CC","CY","CN","CQ","CD","CE","CK","CR","CH","YG","YA","YV","YL","YI","YM","YF","YW","YP","YS","YT","YC","YY","YN","YQ","YD","YE","YK","YR","YH","NG","NA","NV","NL","NI","NM","NF","NW","NP","NS","NT","NC","NY","NN","NQ","ND","NE","NK","NR","NH","QG","QA","QV","QL","QI","QM","QF","QW","QP","QS","QT","QC","QY","QN","QQ","QD","QE","QK","QR","QH","DG","DA","DV","DL","DI","DM","DF","DW","DP","DS","DT","DC","DY","DN","DQ","DD","DE","DK","DR","DH","EG","EA","EV","EL","EI","EM","EF","EW","EP","ES","ET","EC","EY","EN","EQ","ED","EE","EK","ER","EH","KG","KA","KV","KL","KI","KM","KF","KW","KP","KS","KT","KC","KY","KN","KQ","KD","KE","KK","KR","KH","RG","RA","RV","RL","RI","RM","RF","RW","RP","RS","RT","RC","RY","RN","RQ","RD","RE","RK","RR","RH","HG","HA","HV","HL","HI","HM","HF","HW","HP","HS","HT","HC","HY","HN","HQ","HD","HE","HK","HR","HH")
for(l in 1:length(s))
{
for(j in 1:50)
{
if(file.info(paste("/media/_1/BootStrap/All_alpha_proteins/",s[l],"/alp",s[l],"mean",j,".bootstrap",sep=""))$size>1)
{
a<-read.table(paste("/media/_1/BootStrap/All_alpha_proteins/",s[l],"/alp",s[l],"mean",j,".bootstrap",sep=""))
A<-a[,1]
Amax<-max(A)
Amin<-min(A)
Salpha<-summary(A)
df<-data.frame(P=Salpha)
color<-"aquamarine1"
flg<-"All Alpha"
}
if(file.info(paste("/media/_1/BootStrap/All_beta_proteins/",s[l],"/bet",s[l],"mean",j,".bootstrap",sep=""))$size>1)
{
b<-read.table(paste("/media/_1/BootStrap/All_beta_proteins/",s[l],"/bet",s[l],"mean",j,".bootstrap",sep=""))
B<-b[,1]
Bmax<-max(B)
Bmin<-min(B)
Sbeta<-summary(B)
df<-data.frame(df,Q=Sbeta)
color<-c(color,"brown1")
flg<-c(flg,"All Beta")
}
if(file.info(paste("/media/_1/BootStrap/Alpha_and_beta_proteins_aorb/",s[l],"/aob",s[l],"mean",j,".bootstrap",sep=""))$size>1)
{
c<-read.table(paste("/media/_1/BootStrap/Alpha_and_beta_proteins_aorb/",s[l],"/aob",s[l],"mean",j,".bootstrap",sep=""))
C<-c[,1]
Cmax<-max(C)
Cmin<-min(C)
Saob<-summary(C)
df<-data.frame(df,R=Saob)
color<-c(color,"cadetblue1")
flg<-c(flg,"Alpha or Beta")
}
if(file.info(paste("/media/_1/BootStrap/Alpha_and_beta_proteins_a+b/",s[l],"/apb",s[l],"mean",j,".bootstrap",sep=""))$size>1)
{
d<-read.table(paste("/media/_1/BootStrap/Alpha_and_beta_proteins_a+b/",s[l],"/apb",s[l],"mean",j,".bootstrap",sep=""))
D<-d[,1]
Dmax<-max(D)
Dmin<-min(D)
Sapb<-summary(D)
df<-data.frame(df,S=Sapb)
color<-c(color,"darkorchid1")
flg<-c(flg,"Alpha plus Beta")
}
Tmax<-c(Amax,Bmax,Cmax,Dmax)
Tmin<-c(Amin,Bmin,Cmin,Dmin)
pdf
(paste(s[l],"mean",j,"_bootstrap.pdf",sep=""))
boxplot(df,names=flg,ylim=c(min(Tmin),max(Tmax)),col=color,style =
"tukey")
title(paste("Residue Pair Distance
Statistics","\n",s[l]," Paired
Distance of ",j,sep=""),ylab="Distance in
Angstrom",xlab="Protein
Classes")
dev.off()
}
}
*********************************************************************************************
--
******************************************
Deepak.M.R
Research Scholar
Biocomputing Group,
Department of Computer Science,
University of Bologna (Italy).
http://deepak.web.cs.unibo.it/.