Displaying 3 results from an estimated 3 matches for "pred_pres_a_indpdt".
2009 Jul 31
1
automation question
I have columns in two tables that need to be compared in order to get a
statistic. The function calculating the statistic is somers2. Here are the
two columns I must calculate, row by row.
name1 name2
1 Pred_pres_a_indpdt[,3,,] population[,23]
2 Pred_pres_b_indpdt[,3,,] population[,24]
3 Pred_pres_c_indpdt[,3,,] population[,25]
4 Pred_pres_d_indpdt[,3,,] population[,26]
5 Pred_pres_e_indpdt[,3,,] population[,27]
6 Pred_pres_f_indpdt[,3,,] population[,28]
7 Pred_pres_a_indpdt[,4,,] population[,23]
8 Pred_pres_...
2009 Jul 30
3
for loop for file names
I am trying to load binary files in the following fashion
load("pred/Pred_pres_a_indpdt")
load("pred/Pred_pres_b_indpdt")
load("pred/Pred_pres_c_indpdt")
load("pred/Pred_pres_d_indpdt")
load("pred/Pred_pres_e_indpdt")
load("pred/Pred_pres_f_indpdt")
but I would like to set up a for loop to replace the letters a:f
Here is what I...
2009 Jul 31
2
another automation question
This code works:
x<-letters[1:6]
ycols<-23:28
xcols<-rep(c(3,4,5,8),each=length(ycols))
somertime<-function(i,j)somers2(Pred_pres_a_indpdt[,i,,], population[,j])
results<-mapply(somertime,xcols,ycols)
How can I make variable "h" work?
x<-letters[1:6]
ycols<-23:28
xcols<-rep(c(3,4,5,8),each=length(ycols))
somertime<-function(h,i,j)somers2(Pred_pres_h_indpdt[,i,,], population[,j])
results<-mapply(somerti...