Displaying 1 result from an estimated 1 matches for "p_main".
Did you mean:
_main
2009 Oct 28
2
regression on large file
...max){
line<-scan(file=paste(file),nlines=1,skip=(i-1),what="integer",sep=",")
d<-as.numeric(line[-1])
name<-line[1]
modela <- lm(s1~a+a2+b+s+M+W)
modelb <- lm(s2~a+a2+b+s+M+W+d)
modelc <- lm(s3~a+2+b+s+M+W+d+d*s)
p_main <- anova(modela,modelb)$P[2]
p_main_i <- anova(modela,modelc)$P[2]
p_i <- anova(modelb,modelc)$P[2]
cat(c(name,p_main,p_main_i,p_i),file=paste("out",".txt",sep=""),append=T)
cat("\n",file=paste("out",".txt"...