Difference between revisions of "ISSS608 2017-18 T3 Assign Chan En Ying Grace Did Rose Pipit kick the bucket"
(147 intermediate revisions by the same user not shown) | |||
Line 15: | Line 15: | ||
| style="font-family:Arial; font-size:100%; solid #1B338F; background:#E6E6FA; text-align:center;" width="25%" | | | style="font-family:Arial; font-size:100%; solid #1B338F; background:#E6E6FA; text-align:center;" width="25%" | | ||
; | ; | ||
− | [[ISSS608_2017-18_T3_Assign_Chan_En_Ying_Grace_Did Rose Pipit | + | [[ISSS608_2017-18_T3_Assign_Chan_En_Ying_Grace_Did Rose Pipit kicketh the bucket?|<b><font size="3"><font color="#000000">Did Rose Pipit kicketh the bucket?</font></font></b>]] |
| style="font-family:Arial; font-size:100%; solid #1B338F; background:#E6E6FA; text-align:center;" width="25%" | | | style="font-family:Arial; font-size:100%; solid #1B338F; background:#E6E6FA; text-align:center;" width="25%" | | ||
Line 32: | Line 32: | ||
===1. Data Overview=== | ===1. Data Overview=== | ||
− | Different bird species were found throughout the preserve, but Pipits were found typically in the northeastern side of the preserve – that is, closest to the alleged Dumping Site marked by the triangle below. | + | Different bird species were found throughout the preserve (left), but Pipits were found typically in the northeastern side of the preserve – that is, closest to the alleged Dumping Site marked by the triangle below. |
− | We can also see a visible Pipit cluster | + | We can also see a visible Pipit cluster around the Dumping Site (right). This is likely why the alleged dumping by Kasios was “accused” to cause the death of the Pipits since the Pipits were found at that location. |
[[File:Part_1.png|500px]] | [[File:Part_1.png|500px]] | ||
Line 49: | Line 49: | ||
<hr/> | <hr/> | ||
===3. Assignment of Treatment & Control Groups=== | ===3. Assignment of Treatment & Control Groups=== | ||
− | If the dumping indeed took place, one would expect that <i> other birds in this dumping area would also decrease in population </i> as they would migrate elsewhere or die from the dumping substance, not just affecting the Rose Pipits. | + | If the dumping indeed took place, one would expect that <i> other birds species in this dumping area would also decrease in population </i> as they would migrate elsewhere or die from the dumping substance, not just affecting the Rose Pipits. |
− | Other species that are dominantly clustered around the dumping site include the Ordinary Snape and Lesser Birchbeere. | + | Other species that are dominantly clustered around the dumping site include the Ordinary Snape and Lesser Birchbeere. As such, these two species should also be affected by the dumping, if dumping indeed happened. |
+ | |||
+ | We will thus use <b> Ordinary Snape (OS) and Lesser Birchbeere (LB)</b> as our Control Group, and <b> Rose-Crested Blue Pipits (Rose Pipit) </b> as our Treatment Group. | ||
<hr/> | <hr/> | ||
Line 57: | Line 59: | ||
===4. Rose Pipits Only – Across the 4 Seasons=== | ===4. Rose Pipits Only – Across the 4 Seasons=== | ||
Just like other bird species, Rose Pipits population peaked in the Spring. When we analyse its population across time, we should bear in mind that seasonal or monthly population changes would not help us determine the true population patterns accurately as change in population is likely due to seasonal migration, than dumping. | Just like other bird species, Rose Pipits population peaked in the Spring. When we analyse its population across time, we should bear in mind that seasonal or monthly population changes would not help us determine the true population patterns accurately as change in population is likely due to seasonal migration, than dumping. | ||
− | Thus, we will track its population changes <b>across Year</b>, than across month or season. | + | |
+ | Thus, we will track its population changes <b>across Year</b>, than across month or season. (However, if we later wish to zoom into a particular year, we can look at month/season but bear in mind these effects.) | ||
[[File:Pic3.png|800px]] | [[File:Pic3.png|800px]] | ||
<hr/> | <hr/> | ||
− | ===5. Rose Pipits Across Time=== | + | ===5. Rose Pipits Only – Across Time=== |
There are two observations here. | There are two observations here. | ||
− | + | ====a. First, the Pipit population increased in 2013 to 2017.==== | |
− | + | In 2013, the cluster appeared right at the alleged dumping site. We call this the Pipit home range, where the Pipits tended to live and thrive. Population soared in 2015 and then declined gradually till 2017. | |
+ | |||
However, it is important to note that this is a short 5-year period and prior to 2013, Pipits were few and far between. One possible interpretation could be that Pipits only found the Boonseng Preserve as their home in 2013, and lived elsewhere prior to this. | However, it is important to note that this is a short 5-year period and prior to 2013, Pipits were few and far between. One possible interpretation could be that Pipits only found the Boonseng Preserve as their home in 2013, and lived elsewhere prior to this. | ||
− | Another possible interpretation could be due to data quality issue. Although it is mentioned that the data is well represented, an overview of the data showed signs of data quality, especially in the earlier years ( | + | Another possible interpretation could be due to data quality issue. Although it is mentioned that the data is well represented, an overview of the data showed signs of data quality, especially in the earlier years (1990s-2000s) where birds across all species were reported to be very few. Due to the incompleteness of data, we will focus on years 2012-2017 for a more reliable study. We will also omit year 2018 since the data provided is till end of Q1 only. |
[[File:Pic4.png|800px]] | [[File:Pic4.png|800px]] | ||
− | + | ====b. Second, the home range moved away from the dumping site, from 2015 onwards even though that was the year of peak. ==== | |
− | + | This could be signs of dumping. It could be that the dumping caused Pipits to fly slightly further away, but the long-run effect was that the birds population dwindled as they struggled to survive given the decline in population size (as the cluster got smaller). It could be that only those that flew away from the site survived. | |
− | But given that the Kasios company provided some Pipit bird calls, we should also analyse to see if it really belongs to Pipits, and whether these calls were obtained around the dumping site or far away from it. If far away, then it could support the dumping claim | + | But given that the Kasios company provided some Pipit bird calls, we should also analyse to see if it really belongs to Pipits, and whether these calls were obtained around the dumping site or far away from it. If far away, then it could support the dumping claim. |
<hr/> | <hr/> | ||
===6. Kernel Density Plot to visualize clusters (All Birds & Rose Pipit) === | ===6. Kernel Density Plot to visualize clusters (All Birds & Rose Pipit) === | ||
− | To visualize the above clusters more clearly, we plot the Kernel Density to visualize the | + | To visualize the above clusters more clearly, we plot the Kernel Density to visualize the cluster intensity, a method to compute the intensity of a point distribution. |
− | If we look at Rose Pipits (right), the Kernel Density plot gives us clearer visualization that there were 2 (not 1!) clusters formed. This could possibly indicate some sort of migration from one to the other. As the bottom cluster (away from dumping site) looks bigger than the top cluster (at dumping site), it could also be a sign of dumping, as dumping chemicals could cause the population to be smaller and | + | If we look at Rose Pipits (right), the Kernel Density plot gives us clearer visualization that there were 2 (not 1!) clusters formed. This could possibly indicate some sort of migration from one to the other. As the bottom cluster (the one located away from dumping site) looks bigger than the top cluster (at dumping site), it could also be a sign of dumping, as dumping chemicals could cause the population to be smaller and contributed to the Pipits' migration from its original homerange to a new homerange. |
− | To verify this, let us do a Kernel Density plot for Pipits across time for a six-year period from 2012-2017. | + | To verify this, let us do a Kernel Density plot for Rose Pipits across time for a six-year period from 2012-2017. |
[[File:Pic5.png|800px]] | [[File:Pic5.png|800px]] | ||
Line 90: | Line 94: | ||
===7. Kernel Density Plot of Rose Pipits Across Time (2012-2017)=== | ===7. Kernel Density Plot of Rose Pipits Across Time (2012-2017)=== | ||
− | Interestingly, when we plot | + | ====Two Clusters Become One; New Cluster "overtook" Old ==== |
+ | Interestingly, when we plot the Kernel Density across time, we can see that there were originally two clusters in 2012 – one at the dumping site, and one slightly away from the dumping site. The cluster at the dumping site grew in size <b>until</b> 2014, while the cluster away from the dumping site "took over" the growth and started to grow <b>from</b> 2014. | ||
+ | ====Old Cluster Vanished in 2015==== | ||
More importantly, the cluster at the dumping site vanished totally from 2015 onwards, leaving only the cluster away from the dumping site, which grew in size. This could indicate a movement of birds from the dumping site cluster towards the cluster away from the dumping site, possibly indicating a trend of dumping. | More importantly, the cluster at the dumping site vanished totally from 2015 onwards, leaving only the cluster away from the dumping site, which grew in size. This could indicate a movement of birds from the dumping site cluster towards the cluster away from the dumping site, possibly indicating a trend of dumping. | ||
Line 99: | Line 105: | ||
===8. Spatial Point Pattern Analysis (Across Time)=== | ===8. Spatial Point Pattern Analysis (Across Time)=== | ||
− | To more objectively determine the strength and significance of the clusters, we apply distance-based measures to test for cluster statistical significance, that could shed light towards the dumping allegation – the Quadrat Test, K-Nearest Neighbour as well as the K-Function. | + | To more objectively determine the strength and significance of the clusters, we apply density-based followed by distance-based measures to test for cluster statistical significance, that could shed light towards the dumping allegation – namely the Quadrat Test, K-Nearest Neighbour as well as the K-Function. |
+ | |||
+ | ====a. Quadrat Analysis==== | ||
− | + | =====Clusters still significant, but became less compact from 2015===== | |
− | Rose Pipit clusters were statistically significant at the 1% level of significance across 2012 to 2017, when a quadrat analysis was applied with the alternative hypothesis set as clustered. | + | Rose Pipit clusters were statistically significant at the 1% level of significance across 2012 to 2017, when a quadrat analysis was applied with the alternative hypothesis set as clustered. The Pearson X2 statistic peaked at 2015 and reduced after that, showing signs of reduction in cluster compactness (though still statistically significant). |
[[File:Table1.png|500px]] | [[File:Table1.png|500px]] | ||
+ | ====b. Nearest Neighbour, followed by Clarke-Evans Test==== | ||
+ | |||
+ | =====Clusters still significant, but started to disperse after 2015===== | ||
+ | |||
+ | Clusters were most compact in 2015, where majority of each Pipit’s neighbours were within 10units distance away. From 2016 however, the neighbouring distance expanded to beyond 10 towards 40. Though clusters were statistically significant from 2012 to 2017 (applying the Clarke-Evans test), the clusters started to disperse after 2015 which was the suspected year of dumping. This supports the same findings above which saw that year 2015 had the strongest/most compact cluster. | ||
− | + | [[File:Pic9.png|1000px]] | |
− | + | Note: When we conduct the Clarke-Evans test, though all clusters were significant, we see that the coefficient was highest in 2016 (not 2015). Nevertheless, 2015 saw the most compact cluster (smallest distance between neighbours). | |
[[File:Table2.png|500px]] | [[File:Table2.png|500px]] | ||
− | [[File: | + | ====c. K-Function==== |
+ | |||
+ | =====Cluster lost its significance in 2017, for radius <5===== | ||
+ | |||
+ | Likewise, when a K-Function was applied, the size/strength of clusters reduced overtime after year 2015 as the grey shaded envelope increased. | ||
+ | |||
+ | Eventually in 2017, the cluster lost its significance when radius was below 5. Clusters were significant only when the radius was 5 or greater. In other words, the cluster became less compact in 2017. | ||
+ | |||
+ | That said, we note that in year 2013 and 2014, the clusters were not as strongly significant for small values of radius (i.e. compact clusters not significant). Ceteris paribus, this is because there were two clusters if we refer to the earlier Kernel Density plot, rather than a decline in Pipit population, so this is not a cause for concern. | ||
+ | |||
+ | [[File:Pic10.png|1500px]] | ||
+ | |||
+ | <hr/> | ||
+ | <div align="center">[[File:Mid1.png|800px]]</div> | ||
+ | |||
+ | <hr/> | ||
+ | |||
+ | ===9. Did Rose Pipits die from the dumping chemical? If not, where did they flock to?=== | ||
+ | Rose Pipits could either have left the Preserve as a flight mechanism, or that they were defeated and simply died and dwindled in population due to the chemical. | ||
+ | |||
+ | I don’t think it had left for another place, but rather, I think it slowly dwindled in population (i.e. Rose Pipits kicked the bucket). Because if it left, we should see a trend of them flying away from the Dumping Site and moving towards an edge of the Preserve to migrate to another area outside of the Preserve, for example. It was unlikely for the flight to take place immediately and we should see a staggered flight pattern across the months as birds migrate elsewhere. | ||
+ | |||
+ | ====Facet Grid Plot==== | ||
+ | To test my hypothesis, it'd be ideal to do a Kernel Density plot by month for 2016 and 2017. Unfortunately, the plot is not possible for small data points (as some months only have 1-2 birds of data). So let’s start off with a Facet Grid Plot to see if there’s a trend where the birds are migrating. Later on, we will do the Kernel Density Plot but by quarter. | ||
+ | |||
+ | [[File:Pic11.png|700px]] | ||
+ | |||
+ | The pattern is not so clear in 2016. Birds moved away from the dumping site but also dwindled at that location. | ||
+ | |||
+ | However, in 2017, this pattern is clearer. The remaining birds were flocking to the South East of Kasios Island, even in Spring (April). This could be because their alternate homerange is in another region of the country, southeast of Kasios. | ||
+ | |||
+ | [[File:Pic12.png|700px]] | ||
+ | |||
+ | ====Kernel Density Plot – by Quarter==== | ||
+ | |||
+ | However, this is not that convincing. Let’s attempt to run a Kernel Density Plot - instead of by month (which we lack data of some months especially in the 2nd half of the year), we plot the kernel density for each quarter in 2016 and 2017. | ||
+ | |||
+ | Note that Q3 2016 and Q3 2017's data is too few for a Kernel Density plot, so we have no choice but to omit Q3 from our analysis. | ||
+ | |||
+ | [[File:Pic13.png|1000px]] | ||
+ | |||
+ | ====Why did Pipits fly back to the dumping site every Winter?==== | ||
+ | There is no pattern of the birds flying to the Southeast of the Preserve. | ||
+ | |||
+ | But interestingly, the birds actually flew back <b>closer</b> to the dumping site every <b>end of the year (winter)</b>. | ||
+ | |||
+ | Could it be that the dumping was not once-off, but during certain parts of the year e.g. did dumping took place in Q3, explaining the absence of birds? Or was the absence of birds due to seasonal changes? | ||
+ | |||
+ | Well, when we look at the Pipit distribution by season in the earlier graph, we can see that Summer should be the 2nd highest peak season of Pipit population. There is no reason why it should be 0 in population, which is the case now since Q3 lacked data. Just possibly, this could be signs of dumping takes place every year in Spring/Summer from 2015 onwards. | ||
+ | |||
+ | However, I must note that the data points are too few to analyse the birds in 2016 and 2017 even by quarter and make such concusions. Even by visualizing the facet plot, it is clear that the KD plot can be not representative because of the few data points. Instead, let us focus on our Control Groups to shed greater light onto the fate of our Rose Pipits. | ||
+ | |||
+ | |||
+ | <hr/> | ||
+ | |||
+ | ===10. Were Rose Pipits the only ones affected? === | ||
+ | Rose Pipits could either have left the Preserve as a flight mechanism, or that they were defeated and simply died and dwindled in population due to the chemical. | ||
+ | |||
+ | ====Treatment vs Control Groups==== | ||
+ | If there were dumping, even our Control Groups should also be affected unless there were biological reasons why they were “immune” to the chemical. | ||
+ | Interestingly, our Control Groups showed that they were <b>not</b> as affected by the Dumping Site. The OS & LB were chosen as control groups because they were found to be clustered around or near the Dumping Site. | ||
+ | |||
+ | ====a) Ordinary Snape (Control Group 1)==== | ||
+ | 1. For the Ordinary Snape birds, their homerange remained near the Dumping Site, unlike the Rose Pipits whose homerange moved away from the dumping site. | ||
+ | |||
+ | 2. Moreover, the Ordinary Snape increased in population, unlike the Rose Pipits that decreased. | ||
+ | |||
+ | ====Facet Plot==== | ||
+ | This could be a sign that dumping did not take place since the OS were not affected by it. However, we cannot rule out dumping. Because this could be also demonstrate Darwin's evolution theory of Survival of the Fittest, where the OS were less affected by the chemical of the dumping site, and were the stronger species that survived since it did not seem to be as affected in the same period. | ||
+ | |||
+ | [[File:Pic14.png|800px]] | ||
+ | |||
+ | ====Kernel Density Plot==== | ||
+ | To objectively confirm this, we apply the Kernel Density Plot across 2012-2017 to confirm the significance and size of the clusters. True enough, the homerange did not move and the clusters did increase very slightly in 2016-2017 as can be seen by the smudging of the kernel density plot towards the dumping site. | ||
+ | |||
+ | [[File:Pic15.png|1000px]] | ||
+ | |||
+ | ====K-Function==== | ||
+ | The significance of clusters also grew from 2012 to 2017, thus confirming our hypothesis that the <b> Ordinary Snapes were not affected by the dumping</b>, if dumping did take place. | ||
+ | |||
+ | [[File:KCSR1.png|1250px]] | ||
+ | |||
+ | ====K-Cross==== | ||
+ | For more rigorous analysis, we will also perform the K-Cross between the Rose Pipit and Ordinary Snape. | ||
+ | |||
+ | The aim of the bivariate analysis is to identify if the spatial distribution of Ordinary Snape is independent of the distribution of Rose Pipit i.e. our null hypothesis. An estimate of <b>K</b>ordinarysnape,rosepipit(r) is a useful summary statistic in exploratory data analysis of a multitype point pattern. If the migratory behaviour of Ordinary Snape were independent of the migratory behaviour of Rose Pipits, then <b>K</b>ordinarysnape,rosepipit(r) would equal πr^2. Deviations between the empirical <b>K</b>ordinarysnape,rosepipit(r) curve and the theoretical curve πr^2 may suggest dependence between the points of Ordinary Snape and Pipits. | ||
+ | |||
+ | [[File:KC_OS.png|700px]] | ||
+ | |||
+ | ====="Close" Spatial Dependence found between Pipits & OS===== | ||
+ | |||
+ | From above, our observed K-cross curve above the theoretical curve means that type j points (Ordinary Snapes) prefer to cluster around type i points (Rose Pipits). That is, Ordinary Snapes are found closer to Rose Pipits than would be expected under complete spatial randomness. | ||
+ | |||
+ | It looks like there is strong spatial dependence between the 2 species (i.e. they are located close together), but let's look at it by Year to better understand how this dependence changes across time, in relation to the suspected dumping year of 2015. | ||
+ | |||
+ | ====K-Cross Across Years (2012-2017)==== | ||
+ | |||
+ | [[File:KCROSSOS.png|1250px]] | ||
+ | |||
+ | When we compare across years from 2012-2017, we can see that the spatial dependence between Pipits & OS was strongest in 2015. If we recall, 2015 was the year of population peak for the Rose Pipits. | ||
+ | |||
+ | In 2016, these two species were still spatially dependent. But in 2017, they were spatially dependent only up till radius = 40. | ||
+ | |||
+ | ====Evidence Against Dumping==== | ||
+ | What does this imply? One implication is that the 2 species were spatially close, but yet, Ordinary Snapes thrived and Pipits started to dwindle. If dumping too place, then the OS must be affected too, especially in 2015 since the 2 species were spatially dependent! However, the OS were not affected in population and homerange. Thus, we cannot say for sure that dumping too place. Well, possibly, the chemical (released in 2015) only affected the Rose Pipits. | ||
+ | |||
+ | =====Evidence For Dumping===== | ||
+ | Another observation is as follows. Since 2015 was the year when the 2 species were spatially closest and also the year where Rose Pipits moved away from the dumping site, coupled with the fact that the home range of the OS was not as close to the dumping site as compared to the Rose Pipits originally, it is perhaps only normal that the spatial dependence is closest in 2015 if dumping indeed affected the Pipits. This is because the Pipits moved from the dumping site and hence is now closer to the OS homerange, thus a stronger K-Cross significance. | ||
+ | |||
+ | We also see that the strength of the spatial dependence fell after 2015 (though still statistically significant) which meant that perhaps there was some slow recovery of the Rose Pipits, after being affected by the 2015 dumping. | ||
+ | |||
+ | <hr/> | ||
+ | |||
+ | ====b) Lesser Birchbeere (Control Group 2)==== | ||
+ | 1. As for the Lesser Birchbeere, its population dwindled slightly. | ||
+ | <br>2. But its homerange moved even closer towards the Dumping Site, as can be seen by the Kernel Density plot. | ||
+ | |||
+ | Since the deaths were only limited to Rose Pipits, one hypothesis is that the dumping did not take place. | ||
+ | |||
+ | Another hypothesis is that the dumping did take place, but the chemicals being released were only affecting the Pipits, where the other species survived and some even thrived. Possibly because Pipit population was reduced, there could be less competition for food so other species thrived over there. | ||
+ | |||
+ | ====Facet Plot==== | ||
+ | |||
+ | [[File:Pic17.png|800px]] | ||
+ | |||
+ | ====Kernel Density Plot==== | ||
+ | |||
+ | To objectively confirm this, we apply the Kernel Density Plot to confirm the significance and size of the clusters. True enough, the homerange grew closer to the dumping site, whereby a new cluster was formed near the dumping site in 2016 and grew almost as large as the former cluster, in 2017. | ||
+ | |||
+ | [[File:Pic18.png|800px]] | ||
+ | |||
+ | ====K Function==== | ||
+ | However, the clusters were not statistically significant in 2016 (until radius > 15). This could be because of the emergence of a new cluster, as can be seen in the plot. That said, the clusters became significant in 2017 when the new cluster grew in size. We thus still stick to our hypothesis that the <b> Lesser Birchbeeres were not affected by the dumping </b>, if dumping did take place. | ||
+ | |||
+ | [[File:KCSR2.png|1250px]] | ||
+ | |||
+ | ====K-Cross==== | ||
+ | |||
+ | For more rigorous analysis, we will also perform the K-Cross between the Rose Pipit and Lesser Birchbeere, to identify if the spatial distribution of both species are independent. The below graph shows complete spatial randomness between the Pipits and LBs i.e. they are <i>not</i> located close together. But let's analyse this across time to see if there is any change. | ||
+ | |||
+ | [[File:KC_LB.png|700px]] | ||
+ | |||
+ | =====Complete Spatial Randomness between Pipit & LB===== | ||
+ | Contrary to Ordinary Snapes, the observed Kcross-function for Lesser Birchbeere is within the confidence envelope (for r<40). This means that there is complete spatial randomness between Pipits & LBs for small distances, within radius of 40. | ||
+ | |||
+ | Note: However, for radius > 40, there is some statistical significance where LBs and Pipits were spatially dispersed. I.e. they are located far away from one together. Does this mean they had some 'repelling' effect and distanced themselves from one another? But, since the line is almost touching the grey envelope, this is not a strong finding and can be taken lightly. | ||
+ | |||
+ | ====K-Cross Across Years (2012-2017)==== | ||
+ | Even across time, the LBs displayed complete spatial randomness in relation to the Rose Pipits, across all 6 years. That is, the Rose Pipits and LBs were never "spatially close" together in statistical terms. | ||
+ | |||
+ | This sheds light to whether LBs were a good control to select in the first place. We see that in 2012, there was already spatial randomness between LBs and Rose Pipits. By visualising the facet plot of species across time in the earlier part of this analysis, we had earlier selected the Lesser Birchbeere (and the Ordinary Snapes) as our Control Groups - as we assumed that they were located close to the dumping site based on visualisation. In other words, if that assumption holds, LBs actually should display some spatial <i>dependence</i> with the Rose Pipits (at least in the earlier years), since the Pipits were also located at the dumping site. | ||
+ | |||
+ | However, the K-Cross showed us that the LBs did not have spatial dependence with the Rose Pipits across <i>all</i> 6 years, including 2012, and hence LB could be a poor control. Nevertheless, our OSs still serve as good controls so we will place greater emphasis on the results arising from the OSs. | ||
+ | |||
+ | |||
+ | [[File:KCROSSLB.png|1250px]] | ||
+ | |||
+ | <hr/> | ||
+ | |||
+ | ====c) All Species (Control 3 – External Factors)==== | ||
+ | To see whether there was another external event happening on Kasios Island that was causing any decline of bird species, I also plotted the population across all birds across time. Since overall population of the birds increased, there is unlikely to be any external event happening. | ||
+ | |||
+ | [[File:Pic20.png|800px]] | ||
+ | |||
+ | <hr/> | ||
+ | |||
+ | ===11. Type of Bird Collection: Perhaps Bird Song is indicative of thriving population, while Bird Call is a sign of distress?=== | ||
+ | |||
+ | Some research suggests that bird songs could be indicative of happier birds while calls are signs of distress. Let’s plot the composition of calls and songs for Rose Pipits, Ordinary Snape as well as Lesser Birchbeere. | ||
+ | |||
+ | True enough, Rose Pipits stopped singing and there were more signals of distress (calls) from 2015 to 2017. This could be because of the dumping. When there was dumping, the songs turned to calls, especially since 2015 was the year where we see the most calls. | ||
+ | |||
+ | [[File:PipitStopSinging.png|800px]] | ||
+ | |||
+ | ====K-Cross between Happy & Distressed Pipits==== | ||
+ | To confirm this, let us do a K-Cross between <b>Pipits who called</b> and <b>Pipits who sang</b>. As hypothesised, year 2015 which was the suspected year of dumping, saw the most statistically significant spatial dependence between Pipits who called and Pipits who sang i.e. calling and singing Pipits are located closely together. This supports our hypothesis as during the year of dumping, the songs turned to calls. | ||
+ | |||
+ | [[File:SONG1.png|1250px]] | ||
− | + | For completeness, we also look at the song-call composition of our controls. Though, the LB birds and OS birds would not be good controls for bird call because LB birds never sang and OS birds did not sing much to begin with in the earlier years, so we can't see whether they were more distressed than before. | |
− | + | [[File:LSneversang.png|800px]] | |
− | + | [[File:OSdontsingmuch.png|800px]] | |
− | + | ==<b>Concluding Hypothesis: Pipits surviving, but not thriving on the Preserve.</b>== | |
+ | ====Key Observations==== | ||
+ | 1. Pipit clusters were significant in 2012 to 2017. | ||
+ | <br>2. Pipit population peaked in 2015, but fell from 2016. | ||
+ | <br>3. Pipit homerange moved away from dumping site, from 2015 | ||
+ | <br>4. Pipit clusters became less compact from 2015 and lost its significance in 2017 for radius < 5 | ||
+ | <br>5. Pipits stopped singing after 2015. Songs turned into Calls - a sign of distress. | ||
+ | <br>6. Pipits were the only species affected (i.e. their homerange & population). | ||
+ | <br>7. Control Groups (Ordinary Snape & Lesser Birchbeere) thrived and even have their homerange move closer to the dumping site | ||
+ | <br>8. Pipits were spatially dependent to the Ordinary Snapes in 2015, so they should be both exposed and given the same "treatment" if dumping did actually occur, however the Ordinary Snapes were not affected. | ||
− | + | ====Hypothesis: Signs of Dumping But Affected Pipits Only==== | |
+ | Rose Pipits were still surviving as its clusters still exist and were significant, until 2017. However, they were not thriving at the dumping site and had to move away from it. Moreover, its population had fallen. This was especially since 2015 was the year that songs turned to calls, moreover, the control groups – OS and LB – did not experience a fall in population and in fact even increased in population and moved closer to the dumping site from 2015, respectively. | ||
− | + | So, I conclude that there were signs of Dumping and this was likely to take place in 2015, but the Dumping most likely consisted of chemicals that affected mainly the Rose Pipits, and not the other species. |
Latest revision as of 01:06, 1 July 2018
|
|
|
|
|
Contents
- 1 Cluster Visualisation: "Did Rose Pipit actually kicketh the bucket?"
- 1.1 1. Data Overview
- 1.2 2. Clustering by Species
- 1.3 3. Assignment of Treatment & Control Groups
- 1.4 4. Rose Pipits Only – Across the 4 Seasons
- 1.5 5. Rose Pipits Only – Across Time
- 1.6 6. Kernel Density Plot to visualize clusters (All Birds & Rose Pipit)
- 1.7 7. Kernel Density Plot of Rose Pipits Across Time (2012-2017)
- 1.8 8. Spatial Point Pattern Analysis (Across Time)
- 1.9 9. Did Rose Pipits die from the dumping chemical? If not, where did they flock to?
- 1.10 10. Were Rose Pipits the only ones affected?
- 1.10.1 Treatment vs Control Groups
- 1.10.2 a) Ordinary Snape (Control Group 1)
- 1.10.3 Facet Plot
- 1.10.4 Kernel Density Plot
- 1.10.5 K-Function
- 1.10.6 K-Cross
- 1.10.7 K-Cross Across Years (2012-2017)
- 1.10.8 Evidence Against Dumping
- 1.10.9 b) Lesser Birchbeere (Control Group 2)
- 1.10.10 Facet Plot
- 1.10.11 Kernel Density Plot
- 1.10.12 K Function
- 1.10.13 K-Cross
- 1.10.14 K-Cross Across Years (2012-2017)
- 1.10.15 c) All Species (Control 3 – External Factors)
- 1.11 11. Type of Bird Collection: Perhaps Bird Song is indicative of thriving population, while Bird Call is a sign of distress?
- 2 Concluding Hypothesis: Pipits surviving, but not thriving on the Preserve.
Cluster Visualisation: "Did Rose Pipit actually kicketh the bucket?"
1. Data Overview
Different bird species were found throughout the preserve (left), but Pipits were found typically in the northeastern side of the preserve – that is, closest to the alleged Dumping Site marked by the triangle below.
We can also see a visible Pipit cluster around the Dumping Site (right). This is likely why the alleged dumping by Kasios was “accused” to cause the death of the Pipits since the Pipits were found at that location.
2. Clustering by Species
To characterize the pattern relative to the alleged dumping site (blue triangle), the following facet plot shows the clustering patterns by Bird Species, independent of time.
From the cluster locations, we can see that the Rose-Crested Blue Pipits are clustered around the alleged dumping site, explaining why the 15 test birds provided are claimed to be Rose-Crested Blue Pipits to prove the allegation wrong. We can use this facet plot to have a good overview of other bird species’ cluster locations, such as the Ordinary Snape and Lesser Birchbeere being also clustered near to the Dumping Site.
3. Assignment of Treatment & Control Groups
If the dumping indeed took place, one would expect that other birds species in this dumping area would also decrease in population as they would migrate elsewhere or die from the dumping substance, not just affecting the Rose Pipits.
Other species that are dominantly clustered around the dumping site include the Ordinary Snape and Lesser Birchbeere. As such, these two species should also be affected by the dumping, if dumping indeed happened.
We will thus use Ordinary Snape (OS) and Lesser Birchbeere (LB) as our Control Group, and Rose-Crested Blue Pipits (Rose Pipit) as our Treatment Group.
4. Rose Pipits Only – Across the 4 Seasons
Just like other bird species, Rose Pipits population peaked in the Spring. When we analyse its population across time, we should bear in mind that seasonal or monthly population changes would not help us determine the true population patterns accurately as change in population is likely due to seasonal migration, than dumping.
Thus, we will track its population changes across Year, than across month or season. (However, if we later wish to zoom into a particular year, we can look at month/season but bear in mind these effects.)
5. Rose Pipits Only – Across Time
There are two observations here.
a. First, the Pipit population increased in 2013 to 2017.
In 2013, the cluster appeared right at the alleged dumping site. We call this the Pipit home range, where the Pipits tended to live and thrive. Population soared in 2015 and then declined gradually till 2017.
However, it is important to note that this is a short 5-year period and prior to 2013, Pipits were few and far between. One possible interpretation could be that Pipits only found the Boonseng Preserve as their home in 2013, and lived elsewhere prior to this. Another possible interpretation could be due to data quality issue. Although it is mentioned that the data is well represented, an overview of the data showed signs of data quality, especially in the earlier years (1990s-2000s) where birds across all species were reported to be very few. Due to the incompleteness of data, we will focus on years 2012-2017 for a more reliable study. We will also omit year 2018 since the data provided is till end of Q1 only.
b. Second, the home range moved away from the dumping site, from 2015 onwards even though that was the year of peak.
This could be signs of dumping. It could be that the dumping caused Pipits to fly slightly further away, but the long-run effect was that the birds population dwindled as they struggled to survive given the decline in population size (as the cluster got smaller). It could be that only those that flew away from the site survived. But given that the Kasios company provided some Pipit bird calls, we should also analyse to see if it really belongs to Pipits, and whether these calls were obtained around the dumping site or far away from it. If far away, then it could support the dumping claim.
6. Kernel Density Plot to visualize clusters (All Birds & Rose Pipit)
To visualize the above clusters more clearly, we plot the Kernel Density to visualize the cluster intensity, a method to compute the intensity of a point distribution.
If we look at Rose Pipits (right), the Kernel Density plot gives us clearer visualization that there were 2 (not 1!) clusters formed. This could possibly indicate some sort of migration from one to the other. As the bottom cluster (the one located away from dumping site) looks bigger than the top cluster (at dumping site), it could also be a sign of dumping, as dumping chemicals could cause the population to be smaller and contributed to the Pipits' migration from its original homerange to a new homerange.
To verify this, let us do a Kernel Density plot for Rose Pipits across time for a six-year period from 2012-2017.
7. Kernel Density Plot of Rose Pipits Across Time (2012-2017)
Two Clusters Become One; New Cluster "overtook" Old
Interestingly, when we plot the Kernel Density across time, we can see that there were originally two clusters in 2012 – one at the dumping site, and one slightly away from the dumping site. The cluster at the dumping site grew in size until 2014, while the cluster away from the dumping site "took over" the growth and started to grow from 2014.
Old Cluster Vanished in 2015
More importantly, the cluster at the dumping site vanished totally from 2015 onwards, leaving only the cluster away from the dumping site, which grew in size. This could indicate a movement of birds from the dumping site cluster towards the cluster away from the dumping site, possibly indicating a trend of dumping.
8. Spatial Point Pattern Analysis (Across Time)
To more objectively determine the strength and significance of the clusters, we apply density-based followed by distance-based measures to test for cluster statistical significance, that could shed light towards the dumping allegation – namely the Quadrat Test, K-Nearest Neighbour as well as the K-Function.
a. Quadrat Analysis
Clusters still significant, but became less compact from 2015
Rose Pipit clusters were statistically significant at the 1% level of significance across 2012 to 2017, when a quadrat analysis was applied with the alternative hypothesis set as clustered. The Pearson X2 statistic peaked at 2015 and reduced after that, showing signs of reduction in cluster compactness (though still statistically significant).
b. Nearest Neighbour, followed by Clarke-Evans Test
Clusters still significant, but started to disperse after 2015
Clusters were most compact in 2015, where majority of each Pipit’s neighbours were within 10units distance away. From 2016 however, the neighbouring distance expanded to beyond 10 towards 40. Though clusters were statistically significant from 2012 to 2017 (applying the Clarke-Evans test), the clusters started to disperse after 2015 which was the suspected year of dumping. This supports the same findings above which saw that year 2015 had the strongest/most compact cluster.
Note: When we conduct the Clarke-Evans test, though all clusters were significant, we see that the coefficient was highest in 2016 (not 2015). Nevertheless, 2015 saw the most compact cluster (smallest distance between neighbours).
c. K-Function
Cluster lost its significance in 2017, for radius <5
Likewise, when a K-Function was applied, the size/strength of clusters reduced overtime after year 2015 as the grey shaded envelope increased.
Eventually in 2017, the cluster lost its significance when radius was below 5. Clusters were significant only when the radius was 5 or greater. In other words, the cluster became less compact in 2017.
That said, we note that in year 2013 and 2014, the clusters were not as strongly significant for small values of radius (i.e. compact clusters not significant). Ceteris paribus, this is because there were two clusters if we refer to the earlier Kernel Density plot, rather than a decline in Pipit population, so this is not a cause for concern.
9. Did Rose Pipits die from the dumping chemical? If not, where did they flock to?
Rose Pipits could either have left the Preserve as a flight mechanism, or that they were defeated and simply died and dwindled in population due to the chemical.
I don’t think it had left for another place, but rather, I think it slowly dwindled in population (i.e. Rose Pipits kicked the bucket). Because if it left, we should see a trend of them flying away from the Dumping Site and moving towards an edge of the Preserve to migrate to another area outside of the Preserve, for example. It was unlikely for the flight to take place immediately and we should see a staggered flight pattern across the months as birds migrate elsewhere.
Facet Grid Plot
To test my hypothesis, it'd be ideal to do a Kernel Density plot by month for 2016 and 2017. Unfortunately, the plot is not possible for small data points (as some months only have 1-2 birds of data). So let’s start off with a Facet Grid Plot to see if there’s a trend where the birds are migrating. Later on, we will do the Kernel Density Plot but by quarter.
The pattern is not so clear in 2016. Birds moved away from the dumping site but also dwindled at that location.
However, in 2017, this pattern is clearer. The remaining birds were flocking to the South East of Kasios Island, even in Spring (April). This could be because their alternate homerange is in another region of the country, southeast of Kasios.
Kernel Density Plot – by Quarter
However, this is not that convincing. Let’s attempt to run a Kernel Density Plot - instead of by month (which we lack data of some months especially in the 2nd half of the year), we plot the kernel density for each quarter in 2016 and 2017.
Note that Q3 2016 and Q3 2017's data is too few for a Kernel Density plot, so we have no choice but to omit Q3 from our analysis.
Why did Pipits fly back to the dumping site every Winter?
There is no pattern of the birds flying to the Southeast of the Preserve.
But interestingly, the birds actually flew back closer to the dumping site every end of the year (winter).
Could it be that the dumping was not once-off, but during certain parts of the year e.g. did dumping took place in Q3, explaining the absence of birds? Or was the absence of birds due to seasonal changes?
Well, when we look at the Pipit distribution by season in the earlier graph, we can see that Summer should be the 2nd highest peak season of Pipit population. There is no reason why it should be 0 in population, which is the case now since Q3 lacked data. Just possibly, this could be signs of dumping takes place every year in Spring/Summer from 2015 onwards.
However, I must note that the data points are too few to analyse the birds in 2016 and 2017 even by quarter and make such concusions. Even by visualizing the facet plot, it is clear that the KD plot can be not representative because of the few data points. Instead, let us focus on our Control Groups to shed greater light onto the fate of our Rose Pipits.
10. Were Rose Pipits the only ones affected?
Rose Pipits could either have left the Preserve as a flight mechanism, or that they were defeated and simply died and dwindled in population due to the chemical.
Treatment vs Control Groups
If there were dumping, even our Control Groups should also be affected unless there were biological reasons why they were “immune” to the chemical. Interestingly, our Control Groups showed that they were not as affected by the Dumping Site. The OS & LB were chosen as control groups because they were found to be clustered around or near the Dumping Site.
a) Ordinary Snape (Control Group 1)
1. For the Ordinary Snape birds, their homerange remained near the Dumping Site, unlike the Rose Pipits whose homerange moved away from the dumping site.
2. Moreover, the Ordinary Snape increased in population, unlike the Rose Pipits that decreased.
Facet Plot
This could be a sign that dumping did not take place since the OS were not affected by it. However, we cannot rule out dumping. Because this could be also demonstrate Darwin's evolution theory of Survival of the Fittest, where the OS were less affected by the chemical of the dumping site, and were the stronger species that survived since it did not seem to be as affected in the same period.
Kernel Density Plot
To objectively confirm this, we apply the Kernel Density Plot across 2012-2017 to confirm the significance and size of the clusters. True enough, the homerange did not move and the clusters did increase very slightly in 2016-2017 as can be seen by the smudging of the kernel density plot towards the dumping site.
K-Function
The significance of clusters also grew from 2012 to 2017, thus confirming our hypothesis that the Ordinary Snapes were not affected by the dumping, if dumping did take place.
K-Cross
For more rigorous analysis, we will also perform the K-Cross between the Rose Pipit and Ordinary Snape.
The aim of the bivariate analysis is to identify if the spatial distribution of Ordinary Snape is independent of the distribution of Rose Pipit i.e. our null hypothesis. An estimate of Kordinarysnape,rosepipit(r) is a useful summary statistic in exploratory data analysis of a multitype point pattern. If the migratory behaviour of Ordinary Snape were independent of the migratory behaviour of Rose Pipits, then Kordinarysnape,rosepipit(r) would equal πr^2. Deviations between the empirical Kordinarysnape,rosepipit(r) curve and the theoretical curve πr^2 may suggest dependence between the points of Ordinary Snape and Pipits.
"Close" Spatial Dependence found between Pipits & OS
From above, our observed K-cross curve above the theoretical curve means that type j points (Ordinary Snapes) prefer to cluster around type i points (Rose Pipits). That is, Ordinary Snapes are found closer to Rose Pipits than would be expected under complete spatial randomness.
It looks like there is strong spatial dependence between the 2 species (i.e. they are located close together), but let's look at it by Year to better understand how this dependence changes across time, in relation to the suspected dumping year of 2015.
K-Cross Across Years (2012-2017)
When we compare across years from 2012-2017, we can see that the spatial dependence between Pipits & OS was strongest in 2015. If we recall, 2015 was the year of population peak for the Rose Pipits.
In 2016, these two species were still spatially dependent. But in 2017, they were spatially dependent only up till radius = 40.
Evidence Against Dumping
What does this imply? One implication is that the 2 species were spatially close, but yet, Ordinary Snapes thrived and Pipits started to dwindle. If dumping too place, then the OS must be affected too, especially in 2015 since the 2 species were spatially dependent! However, the OS were not affected in population and homerange. Thus, we cannot say for sure that dumping too place. Well, possibly, the chemical (released in 2015) only affected the Rose Pipits.
Evidence For Dumping
Another observation is as follows. Since 2015 was the year when the 2 species were spatially closest and also the year where Rose Pipits moved away from the dumping site, coupled with the fact that the home range of the OS was not as close to the dumping site as compared to the Rose Pipits originally, it is perhaps only normal that the spatial dependence is closest in 2015 if dumping indeed affected the Pipits. This is because the Pipits moved from the dumping site and hence is now closer to the OS homerange, thus a stronger K-Cross significance.
We also see that the strength of the spatial dependence fell after 2015 (though still statistically significant) which meant that perhaps there was some slow recovery of the Rose Pipits, after being affected by the 2015 dumping.
b) Lesser Birchbeere (Control Group 2)
1. As for the Lesser Birchbeere, its population dwindled slightly.
2. But its homerange moved even closer towards the Dumping Site, as can be seen by the Kernel Density plot.
Since the deaths were only limited to Rose Pipits, one hypothesis is that the dumping did not take place.
Another hypothesis is that the dumping did take place, but the chemicals being released were only affecting the Pipits, where the other species survived and some even thrived. Possibly because Pipit population was reduced, there could be less competition for food so other species thrived over there.
Facet Plot
Kernel Density Plot
To objectively confirm this, we apply the Kernel Density Plot to confirm the significance and size of the clusters. True enough, the homerange grew closer to the dumping site, whereby a new cluster was formed near the dumping site in 2016 and grew almost as large as the former cluster, in 2017.
K Function
However, the clusters were not statistically significant in 2016 (until radius > 15). This could be because of the emergence of a new cluster, as can be seen in the plot. That said, the clusters became significant in 2017 when the new cluster grew in size. We thus still stick to our hypothesis that the Lesser Birchbeeres were not affected by the dumping , if dumping did take place.
K-Cross
For more rigorous analysis, we will also perform the K-Cross between the Rose Pipit and Lesser Birchbeere, to identify if the spatial distribution of both species are independent. The below graph shows complete spatial randomness between the Pipits and LBs i.e. they are not located close together. But let's analyse this across time to see if there is any change.
Complete Spatial Randomness between Pipit & LB
Contrary to Ordinary Snapes, the observed Kcross-function for Lesser Birchbeere is within the confidence envelope (for r<40). This means that there is complete spatial randomness between Pipits & LBs for small distances, within radius of 40.
Note: However, for radius > 40, there is some statistical significance where LBs and Pipits were spatially dispersed. I.e. they are located far away from one together. Does this mean they had some 'repelling' effect and distanced themselves from one another? But, since the line is almost touching the grey envelope, this is not a strong finding and can be taken lightly.
K-Cross Across Years (2012-2017)
Even across time, the LBs displayed complete spatial randomness in relation to the Rose Pipits, across all 6 years. That is, the Rose Pipits and LBs were never "spatially close" together in statistical terms.
This sheds light to whether LBs were a good control to select in the first place. We see that in 2012, there was already spatial randomness between LBs and Rose Pipits. By visualising the facet plot of species across time in the earlier part of this analysis, we had earlier selected the Lesser Birchbeere (and the Ordinary Snapes) as our Control Groups - as we assumed that they were located close to the dumping site based on visualisation. In other words, if that assumption holds, LBs actually should display some spatial dependence with the Rose Pipits (at least in the earlier years), since the Pipits were also located at the dumping site.
However, the K-Cross showed us that the LBs did not have spatial dependence with the Rose Pipits across all 6 years, including 2012, and hence LB could be a poor control. Nevertheless, our OSs still serve as good controls so we will place greater emphasis on the results arising from the OSs.
c) All Species (Control 3 – External Factors)
To see whether there was another external event happening on Kasios Island that was causing any decline of bird species, I also plotted the population across all birds across time. Since overall population of the birds increased, there is unlikely to be any external event happening.
11. Type of Bird Collection: Perhaps Bird Song is indicative of thriving population, while Bird Call is a sign of distress?
Some research suggests that bird songs could be indicative of happier birds while calls are signs of distress. Let’s plot the composition of calls and songs for Rose Pipits, Ordinary Snape as well as Lesser Birchbeere.
True enough, Rose Pipits stopped singing and there were more signals of distress (calls) from 2015 to 2017. This could be because of the dumping. When there was dumping, the songs turned to calls, especially since 2015 was the year where we see the most calls.
K-Cross between Happy & Distressed Pipits
To confirm this, let us do a K-Cross between Pipits who called and Pipits who sang. As hypothesised, year 2015 which was the suspected year of dumping, saw the most statistically significant spatial dependence between Pipits who called and Pipits who sang i.e. calling and singing Pipits are located closely together. This supports our hypothesis as during the year of dumping, the songs turned to calls.
For completeness, we also look at the song-call composition of our controls. Though, the LB birds and OS birds would not be good controls for bird call because LB birds never sang and OS birds did not sing much to begin with in the earlier years, so we can't see whether they were more distressed than before.
Concluding Hypothesis: Pipits surviving, but not thriving on the Preserve.
Key Observations
1. Pipit clusters were significant in 2012 to 2017.
2. Pipit population peaked in 2015, but fell from 2016.
3. Pipit homerange moved away from dumping site, from 2015
4. Pipit clusters became less compact from 2015 and lost its significance in 2017 for radius < 5
5. Pipits stopped singing after 2015. Songs turned into Calls - a sign of distress.
6. Pipits were the only species affected (i.e. their homerange & population).
7. Control Groups (Ordinary Snape & Lesser Birchbeere) thrived and even have their homerange move closer to the dumping site
8. Pipits were spatially dependent to the Ordinary Snapes in 2015, so they should be both exposed and given the same "treatment" if dumping did actually occur, however the Ordinary Snapes were not affected.
Hypothesis: Signs of Dumping But Affected Pipits Only
Rose Pipits were still surviving as its clusters still exist and were significant, until 2017. However, they were not thriving at the dumping site and had to move away from it. Moreover, its population had fallen. This was especially since 2015 was the year that songs turned to calls, moreover, the control groups – OS and LB – did not experience a fall in population and in fact even increased in population and moved closer to the dumping site from 2015, respectively.
So, I conclude that there were signs of Dumping and this was likely to take place in 2015, but the Dumping most likely consisted of chemicals that affected mainly the Rose Pipits, and not the other species.