diff --git a/DataCleaningScripts/get_nest_success.R b/DataCleaningScripts/get_nest_success.R index 55cc2ce..4f13c61 100644 --- a/DataCleaningScripts/get_nest_success.R +++ b/DataCleaningScripts/get_nest_success.R @@ -73,6 +73,35 @@ unique(success$species[which(!(success$species %in% species$species))]) success <- success %>% dplyr::arrange(year,colony,species) write.table(success, "Nesting/nest_success.csv", row.names = FALSE, na = "", sep = ",", quote = 18) + +success_summary <- read.csv("Nesting/nest_success_summary.csv") +success_summary_new <- read.csv("../Desktop/nest_success_summary.csv") %>% + dplyr::rename_with(tolower) %>% + dplyr::mutate(metric = tolower(metric), + metric = gsub(paste(c("[(]", "[)]"), collapse = "|"), "", metric), + metric = gsub("overall", "", metric), + metric = gsub("success", "", metric), + metric = gsub(" ", "", metric), + colony = replace(colony, colony=="cuthbert", "cuthbert_lake"), + colony = replace(colony, colony=="paurotis", "paurotis_pond")) %>% + tidyr::pivot_longer(cols = !c(year,colony,type,metric), + names_to = "species", + values_to = "value") %>% + dplyr::filter(!is.na(value)) %>% + dplyr::mutate(variable = paste(type,metric,sep="_")) %>% + tidyr::pivot_wider(id_cols = c(year,colony,species), names_from = variable, values_from = value, + values_fill = NA) %>% + dplyr::mutate(year=as.integer(year)) %>% + dplyr::arrange(year,colony,species) + +unique(success_summary_new$colony[which(!(success_summary_new$colony %in% colonies$colony))]) +unique(success_summary_new$species[which(!(success_summary_new$species %in% species$species))]) + +success_summary <- success_summary %>% + dplyr::bind_rows(success_summary_new) %>% + dplyr::arrange(year,colony,species) +write.table(success_summary, "Nesting/nest_success_summary.csv", row.names = FALSE, na = "", sep = ",") + # make metadata success_metadata <- data.frame(name = as.vector(colnames(success))) %>% dplyr::mutate(class=sapply(success,class), diff --git a/DataCleaningScripts/old_counts.R b/DataCleaningScripts/old_counts.R index bbb5c5e..e6add37 100644 --- a/DataCleaningScripts/old_counts.R +++ b/DataCleaningScripts/old_counts.R @@ -59,7 +59,7 @@ counts <- counts %>% dplyr::arrange(counts) write.csv(counts, "Counts/maxcounts.csv", row.names = FALSE, na = "", quote = 9) species <- species %>% dplyr::arrange(species) -write.csv(species, "SiteandMethods/species_list.csv", row.names = FALSE, na = "", quote = 5:25) +write.csv(species, "SiteandMethods/species_list.csv", row.names = FALSE, na = "", quote = 7:27) colonies <- colonies %>% dplyr::arrange(colony) write.csv(colonies, "SiteandMethods/colonies.csv", row.names = FALSE, na = "", quote = c(7,8)) diff --git a/Nesting/README.md b/Nesting/README.md index 48cff8c..0fbf13e 100644 --- a/Nesting/README.md +++ b/Nesting/README.md @@ -3,12 +3,49 @@ Please see the [methods](../SiteandMethods/methods.md) for details about collection and how to use the data. ### Nest Checks + year + colony + nest - nest id + species + date + eggs - number of eggs + chicks - number of chicks + stage - nesting stage -notes \ No newline at end of file + +notes + +### Nest Success + +We use the Mayfield method to calculate nest success. The nest-level metrics required for calculations are stored in this table. See the [metadata](../Nesting/nest_success_metadata.csv) for data details. + +### Nest Success Summary + +This is colony-level nest success, based on the data from the [nest_success table](../Nesting/nest_success.csv). Success is summarized for the incubation stage, nestling stage, and overall for the whole nesting period. + +k - total number of nests + +sumy - total number of nests that hatched (s(i)) + +e - total number of days monitored (n(i)) + +p - 1-((k-sumy)/e), success + +j - incubation or nestling period, in days (constant by species, can be found in the [species table](../SiteandMethods/species.csv)) + +pj - p\^j, hatch success + +varp - (p\*(1-p))/e, variation of p + +varpj - varp\*((j\*(p\^(j-1)))\^2), variation of pj + +sdp - standard deviation of p + +sdpj - standard deviation of pj diff --git a/Nesting/nest_success_summary.csv b/Nesting/nest_success_summary.csv new file mode 100644 index 0000000..f16bd6f --- /dev/null +++ b/Nesting/nest_success_summary.csv @@ -0,0 +1,174 @@ +"year","colony","species","incubation_k","incubation_sumy","incubation_e","incubation_p","incubation_j","incubation_pj","incubation_varp","incubation_varpj","incubation_sdp","incubation_sdpj","nestling_k","nestling_sumy","nestling_e","nestling_p","nestling_j","nestling_pj","nestling_varp","nestling_varpj","nestling_sdp","nestling_sdpj","overall_p","overall_varp","overall_sd" +2003,"alley_north","greg",44,9,585.5,0.940222032,,0.27405574,9.59941e-05,0.003596667,,0.059972218,5,8,171,0.976608187,,0.717934875,0.000133594,0.014150559,,0.11895612,0.196754173,0.002916631,0.054005846 +2003,"alley_north","whib",50,12,527.5,0.927962085,,0.208034266,0.000126727,0.002808774,,0.052997867,10,16,151.5,0.98679868,,0.830232705,8.59871e-05,0.011929766,,0.109223467,0.172716852,0.002452349,0.049521197 +2003,"hidden","greg",33,27,802.5,0.992523364,,0.85419204,9.24702e-06,0.003020439,,0.05495852,17,29,491,0.979633401,,0.749704607,4.0635e-05,0.004664551,,0.06829752,0.640391708,0.005101121,0.071422131 +2003,"mud_east","greg",34,26,827,0.990326481,,0.815354057,1.1584e-05,0.003462834,,0.058845847,20,34,494,0.987854251,,0.842752837,2.42879e-05,0.003464663,,0.058861388,0.687141944,0.004762732,0.069012551 +2004,"alley_north","greg",65,54,1152.5,0.990455531,28,0.76450352,8.20249e-06,0.003831342,0.002863999,0.061897832,81,74,1615.5,0.995666976,21,0.912843149,2.67053e-06,0.000989922,0.001634177,0.031463024,0.697871801,0.003774958,0.06144069 +2004,"alley_north","smhe",26,12,363.5,0.961485557,22,0.421444924,0.000101874,0.009473347,0.010093249,0.097331122,,,,,,,,,,,0.040770173,0.003275985,0.05723622 +2004,"alley_north","smwh",,,,,,,,,,,3,0,28.5,0.894736842,21,0.09673903,0.003304661,0.017036417,0.057486177,0.130523626,,, +2004,"alley_north","whib",90,58,964.5,0.966822188,21,0.492355483,3.32577e-05,0.003803602,0.005766948,0.061673348,66,65,710,0.998591549,14,0.980461197,1.98094e-06,0.000374294,0.001407458,0.019346667,0.482735446,0.003748576,0.061225615 +2004,"false_l67","greg",24,18,450.5,0.986681465,28,0.686997505,2.91701e-05,0.011086924,0.00540094,0.105294463,30,29,593,0.998313659,21,0.965177703,2.83895e-06,0.001170245,0.001684918,0.034208842,0.663074674,0.010893513,0.104371995 +2004,"false_l67","smhe",0,0,0,,22,,,,,,0,0,0,,21,,,,,,,, +2004,"false_l67","whib",0,0,0,,21,,,,,,0,0,0,,14,,,,,,,, +2004,"hidden","greg",31,27,489,0.991820041,28,0.79454889,1.65911e-05,0.008347698,0.004073217,0.091365735,33,29,663,0.993966817,21,0.880662715,9.04492e-06,0.003131258,0.003007478,0.055957645,0.699729583,0.008477124,0.092071299 +2004,"hidden","smhe",10,0,80,0.875,28,0.023780747,0.001367188,0.000791733,0.036975499,0.028137759,0,0,20,1,50,1,0,0,0,0,0.023780747,0.000791733,0.028137759 +2004,"hidden","whib",23,21,375,0.994666667,21,0.893776323,1.41464e-05,0.005037165,0.003761166,0.070972988,36,36,498,1,14,1,0,0,0,0,0.893776323,0.005037165,0.070972988 +2004,"tamiami_west","wost",10,0,80,0.875,28,0.023780747,0.001367188,0.000791733,0.036975499,0.028137759,0,0,1,1,50,1,0,0,0,0,0.023780747,0.000791733,0.028137759 +2011,"joule","gbhe",1,1,4.5,1,28,1,0,0,,,2,2,53.5,1,60,1,0,0,,,1,0,0 +2011,"joule","greg",54,50,856.5,0.995329831,28,0.877155458,5.42716e-06,0.003304513,,,53,35.5,836.5,0.979079498,21,0.641471086,2.44864e-05,0.004635333,,,0.562669865,0.004941509,0.070295866 +2011,"joule","rosp",3,2,49,0.979591837,28,0.561389655,0.000407993,0.105052782,,,2,2,20,1,21,1,0,0,,,0.561389655,0.105052782,0.324118469 +2011,"joule","smhe",2,2,48,1,22,1,0,0,,,2,2,10,1,14,1,0,0,,,1,0,0 +2011,"tamiami_west","bcnh",65,43,716.5,0.969295185,28,0.417608711,4.15381e-05,0.006044904,0.006445004,0.07774898,66,47,950.5,0.980010521,21,0.654403327,2.06101e-05,0.004052735,0.004539835,0.063661095,0.27328453,0.003319975,0.057619228 +2011,"tamiami_west","greg",15,11,257,0.984435798,28,0.644535941,5.96185e-05,0.020036261,0.007721303,0.141549502,11,10,217,0.995391705,21,0.907557946,2.11385e-05,0.007749503,0.004597664,0.08803126,0.584953714,0.019877716,0.140988354 +2011,"tamiami_west","lbhe",23,7,308,0.948051948,22,0.309248049,0.000159901,0.008234661,0.01264519,0.090745032,9,8,123.5,0.991902834,14,0.892417002,6.50332e-05,0.010317838,0.008064317,0.101576758,0.275978217,0.007629854,0.087349037 +2011,"tamiami_west","sneg",28,12,412.5,0.961212121,22,0.418815989,9.03839e-05,0.008305097,0.009507047,0.091132305,13,12,179.5,0.994428969,14,0.924767898,3.08635e-05,0.005231412,0.005555491,0.072328498,0.387307582,0.008063556,0.089797303 +2011,"tamiami_west","trhe",24,9,327,0.95412844,22,0.355918457,0.000133845,0.009014368,0.011569146,0.094944027,10,9,137.5,0.992727273,14,0.90285778,5.25079e-05,0.008512545,0.007246233,0.092263453,0.321343748,0.008503169,0.092212628 +2011,"tamiami_west","whib",66,40,931,0.97207304,21,0.551666818,2.9159e-05,0.004141593,0.005399909,0.064355208,40,33,515.5,0.986420951,14,0.825794578,2.59838e-05,0.003569263,0.005097433,0.059743308,0.455563467,0.003925343,0.062652556 +2011,"tamiami_west","wost",19,14,374,0.986631016,28,0.686014651,3.52681e-05,0.013367644,0.005938692,0.115618529,14,10,428.5,0.990665111,50,0.62566724,2.15817e-05,0.021520795,0.004645608,0.146699676,0.429216893,0.015648606,0.125094389 +2011,"vacation","gbhe",6,6,104.5,1,28,1,0,0,,,6,6,183.5,1,60,1,0,0,,,1,0, +2011,"vacation","greg",39,35,773.5,0.994828701,28,0.864873463,6.65101e-06,0.003941053,,,34,32,562.5,0.996444444,21,0.927929311,6.29851e-06,0.002408797,,,0.802541436,0.005204743, +2011,"vacation","smhe",23,18,438,0.988584475,22,0.776789397,2.57653e-05,0.007699454,,,13,13,91.5,1,14,1,0,0,,,0.776789397,0.007699454, +2011,"yonteau","glib",1,1,20.5,1,21,1,0,0,,,1,0,11.5,0.913043478,21,0.148020079,0.00690392,0.080018926,,,0.148020079,0.080018926, +2011,"yonteau","greg",2,2,12,1,28,1,0,0,,,2,2,42,1,21,1,0,0,,,1,0, +2011,"yonteau","lbhe",24,9,313.5,0.95215311,22,0.340055174,0.000145319,0.008971261,,,9,8,122,0.991803279,14,0.891163838,6.66355e-05,0.0105445,,,0.303044873,0.008438671, +2011,"yonteau","sneg",27,12,450,0.966666667,22,0.474338472,7.16049e-05,0.008344716,,,14,13,164,0.993902439,14,0.917936394,3.69535e-05,0.006178027,,,0.435412546,0.008472909, +2011,"yonteau","trhe",43,28,813.5,0.981561156,22,0.664021269,2.22481e-05,0.004927971,,,28,26,384,0.994791667,14,0.929501168,1.34927e-05,0.002308827,,,0.617208545,0.005287026, +2011,"yonteau","whib",60,40,887.5,0.977464789,21,0.619617331,2.48196e-05,0.004398229,,,40,25,494.5,0.96966633,14,0.649699329,5.94814e-05,0.005233802,,,0.402564964,0.003888943, +2012,"6th_bridge","bcnh",2,2,32,1,28,1,0,0,0,0,2,1,39,0.974358974,21,0.579560504,0.000640604,0.099951185,0.02531016,0.316150574,0.579560504,0.099951185,0.316150574 +2012,"6th_bridge","greg",12,6,189.5,0.968337731,28,0.406211225,0.000161793,0.022321628,0.012719787,0.149404244,6,3,81.5,0.963190184,21,0.454938832,0.000435029,0.04279949,0.020857346,0.206880376,0.18480126,0.012637486,0.112416575 +2012,"6th_bridge","rosp",3,2,45.5,0.978021978,28,0.536736588,0.000472417,0.111549353,0.021735163,0.33399005,2,1,38,0.973684211,21,0.571190106,0.000674297,0.102332744,0.025967221,0.319894896,0.306578629,0.077289684,0.278010222 +2012,"6th_bridge","smhe",26,12,363.5,0.961485557,22,0.421444924,0.000101874,0.009473347,0.010093249,0.097331122,12,5,120,0.941666667,14,0.431082869,0.000457755,0.018802515,0.021395201,0.137122264,0.181677687,0.005352583,0.073161347 +2012,"6th_bridge","trhe",,,,,,,,,,,5,4,68,0.985294118,14,0.812685404,0.000213083,0.028413012,0.01459735,0.168561597,,, +2012,"6th_bridge","whib",16,15,291.5,0.996569468,21,0.930377353,1.17282e-05,0.004507878,0.003424642,0.067140733,15,4,99,0.888888889,14,0.192249264,0.000997631,0.009146612,0.031585291,0.095637922,0.178864361,0.008125168,0.090139713 +2012,"paurotis_pond","wost",34,19,479.5,0.968717414,28,0.410694599,6.31991e-05,0.00890576,0.007949788,0.094370333,31,0,483.5,0.935884178,50,0.036400325,0.000124105,0.00046935,0.011140262,0.021664495,0.014949417,9.51453e-05,0.009754243 +2012,"tamiami_west","bcnh",3,3,22,1,28,1,0,0,0,0,3,3,63,1,21,1,0,0,0,0,1,0,0 +2012,"tamiami_west","greg",1,1,13,1,28,1,0,0,0,0,1,1,21,1,21,1,0,0,0,0,1,0,0 +2012,"tamiami_west","smhe",28,23,411.5,0.987849332,22,0.764179932,2.9169e-05,0.008448434,0.00540083,0.091915365,23,17,201.5,0.970223325,14,0.654943683,0.000143375,0.012805408,0.011973922,0.113160983,0.500494819,0.011104707,0.105378873 +2012,"tamiami_west","trhe",,,,,,,,,,,3,3,63,1,21,1,0,0,0,0,,, +2012,"tamiami_west","whib",50,42,790.5,0.989879823,21,0.807666205,1.26727e-05,0.003720542,0.003559872,0.060996243,40,32,512.5,0.984390244,14,0.802310012,2.99826e-05,0.003903691,0.005475639,0.062479525,0.647998682,0.004955916,0.070398265 +2012,"tamiami_west","wost",12,11,183,0.994535519,28,0.857765059,2.96974e-05,0.017319298,0.00544953,0.131602804,11,0,252.5,0.956435644,50,0.107841925,0.000165016,0.005244806,0.01284585,0.072421034,0.092503035,0.004151181,0.064429659 +2012,"vacation","greg",23,12,303,0.96369637,28,0.355080026,0.000115464,0.012289544,0.010745431,0.110858214,12,11,245.5,0.99592668,21,0.917856329,1.65243e-05,0.00618953,0.004065015,0.078673563,0.325912449,0.011209905,0.105876841 +2013,"6th_bridge","greg",32,16,459,0.965141612,28,0.370296101,7.32969e-05,0.008458978,0.008561361,0.091972704,21,14,372.5,0.981208054,21,0.671402862,4.95002e-05,0.010220905,0.007035635,0.101098491,0.248617862,0.005301094,0.072808614 +2013,"6th_bridge","rosp",3,3,33.5,1,22,1,0,0,0,0,3,3,63,1,21,1,0,0,0,0,1,0,0 +2013,"6th_bridge","smhe",28,14,398,0.964824121,22,0.454840733,8.52727e-05,0.009172291,0.009234322,0.095772076,16,4,151.5,0.920792079,14,0.314964819,0.000481413,0.01104014,0.021941119,0.105072069,0.143258829,0.003295166,0.057403534 +2013,"6th_bridge","trhe",,,,,22,,,,,,10,3,124,0.943548387,14,0.443300756,0.000429555,0.018584185,0.02072571,0.136323822,,, +2013,"joule","greg",42,11,298.5,0.896147404,28,0.046411395,0.000311783,0.00065563,0.017657379,0.025605275,13,11,240.5,0.991683992,21,0.839149531,3.42904e-05,0.01082789,0.005855804,0.104057148,0.038946101,0.000485,0.022022711 +2013,"joule","sneg",2,2,43,1,28,1,0,0,0,0,2,0,16.5,0.878787879,21,0.066307892,0.006455742,0.016208659,0.080347632,0.12731323,0.066307892,0.016208659,0.12731323 +2013,"paurotis_pond","wost",89,46,719.5,0.940236275,28,0.178088581,7.80987e-05,0.002196636,0.008837347,0.046868282,18,17,119.5,0.991631799,55,0.629904163,6.94408e-05,0.08475944,0.008333114,0.291134746,0.112178739,0.003745957,0.061204222 +2013,"tamiami_west","bcnh",16,14,67.5,0.97037037,28,0.430775267,0.000425951,0.065811609,0.020638588,0.256537733,25,4,296.5,0.929173693,21,0.213813447,0.000221956,0.005183005,0.014898186,0.071993089,0.092105545,0.004200014,0.064807518 +2013,"tamiami_west","greg",10,9,107.5,0.990697674,28,0.769754114,8.57283e-05,0.040575327,0.009258958,0.201433182,17,15,320,0.99375,21,0.876637359,1.94092e-05,0.006660877,0.004405585,0.081614198,0.674795214,0.035398837,0.188145786 +2013,"tamiami_west","smhe",15,7,168,0.952380952,22,0.341849871,0.000269949,0.016833597,0.016430132,0.129744352,10,8,113.5,0.982378855,14,0.779661487,0.000152517,0.018829011,0.012349763,0.137218842,0.266527179,0.012750016,0.112915969 +2013,"tamiami_west","trhe",,,,,22,,,,,,7,5,78.5,0.974522293,14,0.696762707,0.000316288,0.031690185,0.017784482,0.178017372,,, +2013,"tamiami_west","whib",50,36,353,0.960339943,21,0.427488907,0.000107896,0.009428485,0.010387279,0.097100387,58,21,668.5,0.944652206,14,0.45061662,7.82115e-05,0.003488168,0.008843729,0.059060718,0.192633606,0.002584844,0.05084136 +2013,"tamiami_west","wost",30,23,339.5,0.979381443,28,0.558023371,5.94799e-05,0.015138661,0.007712322,0.123039264,29,26,1535,0.998045603,55,0.897989382,1.27073e-06,0.003111873,0.001127269,0.055784166,0.501099062,0.013176594,0.114789347 +2015,"alley_north","greg",42,17,618,0.959546926,28,0.31466911,6.28101e-05,0.005295673,0.00792528,0.072771372,16,16,320,1,21,1,0,0,0,0,0.31466911,0.005295673,0.072771372 +2015,"alley_north","rosp",6,6,67,1,22,1,0,0,0,0,11,11,231,1,21,1,0,0,0,0,1,0,0 +2015,"alley_north","whib",78,41,894,0.958612975,22,0.394595991,4.43782e-05,0.003639436,0.006661699,0.060327738,64,64,147,1,14,1,0,0,0,0,0.394595991,0.003639436,0.060327738 +2015,"cuthbert_lake","greg",53,51,1084,0.998154982,28,0.94960589,1.69891e-06,0.001205525,0.001303421,0.034720675,51,48,1022.5,0.997066015,21,0.940160893,2.861e-06,0.001121795,0.00169145,0.033493215,0.892782322,0.0020785,0.045590566 +2015,"jerrod","greg",20,16,326,0.987730061,28,0.707736519,3.7176e-05,0.014963934,0.006097215,0.122327161,15,15,315,1,21,1,0,0,0,0,0.707736519,0.014963934,0.122327161 +2015,"joule","greg",20,9,243.5,0.954825462,28,0.274076479,0.000177141,0.011442749,0.013309427,0.106970784,9,9,189,1,21,1,0,0,0,0,0.274076479,0.011442749,0.106970784 +2015,"paurotis_pond","greg",71,58,1062,0.987758945,28,0.708316242,1.13853e-05,0.004590012,0.003374214,0.067749627,58,56,1172.5,0.998294243,21,0.964783571,1.45232e-06,0.000598197,0.001205123,0.024458063,0.683371874,0.004575285,0.067640853 +2015,"paurotis_pond","rosp",2,2,27,1,22,1,0,0,0,0,2,2,42,1,21,1,0,0,0,0,1,0,0 +2015,"rookery_branch","greg",28,23,427.5,0.988304094,28,0.719344021,2.70389e-05,0.011230433,0.00519989,0.105973739,23,22,315,0.996825397,21,0.935407785,1.00461e-05,0.003901219,0.00316956,0.06245974,0.672879997,0.011889012,0.109036747 +2015,"rookery_branch","smhe",8,8,163,1,22,1,0,0,0,0,8,8,112,1,14,1,0,0,0,0,1,0,0 +2015,"tamiami_west","bcnh",24,14,313.5,0.968102073,28,0.40345231,9.85022e-05,0.01341233,0.009924829,0.115811615,15,13,276,0.992753623,21,0.858363211,2.60647e-05,0.008593125,0.005105363,0.09269911,0.346308621,0.011396025,0.106752167 +2015,"tamiami_west","greg",22,21,478,0.99790795,28,0.943047369,4.36752e-06,0.003057995,0.002089861,0.055299143,22,22,456,1,21,1,0,0,0,0,0.943047369,0.003057995,0.055299143 +2015,"tamiami_west","smhe",35,31,422,0.990521327,22,0.810969046,2.22484e-05,0.007218138,0.004716822,0.084959625,33,31,450.5,0.995560488,14,0.939608915,9.81088e-06,0.001712865,0.003132233,0.041386771,0.761993746,0.007511506,0.086668942 +2015,"tamiami_west","whib",5,4,58.5,0.982905983,21,0.696228105,0.00028721,0.063550355,0.016947285,0.252091957,4,2,43.5,0.954022989,14,0.517396822,0.001008348,0.058129452,0.031754491,0.241100501,0.360226209,0.048883841,0.221096903 +2015,"tamiami_west","wost",6,2,81.5,0.950920245,28,0.244362537,0.000572649,0.029647358,0.023930096,0.172184081,2,1,71.5,0.986013986,50,0.494485726,0.000192873,0.121269718,0.013887865,0.34823802,0.120833787,0.018085969,0.134484083 +2016,"6th_bridge","greg",27,18,492,0.981707317,28,0.596343903,3.65001e-05,0.010559416,0.006041533,0.10275902,18,17,374,0.997326203,21,0.945326477,7.13007e-06,0.002825022,0.00267022,0.053150934,0.563739681,0.010470821,0.102327031 +2016,"6th_bridge","rosp",3,2,42,0.976190476,22,0.58851896,0.000553396,0.097349367,0.02352437,0.312008601,2,2,42,1,21,1,0,0,0,0,0.58851896,0.097349367,0.312008601 +2016,"6th_bridge","smhe",4,2,47,0.957446809,22,0.384169098,0.00086686,0.067547587,0.029442485,0.259899187,1,1,26,1,14,1,0,0,0,0,0.384169098,0.067547587,0.259899187 +2016,"6th_bridge","whib",25,14,182,0.93956044,21,0.270034448,0.000312014,0.011365855,0.017663929,0.106610765,15,15,208.5,1,14,1,0,0,0,0,0.270034448,0.011365855,0.106610765 +2016,"alley_north","greg",24,11,299,0.956521739,28,0.288041892,0.00013909,0.00988857,0.011793641,0.09944129,11,11,222,1,21,1,0,0,0,0,0.288041892,0.00988857,0.09944129 +2016,"alley_north","smhe",5,1,30.5,0.868852459,22,0.045374595,0.003735996,0.00493157,0.061122791,0.070225135,1,1,14,1,14,1,0,0,0,0,0.045374595,0.00493157,0.070225135 +2016,"alley_north","whib",7,5,73.5,0.972789116,21,0.560264051,0.000360142,0.052681814,0.018977413,0.229525192,6,6,72,1,14,1,0,0,0,0,0.560264051,0.052681814,0.229525192 +2016,"cuthbert_lake","greg",44,44,480,1,28,1,0,0,0,0,46,44,943.5,0.997880233,21,0.956415965,2.24194e-06,0.000908239,0.001497312,0.030136998,0.956415965,0.000908239,0.030136998 +2016,"paurotis_pond","greg",12,12,103,1,28,1,0,0,0,0,19,19,376,1,21,1,0,0,0,0,1,0,0 +2016,"paurotis_pond","smhe",22,21,234.5,0.995735608,22,0.910267072,1.81075e-05,0.007324084,0.00425529,0.08558086,21,18,268.5,0.988826816,14,0.854443392,4.11484e-05,0.006021922,0.006414702,0.077601042,0.777771684,0.010380905,0.101886728 +2016,"paurotis_pond","whib",12,10,188.5,0.98938992,21,0.799313437,5.56897e-05,0.016029229,0.007462552,0.12660659,9,9,102.5,1,14,1,0,0,0,0,0.799313437,0.016029229,0.12660659 +2016,"tamiami_west","greg",36,16,504,0.96031746,28,0.321821538,7.56108e-05,0.006657334,0.008695446,0.08159249,11,10,194.5,0.994858612,21,0.897405315,2.6298e-05,0.009436594,0.005128154,0.097142135,0.288804359,0.006401555,0.08000972 +2016,"tamiami_west","smhe",19,10,240.5,0.962577963,22,0.432105826,0.000149778,0.014608379,0.012238386,0.120865129,10,8,135.5,0.985239852,14,0.812059005,0.000107323,0.014290278,0.010359688,0.119541951,0.350895428,0.012510321,0.111849545 +2018,"6th_bridge","bcnh",2,2,29.5,1,28,1,0,0,0,0,2,2,42,1,21,1,0,0,0,0,1,0,0 +2018,"6th_bridge","glib",57,40,735.5,0.976886472,21,0.611964182,3.06992e-05,0.005312881,0.00554069,0.072889511,42,22,436.5,0.954180985,14,0.518597726,0.00010016,0.005798932,0.010007973,0.07615072,0.317363233,0.003631375,0.060260894 +2018,"6th_bridge","greg",30,27,297.5,0.989915966,28,0.752927594,3.35541e-05,0.015218487,0.00579259,0.123363232,28,21,524,0.986641221,21,0.753953182,2.51533e-05,0.006477442,0.005015305,0.080482559,0.567672155,0.012421518,0.111451862 +2018,"6th_bridge","rosp",21,17,206,0.980582524,22,0.649607873,9.24293e-05,0.01963312,0.009614016,0.140118236,15,12,275,0.989090909,21,0.794255835,3.92367e-05,0.011157813,0.006263918,0.105630549,0.515954844,0.017312956,0.131578705 +2018,"6th_bridge","smhe",74,58,1124.5,0.985771454,22,0.729587238,1.24732e-05,0.003306927,0.003531739,0.057505889,59,49,756.5,0.986781229,14,0.830027184,1.72426e-05,0.002391123,0.004152422,0.048899113,0.605577241,0.003558988,0.059657252 +2018,"6th_bridge","whib",2,0,24.5,0.918367347,21,0.167241906,0.00305995,0.044751768,0.05531681,0.211546137,,,,,14,,,,,,0,0,0 +2018,"henry","bcnh",50,47,676.5,0.99556541,28,0.882987134,6.52613e-06,0.004024766,0.002554628,0.063441046,48,45,980,0.996938776,21,0.937644582,3.11414e-06,0.001214831,0.001764692,0.034854419,0.827928102,0.004490535,0.067011456 +2018,"henry","greg",59,58,957.5,0.998955614,28,0.971165773,1.0896e-06,0.000807382,0.001043841,0.028414467,58,56,1199.5,0.998332639,21,0.965563111,1.38773e-06,0.000572472,0.001178019,0.023926384,0.937721846,0.001293128,0.035960095 +2018,"henry","trhe",1,1,21,1,22,1,0,0,0,0,1,,6,,14,,,,,,0,, +2018,"jetport_south","wost",79,64,1251.5,0.988014383,28,0.713463022,9.46222e-06,0.003868346,0.003076071,0.062196027,72,58,3113,0.99550273,50,0.798220317,1.43818e-06,0.002311601,0.00119924,0.048079114,0.56950068,0.003650354,0.060418157 +2018,"joule","greg",46,44,630,0.996825397,28,0.914817871,5.02306e-06,0.003316768,0.002241217,0.057591392,45,45,945,1,21,1,0,0,0,0,0.914817871,0.003316768,0.057591392 +2018,"tamiami_west","bcnh",43,39,660.5,0.993943982,28,0.843593913,9.11331e-06,0.005146777,0.003018826,0.071741043,39,35,788.5,0.994927077,21,0.898703132,6.401e-06,0.002303225,0.00253002,0.047991931,0.758140492,0.00580783,0.07620912 +2018,"tamiami_west","greg",25,24,344.5,0.997097242,28,0.921829122,8.40154e-06,0.005629903,0.002898542,0.075032682,25,24,510,0.998039216,21,0.959620973,3.83714e-06,0.001564408,0.001958861,0.039552601,0.884606559,0.006522616,0.080762712 +2018,"tamiami_west","rosp",1,1,1,1,22,1,0,0,0,0,1,0,2.5,0.6,21,2.1937e-05,0.096,5.65926e-08,0.309838668,0.000237892,2.1937e-05,5.65926e-08,0.000237892 +2018,"tamiami_west","smhe",11,8,94.5,0.968253968,22,0.491772549,0.000325272,0.040610869,0.018035304,0.201521385,8,4,85.5,0.953216374,14,0.511306032,0.000521578,0.029414006,0.022838082,0.171505118,0.25144627,0.018925075,0.137568436 +2018,"tamiami_west","whib",56,45,716,0.984636872,21,0.722433686,2.11272e-05,0.005015626,0.004596438,0.070821082,48,32,611,0.973813421,14,0.689700557,4.17362e-05,0.004103354,0.006460359,0.064057431,0.498262916,0.004548032,0.067439095 +2018,"vacation","bcnh",2,0,15,0.866666667,28,0.018191085,0.007703704,0.002660901,0.087770745,0.051583919,,,,,,,,,,,,, +2018,"vacation","greg",46,33,790,0.983544304,28,0.62839102,2.04872e-05,0.006556488,0.004526281,0.080972141,33,29,645,0.99379845,21,0.877535334,9.55518e-06,0.003285565,0.003091145,0.057319848,0.551435324,0.006367873,0.079798955 +2019,"alley_north","whib",150,65,1294,0.93431221,21,0.240067028,4.74288e-05,0.001380899,0.006886859,0.037160452,62,49,767,0.983050847,14,0.787161323,2.17234e-05,0.002729985,0.004660841,0.052249255,0.188971479,0.001012972,0.031827218 +2020,"6th_bridge","bcnh",11,6,128.5,0.961089494,28,0.329144982,0.000291023,0.026760268,0.017059401,0.163585659,6,4,105,0.980952381,21,0.667738532,0.000177951,0.036362601,0.013339811,0.190689803,0.219782787,0.016844195,0.129785186 +2020,"6th_bridge","greg",4,4,56,1,28,1,0,0,0,0,4,4,84,1,21,1,0,0,0,0,1,0,0 +2020,"6th_bridge","rosp",16,16,119,1,22,1,0,0,0,0,16,16,420,1,21,1,0,0,0,0,1,0,0 +2020,"6th_bridge","smhe",101,92,1716,0.994755245,22,0.890752864,3.04035e-06,0.001179917,0.001743661,0.034349922,92,82,1226.5,0.991846718,14,0.891710438,6.5934e-06,0.001044538,0.002567762,0.032319308,0.794293626,0.001768219,0.0420502 +2020,"6th_bridge","whib",75,62,1257,0.989657916,21,0.803872471,8.1425e-06,0.002369194,0.002853507,0.048674365,58,53,806,0.993796526,14,0.916567893,7.64887e-06,0.001275228,0.002765659,0.035710333,0.736803696,0.002817439,0.053079557 +2020,"henry","greg",22,20,271.5,0.992633517,28,0.812999347,2.69327e-05,0.014164389,0.005189669,0.119014239,20,19,420,0.997619048,21,0.951172715,5.65544e-06,0.002267222,0.002378116,0.047615354,0.773302796,0.014345616,0.119773186 +2020,"henry","rosp",15,13,200,0.99,22,0.80163059,4.95e-05,0.015708284,0.007035624,0.125332692,13,12,269.5,0.996289425,21,0.924902452,1.37173e-05,0.005213476,0.003703685,0.072204405,0.741430098,0.016869701,0.129883411 +2020,"jerrod","greg",44,36,913,0.991237678,28,0.781589015,9.51319e-06,0.004637069,0.003084346,0.06809603,35,34,735,0.998639456,21,0.971813969,1.84856e-06,0.000772007,0.001359618,0.027785021,0.759559123,0.004854537,0.069674507 +2020,"joule","greg",48,39,817,0.988984088,28,0.733331801,1.33348e-05,0.005748125,0.00365169,0.075816389,35,32,695.5,0.995686556,21,0.913220205,6.17518e-06,0.002290836,0.002484991,0.047862676,0.669693418,0.006038894,0.077710319 +2020,"start_mel","greg",36,32,636,0.993710692,28,0.838067412,9.82666e-06,0.005479738,0.00313475,0.07402525,32,32,672,1,21,1,0,0,0,0,0.838067412,0.005479738,0.07402525 +2020,"tamiami_west","bcnh",17,13,208,0.980769231,28,0.58059235,9.06776e-05,0.024912951,0.00952248,0.157838369,13,10,239.5,0.987473904,21,0.767428942,5.16459e-05,0.013756233,0.007186508,0.117286969,0.445563373,0.019652174,0.140186213 +2020,"tamiami_west","greg",24,20,384.5,0.989596879,28,0.746161566,2.67748e-05,0.011934131,0.005174434,0.109243448,19,18,398.5,0.99749059,21,0.948604002,6.28134e-06,0.002505204,0.00250626,0.050052009,0.707811847,0.01216361,0.110288756 +2020,"tamiami_west","rosp",2,1,21.5,0.953488372,22,0.35070247,0.002062711,0.135061227,0.045417083,0.367506771,1,1,21,1,21,1,0,0,0,0,0.35070247,0.135061227,0.367506771 +2020,"tamiami_west","smhe",15,13,211.5,0.990543735,22,0.81137276,4.42877e-05,0.014382083,0.006654899,0.119925324,12,11,156.5,0.993610224,14,0.914165276,4.05684e-05,0.006730708,0.006369329,0.082040895,0.741728803,0.016546881,0.12863468 +2020,"tamiami_west","whib",108,59,1382.5,0.964556962,21,0.468689777,2.47283e-05,0.002574824,0.004972752,0.050742719,59,44,711.5,0.978917779,14,0.742073681,2.9006e-05,0.003266972,0.005385721,0.057157428,0.347802348,0.002143955,0.046302857 +2020,"tamiami_west","wost",41,37,822.5,0.995136778,28,0.872404224,5.88398e-06,0.003545333,0.002425691,0.059542696,37,18,1534.5,0.987618117,50,0.536355078,7.96909e-06,0.005875916,0.002822958,0.076654523,0.467918436,0.005512838,0.074248486 +2020,"vacation","greg",27,26,578,0.998269896,28,0.952671762,2.98808e-06,0.002133531,0.001728607,0.046190163,27,27,1199.5,1,21,1,0,0,0,0,0.952671762,0.002133531,0.046190163 +2021,"6th_bridge","greg",32,22,639.5,0.984362783,28,0.643198759,2.40699e-05,0.008056941,0.004906107,0.089760464,21,18,424,0.992924528,21,0.861471718,1.65694e-05,0.005500396,0.004070548,0.074164653,0.55409754,0.008299182,0.091099844 +2021,"6th_bridge","rosp",10,10,175,1,22,1,0,0,0,0,9,9,189,1,21,1,0,0,0,0,1,0,0 +2021,"alley_north","greg",21,10,239.5,0.954070981,28,0.268076797,0.000182963,0.011324929,0.013526367,0.106418649,10,9,199.5,0.994987469,21,0.899849402,2.49995e-05,0.009017265,0.004999953,0.094959282,0.241228746,0.00992027,0.099600551 +2021,"alley_north","rosp",3,3,20.5,1,22,1,0,0,0,0,4,4,84,1,21,1,0,0,0,0,1,0,0 +2021,"alley_north","whib",160,102,2260,0.974336283,21,0.579277133,1.10642e-05,0.0017247,0.003326289,0.041529509,102,97,1408.5,0.996450124,14,0.951432361,2.51138e-06,0.000448758,0.001584732,0.02118392,0.551143011,0.001712599,0.041383565 +2021,"cypress_city","bcnh",1,0,24.5,0.959183673,28,0.311350653,0.001597974,0.132002355,0.039974662,0.363321283,0,0,0,,21,,,,,,,, +2021,"cypress_city","greg",43,31,628.5,0.980906921,28,0.582878947,2.97988e-05,0.008249273,0.005458826,0.090825509,31,29,645,0.996899225,21,0.936863726,4.7925e-06,0.001866598,0.002189177,0.043204141,0.546078142,0.00789007,0.088826068 +2021,"cypress_city","rosp",9,8,170,0.994117647,22,0.878276437,3.43985e-05,0.012994878,0.005865026,0.113995078,8,6,166,0.987951807,21,0.775266397,7.1705e-05,0.019472359,0.008467882,0.139543395,0.680898208,0.02308384,0.15193367 +2021,"henry","bcnh",3,3,45,1,28,1,0,0,0,0,3,3,63,1,21,1,0,0,0,0,1,0,0 +2021,"henry","greg",84,74,1661,0.99397953,28,0.844439117,3.60278e-06,0.002038619,0.0018981,0.045151067,74,73,1551,0.999355255,21,0.986547289,4.15429e-07,0.000178538,0.000644537,0.013361815,0.833079121,0.002111813,0.045954471 +2021,"henry","rosp",3,3,58,1,22,1,0,0,0,0,3,3,63,1,21,1,0,0,0,0,1,0,0 +2021,"hidden","whib",112,101,1647.5,0.993323217,21,0.868764996,4.02562e-06,0.001357983,0.002006394,0.036850827,101,79,2567,0.991429684,14,0.886475734,3.31004e-06,0.000518679,0.001819351,0.022774529,0.770139087,0.001459335,0.038201249 +2021,"jerrod","greg",34,17,371.5,0.954239569,28,0.269406333,0.000117541,0.007345232,0.010841624,0.08570433,16,15,334,0.997005988,21,0.938972992,8.93727e-06,0.00349586,0.002989527,0.059125802,0.25296527,0.00675548,0.082191724 +2021,"jetport_south","rosp",1,1,15,1,22,1,0,0,0,0,1,1,21,1,21,1,0,0,0,0,1,0,0 +2021,"jetport_south","wost",91,81,1785.5,0.994399328,28,0.854482198,3.11918e-06,0.001805683,0.001766121,0.042493332,81,67,3677,0.996192548,50,0.826350777,1.03154e-06,0.001774461,0.001015645,0.042124348,0.706102028,0.002531829,0.050317288 +2021,"joule","greg",41,25,760,0.978947368,28,0.551139593,2.71177e-05,0.006738648,0.005207462,0.082089266,25,25,525,1,21,1,0,0,0,0,0.551139593,0.006738648,0.082089266 +2022,"6th_bridge","bcnh",16,5,148,0.925675676,28,0.115039729,0.000464866,0.005628882,0.021560759,0.075025877,5,3,90,0.977777778,21,0.623797194,0.000241427,0.043334146,0.015537909,0.208168551,0.07176146,0.00300774,0.054842868 +2022,"6th_bridge","greg",26,20,428.5,0.985997666,28,0.673790428,3.222e-05,0.011796127,0.005676266,0.108609977,15,12,299,0.989966555,21,0.809153613,3.322e-05,0.00978721,0.005763678,0.098930328,0.545199959,0.012282054,0.110824431 +2022,"6th_bridge","rosp",10,7,160,0.98125,22,0.659405757,0.00011499,0.02513347,0.01072335,0.15853539,3,3,63,1,21,1,0,0,0,0,0.659405757,0.02513347,0.15853539 +2022,"6th_bridge","smhe",59,51,896,0.991071429,22,0.820935496,9.87595e-06,0.003279687,0.003142603,0.057268554,,,,,14,,,,,,,, +2022,"6th_bridge","whib",63,43,990.5,0.979808178,21,0.651571766,1.99739e-05,0.003895319,0.004469213,0.062412494,,,,,14,,,,,,,, +2022,"cypress_city","greg",38,31,765,0.990849673,28,0.773067781,1.18518e-05,0.00565612,0.003442639,0.07520718,31,30,650,0.998461538,21,0.968184539,2.36322e-06,0.000979934,0.001537278,0.031303897,0.748472273,0.005893125,0.076766694 +2022,"cypress_city","rosp",19,15,238,0.983193277,22,0.688740795,6.94297e-05,0.016490133,0.008332446,0.128413912,15,14,292.5,0.996581197,21,0.930607316,1.16483e-05,0.004479268,0.003412954,0.066927337,0.640947223,0.016479616,0.128372957 +2022,"hidden","greg",60,17,575,0.925217391,28,0.11345563,0.000120331,0.001418591,0.010969536,0.037664185,17,13,316,0.987341772,21,0.765275376,3.95506e-05,0.010478338,0.006288929,0.102363751,0.0868248,0.000980536,0.031313514 +2022,"hidden","rosp",21,19,126,0.984126984,22,0.703274841,0.000123977,0.030643142,0.011134482,0.175051826,19,17,379,0.994722955,21,0.894839085,1.38501e-05,0.004942856,0.003721576,0.070305444,0.629317815,0.027133276,0.164721813 +2022,"jerrod","greg",30,11,384.5,0.950585176,28,0.241963045,0.000122166,0.006205603,0.011052893,0.078775653,11,9,202.5,0.990123457,21,0.811851014,4.82913e-05,0.014317983,0.006949197,0.119657776,0.196437943,0.00501724,0.070832479 +2022,"joule","greg",28,11,392,0.956632653,28,0.288978558,0.000105833,0.007571447,0.010287527,0.087014062,11,9,198,0.98989899,21,0.807994683,5.04999e-05,0.014837623,0.007106328,0.121809782,0.233493139,0.006294472,0.079337706 +2022,"start_mel","greg",21,2,206.5,0.907990315,28,0.067030746,0.000404571,0.001728608,0.020113949,0.041576528,2,2,42,1,21,1,0,0,0,0,0.067030746,0.001728608,0.041576528 +2022,"vacation","greg",64,43,1112,0.981115108,28,0.58635276,1.66621e-05,0.004665775,0.004081923,0.068306476,42,38,861,0.995354239,21,0.906840859,5.37071e-06,0.001965969,0.002317478,0.044339254,0.531728641,0.00452204,0.067246115 +2023,"6th_bridge","bcnh",30,24,407,0.985,24,0.7,0,0.01,0.006,0.102,22,13,365,0.975,21,0.592,0,0.011,0.008,0.103,0.414292523,0.008997512,0.094855214 +2023,"6th_bridge","glib",57,19,702.5,0.946,21,0.311,0,0.003,0.009,0.059,19,0,139,0.863,14,0.128,0.001,0.004,0.029,0.06,0.039733062,0.000421794,0.020537616 +2023,"6th_bridge","greg",30,22,370,0.978,28,0.542,0,0.014,0.008,0.117,22,21,462,0.998,21,0.956,0,0.002,0.002,0.043,0.517685494,0.013164874,0.114738284 +2023,"6th_bridge","rosp",4,4,38,1,22,1,0,0,0,0,4,3,69,0.985,21,0.734,0,0.051,0.014,0.227,0.734306256,0.051427889,0.22677718 +2023,"6th_bridge","smhe",92,82,1502,0.993,22,0.863,0,0.002,0.002,0.04,75,41,892.5,0.962,14,0.581,0,0.003,0.006,0.054,0.501215665,0.00273221,0.052270547 +2023,"6th_bridge","smwh",6,6,100.5,1,22,1,0,0,0,0,62,35,752,0.964,14,0.599,0,0.003,0.007,0.059,0.599352503,0.003486813,0.059049241 +2023,"6th_bridge","trhe",1,1,20,1,22,1,0,0,0,0,8,5,101.5,0.97,14,0.657,0,0.025,0.017,0.159,0.65702612,0.025388692,0.159338294 +2023,"6th_bridge","whib",136,82,1800.5,0.97,21,0.528,0,0.002,0.004,0.046,72,54,902.5,0.98,14,0.754,0,0.003,0.005,0.05,0.397917943,0.0019042,0.043637145 +2023,"hidden","bcnh",16,13,247,0.987854251,24,0.745811291,4.85758e-05,0.015948317,0.006969636,0.126286647,13,5,90,0.911111111,21,0.141578792,0.000899863,0.009582275,0.029997714,0.097889094,0.105591061,0.00580249,0.076174075 +2023,"hidden","greg",87,54,926,0.964362851,28,0.362020548,3.71135e-05,0.004100473,0.006092089,0.064034939,53,41,1044,0.988505747,21,0.784446214,1.08833e-05,0.003022498,0.00329898,0.054977247,0.283985648,0.002931769,0.054145814 +2023,"hidden","rosp",6,4,85.5,0.976608187,22,0.59408411,0.000267189,0.047854044,0.016345908,0.218755671,4,3,74.5,0.986577181,21,0.752926171,0.000177754,0.0456562,0.013332428,0.213673115,0.447301474,0.0454269,0.213135873 +2023,"hidden","smhe",9,8,141.5,0.992932862,22,0.85553436,4.95915e-05,0.017819172,0.007042121,0.13348847,8,2,89,0.93258427,14,0.376385774,0.000706414,0.022553068,0.026578451,0.150176788,0.322010962,0.019433723,0.139404889 +2023,"hidden","smwh",6,6,100.5,1,22,1,0,0,0,0,6,2,72,0.944444444,14,0.449231114,0.000728738,0.032315755,0.026995148,0.179765834,0.449231114,0.032315755,0.179765834 +2023,"hidden","trhe",1,1,20,1,22,1,0,0,0,0,1,0,12,0.916666667,14,0.295774104,0.006365741,0.129897991,0.079785592,0.360413639,0.295774104,0.129897991,0.360413639 +2023,"hidden","whib",50,13,650,0.943076923,21,0.292071698,8.2589e-05,0.003493377,0.009087848,0.059104795,13,1,142.5,0.915789474,14,0.291836133,0.000541187,0.010771844,0.023263421,0.103787496,0.085237075,0.001254057,0.035412664 +2023,"jerrod","greg",12,4,112.5,0.928888889,28,0.126760955,0.00058715,0.008572515,0.024231168,0.092587876,4,3,70.5,0.985815603,21,0.740814441,0.000198343,0.049395063,0.014083439,0.222250002,0.093906346,0.005921784,0.076953132 +2023,"little_a","greg",74,54,1215.5,0.983545866,28,0.628418966,1.33142e-05,0.00426128,0.003648861,0.065278481,54,38,1004,0.984063745,21,0.713654268,1.56198e-05,0.003622794,0.003952191,0.060189651,0.448473877,0.003616397,0.060136487 +2023,"vacation","greg",47,15,753,0.95750332,28,0.296434044,5.40381e-05,0.004060616,0.007351063,0.063722966,15,9,276.5,0.978300181,21,0.63083359,7.67773e-05,0.014078522,0.008762268,0.118652946,0.187000552,0.002910218,0.053946435 diff --git a/Reports/sfwmd_report_tables.R b/Reports/sfwmd_report_tables.R new file mode 100644 index 0000000..5bf77dc --- /dev/null +++ b/Reports/sfwmd_report_tables.R @@ -0,0 +1,388 @@ +# SFWMD Report Table +report_year <- 2023 + +library(wader) +library(xlsx) + +wb <- createWorkbook() +#-------------------- +# DEFINE CELL STYLES +#-------------------- +# title and subtitle styles +title_style <- CellStyle(wb) + + Font(wb, heightInPoints = 12, + name="Calibri", + underline = 1) + +subtitle_style <- CellStyle(wb) + + Font(wb, + name="Calibri", + heightInPoints = 10) + +nesting_style <- CellStyle(wb) + + Font(wb, heightInPoints = 11, + name="Calibri", + color = "red", + isItalic = FALSE, + isBold = FALSE) + +# data table styles +cs <- CellStyle(wb) + + Font(wb, heightInPoints=11, name="Calibri") + + Alignment(horizontal = "ALIGN_CENTER") + +rowname_style <- CellStyle(wb) + + Font(wb, name = "Calibri") + + Alignment(horizontal = "ALIGN_CENTER") + + Border(position = c("RIGHT"), + pen = c("BORDER_THIN")) + +colname_style <- CellStyle(wb) + + Font(wb, name = "Calibri") + + Alignment(horizontal = "ALIGN_CENTER") + + Border(position = c("TOP", "BOTTOM"), + pen = c("BORDER_THIN")) + +xlsx.addTitle <- function(sheet, rowIndex, title, titleStyle) { + rows <- createRow(sheet, rowIndex = rowIndex) + sheetTitle <- createCell(rows, colIndex = 1) + setCellValue(sheetTitle[[1,1]], title) + setCellStyle(sheetTitle[[1,1]], titleStyle) +} + +#### Get data ### +# counts +main <- max_counts() %>% + dplyr::ungroup() %>% + dplyr::filter(year==report_year) %>% + dplyr::select(-year) %>% + tidyr::pivot_wider(id_cols = colony, names_from = species, values_from = count) %>% + dplyr::mutate(total = rowSums(dplyr::across(where(is.numeric)), na.rm=TRUE)) %>% + dplyr::full_join(load_datafile("SiteandMethods/colonies.csv", path = get_default_data_path()),by = dplyr::join_by(colony)) %>% + dplyr::select(-c(group_id,subregion,colony,aka)) %>% + dplyr::rename(wca=region, colony=display_name) +under40 <- load_datafile("Counts/maxcounts_under40.csv", path = get_default_data_path()) %>% + dplyr::filter(year==report_year) %>% + dplyr::select(-c(year,notes,colony_old,group_id)) + +# by WCA +lox <- main %>% + dplyr::select(colony,wca,latitude, longitude,everything()) %>% + dplyr::relocate(total, .after = last_col()) %>% + dplyr::filter(total >= 40, wca=="1") %>% + dplyr::arrange(wca, desc(total)) %>% + dplyr::mutate(latitude = as.character(latitude), longitude = as.character(longitude)) %>% + dplyr::bind_rows(dplyr::summarise_all(., ~if(is.numeric(.)) sum(., na.rm=TRUE) else "Total")) +under40_lox <- under40 %>% + dplyr::filter(wca=="1") %>% + dplyr::select(colnames(lox), + -c(colony,latitude,longitude,wca)) %>% + dplyr::summarise_all(~ sum(., na.rm=TRUE)) %>% + dplyr::mutate(colony=c("Colonies < 40 nests**")) +totals_lox <- tail(lox[,-c(2:4)],1) %>% + dplyr::bind_rows(under40_lox) %>% + dplyr::select(-colony) %>% + dplyr::summarise_all(~ sum(., na.rm=TRUE)) %>% + dplyr::mutate(colony=c("Total nests by species")) +lox <- lox %>% + dplyr::bind_rows(under40_lox) %>% + dplyr::bind_rows(totals_lox) %>% + dplyr::add_row(colony="Total nests excluding ANHI",total=tail(lox$total,1)-tail(lox$anhi,1)) %>% + dplyr::mutate(colony = replace(colony,colony=="Total","Colonies > 40 nests"), + wca = replace(wca,wca=="Total",NA), + latitude = replace(latitude,latitude=="Total",NA), + longitude = replace(longitude,longitude=="Total",NA), + ycnh = NA, + smda = NA) %>% + dplyr::rename_with(toupper, .cols = -colony) %>% + dplyr::select(Colony=colony,WCA,Latitude=LATITUDE,Longitude=LONGITUDE,GREG,WHIB,WOST,ROSP,SNEG,GBHE,LBHE, + TRHE,GLIB,BCNH,CAEG,YCNH,SMDA,SMWH,ANHI,Total=TOTAL) %>% + `row.names<-`(., NULL) %>% + tibble::column_to_rownames(var = "Colony") + +wcas <- main %>% + dplyr::select(colony,wca,latitude, longitude,everything()) %>% + dplyr::relocate(total, .after = last_col()) %>% + dplyr::filter(total >= 40, wca %in% c("2","3")) %>% + dplyr::arrange(wca, desc(total)) %>% + dplyr::mutate(latitude = as.character(latitude), longitude = as.character(longitude)) %>% + dplyr::bind_rows(dplyr::summarise_all(., ~if(is.numeric(.)) sum(., na.rm=TRUE) else "Total")) +under40_wcas <- under40 %>% + dplyr::filter(wca %in% c("2","3")) %>% + dplyr::select(colnames(wcas), + -c(colony,latitude,longitude,wca)) %>% + dplyr::summarise_all(~ sum(., na.rm=TRUE)) %>% + dplyr::mutate(colony="Colonies < 40 nests**") +totals_wcas <- tail(wcas[,-c(2:4)],1) %>% + dplyr::bind_rows(under40_wcas) %>% + dplyr::select(-colony) %>% + dplyr::summarise_all(~ sum(., na.rm=TRUE)) %>% + dplyr::mutate(colony=c("Total nests by species")) +wcas <- wcas %>% + dplyr::bind_rows(under40_wcas) %>% + dplyr::bind_rows(totals_wcas) %>% + dplyr::add_row(colony="Total nests excluding ANHI",total=tail(wcas$total,1)-tail(wcas$anhi,1)) %>% + dplyr::mutate(colony = replace(colony,colony=="Total","Colonies > 40 nests"), + wca = replace(wca,wca=="Total",NA), + latitude = replace(latitude,latitude=="Total",NA), + longitude = replace(longitude,longitude=="Total",NA), + ycnh = NA, + smda = NA) %>% + dplyr::rename_with(toupper, .cols = -colony) %>% + dplyr::select(Colony=colony,WCA,Latitude=LATITUDE,Longitude=LONGITUDE,GREG,WHIB,WOST,ROSP,SNEG,GBHE,LBHE, + TRHE,GLIB,BCNH,CAEG,YCNH,SMDA,SMWH,ANHI,Total=TOTAL) %>% + `row.names<-`(., NULL) %>% + tibble::column_to_rownames(var = "Colony") + +# all wcas +all <- main %>% + dplyr::select(colony,wca,latitude, longitude,everything()) %>% + dplyr::relocate(total, .after = last_col()) %>% + dplyr::filter(total >= 40) %>% + dplyr::arrange(wca, desc(total)) %>% + dplyr::mutate(latitude = as.character(latitude), longitude = as.character(longitude)) %>% + dplyr::bind_rows(dplyr::summarise_all(., ~if(is.numeric(.)) sum(., na.rm=TRUE) else "Total")) +under40_all <- under40 %>% + dplyr::select(colnames(all), + -c(colony,latitude,longitude,wca)) %>% + dplyr::summarise_all(~ sum(., na.rm=TRUE)) %>% + dplyr::mutate(colony="Colonies < 40 nests**") +totals_all <- tail(all[,-c(2:4)],1) %>% + dplyr::bind_rows(under40_all) %>% + dplyr::select(-colony) %>% + dplyr::summarise_all(~ sum(., na.rm=TRUE)) %>% + dplyr::mutate(colony=c("Total nests by species")) +all <- all %>% + dplyr::bind_rows(under40_all) %>% + dplyr::bind_rows(totals_all) %>% + dplyr::add_row(colony="Total nests excluding ANHI",total=tail(all$total,1)-tail(all$anhi,1)) %>% + dplyr::mutate(colony = replace(colony,colony=="Total","Colonies > 40 nests"), + wca = replace(wca,wca=="Total",NA), + latitude = replace(latitude,latitude=="Total",NA), + longitude = replace(longitude,longitude=="Total",NA), + ycnh=NA, + smda=NA) %>% + dplyr::rename_with(toupper, .cols = -colony) %>% + dplyr::select(Colony=colony,WCA,Latitude=LATITUDE,Longitude=LONGITUDE,GREG,WHIB,WOST,ROSP,SNEG,GBHE,LBHE, + TRHE,GLIB,BCNH,CAEG,YCNH,SMDA,SMWH,ANHI,Total=TOTAL) %>% + `row.names<-`(., NULL) %>% + tibble::column_to_rownames(var = "Colony") + +# appendix +append <- merge(main, under40, all = TRUE) %>% + dplyr::select(colony,wca,latitude, longitude,everything()) %>% + dplyr::relocate(total, .after = last_col()) %>% + dplyr::filter(total > 0) %>% + dplyr::arrange(wca, desc(total)) %>% + dplyr::mutate(latitude = as.character(latitude), longitude = as.character(longitude)) %>% + dplyr::bind_rows(dplyr::summarise_all(., ~if(is.numeric(.)) sum(., na.rm=TRUE) else "Total")) %>% + dplyr::rename_with(toupper, .cols = -colony) +append <- append %>% + dplyr::add_row(colony="Total nests excluding ANHI",TOTAL=tail(append$TOTAL,1)-tail(append$ANHI,1)) %>% + dplyr::mutate(colony = replace(colony,colony=="Total","Total nests by species"), + WCA = replace(WCA,WCA=="Total",NA), + LATITUDE = replace(LATITUDE,LATITUDE=="Total",NA), + LONGITUDE = replace(LONGITUDE,LONGITUDE=="Total",NA), + SMDA=NA) %>% + dplyr::select(Colony=colony,WCA,Latitude=LATITUDE,Longitude=LONGITUDE,GREG,WHIB,WOST,ROSP,SNEG,GBHE,LBHE, + TRHE,GLIB,BCNH,CAEG,YCNH,SMDA,SMWH,ANHI,Total=TOTAL) + +# nesting +success <- read.csv("Nesting/nest_success_summary.csv") %>% + dplyr::group_by(year, species) %>% + dplyr::mutate(species = toupper(species)) %>% + dplyr::filter(year==report_year) + +colony_number <- length(unique(success$colony)) +species_list <- unique(success$species) + +success <- success %>% + dplyr::summarise(incubation_N=sum(incubation_k), incubation_sumy=sum(incubation_sumy), + incubation_e=sum(incubation_e), incubation_j=mean(incubation_j), + nestling_N=sum(nestling_k), nestling_sumy=sum(nestling_sumy), nestling_e=sum(nestling_e), + nestling_j=mean(nestling_j)) %>% + dplyr::mutate(incubation_p = 1-((incubation_N-incubation_sumy)/incubation_e), + incubation_Success = incubation_p^incubation_j, + incubation_varp=(incubation_p*(1-incubation_p))/incubation_e, + incubation_var = incubation_varp*((incubation_j*(incubation_p^(incubation_j-1)))^2), + incubation_SD = sqrt(incubation_var), + nestling_p = 1-((nestling_N-nestling_sumy)/nestling_e), + nestling_Success = nestling_p^nestling_j, + nestling_varp=(nestling_p*(1-nestling_p))/nestling_e, + nestling_var = nestling_varp*((nestling_j*(nestling_p^(nestling_j-1)))^2), + nestling_SD = sqrt(nestling_var), + overall_Success = (incubation_p^incubation_j)*(nestling_p^nestling_j), + overall_var = ((incubation_Success^2)*nestling_var)+((nestling_Success^2)*incubation_var)+(incubation_var*nestling_var), + overall_SD = sqrt(overall_var)) + + incubation <- success %>% + dplyr::ungroup() %>% + dplyr::select(species, N = incubation_N , Success = incubation_Success, SD = incubation_SD) %>% + t() + colnames(incubation) <- incubation[1,] + incubation <- incubation[-1,] + + nestling <- success %>% + dplyr::ungroup() %>% + dplyr::select(species, N = nestling_N, Success = nestling_Success, SD = nestling_SD) %>% + t() + colnames(nestling) <- nestling[1,] + nestling <- nestling[-1,] + + overall <- success %>% + dplyr::ungroup() %>% + dplyr::select(species, Success = overall_Success, SD = overall_SD) %>% + t() + colnames(overall) <- overall[1,] + overall <- overall[-1,] + + clutch_fledge <- read.csv("Nesting/nest_success.csv") %>% + dplyr::mutate(species = toupper(species)) %>% + dplyr::group_by(year, species) %>% + dplyr::filter(year==report_year, species %in% species_list) %>% + dplyr::group_by(species) %>% + dplyr::summarise(N_clutch = sum(!is.na(clutch)), N_fledge = sum(!is.na(fledged)), + "Mean Clutch Size" = mean(clutch, na.rm = TRUE), + "Mean Number of Chicks" = mean(fledged, na.rm = TRUE), + SD_clutch = sd(clutch, na.rm = TRUE), + SD_fledge = sd(fledged, na.rm = TRUE)) + + clutch <- clutch_fledge %>% + dplyr::ungroup() %>% + dplyr::select(species, N = N_clutch, "Mean Clutch Size", SD = SD_clutch) %>% + t() + colnames(clutch) <- clutch[1,] + clutch <- clutch[-1,] + + fledge <- clutch_fledge %>% + dplyr::ungroup() %>% + dplyr::select(species, N = N_fledge, "Mean Number of Chicks", SD = SD_fledge) %>% + t() + colnames(fledge) <- fledge[1,] + fledge <- fledge[-1,] + +#-------------------- +# BUILD WORKBOOK +#-------------------- + +#### Build WCAS sheet #### +WCAS <- createSheet(wb, sheetName = "WCAs") + +xlsx.addTitle(sheet = WCAS, rowIndex = 1, + title = paste("Number of nesting pairs found in Loxahatchee NWR during aerial surveys, January through June of ", report_year,".", sep=""), + titleStyle = title_style) + +addDataFrame(lox, sheet = WCAS, startRow = 3, startColumn = 1, + colnamesStyle = colname_style, colStyle = cs, + rownamesStyle = rowname_style) +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+4, + title = "*Totals by species including ANHI", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+5, + title = "** Includes count of wading bird nesting pairs from ground surveys", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+6, + title = "***present but not counted", + titleStyle = subtitle_style) + +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+8, + title = paste("Number of nesting pairs found in WCAs 2 and 3 during aerial surveys, January through June of ", report_year,".", sep=""), + titleStyle = title_style) + +addDataFrame(wcas, sheet = WCAS, startRow = dim(lox)[1]+10, startColumn = 1, + colnamesStyle = colname_style, colStyle = cs, + rownamesStyle = rowname_style) +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+10+dim(wcas)[1]+1, + title = "*Totals by species including ANHI", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+10+dim(wcas)[1]+2, + title = "** Includes count of wading bird nesting pairs from ground surveys", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = WCAS, rowIndex = dim(lox)[1]+10+dim(wcas)[1]+3, + title = "***present but not counted", + titleStyle = subtitle_style) +setColumnWidth(WCAS, 1, 22) +autoSizeColumn(WCAS, 2:26) + +#### Build Colonies sheet #### +all_colonies <- createSheet(wb, sheetName = "All Colonies") +xlsx.addTitle(sheet = all_colonies, rowIndex = 1, + title = paste("Number of nesting pairs found in Loxahatchee NWR and WCA's 2 and 3 during aerial surveys, January through June of ", report_year,".", sep=""), + titleStyle = title_style) + +addDataFrame(all, sheet = all_colonies, startRow = 3, + colnamesStyle = colname_style, colStyle = cs, + rownamesStyle = rowname_style) +xlsx.addTitle(sheet = all_colonies, rowIndex = dim(all)[1]+4, + title = "*Totals by species including ANHI", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = all_colonies, rowIndex = dim(all)[1]+5, + title = "** Includes count of wading bird nesting pairs from ground surveys", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = all_colonies, rowIndex = dim(all)[1]+6, + title = "***present but not counted", + titleStyle = subtitle_style) +setColumnWidth(all_colonies, 1, 22) +autoSizeColumn(all_colonies, 2:26) + +#### Build Appendix sheet #### +Appendix <- createSheet(wb, sheetName = "Appendix") +xlsx.addTitle(sheet = Appendix, rowIndex = 1, + title = paste("Number of nesting pairs found in Loxahatchee NWR and WCA's 2 and 3 during aerial surveys, January through June of ", report_year,".", sep=""), + titleStyle = title_style) +addDataFrame(append, sheet = Appendix, startRow = 3, + colnamesStyle = colname_style, colStyle = cs, + rownamesStyle = rowname_style, row.names = FALSE) +xlsx.addTitle(sheet = Appendix, rowIndex = dim(append)[1]+4, + title = "*Totals by species including ANHI", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = Appendix, rowIndex = dim(append)[1]+5, + title = "** Includes count of wading bird nesting pairs from ground surveys", + titleStyle = subtitle_style) +xlsx.addTitle(sheet = Appendix, rowIndex = dim(append)[1]+6, + title = "***present but not counted", + titleStyle = subtitle_style) +setColumnWidth(Appendix, 1, 22) +autoSizeColumn(Appendix, 2:26) + +#### Build nesting sheet #### +Nest_Success <- createSheet(wb, sheetName = "Nest Success") +xlsx.addTitle(sheet = Nest_Success, rowIndex = 1, + title = paste("Nest success metrics of wading birds monitored in ", colony_number, " colonies throughout WCA 3 in ", report_year,".", sep=""), + titleStyle = title_style) +xlsx.addTitle(sheet = Nest_Success, rowIndex = 3, + title = "Incubation", + titleStyle = nesting_style) +addDataFrame(incubation, sheet = Nest_Success, startRow = 4, + colnamesStyle = colname_style, + rownamesStyle = rowname_style, colStyle = cs) +xlsx.addTitle(sheet = Nest_Success, rowIndex = 9, + title = "Nestling", + titleStyle = nesting_style) +addDataFrame(nestling, sheet = Nest_Success, startRow = 10, + colnamesStyle = colname_style, + rownamesStyle = rowname_style, colStyle = cs) +xlsx.addTitle(sheet = Nest_Success, rowIndex = 15, + title = "Overall Success (P)", + titleStyle = nesting_style) +addDataFrame(overall, sheet = Nest_Success, startRow = 16, + colnamesStyle = colname_style, + rownamesStyle = rowname_style, colStyle = cs) +xlsx.addTitle(sheet = Nest_Success, rowIndex = 20, + title = "Mean Clutch Size", + titleStyle = nesting_style) +addDataFrame(clutch, sheet = Nest_Success, startRow = 21, + colnamesStyle = colname_style, + rownamesStyle = rowname_style, colStyle = cs) +xlsx.addTitle(sheet = Nest_Success, rowIndex = 26, + title = "Mean Number of Chicks Fledged", + titleStyle = nesting_style) +addDataFrame(fledge, sheet = Nest_Success, startRow = 27, + colnamesStyle = colname_style, + rownamesStyle = rowname_style, colStyle = cs) +setColumnWidth(Nest_Success, 1, 27) +autoSizeColumn(Nest_Success, 2:9) + +#### Write #### +saveWorkbook(wb, file = 'Reports/sfwmd_report_table_2023.xlsx') + diff --git a/SiteandMethods/species_list.csv b/SiteandMethods/species_list.csv index 3c66b2a..964b316 100644 --- a/SiteandMethods/species_list.csv +++ b/SiteandMethods/species_list.csv @@ -1,24 +1,24 @@ -"species","commonname","scientificname","target_species","clutch_size","egg_color","nest_size","nest_materials","nest_microhabitat","colony_habitat","brood_size","nest_success","chick_description","typical_timing_range","courtship_period","nestbuilding_period","reproductive_period","nestling_period","incubation_period","branchling_period","feeding_methods","foraging_depth","foraging_type","niche_description","prey" -anhi,Anhinga,Anhinga anhinga,no,"5-Feb","Conspicuously pointed at one end, pale bluish green, and overlaid with a chalky coating","","","","shallow, slow-moving, sheltered waters","1","","Naked, with eyes open","","","","","14 - 21 days","26 - 30 days","","swims slowly underwater, stalking fish around submerged vegetation, spear fish","","","","small- to medium-sized wetland fishes, with very small amounts of crustaceans and invertebrates" -bcnh,Black Crowned Night Heron,Nycticorax nycticorax,no,"","","","","","","","","","","","","","","","","","","","","" -caeg,Cattle Egret,Bubulcus ibis,no,"","","","","","","","","","","","","","","","","","","","","" -coga,Common Gallinule,Gallinula galeata,no,"","","","","","","","","","","","","","","","","","","","","" -dcco,Double-crested Cormorant,Phalacrocorax auritus,no,"","","","","","","","","","","","","","","","","","","","","" -gbhe,Great Blue Heron,Ardea herodias,yes,"6-Feb","Pale blue, fading slightly with age","50 cm - 122 cm","sticks, lined with grass, reeds, moss, cupped","mainly in trees, but will also nest on the ground, on bushes, in mangroves, and on structures such as duck blinds, channel markers, or artificial nest platforms","","2-Jan","","bluish eyes open, covered in pale gray down, able to vocalize","","","","","49 - 81 days","27 - 29 days","","grab smaller prey in their strong mandibles or use their dagger-like bills to impale larger fish, often shaking them to break or relax the sharp spines before gulping them down","","exploiter","","anything within striking distance, including fish, amphibians, reptiles, small mammals, insects, and other birds" -glib,Glossy Ibis,Plegadis falcinellus,no,"","","","","","","","","","","","","","","","","","","","","" -greg,Great Egret,Ardea alba,yes,"1 - 6, 4 normal","light blue green","widths of ground nests (49.2 cm ± 5.9 SD) versus tree nests (42.5 cm ± 3.7 SD) were significantly different, as were nest depths (21.2 cm ± 8.7 SD and 10.5 cm ± 2.2 SD),","sticks, vary in size","open canopy","","3-Jan","","white down, body dak gray, wram pale gray bill balck tip, change to yellow by 7 days.","Jan - May","7 - 11 days","3 - 5days","80","21","27","50+ days","walking slowly, standing and waiting, second to SNEG in number of different foraging behaviors; kleptoparasitism from other species","average 28 cm, GE will increase to 20 - 40 cm then decline","exploiter","depth more important than salinity, but will stay in open areas","98% fish, but opportunistic, size larger than other herons except for GBHE" -grhe,Green Heron,Butorides virescens,no,"","","","","","","","","","","","","","","","","","","","","" -lada,large dark unidentified bird,unknown,no,"","","","","","","","","","","","","","","","","","","","","" -lawh,large white unidentified bird,unknown,no,"","","","","","","","","","","","","","","","","","","","","" -lbhe,Little Blue Heron,Egretta caerulea,yes,"","","","","","","","","","","","","","","","","","","","","" -rosp,Roseate Spoonbill,Platalea ajaja,yes,"1 - 5, 4 typical","dull white, brown spots","","twigs and large sticks","subcanopy","","","","pinkish down, tubular to flattened bills, exceptionally cute","Nov - Feb","","20 days","80","5-6weeks","22 days","60 days","tactolocation, head sweeping, will use grabbing or probing, and will run after abundant or active prey.","<12 cm, max 20 cm based on leg length","high grader","fresh to hypersaline, will feed at night; prefers coastal habitat","" -rsha,Red-shouldered Hawk,Buteo lineatus,no,"","","","","","","","","","","","","","","","","","","","","" -smda,small dark unidentified bird,unknown,no,"","","","","","","","","","","","","","","","","","","","","" -smhe,small heron,unknown,no,"","","","","","","","","","","","","","","","","","","","","" -smwh,small white unidentified bird,unknown,no,"","","","","","","","","","","","","","","","","","","","","" -sneg,Snowy Egret,Egretta thula,yes,"5-Mar","pale greenish blue","","fine twigs , cupped","head to chest height, never canopy","","","","white down, leg yellowish to grayish to black, pink mouth, pale gray iris","mid - late March","","4.4 days","","15 days","22 days","","broadest repertoire of foraging of any n. american heron. Foot stirring, fly-stabbing, foot dragging, lots of habitats,","","high grader","brackish and marine, shallow water","worms, insects, crustaceans, fishes, anurans 75% fish, mosquitofish, needlefish" -trhe,Tricolored Heron,Egretta tricolor,no,"4-Mar","pale bluish green","","small twigs, woody vege, no green material, will nest in needlerush","subcanopy, dense, well shaded 0.15 - 3.7 m above ground","","","","brownish reddish crown, bronish gray on back, white down abdoment,","Feb - March","5 - 6 days","","","83 days","21 days","50 - 56 days","broad - striking, running, wing shading, foot raking,. Slow walking, usually solitary foraging or on edge of groups","deeper than snowy, to belly deep, 18cm max","exploiter","estuarine, coastal primary habitat","primarily fish, opportunistic, fundulus, killifishes" -unkn,Unknown,Unknown,no,"","","","","","","","","","","","","","","","","","","","","" -whib,White Ibis,Eudocimus albus,yes,"3-Feb","creamy white to greenish, brown splotches","25 cm average","various, twigs to grass","various height, very dense","","3-Jan","5 - 70%","dark downy, pink bills with dark markings,","Feb - March in Everglades","5 - 6 days","","","15 - 18 days","22 days","61 days mean, as early as 40 days","walk and probe quickly, will use visual foraging, foot raking, pecking at surface,","14 cm on average","high grader","prefer freshwater in breeding season, will feed in coastal","crustaceans, insects, small fishes, gambusia, crayfish, fiddler crabs" -wost,Wood Stork,Mycteria americana,yes,"5-Jan","creamy white","","twigs and sticks","open at top","historically cypress swamps, now many types including impoundments","","","white, sparse down, eyes closed, 60 gm at hatching","historically november, now Feb - March","","2 - 3 days","99 - 118 days","25 days","28 days","50 - 55 days","grope foraging, walking slowly with bill open in water, side to side motion, snap reflex, unitary foraging technique. Foot stirring and wing flicking","15 - 50 cm","high grader","wetlands,","fish>1 yr old, sunfishes, catfishes, killifishes, mollies," -ycnh,Yellow-crowned Night-Heron,Nyctanassa violacea,no,"","","","","","","","","","","","","","","","","","","","","" +"species","commonname","scientificname","target_species","incubation_j","nestling_j","clutch_size","egg_color","nest_size","nest_materials","nest_microhabitat","colony_habitat","brood_size","nest_success","chick_description","typical_timing_range","courtship_period","nestbuilding_period","reproductive_period","nestling_period","incubation_period","branchling_period","feeding_methods","foraging_depth","foraging_type","niche_description","prey" +anhi,Anhinga,Anhinga anhinga,no,,,"2 - 5","Conspicuously pointed at one end, pale bluish green, and overlaid with a chalky coating","","","","shallow, slow-moving, sheltered waters","1","","Naked, with eyes open","","","","","14 - 21 days","26 - 30 days","","swims slowly underwater, stalking fish around submerged vegetation, spear fish","","","","small- to medium-sized wetland fishes, with very small amounts of crustaceans and invertebrates" +bcnh,Black Crowned Night Heron,Nycticorax nycticorax,no,28,21,"","","","","","","","","","","","","","","","","","","","","" +caeg,Cattle Egret,Bubulcus ibis,no,,,"","","","","","","","","","","","","","","","","","","","","" +coga,Common Gallinule,Gallinula galeata,no,,,"","","","","","","","","","","","","","","","","","","","","" +dcco,Double-crested Cormorant,Phalacrocorax auritus,no,,,"","","","","","","","","","","","","","","","","","","","","" +gbhe,Great Blue Heron,Ardea herodias,yes,28,60,"2 - 6","Pale blue, fading slightly with age","50 cm - 122 cm","sticks, lined with grass, reeds, moss, cupped","mainly in trees, but will also nest on the ground, on bushes, in mangroves, and on structures such as duck blinds, channel markers, or artificial nest platforms","","1 - 2","","bluish eyes open, covered in pale gray down, able to vocalize","","","","","49 - 81 days","27 - 29 days","","grab smaller prey in their strong mandibles or use their dagger-like bills to impale larger fish, often shaking them to break or relax the sharp spines before gulping them down","","exploiter","","anything within striking distance, including fish, amphibians, reptiles, small mammals, insects, and other birds" +glib,Glossy Ibis,Plegadis falcinellus,no,,,"","","","","","","","","","","","","","","","","","","","","" +greg,Great Egret,Ardea alba,yes,28,21,"1 - 6, 4 normal","light blue green","widths of ground nests (49.2 cm ± 5.9 SD) versus tree nests (42.5 cm ± 3.7 SD) were significantly different, as were nest depths (21.2 cm ± 8.7 SD and 10.5 cm ± 2.2 SD),","sticks, vary in size","open canopy","","1 - 3","","white down, body dak gray, wram pale gray bill balck tip, change to yellow by 7 days.","Jan - May","7 - 11 days","3 - 5days","80","21","27","50+ days","walking slowly, standing and waiting, second to SNEG in number of different foraging behaviors; kleptoparasitism from other species","average 28 cm, GE will increase to 20 - 40 cm then decline","exploiter","depth more important than salinity, but will stay in open areas","98% fish, but opportunistic, size larger than other herons except for GBHE" +grhe,Green Heron,Butorides virescens,no,,,"","","","","","","","","","","","","","","","","","","","","" +lada,large dark unidentified bird,unknown,no,,,"","","","","","","","","","","","","","","","","","","","","" +lawh,large white unidentified bird,unknown,no,,,"","","","","","","","","","","","","","","","","","","","","" +lbhe,Little Blue Heron,Egretta caerulea,yes,,,"","","","","","","","","","","","","","","","","","","","","" +rosp,Roseate Spoonbill,Platalea ajaja,yes,22,21,"1 - 5, 4 typical","dull white, brown spots","","twigs and large sticks","subcanopy","","","","pinkish down, tubular to flattened bills, exceptionally cute","Nov - Feb","","20 days","80","5-6weeks","22 days","60 days","tactolocation, head sweeping, will use grabbing or probing, and will run after abundant or active prey.","<12 cm, max 20 cm based on leg length","high grader","fresh to hypersaline, will feed at night; prefers coastal habitat","" +rsha,Red-shouldered Hawk,Buteo lineatus,no,,,"","","","","","","","","","","","","","","","","","","","","" +smda,small dark unidentified bird,unknown,no,,,"","","","","","","","","","","","","","","","","","","","","" +smhe,small heron,unknown,no,22,14,"","","","","","","","","","","","","","","","","","","","","" +smwh,small white unidentified bird,unknown,no,,,"","","","","","","","","","","","","","","","","","","","","" +sneg,Snowy Egret,Egretta thula,yes,,,"3 - 5","pale greenish blue","","fine twigs , cupped","head to chest height, never canopy","","","","white down, leg yellowish to grayish to black, pink mouth, pale gray iris","mid - late March","","4.4 days","","15 days","22 days","","broadest repertoire of foraging of any n. american heron. Foot stirring, fly-stabbing, foot dragging, lots of habitats,","","high grader","brackish and marine, shallow water","worms, insects, crustaceans, fishes, anurans 75% fish, mosquitofish, needlefish" +trhe,Tricolored Heron,Egretta tricolor,no,,,"3 - 4","pale bluish green","","small twigs, woody vege, no green material, will nest in needlerush","subcanopy, dense, well shaded 0.15 - 3.7 m above ground","","","","brownish reddish crown, bronish gray on back, white down abdoment,","Feb - March","5 - 6 days","","","83 days","21 days","50 - 56 days","broad - striking, running, wing shading, foot raking,. Slow walking, usually solitary foraging or on edge of groups","deeper than snowy, to belly deep, 18cm max","exploiter","estuarine, coastal primary habitat","primarily fish, opportunistic, fundulus, killifishes" +unkn,Unknown,Unknown,no,,,"","","","","","","","","","","","","","","","","","","","","" +whib,White Ibis,Eudocimus albus,yes,21,14,"2 - 3","creamy white to greenish, brown splotches","25 cm average","various, twigs to grass","various height, very dense","","1 - 3","5 - 70%","dark downy, pink bills with dark markings,","Feb - March in Everglades","5 - 6 days","","","15 - 18 days","22 days","61 days mean, as early as 40 days","walk and probe quickly, will use visual foraging, foot raking, pecking at surface,","14 cm on average","high grader","prefer freshwater in breeding season, will feed in coastal","crustaceans, insects, small fishes, gambusia, crayfish, fiddler crabs" +wost,Wood Stork,Mycteria americana,yes,28,50,"1 - 5","creamy white","","twigs and sticks","open at top","historically cypress swamps, now many types including impoundments","","","white, sparse down, eyes closed, 60 gm at hatching","historically november, now Feb - March","","2 - 3 days","99 - 118 days","25 days","28 days","50 - 55 days","grope foraging, walking slowly with bill open in water, side to side motion, snap reflex, unitary foraging technique. Foot stirring and wing flicking","15 - 50 cm","high grader","wetlands,","fish>1 yr old, sunfishes, catfishes, killifishes, mollies," +ycnh,Yellow-crowned Night-Heron,Nyctanassa violacea,no,,,"","","","","","","","","","","","","","","","","","","","","" diff --git a/testthat/test-nests.R b/testthat/test-nests.R index aff1dd7..08e11b1 100644 --- a/testthat/test-nests.R +++ b/testthat/test-nests.R @@ -8,6 +8,8 @@ nests <- read.csv("../Nesting/nest_checks.csv", stringsAsFactors = FALSE) %>% nest_cols <- colnames(nests) success <- read.csv("../Nesting/nest_success.csv") success_cols <- colnames(success) +success_summary <- read.csv("../Nesting/nest_success_summary.csv") +success_summary_cols <- colnames(success_summary) test_that("required column names in nests df", { @@ -16,7 +18,13 @@ test_that("required column names in nests df", { expect_identical(success_cols, c("year","colony","nest_number","species","n_days_incubation","incubation_success", "n_days_nestling","nestling_success","clutch","brood","fledged","clutch_type", - "young_lost","real_success","real_failure","start_date","end_date","notes")) + "young_lost","real_success","real_failure","start_date","end_date","notes")) + expect_identical(success_summary_cols, + c("year","colony","species","incubation_k","incubation_sumy","incubation_e","incubation_p", + "incubation_j","incubation_pj","incubation_varp","incubation_varpj","incubation_sdp", + "incubation_sdpj","nestling_k","nestling_sumy","nestling_e","nestling_p","nestling_j", + "nestling_pj","nestling_varp","nestling_varpj","nestling_sdp","nestling_sdpj","overall_p", + "overall_varp","overall_sd")) }) test_that("Colony and species valid", { @@ -25,6 +33,8 @@ test_that("Colony and species valid", { expect_true(all(nests$species %in% species$species)) expect_true(all(success$colony %in% colonies$colony)) expect_true(all(success$species %in% species$species)) + expect_true(all(success_summary$colony %in% colonies$colony)) + expect_true(all(success_summary$species %in% species$species)) }) test_that("Dates valid", { @@ -32,6 +42,8 @@ test_that("Dates valid", { expect_false(any(is.na(nests$date))) expect_false(any(is.na(nests$year))) expect_false(any(is.na(success$year))) + expect_false(any(is.na(success_summary$year))) + expect_false(any(is.na(success_summary$colony))) }) test_that("Egg and chick counts valid", { @@ -40,7 +52,7 @@ test_that("Egg and chick counts valid", { expect_true(all(nests$chicks %in% c(0,1,2,3,4,5,6,7, NA))) }) -# test_that("no duplicated rows", { -# -# expect_false(any(duplicated(nests))) -# }) +test_that("no duplicated rows", { + + expect_false(any(duplicated(success_summary))) +})