Troubleshooting RSF results

[Ryan Bourbour]

Hi Group,

We have a couple questions regarding RSF mean outputs. Quick background, we classified hunting locations for 11 individual owls and ran an RSF using 9 land cover groupings. 

When we plot the individuals, all land cover groups (besides reference) are summarized. However, when we calculate and plot the mean for all owls, only 6 of the groups are reported in the results. When individuals are plotted, only 7 are shown because the null RSF model was selected for the others. 

Is there anything obvious from the outputs that we are missing that may help us adjust our approach? Any insights would be much appreciated on how to investigate why this may be happening (zeros, NAs, etc.?) or if there is a way to correct this so that all results can be reported. Below are the summaries of the 11 individuals and the RSF mean summary output. 

Another question, for the RSF mean summary, are there any special considerations for interpreting the SD?

Thanks in advance,

Ryan

> ## Individual summaries
> lapply(rsfs, summary)
[[1]]
[[1]]$name
[1] "OUF error"

[[1]]$DOF
     mean      area diffusion     speed
 31.36734  29.42258 162.60563 500.81576

[[1]]$CI
                                        low       est      high
area (square kilometers)           5.162805  7.684476 10.699671
τ[position] (hours)                1.179456  1.807255  2.769218
τ[velocity] (minutes)              2.489771  2.862069  3.290038
speed (kilometers/day)            70.644207 73.879874 77.113619
diffusion (square kilometers/day)  8.373433  9.825635 11.392250


[[2]]
[[2]]$name
[1] "OUF error"

[[2]]$DOF
     mean      area diffusion     speed
119.46328  95.09445 167.40093 459.18437

[[2]]$CI
                                         low        est      high
area (square kilometers)           26.710979  33.011182  39.96848
τ[position] (minutes)              33.110135  44.061965  58.63633
τ[velocity] (minutes)               7.482717   8.788679  10.32257
speed (kilometers/day)            130.815277 137.085363 143.35156
diffusion (square kilometers/day)  65.530751  76.711442  88.75983


[[3]]
[[3]]$name
[1] "OUF error"

[[3]]$DOF
      mean       area  diffusion      speed
 11.223394   9.950299 108.997746 396.531766

[[3]]$CI
                                        low       est      high
area (square kilometers)           3.468499  7.248719 12.398593
τ[position] (hours)                1.818217  3.956823  8.610881
τ[velocity] (minutes)              2.623711  3.109948  3.686297
speed (kilometers/day)            44.234229 46.524109 48.812465
diffusion (square kilometers/day)  3.621428  4.410447  5.276065


[[4]]
[[4]]$name
[1] "OUF error"

[[4]]$DOF
     mean      area diffusion     speed
 121.7640  113.9701  145.2447  135.9765

[[4]]$CI
                                            low         est        high
habitat.9_4 (1/habitat.9_4)       -6.748307e+02  -9.1713648  656.487931
habitat.8_4 (1/habitat.8_4)       -2.724741e+00  -0.8305975    1.063546
habitat.7_4 (1/habitat.7_4)        5.723566e-02   1.6471479    3.237060
habitat.6_4 (1/habitat.6_4)        7.900349e-01   2.3630551    3.936075
habitat.5_4 (1/habitat.5_4)       -1.816628e-01   1.7784012    3.738465
habitat.3_4 (1/habitat.3_4)       -4.394110e-01   1.6076081    3.654627
habitat.2_4 (1/habitat.2_4)       -1.920386e+03  -2.4735999 1915.439163
habitat.1_4 (1/habitat.1_4)       -7.670077e+02  -7.6326296  751.742450
area (square kilometers)           6.303925e+00   7.6424660    9.107952
τ[position] (minutes)              1.805447e+01  24.1711014   32.359978
τ[velocity] (minutes)              4.612397e+00   5.5379161    6.649150
speed (kilometers/day)             1.028275e+02 112.2614519  121.685076
diffusion (square kilometers/day)  2.635025e+01  31.2211549   36.499051


[[5]]
[[5]]$name
[1] "OUf error"

[[5]]$DOF
     mean      area diffusion     speed
 72.48112  58.20826 117.12255 272.24717

[[5]]$CI
                                        low       est      high
area (square kilometers)           25.63138  33.73694  42.93906
τ (minutes)                        17.47951  19.06928  20.80364
speed (kilometers/day)            134.54250 143.03888 151.52849
diffusion (square kilometers/day)  82.36671  99.58333 118.40960


[[6]]
[[6]]$name
[1] "OUF error"

[[6]]$DOF
     mean      area diffusion     speed
 66.01991  73.62575  88.64723  84.26881

[[6]]$CI
                                           low        est       high
