Displaying 1 result from an estimated 1 matches for "crossval2".
Did you mean:
crossval
2011 Nov 07
0
Adding lines to scatterplot odd result when creating multiple plots
...ss (mE)"
Y<-"True thickness (mE)"
#load data
#read 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= "Mi...