search for: aveest

...the data skip then read header crossval <- read.table(Infile,skip = 9,sep = "") head<-readLines(Infile,9) #decode the header lines wanted head2<-head[3:9] colnames(crossval)=head2 #Filter out non-estimated crossval2<-crossval[crossval$Estimate>0,] # Compute the means AveEst<-mean(crossval2$Estimate) AveTru<-mean(crossval2$True) #Fit the regression line Fit<-lm(crossval2$True~crossval2$Estimate ) #create plots scatterplot(True ~ Estimate, data=crossval2, xlab= X, ylab= Y, main= "Min 4 Max 8", grid=FALSE, xlim=c(0,8), ylim=c(0,8), pch=21,...