library(R.utils)
mvec<-intToBin(0:15)
ModelMat<-matrix(as.numeric(unlist(strsplit(mvec,split=""))),byrow=T,nrow=16,ncol=4)

xnames<-c("1","x1","x2","x3")
ModelForm<-character(0)
for(j in 1:16){
	if(ModelMat[j,4] == 0){
		mn<-"y~-1"
	}else{
		mn<-"y~ 1"
	}
	if(ModelMat[j,3] == 1){
		mn<-paste(mn,"+x1",sep="")
	}
	if(ModelMat[j,2] == 1){
		mn<-paste(mn,"+x2",sep="")
	}
	if(ModelMat[j,1] == 1){
		mn<-paste(mn,"+x3",sep="")
	}
	ModelForm<-c(ModelForm,mn)

}

#AIC/BICs
dta<-data.frame(x1,x2,x3,y)
AICv<-BICv<-rep(0,16)
for(j in 1:16){
	fj<-lm(formula(ModelForm[j],data=dta))
	AICv[j]<-AIC(fj)	
	BICv[j]<-BIC(fj)
}
astar<-bstar<-rep(" ",16)
astar[which.min(AICv)]<-"*"
bstar[which.min(BICv)]<-"*"
	
result<-data.frame(form=ModelForm,AIC=round(AICv,6),ASel=astar,BIC=round(BICv,6),BSel=bstar)

write(' ',file='Assign3-BIC-Res.txt',append=T)
write(paste(ID.data$ID[i],as.character(ID.data$Name[i])),file='Assign3-BIC-Res.txt',append=T)	
write(' ',file='Assign3-BIC-Res.txt',append=T)

	
capture.output(print(result),file='Assign3-BIC-Res.txt',append=T)
write(' ',file='Assign3-BIC-Res.txt',append=T)
write(' ',file='Assign3-BIC-Res.txt',append=T)
write(' ',file='Assign3-BIC-Res.txt',append=T)
write('-----------------------------',file='Assign3-BIC-Res.txt',append=T)