Displaying 5 results from an estimated 5 matches for "mytrain".
2006 Apr 27
1
scope of variable/object ?
...ndalone, but doesn't work if it's divided into two functions.
The code that works is:
################### WORKS ###############
library(pamr)
set.seed(120)
x <- matrix(rnorm(1000*20),ncol=20)
y <- sample(c(1:4),size=20,replace=TRUE)
mydata <- list(x=x,y=y)
mytrain <- pamr.train(mydata)
new.scales <- pamr.adaptthresh(mytrain,ntries = 1)
########################################
But if I split the lines into two functions, then I get an error message that reads :
'Error in pamr.train(data = myd2, threshold = threshold, threshold.scale...
2006 Apr 27
0
pamr package: pamr.adaptthresh() error rates
...urn(list(mat,actclass))
}
m <- makeColon()
mat <- m[[1]]
actclass <- m[[2]]
mat <- matrix(as.numeric(mat),nrow(mat),ncol(mat))
geneid = as.character(1:nrow(mat))
gs = as.character(1:nrow(mat))
mydata <- list(x= mat,y=factor(actclass),geneid = geneid ,genenames=gs)
mytrain <- pamr.train(mydata)
new.scales <- pamr.adaptthresh(mytrain,ntries = 10, reduction.factor = 0.9) mytrain2 <- pamr.train(mydata,threshold.scale = new.scales)
mycv <- pamr.cv(mytrain2,mydata,nfold = 10)
res1 <- pamr.confusion(mycv, threshold = mytrain2$th...
2006 Apr 27
0
package pamr: pamr.adaptthresh() error rates
...st(mat,actclass)) }
m <- makeColon() mat <- m[[1]] actclass <- m[[2]] mat <- matrix(as.numeric(mat),nrow(mat),ncol(mat)) geneid = as.character(1:nrow(mat)) gs = as.character(1:nrow(mat)) mydata <- list(x= mat,y=factor(actclass),geneid = geneid ,genenames=gs)
mytrain <- pamr.train(mydata) new.scales <- pamr.adaptthresh(mytrain,ntries = 10, reduction.factor = 0.9)
mytrain2 <- pamr.train(mydata,threshold.scale = new.scales) mycv <- pamr.cv(mytrain2,mydata,nfold = 10) res1 <- pamr.confusion(mycv, threshold = mytrain2$threshold.scale,ext...
2006 Apr 27
0
package pamr: pamr.adapthresh() ---- Take 2!
...urn(list(mat,actclass))
}
m <- makeColon()
mat <- m[[1]]
actclass <- m[[2]]
mat <- matrix(as.numeric(mat),nrow(mat),ncol(mat))
geneid = as.character(1:nrow(mat))
gs = as.character(1:nrow(mat))
mydata <- list(x= mat,y=factor(actclass),geneid = geneid ,genenames=gs)
mytrain <- pamr.train(mydata)
new.scales <- pamr.adaptthresh(mytrain,ntries = 10, reduction.factor = 0.9) mytrain2 <- pamr.train(mydata,threshold.scale = new.scales)
mycv <- pamr.cv(mytrain2,mydata,nfold = 10)
res1 <- pamr.confusion(mycv, threshold = mytrain2$th...
2007 Dec 07
1
low level plotting question on R
Dear List,
I am trying to modify the xlab and ylab for a current figure that was
plotted by a package, I searched through the low level plotting command and
they do not seem to contain how to do this (the only way is to use xlab,
ylab as arguments in "plot" command, which I can not do since the plot is
plotted using some other package, not by my own script). Is there any
command for