# NOT RUN {
data(rland)
data(landscape_change)
data(param1)
#First, using simulate graph, simulate the occupation on a dynamic landscape
#(output of span.graph):
sim1 <- simulate_graph( rl=rland, rlist=landscape_change, simulate.start=TRUE,
method="percentage", parm=50, nsew="none", succ = "none",
param_df=param1, kern="op1", conn="op1", colnz="op1",
ext="op1", beta1=NULL, b=1, c1=NULL, c2=NULL, z=NULL, R=NULL)
#Then evaluate species occupancy through the changes suffered by the landscape:
occ <- list.stats(sim_list=sim1, stat="occupation", plotG=TRUE)
#Checking the percentage of occupation in the 40 first landscapes:
head(occ,40)
#Output:
#[1] 50.000000 65.000000 90.000000 96.666667 93.333333 91.666667
#[7] 91.666667 90.000000 93.333333 90.000000 85.000000 83.333333
#[13] 85.000000 88.333333 83.333333 86.666667 81.666667 68.333333
#[19] 70.000000 75.000000 80.000000 73.333333 63.333333 56.666667
#[25] 55.000000 51.666667 46.666667 41.666667 38.333333 21.666667
#[31] 13.333333 13.333333 10.000000 6.666667 5.000000 3.389831
#[37] 1.694915 1.694915 0.000000 0.000000
# }
Run the code above in your browser using DataLab