Displaying 1 result from an estimated 1 matches for "xnam2".
Did you mean:
xname
2007 Dec 19
0
leaps
...sinmat[,i] <- sin(tt*2*pi*i/T)
sincol <- c(sincol,paste("sin",i,sep=""))
}
colnames(cosmat) <- coscol
colnames(sinmat) <- sincol
xnam1 <- NULL
xnam1 <- paste(sep="","cosmat[,",1:nfr,"]",collapse="+")
xnam2 <- NULL
xnam2 <- paste(sep="","sinmat[,",1:nfr,"]",collapse="+")
fmla <- as.formula(paste("yy ~ ", paste(xnam1,"+",xnam2,sep="")))
FTmod <- lm(fmla)
stepmod <- step(FTmod, direction="both")...