Hi,
I think that at this moment there is not a function to automate the 
"joining" of two envelope objects (although I may be wrong and 
probably Adrian will correct me).

However you could easily join them by hand (at least the simulated 
functions) and then construct the envelopes. The "trick" 
is  setting  the argument savefuns=TRUE. For example:


ENV_OBJ_1 = envelope(X,fun=pcf,simulate=sim_pat_list1, savefuns=TRUE)
ENV_OBJ_2 = envelope(X,fun=pcf,simulate=sim_pat_list2, savefuns=TRUE)

ENV_OBJ_10 = envelope(X,fun=pcf,simulate=sim_pat_list10, savefuns=TRUE)

# joining" the results:
ENV_OBJ <- c("ENV_OBJ_1", ENV_OBJ_2", ......, "ENV_OBJ_10")

simfun <-ENV_OBJ_1$r
for ( i in ENV_OBJ) {
    simfun <- data.frame(simfun, data.frame(attributes(get(i))$simfuns)[,-1])
}

envelopes <-  t(apply(simfun[,-1],1, function(x) range(x)))



HTH,

Marcelino

At 11:03 12/09/2011, Quets Jan wrote:
>Hi,
>
>is it possible to join several envelope (spatstat) objects into a 
>single object?
>
>The reason would be saving computation time with parallel computing.
>
>example:
>
>suppose I already have a list of 120 simulated patterns, called sim_pat_list.
>
>and suppose this is what I want:
>
>ENV_OBJ = envelope(X,fun=pcf,simulate=sim_pat_list)
>
>
>but suppose I have 12 cores available.
>
>and that I split the list sim_pat_list into 10 lists, each 
>comprising 12 patterns, and calculate herefrom 10 envelope objects:
>
>ENV_OBJ_1 = envelope(X,fun=pcf,simulate=sim_pat_list1)
>ENV_OBJ_2 = envelope(X,fun=pcf,simulate=sim_pat_list2)
>
>ENV_OBJ_10 = envelope(X,fun=pcf,simulate=sim_pat_list10)
>
>
>is it then possible to produce something like this (pseudo-code):
>
>ENV_OBJ = join(ENV_OBJ_1, ENV_OBJ_2, ..., ENV_OBJ_10)
>
>to get the same result?
>
>Thank you,
>Jan
>_______________________________________________
>R-sig-Geo mailing list
>R-sig-Geo@r-project.org
>https://stat.ethz.ch/mailman/listinfo/r-sig-geo

____________________________________

Marcelino de la Cruz Rot
Depto. Biologia Vegetal
EUIT Agricola
Universidad Politecnica de Madrid

tel: 34 + 913365435
____________________________________  
	[[alternative HTML version deleted]]

