search for: xnam

...anks in advance, Scott Gilpin ## Example code rm(list=ls()) ## Forward stepwise regression - works fine set.seed(2863) x<-matrix(runif(1000),ncol=10) colnames(x)<-1:10 beta<-matrix(c(1,2,3,4,5,0,0,0,0,0),ncol=1) y<-drop(x %*% beta + rnorm(100)) lsdata<-data.frame(cbind(x=x,y=y)) xnam <- names(lsdata)[names(lsdata) != "y"] (fmla <- as.formula(paste("y ~ ", paste(xnam, collapse= "+")))) fit<-lm(y ~ 1,data=lsdata) step.fit<-step(fit,scope = fmla, direction="forward",data=lsdata) ## start from scratch rm(list=ls()) set.seed(2863)...

...repare .asc-files filenames <- c("R2", "adj.R2", "p", "b0", colnames (x)) for (i in 1:items) { write.table (hd, file = paste (outfolder, filenames [i],".asc",sep =""), quote=F, row.names=F, col.names=F) } rm (hd) # Prepare regression xnam <- paste ("x[,", 1:(ncol(x)),"]", sep="") form <- paste("y ~ ", paste(xnam, collapse="+")) rm (xnam) # Loop through rows for (j in 1:rows) { getgrid <- function (j) { print (paste ("Row",j,"/",rows),quote = F)...

Hi, If I create a formula with say 100 terms and then paste it: xnam <- paste("x", 1:100, sep="") fmla <- as.formula(paste("y ~ ", paste(xnam, collapse= "+"))) paste(fmla) The result seems to cut off everything after the first 500 characters and gives no warning message. I have the most recent version of R from the R...

Hi, I am fitting a weighted least square regression and trying to compute SSE,SST and SSReg but I am not getting SST = SSReg + SSE and I dont know what I am coding wrong. Can you help please? xnam <-colnames(X) # colnames Design Matrix fmla1 <- as.formula(paste("Y ~",paste(xnam, collapse= "+"),"-1",sep="")) fitlm <- lm(formula=fmla1,data = data.frame(cbind(X,Y))) ResiSqr &l...

Hi, Can you please help me with this please? What I am trying to do is call a vector from R function and used in the new function So I create 4 functions with these arguments M11 <- function(TrainData,TestData,mdat,nsam) { ls <- list() I have few statments one of them is vectx <- c(,1,2,3,4,5,6,6) vectz <- c(12,34,5,6,78,9,90) and then................ ls(vectx=vtecx,vectz=vectz)

hi netters i have a vector NAMES containing a series of variable names: NAMES=c(x,r,z,m,st,qr,.....nn). i wanna fit a regression tree by using the code: my.tree<-tree(y~x+r+z+m+....nn,my.dataframe) but i don't want to type out "x+r+z+m+....+nn" one by one, as there are so many variables. besides, sometimes i wanna put the code in a function. so i need to have the

Hi, I'd like to apply a 2^k factorial design with k=10 parameters. Obviously this results in a quite long term for the model equation due to the high number of combinations of parameters. How can I specify the equation for the linear model (lm) without writing all combinations explicitly down by hand? Does a R command exist for this problematic? Thanks for your help in advance, Sven

...k is the first line of function iscube, and it seems that the object r.plan cannot be found here. I have no idea how to fix this. In case it helps: The expression generating the formula in function formula.design is as.formula(paste(response, paste(c(names(factor.names), paste("!iscube(",xnam,")",sep="")), collapse="+"),sep="~")) I am grateful for any pointers that direct me to a fix for this problem and/or to a better understanding of its root cause. Thanks and regards, Ulrike -- View this message in context: http://www.nabble.com/Problem-with...

Hi! I had a database with some variables in sequence. Let me say: TX01, TX02, TX03 and TX04. But I need to run some regressions changing the variables... so: variable <- paste("TX0", 1:4, sep="") for(i in 1:4){ test[i] <- lm(variable[i] ~ INCOME, data=database) } But doesn't work... lm tries to find a variable inside database named variable[i] ... Suggestions?

...t;- function(x, type = "h", ylim = c(0, max(x)), lwd = 2, - xlab = NULL, ylab = NULL, frame.plot = is.num, ...) + xlab = NULL, ylab = NULL, frame.plot = is.num, axes = TRUE, ...) { + localaxis <- function(..., col, bg, pch, cex, lty, lwd) axis(...) xnam <- deparse(substitute(x)) rnk <- length(dim(x)) if(rnk == 0L) stop("invalid table 'x'") @@ -121,11 +122,7 @@ plot.table <- plot(x0, unclass(x), type = type, ylim = ylim, xlab = xlab, ylab = ylab, frame.plot = frame.plot, lwd = lwd, ..., xaxt = &q...