fname<-paste(YOURIDNUMBERHERE,'.csv',sep='') Cigi<-read.csv(fname,header=T) x1<-Cigi$TAR x2<-Cigi$NICOTINE x3<-Cigi$WEIGHT y<-Cigi$CO dti<-data.frame(x1,x2,x3,y) fit.full123<-lm(y~x1+x2+x3,data=dti) a.val<-summary(aov(fit.full123))[[1]][4,2] #(a) b.val<-summary(aov(lm(y~x1+x2,data=dti)))[[1]][3,2] #(b) c.val<-drop1(fit.full123,test="F")[4,5] #(c) fit.full321<-lm(y~x3+x2+x1,data=dti) d.val<-c(sum(anova(fit.full321)[1:3,2]),anova(fit.full321)[1:3,2]) #(d) fit.full12<-lm(y~x1+x2,data=dti) e.val<-c(sum(anova(fit.full12)[1:2,2]),anova(fit.full12)[1:2,2]) #(e) #Using X1+X2 f.val1<-drop1(fit.full12,test="F")[3,5] #(f) #Using X1+X2+X3 f.val2<-anova(fit.full123)[2,4] #(f) #Using X1+X2 g.val1<-summary(fit.full12)$f[1] #(g) #Using X1+X2+X3 SSRes0<-anova(lm(y~1,data=dti))[1,2] SSRes012<-anova(lm(y~x1+x2,data=dti))[3,2] MSRes0123<-anova(fit.full123)[4,3] g.val2<-((SSRes0-SSRes012)/2)/MSRes0123 #(g) write(paste(ID.data$ID[i],as.character(ID.data$Name[i])),file='Assign3-Res.txt',append=T) #write(NamesIDs[i,2],file='Assign3-Res.txt',append=T) write(' ',file='Assign3-Res.txt',append=T) write(paste('(a)',round(a.val,6)),file='Assign3-Res.txt',append=T) write(paste('(b)',round(b.val,6)),file='Assign3-Res.txt',append=T) write(paste('(c)',round(c.val,6)),file='Assign3-Res.txt',append=T) write(paste('(d)',round(d.val,6)),ncol=4,file='Assign3-Res.txt',append=T) write(paste('(e)',round(e.val,6)),ncol=3,file='Assign3-Res.txt',append=T) #write('Using X1+X2',file='Assign3-Res.txt',append=T) write(paste('(f)',round(f.val1,6)),file='Assign3-Res.txt',append=T) #write('Using X1+X2',file='Assign3-Res.txt',append=T) write(paste('(g)',round(g.val1,6)),file='Assign3-Res.txt',append=T) write(' ',file='Assign3-Res.txt',append=T) write('-----------------------------',file='Assign3-Res.txt',append=T) write(' ',file='Assign3-Res.txt',append=T)