Running 'test_SongEvo_method.R' [188s]
Running the tests in 'tests/test_SongEvo_method.R' failed.
Complete output:
  > # Basic tests for SongEvo correctness
  > 
  > library(SongEvo)
  > # glo.parms=glo.parms[ ! names(glo.parms) %in% "global.parms$learning.error.d" ]
  > reg.count=40
  > 
  > # Make starting population for prbo in 1969
  > local({
  +   data(song.data)
  +   starting.trait.obs <- subset(song.data, Population %in% "PRBO" & Year %in% 1969)$Trill.FBW
  +   
  +   # Generate additional individuals to reach number of territories for this vignette (40)
  +   starting.trait.gen <- rnorm(reg.count-length(starting.trait.obs), 
  +                               mean=mean(starting.trait.obs), sd=sd(starting.trait.obs))
  +   starting.trait <- c(starting.trait.obs, starting.trait.gen)
  +   
  +   # Build data frame of territory id, age, and song trait value (in this case trill bandwidth)
  +   init.inds <- data.frame(id = 1:reg.count, age = 2, trait = starting.trait,
  +                           row.names=paste0("A",formatC(1:reg.count, width=2, flag = "0")))
  +   
  +   # Generate lat long coordinates for center? of each territory
  +   init.inds$x1 <-  round(runif(reg.count, min=-122.481858, max=-122.447270), digits=8)
  +   init.inds$y1 <-  round(runif(reg.count, min=37.787768, max=37.805645), digits=8)
  +   init.inds
  + }) -> trill.prbo.1969
  > 
  > # Make starting population for Schooner in 1969
  > local({
  +   data(song.data)
  +   starting.trait.obs <- subset(song.data, Population %in% "Schooner" & Year %in% 1969)$Trill.FBW
  +   
  +   # Generate additional individuals to reach number of territories for this vignette (40)
  +   starting.trait.gen <- rnorm(reg.count-length(starting.trait.obs), 
  +                               mean=mean(starting.trait.obs), sd=sd(starting.trait.obs))
  +   starting.trait <- c(starting.trait.obs, starting.trait.gen)
  +   
  +   # Build data frame of territory id, age, and song trait value (in this case trill bandwidth)
  +   init.inds <- data.frame(id = 1:reg.count, age = 2, trait = starting.trait,
  +                           row.names=paste0("A",formatC(1:reg.count, width=2, flag = "0")))
  +   
  +   # Generate lat long coordinates for center? of each territory
  +   init.inds$x1 <-  round(runif(reg.count, min=-122.481858, max=-122.447270), digits=8)
  +   init.inds$y1 <-  round(runif(reg.count, min=37.787768, max=37.805645), digits=8)
  +   init.inds
  + }) -> trill.sb.1969
  > 
  > # Make starting population for Bear Valley in 1969
  > local({
  +   data(song.data)
  +   starting.trait.obs <- subset(song.data, Population %in% "Bear Valley" & Year %in% 1969)$Trill.FBW
  +   
  +   # Generate additional individuals to reach number of territories for this vignette (40)
  +   starting.trait.gen <- rnorm(reg.count-length(starting.trait.obs), 
  +                               mean=mean(starting.trait.obs), sd=sd(starting.trait.obs))
  +   starting.trait <- c(starting.trait.obs, starting.trait.gen)
  +   
  +   # Build data frame of territory id, age, and song trait value (in this case trill bandwidth)
  +   init.inds <- data.frame(id = 1:reg.count, age = 2, trait = starting.trait,
  +                           row.names=paste0("A",formatC(1:reg.count, width=2, flag = "0")))
  +   
  +   # Generate lat long coordinates for center? of each territory
  +   init.inds$x1 <-  round(runif(reg.count, min=-122.481858, max=-122.447270), digits=8)
  +   init.inds$y1 <-  round(runif(reg.count, min=37.787768, max=37.805645), digits=8)
  +   init.inds
  + }) -> trill.bear.1969
  > 
  > # n.territories = 40, 
  > # lifespan = 2.083333, 
  >      
  > 
  > # Default arguments:
  > fledge_counts=c(1, 1, 2, 1, 0, 2, 2, 2, 2, 1, 1, 2, 1, 
  +                 1, 1, 1, 1, 1, 0, 2, 1, 2, 1, 2, 2, 2, 
  +                 1, 1, 1, 2, 1, 1, 
  +                 1, 2, 2, 1, 1, 2, 1, 1)
  > 
  > default.parms <- list(init.inds = NULL, iteration = 10, steps = NULL,  # years / timestep
  +                       timestep = 1, #steps in years by default
  +                       n.territories = reg.count, terr.turnover = 0.5, 
  +                       learning.method = "integrate", integrate.dist = 0.1, 
  +                       learning.error.d = 0, learning.error.sd = 430, 
  +                       mortality.a = 0.468, mortality.j = 0.5, lifespan = NA, 
  +                       phys.lim.min = 1559, phys.lim.max = 4364, 
  +                       male.fledge.n.mean = 1.35, male.fledge.n.sd = 0.5,
  +                       male.fledge.n = fledge_counts, disp.age = 2, 
  +                       disp.distance.mean = 110, disp.distance.sd = 100,
  +                       mate.comp = FALSE, prin = FALSE, all = TRUE)
  > 
  > str(default.parms)
  List of 24
   $ init.inds         : NULL
   $ iteration         : num 10
   $ steps             : NULL
   $ timestep          : num 1
   $ n.territories     : num 40
   $ terr.turnover     : num 0.5
   $ learning.method   : chr "integrate"
   $ integrate.dist    : num 0.1
   $ learning.error.d  : num 0
   $ learning.error.sd : num 430
   $ mortality.a       : num 0.468
   $ mortality.j       : num 0.5
   $ lifespan          : logi NA
   $ phys.lim.min      : num 1559
   $ phys.lim.max      : num 4364
   $ male.fledge.n.mean: num 1.35
   $ male.fledge.n.sd  : num 0.5
   $ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
   $ disp.age          : num 2
   $ disp.distance.mean: num 110
   $ disp.distance.sd  : num 100
   $ mate.comp         : logi FALSE
   $ prin              : logi FALSE
   $ all               : logi TRUE
  > 
  > # SongEvo args for prbo 1969
  > trill.prbo.test=list(parms=default.parms)
  > trill.prbo.test$parms$init.inds = trill.prbo.1969 
  > trill.prbo.test$parms$steps = 2005-1969
  > 
  > # SongEvo args for Schooner Bay 1969
  > trill.sb.test=list(parms=default.parms)
  > trill.sb.test$parms$init.inds = trill.sb.1969 
  > trill.sb.test$parms$steps = 2005-1969
  > 
  > # SongEvo args for Bear Valley 1969
  > trill.bear.test=list(parms=default.parms)
  > trill.bear.test$parms$init.inds = trill.bear.1969 
  > trill.bear.test$parms$steps = 2005-1969
  > 
  > str(trill.prbo.test)
  List of 1
   $ parms:List of 24
    ..$ init.inds         :'data.frame':	40 obs. of  5 variables:
    .. ..$ id   : int [1:40] 1 2 3 4 5 6 7 8 9 10 ...
    .. ..$ age  : num [1:40] 2 2 2 2 2 2 2 2 2 2 ...
    .. ..$ trait: num [1:40] 4005 3765 3237 3621 3285 ...
    .. ..$ x1   : num [1:40] -122 -122 -122 -122 -122 ...
    .. ..$ y1   : num [1:40] 37.8 37.8 37.8 37.8 37.8 ...
    ..$ iteration         : num 10
    ..$ steps             : num 36
    ..$ timestep          : num 1
    ..$ n.territories     : num 40
    ..$ terr.turnover     : num 0.5
    ..$ learning.method   : chr "integrate"
    ..$ integrate.dist    : num 0.1
    ..$ learning.error.d  : num 0
    ..$ learning.error.sd : num 430
    ..$ mortality.a       : num 0.468
    ..$ mortality.j       : num 0.5
    ..$ lifespan          : logi NA
    ..$ phys.lim.min      : num 1559
    ..$ phys.lim.max      : num 4364
    ..$ male.fledge.n.mean: num 1.35
    ..$ male.fledge.n.sd  : num 0.5
    ..$ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
    ..$ disp.age          : num 2
    ..$ disp.distance.mean: num 110
    ..$ disp.distance.sd  : num 100
    ..$ mate.comp         : logi FALSE
    ..$ prin              : logi FALSE
    ..$ all               : logi TRUE
  > str(trill.sb.test)
  List of 1
   $ parms:List of 24
    ..$ init.inds         :'data.frame':	40 obs. of  5 variables:
    .. ..$ id   : int [1:40] 1 2 3 4 5 6 7 8 9 10 ...
    .. ..$ age  : num [1:40] 2 2 2 2 2 2 2 2 2 2 ...
    .. ..$ trait: num [1:40] 3261 2494 2806 2878 2758 ...
    .. ..$ x1   : num [1:40] -122 -122 -122 -122 -122 ...
    .. ..$ y1   : num [1:40] 37.8 37.8 37.8 37.8 37.8 ...
    ..$ iteration         : num 10
    ..$ steps             : num 36
    ..$ timestep          : num 1
    ..$ n.territories     : num 40
    ..$ terr.turnover     : num 0.5
    ..$ learning.method   : chr "integrate"
    ..$ integrate.dist    : num 0.1
    ..$ learning.error.d  : num 0
    ..$ learning.error.sd : num 430
    ..$ mortality.a       : num 0.468
    ..$ mortality.j       : num 0.5
    ..$ lifespan          : logi NA
    ..$ phys.lim.min      : num 1559
    ..$ phys.lim.max      : num 4364
    ..$ male.fledge.n.mean: num 1.35
    ..$ male.fledge.n.sd  : num 0.5
    ..$ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
    ..$ disp.age          : num 2
    ..$ disp.distance.mean: num 110
    ..$ disp.distance.sd  : num 100
    ..$ mate.comp         : logi FALSE
    ..$ prin              : logi FALSE
    ..$ all               : logi TRUE
  > str(trill.bear.test)
  List of 1
   $ parms:List of 24
    ..$ init.inds         :'data.frame':	40 obs. of  5 variables:
    .. ..$ id   : int [1:40] 1 2 3 4 5 6 7 8 9 10 ...
    .. ..$ age  : num [1:40] 2 2 2 2 2 2 2 2 2 2 ...
    .. ..$ trait: num [1:40] 2278 2182 2926 2614 2686 ...
    .. ..$ x1   : num [1:40] -122 -122 -122 -122 -122 ...
    .. ..$ y1   : num [1:40] 37.8 37.8 37.8 37.8 37.8 ...
    ..$ iteration         : num 10
    ..$ steps             : num 36
    ..$ timestep          : num 1
    ..$ n.territories     : num 40
    ..$ terr.turnover     : num 0.5
    ..$ learning.method   : chr "integrate"
    ..$ integrate.dist    : num 0.1
    ..$ learning.error.d  : num 0
    ..$ learning.error.sd : num 430
    ..$ mortality.a       : num 0.468
    ..$ mortality.j       : num 0.5
    ..$ lifespan          : logi NA
    ..$ phys.lim.min      : num 1559
    ..$ phys.lim.max      : num 4364
    ..$ male.fledge.n.mean: num 1.35
    ..$ male.fledge.n.sd  : num 0.5
    ..$ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
    ..$ disp.age          : num 2
    ..$ disp.distance.mean: num 110
    ..$ disp.distance.sd  : num 100
    ..$ mate.comp         : logi FALSE
    ..$ prin              : logi FALSE
    ..$ all               : logi TRUE
  > 
  > # Execution tests
  > trill.prbo.test$res=do.call(SongEvo,trill.prbo.test$parms)
  > trill.sb.test$res=  do.call(SongEvo,trill.sb.test$parms)
  > trill.bear.test$res=do.call(SongEvo,trill.bear.test$parms)
  > 
  > str(trill.prbo.test)
  List of 2
   $ parms:List of 24
    ..$ init.inds         :'data.frame':	40 obs. of  5 variables:
    .. ..$ id   : int [1:40] 1 2 3 4 5 6 7 8 9 10 ...
    .. ..$ age  : num [1:40] 2 2 2 2 2 2 2 2 2 2 ...
    .. ..$ trait: num [1:40] 4005 3765 3237 3621 3285 ...
    .. ..$ x1   : num [1:40] -122 -122 -122 -122 -122 ...
    .. ..$ y1   : num [1:40] 37.8 37.8 37.8 37.8 37.8 ...
    ..$ iteration         : num 10
    ..$ steps             : num 36
    ..$ timestep          : num 1
    ..$ n.territories     : num 40
    ..$ terr.turnover     : num 0.5
    ..$ learning.method   : chr "integrate"
    ..$ integrate.dist    : num 0.1
    ..$ learning.error.d  : num 0
    ..$ learning.error.sd : num 430
    ..$ mortality.a       : num 0.468
    ..$ mortality.j       : num 0.5
    ..$ lifespan          : logi NA
    ..$ phys.lim.min      : num 1559
    ..$ phys.lim.max      : num 4364
    ..$ male.fledge.n.mean: num 1.35
    ..$ male.fledge.n.sd  : num 0.5
    ..$ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
    ..$ disp.age          : num 2
    ..$ disp.distance.mean: num 110
    ..$ disp.distance.sd  : num 100
    ..$ mate.comp         : logi FALSE
    ..$ prin              : logi FALSE
    ..$ all               : logi TRUE
   $ res  :List of 4
    ..$ summary.results: num [1:10, 1:36, 1:5] 56 53 43 52 50 52 49 46 49 49 ...
    .. ..- attr(*, "dimnames")=List of 3
    .. .. ..$ iteration: chr [1:10] "iteration 1" "iteration 2" "iteration 3" "iteration 4" ...
    .. .. ..$ step     : chr [1:36] "1" "2" "3" "4" ...
    .. .. ..$ feature  : chr [1:5] "sample.n" "trait.pop.mean" "trait.pop.variance" "lci" ...
    ..$ inds           :Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
    .. .. ..@ data       :'data.frame':	50 obs. of  10 variables:
    .. .. .. ..$ id             : int [1:50] 1530 1555 1647 1705 1750 1764 1772 1821 1826 1827 ...
    .. .. .. ..$ age            : num [1:50] 9 9 7 6 5 5 5 4 4 4 ...
    .. .. .. ..$ trait          : num [1:50] 1953 4364 2702 3997 3149 ...
    .. .. .. ..$ x1             : num [1:50] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y1             : num [1:50] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. .. ..$ male.fledglings: num [1:50] 2 0 0 0 1 1 0 1 1 2 ...
    .. .. .. ..$ territory      : num [1:50] 1 0 0 0 1 1 0 1 1 1 ...
    .. .. .. ..$ father         : num [1:50] 1407 1471 1473 1586 1589 ...
    .. .. .. ..$ x0             : num [1:50] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y0             : num [1:50] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. ..@ coords.nrs : int [1:2] 11 12
    .. .. ..@ coords     : num [1:50, 1:2] -122 -122 -122 -122 -122 ...
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:50] "1530" "1555" "1647" "1705" ...
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. ..@ bbox       : num [1:2, 1:2] -122.5 37.8 -122.4 37.8
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. .. .. ..$ : chr [1:2] "min" "max"
    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
    .. .. .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84"
    ..$ all.inds       :Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
    .. .. ..@ data       :'data.frame':	19975 obs. of  12 variables:
    .. .. .. ..$ id             : int [1:19975] 1 2 3 4 5 6 7 8 9 10 ...
    .. .. .. ..$ age            : num [1:19975] 2 2 2 2 2 2 2 2 2 2 ...
    .. .. .. ..$ trait          : num [1:19975] 4005 3765 3237 3621 3285 ...
    .. .. .. ..$ x1             : num [1:19975] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y1             : num [1:19975] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. .. ..$ male.fledglings: num [1:19975] 1 1 2 1 0 2 2 2 2 1 ...
    .. .. .. ..$ territory      : num [1:19975] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. .. ..$ father         : num [1:19975] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ x0             : num [1:19975] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ y0             : num [1:19975] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ timestep       : int [1:19975] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. .. ..$ iteration      : int [1:19975] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. ..@ coords.nrs : int [1:2] 11 12
    .. .. ..@ coords     : num [1:19975, 1:2] -122 -122 -122 -122 -122 ...
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:19975] "I1.T1.A01" "I1.T1.A02" "I1.T1.A03" "I1.T1.A04" ...
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. ..@ bbox       : num [1:2, 1:2] -122.5 37.8 -122.4 37.8
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. .. .. ..$ : chr [1:2] "min" "max"
    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
    .. .. .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84"
    ..$ time           : 'proc_time' Named num [1:5] 3.64 0.1 3.74 NA NA
    .. ..- attr(*, "names")= chr [1:5] "user.self" "sys.self" "elapsed" "user.child" ...
  > str(trill.sb.test)
  List of 2
   $ parms:List of 24
    ..$ init.inds         :'data.frame':	40 obs. of  5 variables:
    .. ..$ id   : int [1:40] 1 2 3 4 5 6 7 8 9 10 ...
    .. ..$ age  : num [1:40] 2 2 2 2 2 2 2 2 2 2 ...
    .. ..$ trait: num [1:40] 3261 2494 2806 2878 2758 ...
    .. ..$ x1   : num [1:40] -122 -122 -122 -122 -122 ...
    .. ..$ y1   : num [1:40] 37.8 37.8 37.8 37.8 37.8 ...
    ..$ iteration         : num 10
    ..$ steps             : num 36
    ..$ timestep          : num 1
    ..$ n.territories     : num 40
    ..$ terr.turnover     : num 0.5
    ..$ learning.method   : chr "integrate"
    ..$ integrate.dist    : num 0.1
    ..$ learning.error.d  : num 0
    ..$ learning.error.sd : num 430
    ..$ mortality.a       : num 0.468
    ..$ mortality.j       : num 0.5
    ..$ lifespan          : logi NA
    ..$ phys.lim.min      : num 1559
    ..$ phys.lim.max      : num 4364
    ..$ male.fledge.n.mean: num 1.35
    ..$ male.fledge.n.sd  : num 0.5
    ..$ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
    ..$ disp.age          : num 2
    ..$ disp.distance.mean: num 110
    ..$ disp.distance.sd  : num 100
    ..$ mate.comp         : logi FALSE
    ..$ prin              : logi FALSE
    ..$ all               : logi TRUE
   $ res  :List of 4
    ..$ summary.results: num [1:10, 1:36, 1:5] 47 51 54 50 56 50 43 48 54 45 ...
    .. ..- attr(*, "dimnames")=List of 3
    .. .. ..$ iteration: chr [1:10] "iteration 1" "iteration 2" "iteration 3" "iteration 4" ...
    .. .. ..$ step     : chr [1:36] "1" "2" "3" "4" ...
    .. .. ..$ feature  : chr [1:5] "sample.n" "trait.pop.mean" "trait.pop.variance" "lci" ...
    ..$ inds           :Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
    .. .. ..@ data       :'data.frame':	67 obs. of  10 variables:
    .. .. .. ..$ id             : int [1:67] 1601 1635 1693 1745 1749 1756 1759 1773 1774 1775 ...
    .. .. .. ..$ age            : num [1:67] 8 7 6 5 5 5 5 4 4 4 ...
    .. .. .. ..$ trait          : num [1:67] 2710 2111 1835 1559 1559 ...
    .. .. .. ..$ x1             : num [1:67] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y1             : num [1:67] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. .. ..$ male.fledglings: num [1:67] 1 0 2 0 0 0 0 0 3 3 ...
    .. .. .. ..$ territory      : num [1:67] 1 0 1 0 0 0 0 0 1 1 ...
    .. .. .. ..$ father         : num [1:67] 1520 1520 1593 1652 1671 ...
    .. .. .. ..$ x0             : num [1:67] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y0             : num [1:67] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. ..@ coords.nrs : int [1:2] 11 12
    .. .. ..@ coords     : num [1:67, 1:2] -122 -122 -122 -122 -122 ...
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:67] "1601" "1635" "1693" "1745" ...
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. ..@ bbox       : num [1:2, 1:2] -122.5 37.8 -122.4 37.8
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. .. .. ..$ : chr [1:2] "min" "max"
    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
    .. .. .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84"
    ..$ all.inds       :Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
    .. .. ..@ data       :'data.frame':	19993 obs. of  12 variables:
    .. .. .. ..$ id             : int [1:19993] 1 2 3 4 5 6 7 8 9 10 ...
    .. .. .. ..$ age            : num [1:19993] 2 2 2 2 2 2 2 2 2 2 ...
    .. .. .. ..$ trait          : num [1:19993] 3261 2494 2806 2878 2758 ...
    .. .. .. ..$ x1             : num [1:19993] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y1             : num [1:19993] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. .. ..$ male.fledglings: num [1:19993] 1 1 2 1 0 2 2 2 2 1 ...
    .. .. .. ..$ territory      : num [1:19993] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. .. ..$ father         : num [1:19993] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ x0             : num [1:19993] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ y0             : num [1:19993] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ timestep       : int [1:19993] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. .. ..$ iteration      : int [1:19993] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. ..@ coords.nrs : int [1:2] 11 12
    .. .. ..@ coords     : num [1:19993, 1:2] -122 -122 -122 -122 -122 ...
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:19993] "I1.T1.A01" "I1.T1.A02" "I1.T1.A03" "I1.T1.A04" ...
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. ..@ bbox       : num [1:2, 1:2] -122.5 37.8 -122.4 37.8
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. .. .. ..$ : chr [1:2] "min" "max"
    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
    .. .. .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84"
    ..$ time           : 'proc_time' Named num [1:5] 3.53 0.09 3.62 NA NA
    .. ..- attr(*, "names")= chr [1:5] "user.self" "sys.self" "elapsed" "user.child" ...
  > str(trill.bear.test)
  List of 2
   $ parms:List of 24
    ..$ init.inds         :'data.frame':	40 obs. of  5 variables:
    .. ..$ id   : int [1:40] 1 2 3 4 5 6 7 8 9 10 ...
    .. ..$ age  : num [1:40] 2 2 2 2 2 2 2 2 2 2 ...
    .. ..$ trait: num [1:40] 2278 2182 2926 2614 2686 ...
    .. ..$ x1   : num [1:40] -122 -122 -122 -122 -122 ...
    .. ..$ y1   : num [1:40] 37.8 37.8 37.8 37.8 37.8 ...
    ..$ iteration         : num 10
    ..$ steps             : num 36
    ..$ timestep          : num 1
    ..$ n.territories     : num 40
    ..$ terr.turnover     : num 0.5
    ..$ learning.method   : chr "integrate"
    ..$ integrate.dist    : num 0.1
    ..$ learning.error.d  : num 0
    ..$ learning.error.sd : num 430
    ..$ mortality.a       : num 0.468
    ..$ mortality.j       : num 0.5
    ..$ lifespan          : logi NA
    ..$ phys.lim.min      : num 1559
    ..$ phys.lim.max      : num 4364
    ..$ male.fledge.n.mean: num 1.35
    ..$ male.fledge.n.sd  : num 0.5
    ..$ male.fledge.n     : num [1:40] 1 1 2 1 0 2 2 2 2 1 ...
    ..$ disp.age          : num 2
    ..$ disp.distance.mean: num 110
    ..$ disp.distance.sd  : num 100
    ..$ mate.comp         : logi FALSE
    ..$ prin              : logi FALSE
    ..$ all               : logi TRUE
   $ res  :List of 4
    ..$ summary.results: num [1:10, 1:36, 1:5] 40 39 41 53 47 46 44 52 54 53 ...
    .. ..- attr(*, "dimnames")=List of 3
    .. .. ..$ iteration: chr [1:10] "iteration 1" "iteration 2" "iteration 3" "iteration 4" ...
    .. .. ..$ step     : chr [1:36] "1" "2" "3" "4" ...
    .. .. ..$ feature  : chr [1:5] "sample.n" "trait.pop.mean" "trait.pop.variance" "lci" ...
    ..$ inds           :Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
    .. .. ..@ data       :'data.frame':	50 obs. of  10 variables:
    .. .. .. ..$ id             : int [1:50] 1547 1600 1699 1726 1769 1778 1782 1784 1811 1812 ...
    .. .. .. ..$ age            : num [1:50] 9 8 6 6 5 5 5 5 4 4 ...
    .. .. .. ..$ trait          : num [1:50] 1559 3767 2832 2343 2354 ...
    .. .. .. ..$ x1             : num [1:50] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y1             : num [1:50] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. .. ..$ male.fledglings: num [1:50] 0 0 0 2 0 1 1 0 0 2 ...
    .. .. .. ..$ territory      : num [1:50] 0 0 0 1 0 1 1 0 0 1 ...
    .. .. .. ..$ father         : num [1:50] 1450 1483 1532 1646 1648 ...
    .. .. .. ..$ x0             : num [1:50] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y0             : num [1:50] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. ..@ coords.nrs : int [1:2] 11 12
    .. .. ..@ coords     : num [1:50, 1:2] -122 -122 -122 -122 -122 ...
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:50] "1547" "1600" "1699" "1726" ...
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. ..@ bbox       : num [1:2, 1:2] -122.5 37.8 -122.4 37.8
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. .. .. ..$ : chr [1:2] "min" "max"
    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
    .. .. .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84"
    ..$ all.inds       :Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
    .. .. ..@ data       :'data.frame':	19894 obs. of  12 variables:
    .. .. .. ..$ id             : int [1:19894] 1 2 3 4 5 6 7 8 9 10 ...
    .. .. .. ..$ age            : num [1:19894] 2 2 2 2 2 2 2 2 2 2 ...
    .. .. .. ..$ trait          : num [1:19894] 2278 2182 2926 2614 2686 ...
    .. .. .. ..$ x1             : num [1:19894] -122 -122 -122 -122 -122 ...
    .. .. .. ..$ y1             : num [1:19894] 37.8 37.8 37.8 37.8 37.8 ...
    .. .. .. ..$ male.fledglings: num [1:19894] 1 1 2 1 0 2 2 2 2 1 ...
    .. .. .. ..$ territory      : num [1:19894] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. .. ..$ father         : num [1:19894] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ x0             : num [1:19894] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ y0             : num [1:19894] 0 0 0 0 0 0 0 0 0 0 ...
    .. .. .. ..$ timestep       : int [1:19894] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. .. ..$ iteration      : int [1:19894] 1 1 1 1 1 1 1 1 1 1 ...
    .. .. ..@ coords.nrs : int [1:2] 11 12
    .. .. ..@ coords     : num [1:19894, 1:2] -122 -122 -122 -122 -122 ...
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:19894] "I1.T1.A01" "I1.T1.A02" "I1.T1.A03" "I1.T1.A04" ...
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. ..@ bbox       : num [1:2, 1:2] -122.5 37.8 -122.4 37.8
    .. .. .. ..- attr(*, "dimnames")=List of 2
    .. .. .. .. ..$ : chr [1:2] "x" "y"
    .. .. .. .. ..$ : chr [1:2] "min" "max"
    .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
    .. .. .. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +datum=WGS84"
    ..$ time           : 'proc_time' Named num [1:5] 3.22 0.11 3.33 NA NA
    .. ..- attr(*, "names")= chr [1:5] "user.self" "sys.self" "elapsed" "user.child" ...
  > 
  > 
  > # Save results if you wish to do explore the results more:
  > # save(trill.prbo.test,trill.sb.test,trill.bear.test,
  > #      file=(most_recent_SongEvo_test<-paste0("test_SongEvo",format(Sys.time(),format="%Y%m%d%H%M"),".RData")))
  > 
  > 
  > 
  > proc.time()
     user  system elapsed 
    11.12    0.50   11.60