[R] automate multiple object creation
waltzmiester
cwalte03 at shepherd.edu
Fri Jul 31 15:40:14 CEST 2009
I appreciate all the help I have received from this list and also not being
flamed because I am new to R. Many of my problems are in automation so far.
I am trying to create multiple objects that are outputs of functions. Here
are the tasks:
aGAM<-somers2(Pred_pres_a_indpdt[,3,,],population[,23])
aGBM<-somers2(Pred_pres_a_indpdt[,4,,],population[,23])
cGLM<-somers2(Pred_pres_a_indpdt[,5,,],population[,23])
cRF<-somers2(Pred_pres_a_indpdt[,8,,],population[,23])
bGAM<-somers2(Pred_pres_b_indpdt[,3,,],population[,24])
bGBM<-somers2(Pred_pres_b_indpdt[,4,,],population[,24])
bGLM<-somers2(Pred_pres_b_indpdt[,5,,],population[,24])
bRF<-somers2(Pred_pres_b_indpdt[,8,,],population[,24])
...and so on through the letter f.
so the variables (x,y,z in this example) are as follows:
xGAM<-somers2(Pred_pres_x_indpdt[,y,,],population[,z])
the GAM, GBM, GLM, and RF are stored in a column y in "Pred_pres_x_indpdt."
x is denoting a species, which corresponds to column z in "population," so
a:f == 23:28
~~~~~~~
This may be pushing it, but it would be great if I could get all of these in
a matrix of dim[1,96]
Here is what I have so far:
x <- c("a","b","c","d","e","f")
nums <- c(23:28)
rank<-function(x){
for(y in nums){
xGAM<-somers2(Pred_pres_x_indpdt[,3,,],population[,y])
xGBM<-somers2(Pred_pres_x_indpdt[,4,,],population[,y])
xGLM<-somers2(Pred_pres_x_indpdt[,5,,],population[,y])
xRF<-somers2(Pred_pres_x_indpdt[,8,,],population[,y])
rank(x)
}}
#for each species, there is 16 columns of output, so a total of 96 columns
is needed
x_matrix<-matrix(c(xGAM[1:2],Evaluation.results.TSS[[1]][1,3:4],xGBM[1:2],Evaluation.results.TSS[[1]][2,3:4],
xGLM[1:2],Evaluation.results.TSS[[1]][3,3:4],xRF[1:2],Evaluation.results.TSS[[1]][4,3:4]),
nrow=1,ncol=16,
dimnames=list(c(""),
c("x_gam_AUC", "x_gam_Dxy","x_gam_TSS","x_gam_Cutoff","x_gbm_AUC",
"x_gbm_Dxy","x_gbm_TSS","x_gbm_Cutoff",
"x_glm_AUC", "x_glm_Dxy","x_glm_TSS","x_glm_Cutoff","x_rf_AUC",
"x_rf_Dxy","x_rf_TSS","x_rf_Cutoff")))
Thanks again for helping me out.
Chris
--
View this message in context: http://www.nabble.com/automate-multiple-object-creation-tp24756632p24756632.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list