habitat.9_4 (1/habitat.9_4)         -3.4475189  0.2200845   3.887688
habitat.8_4 (1/habitat.8_4)         -2.4725135  0.7352740   3.943062
habitat.7_4 (1/habitat.7_4)         -1.9269923  0.9982207   3.923434
habitat.6_4 (1/habitat.6_4)         -1.6222206  1.2718584   4.165937
habitat.5_4 (1/habitat.5_4)          0.9442078  3.8192768   6.694346
habitat.3_4 (1/habitat.3_4)         -1.7475788  1.2072161   4.162011
habitat.2_4 (1/habitat.2_4)       -124.8481042 -3.8684061 117.111292
habitat.1_4 (1/habitat.1_4)         -3.3220117  0.7056885   4.733389
area (square kilometers)             1.5618253  1.9903375   2.470009
τ[position] (minutes)               34.7441506 45.6026914  59.854837
τ[velocity] (minutes)                2.9871454  3.3940923   3.856479
speed (kilometers/day)              47.6181651 53.3072863  58.988710
diffusion (square kilometers/day)    4.3499014  5.4190019   6.603807


[[7]]
[[7]]$name
[1] "OUF error"

[[7]]$DOF
     mean      area diffusion     speed
 93.43592  90.33743  89.76983  98.29601

[[7]]$CI
                                            low         est        high
habitat.9_4 (1/habitat.9_4)          -3.3124826  -0.6797524    1.952978
habitat.8_4 (1/habitat.8_4)          -3.4769642  -0.6324431    2.212078
habitat.7_4 (1/habitat.7_4)           0.3501424   2.2897876    4.229433
habitat.6_4 (1/habitat.6_4)          -0.9896897   0.9887335    2.967157
habitat.5_4 (1/habitat.5_4)          -1.6379717   0.6147999    2.867571
habitat.3_4 (1/habitat.3_4)          -0.6885514   1.2793413    3.247234
habitat.2_4 (1/habitat.2_4)       -1123.3951884 -11.4284692 1100.538250
habitat.1_4 (1/habitat.1_4)        -710.8882016  -9.1425955  692.603011
area (square kilometers)             14.3668450  17.8589052   21.725154
τ[position] (minutes)                37.7719758  52.1807380   72.085967
τ[velocity] (minutes)                 6.5879956   7.8232633    9.290147
speed (kilometers/day)               88.5882103  98.3029786  108.005428
diffusion (square kilometers/day)    30.1229031  37.4718384   45.610861


[[8]]
[[8]]$name
[1] "OUF error"

[[8]]$DOF
     mean      area diffusion     speed
 35.60143  37.52335  28.52272  30.92778

[[8]]$CI
                                           low        est       high
habitat.9_4 (1/habitat.9_4)         -1.2586227  1.8173975   4.893418
habitat.8_4 (1/habitat.8_4)         -1.0630052  1.7071970   4.477399
habitat.7_4 (1/habitat.7_4)         -2.2946282  0.5483409   3.391310
habitat.6_4 (1/habitat.6_4)         -0.0775807  2.8220147   5.721610
habitat.5_4 (1/habitat.5_4)         -1.7825630  1.5428541   4.868271
habitat.3_4 (1/habitat.3_4)         -1.3208155  1.5665118   4.453839
habitat.2_4 (1/habitat.2_4)         -1.8426573  1.4378909   4.718439
habitat.1_4 (1/habitat.1_4)       -229.5334442 -5.6281080 218.277228
area (square kilometers)            13.2239525 18.7314537  25.182753
τ[position] (hours)                  2.8915597  4.4590909   6.876390
τ[velocity] (minutes)                3.5710344  4.1910917   4.918813
speed (kilometers/day)              50.1979936 60.9136526  71.608918
diffusion (square kilometers/day)    6.6903466 10.0266803  14.027225


[[9]]
[[9]]$name
[1] "OUF error"

[[9]]$DOF
     mean      area diffusion     speed
 75.52661  70.66023  81.04194  87.45983

[[9]]$CI
                                         low         est       high
habitat.9_4 (1/habitat.9_4)       -0.8347109 -0.03868025  0.7573504
habitat.8_4 (1/habitat.8_4)       -3.8638996 -2.13911581 -0.4143320
habitat.7_4 (1/habitat.7_4)       -2.3626817  2.27296432  6.9086104
habitat.6_4 (1/habitat.6_4)       -3.1247598  0.24462578  3.6140113
habitat.5_4 (1/habitat.5_4)        1.3109229  2.06483868  2.8187544
habitat.3_4 (1/habitat.3_4)        0.2468971  0.86837113  1.4898451
habitat.2_4 (1/habitat.2_4)       -0.4700699  0.55537013  1.5808102
habitat.1_4 (1/habitat.1_4)       -5.5925836  0.62438270  6.8413490
area (square kilometers)           7.0034108  8.97279810 11.1824881
τ[position] (minutes)             28.0001610 44.49040112 70.6923003
τ[velocity] (minutes)              8.6237864 12.07032565 16.8942916
speed (kilometers/day)            54.3960268 60.76144066 67.1183744
diffusion (square kilometers/day) 15.0787973 18.98627816 23.3371562


[[10]]
[[10]]$name
[1] "OUf error"

[[10]]$DOF
     mean      area diffusion     speed
 82.91700  58.68078 819.59102 204.89822

[[10]]$CI
                                           low         est       high
