Hi,
   
  I must be missing something here...Essentially, a short piece of code works if it's standalone, 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 = all.scales[i+ : object "myd2" not found.'
   
  The code that doesn't work is:
   
  
################### DOESN'T WORK ############
library(pamr)
     
domat <- function(myd){
      myd2 <- myd
    mytrain <-   pamr.train(myd2)
    new.scales <- pamr.adaptthresh(mytrain)
  }
  dom <- function(){
     set.seed(120)
     x <- matrix(rnorm(1000*20),ncol=20)
     y <- sample(c(1:4),size=20,replace=TRUE)
     myda <- list(x=x,y=y)
     domat(myda)
}
  dom()
   
  #########################################
   
  Did I do something really goofy? How can I find out what's happening? 
   
  many thanks.
   
   

		
---------------------------------

	[[alternative HTML version deleted]]


