


Dear all,

I have a
big data frame:


str(data1)

'data.frame':   18272 obs. of  11 variables:

 $ tag  :
int  100001 100002 100003 100005 100007
100008 100009 100011 100012 100014 ...

 $ sp   :
Factor w/ 18 levels "acassp","acocar",..: 13 5 7 14 14 18 3
11 13 10 ...

 $ gx   :
num  20 10 35 68 88 63 123 115 137 136
...

 $ gy   :
num  30 25 24 1 10 40 45 25 23 45 ...

 $ d0   :
int  0 0 0 0 0 0 0 0 0 0 ...

 $ d1   :
int  395 395 395 395 395 395 395 395 395
395 ...

 $ d2   :
int  751 751 751 751 751 751 751 751 751
751 ...

 $ d3   :
int  1515 1515 1515 1515 1515 1515 1515
1515 1515 1515 ...

 $ d4   :
int  2562 2562 2562 2562 2562 2562 2562
2562 2562 2562 ...

 $ block: int 
1 1 1 1 1 1 1 1 1 1 ...

 $ treat: Factor w/ 4 levels
"I","M","N","T": 1 1 1 1 1 1 1 1 1 1
...


And I need to
do multiple selections of gx and gy for all the levels of: sp, block, treat, and when d0!="NA". Then
calculate some spatial functions with the selected gx and gy coordinates, and
save the results with a name according to the selection.
 

One single
selection could be done and named like that:

acocar1I=subset(data1,(treat=="I"
& data1$block==1 & data1$sp=="acocar" & data1$d0!="NA"))


These are some
of the functions I have to calculate:

acocar1I.spp<-spp(x=acocar1I$gx,
y=acocar1I$gy, window=wA)

acocar1I.dp<-dval(acocar1I.spp,25,2.5,18,20)


And I want
to create a 'results' object to access easily all the results:

acocar1I.res<-alist()

acocar1I.res$data1<-acocar1I

acocar1I.res$spp<-acocar1I.spp

acocar1I.res$dp<-acocar1I.dp


Can I do everything
in a single loop or a single function?

Many
thanks!











 		 	   		  
	[[alternative HTML version deleted]]