habitat.9_4 (1/habitat.9_4)          1.0890132   3.3514930   5.613973
habitat.8_4 (1/habitat.8_4)          0.5925491   2.8656071   5.138665
habitat.7_4 (1/habitat.7_4)       -664.1138054  -5.7167227 652.680360
habitat.6_4 (1/habitat.6_4)          1.8946113   4.4954825   7.096354
habitat.5_4 (1/habitat.5_4)         -2.1209454   1.0676611   4.256268
habitat.3_4 (1/habitat.3_4)         -1.3044506   1.1447119   3.593874
habitat.2_4 (1/habitat.2_4)         -3.2892713  -0.3937563   2.501759
habitat.1_4 (1/habitat.1_4)       -635.0402786  -6.0554636 622.929351
area (square kilometers)            38.0101729  49.9710794  63.543141
τ (minutes)                         18.2372518  19.5531603  20.964018
speed (kilometers/day)             158.1765682 169.8022402 181.417318
diffusion (square kilometers/day)  134.1712003 143.8523358 153.865926


[[11]]
[[11]]$name
[1] "OUF error"

[[11]]$DOF
     mean      area diffusion     speed
 58.60818  52.17899  81.71008  78.06974

[[11]]$CI
                                            low        est         high
habitat.9_4 (1/habitat.9_4)       -9.354438e+02  10.031142 9.555061e+02
habitat.8_4 (1/habitat.8_4)       -9.344216e+02  11.052218 9.565260e+02
habitat.7_4 (1/habitat.7_4)       -1.390977e+08 -18.452643 1.390977e+08
habitat.6_4 (1/habitat.6_4)       -1.304758e+08 -17.295379 1.304757e+08
habitat.5_4 (1/habitat.5_4)       -9.325173e+02  12.956507 9.584304e+02
habitat.3_4 (1/habitat.3_4)       -9.335886e+02  11.885187 9.573590e+02
habitat.2_4 (1/habitat.2_4)       -9.347882e+02  10.685851 9.561599e+02
habitat.1_4 (1/habitat.1_4)       -1.159659e+05  -8.770544 1.159483e+05
area (square kilometers)           1.115073e+01  14.922356 1.923497e+01
τ[position] (minutes)              2.903191e+01  41.769700 6.009621e+01
τ[velocity] (minutes)              7.381977e+00   9.157785 1.136078e+01
speed (kilometers/day)             8.256347e+01  92.859048 1.031403e+02
diffusion (square kilometers/day)  2.837736e+01  35.695322 4.384009e+01

> summary(rsf_mean)
$name
[1] "OUF anisotropic-μ"

$DOF
      mean       area  diffusion      speed
142.821282  13.711910   7.252066  15.134674

$CI
                                         low        est       high
habitat.9_4 (1/habitat.9_4)       -0.4325417  0.6494788  1.7314993
SD[habitat.9_4] (1/habitat.9_4)    1.1047805  1.3203281  1.5354806
habitat.8_4 (1/habitat.8_4)       -1.0277182  0.1575992  1.3429166
SD[habitat.8_4] (1/habitat.8_4)    1.4289514  1.4963600  1.5637274
habitat.7_4 (1/habitat.7_4)        0.2624142  1.0574692  1.8525243
SD[habitat.7_4] (1/habitat.7_4)    0.6569067  0.7200556  0.7831323
habitat.6_4 (1/habitat.6_4)        0.9297998  1.9029734  2.8761470
SD[habitat.6_4] (1/habitat.6_4)    1.1947199  1.2431531  1.2915607
habitat.5_4 (1/habitat.5_4)        0.9638450  1.6503354  2.3368258
SD[habitat.5_4] (1/habitat.5_4)    0.7299402  0.7828240  0.8356602
habitat.3_4 (1/habitat.3_4)        0.5359506  0.9290031  1.3220556
SD[habitat.3_4] (1/habitat.3_4)    0.0000000  0.0000000        Inf
area (square kilometers)           4.2190880  7.7736448 12.3963540
CoV[area]                          0.4401657  0.9063138  1.3761633
τ[position] (hours)                0.5867928  1.0765536  1.9750885
CoV[τ[position]]                   0.2571299  0.8966907  1.5620427
τ[velocity] (minutes)              3.3333052  5.2005950  8.1139249
CoV[τ[velocity]]                   0.4162390  0.6922150  0.9684641
speed (kilometers/day)            53.9427618 71.9970498 90.0164790
CoV[speed]                         0.1145702  0.4373971  0.7756178
diffusion (square kilometers/day)  6.0549196 14.7686834 27.2999968
CoV[diffusion]                     0.3850983  1.2684673  2.1835699

[Christen Fleming]

Hi Ryan,

In more recent versions of the package, I disabled feature selection on such parameters in the population model. You should be able to update from GitHub and get those values, but they are probably very insignificant (crossing zero).

With 9 habitat types, it might be a good idea to consolidate them to some degree.

The SD is just the population standard deviation in the individual parameters. If it's (0,Inf) then its not supported in the selected model.

Best,

Chris