Simple data analysis for biologists

This document provides a simple introduction to research methods and analysis tools for biologists or environmental scientists, with particular emphasis on fish biology in devleoping countries.

An evaluation of length-frequency sampling procedures and subsequent data analysis for purse seine-caught yellowfin tuna in the Eastern Pacific Ocean

ENGLISH: A two-stage sampling design is used to estimate the variances of the numbers of yellowfin in different age groups caught in the eastern Pacific Ocean. For purse seiners, the primary sampling unit (n) is a brine well containing fish from a month-area stratum; the number of fish lengths (m) measured from each well are the secondary units. The fish cannot be selected at random from the wells because of practical limitations. The effects of different sampling methods and other factors on the reliability and precision of statistics derived from the length-frequency data were therefore examined. Modifications are recommended where necessary. Lengths of fish measured during the unloading of six test wells revealed two forms of inherent size stratification: 1) short-term disruptions of existing pattern of sizes, and 2) transition zones between long-term trends in sizes. To some degree, all wells exhibited cyclic changes in mean size and variance during unloading. In half of the wells, it was observed that size selection by the unloaders induced a change in mean size. As a result of stratification, the sequence of sizes removed from all wells was non-random, regardless of whether a well contained fish from a single set or from more than one set. The number of modal sizes in a well was not related to the number of sets. In an additional well composed of fish from several sets, an experiment on vertical mixing indicated that a representative sample of the contents may be restricted to the bottom half of the well. The contents of the test wells were used to generate 25 simulated wells and to compare the results of three sampling methods applied to them. The methods were: (1) random sampling (also used as a standard), (2) protracted sampling, in which the selection process was extended over a large portion of a well, and (3) measuring fish consecutively during removal from the well. Repeated sampling by each method and different combinations indicated that, because the principal source of size variation occurred among primary units, increasing n was the most effective way to reduce the variance estimates of both the age-group sizes and the total number of fish in the landings. Protracted sampling largely circumvented the effects of size stratification, and its performance was essentially comparable to that of random sampling. Sampling by this method is recommended. Consecutive-fish sampling produced more biased estimates with greater variances. Analysis of the 1988 length-frequency samples indicated that, for age groups that appear most frequently in the catch, a minimum sampling frequency of one primary unit in six for each month-area stratum would reduce the coefficients of variation (CV) of their size estimates to approximately 10 percent or less. Additional stratification of samples by set type, rather than month-area alone, further reduced the CV's of scarce age groups, such as the recruits, and potentially improved their accuracy. The CV's of recruitment estimates for completely-fished cohorts during the 198184 period were in the vicinity of 3 to 8 percent. Recruitment estimates and their variances were also relatively insensitive to changes in the individual quarterly catches and variances, respectively, of which they were composed. SPANISH: Se usa un diseño de muestreo de dos etapas para estimar las varianzas de los números de aletas amari11as en distintos grupos de edad capturados en el Océano Pacifico oriental. Para barcos cerqueros, la unidad primaria de muestreo (n) es una bodega de salmuera que contenía peces de un estrato de mes-área; el numero de ta11as de peces (m) medidas de cada bodega es la unidad secundaria. Limitaciones de carácter practico impiden la selección aleatoria de peces de las bodegas. Por 10 tanto, fueron examinados los efectos de distintos métodos de muestreo y otros factores sobre la confiabilidad y precisión de las estadísticas derivadas de los datos de frecuencia de ta11a. Se recomiendan modificaciones donde sean necesarias. Las ta11as de peces medidas durante la descarga de seis bodegas de prueba revelaron dos formas de estratificación inherente por ta11a: 1) perturbaciones a corto plazo en la pauta de ta11as existente, y 2) zonas de transición entre las tendencias a largo plazo en las ta11as. En cierto grado, todas las bodegas mostraron cambios cíclicos en ta11a media y varianza durante la descarga. En la mitad de las bodegas, se observo que selección por ta11a por los descargadores indujo un cambio en la ta11a media. Como resultado de la estratificación, la secuencia de ta11as sacadas de todas las bodegas no fue aleatoria, sin considerar si una bodega contenía peces de un solo lance 0 de mas de uno. El numero de ta11as modales en una bodega no estaba relacionado al numero de lances. En una bodega adicional compuesta de peces de varios lances, un experimento de mezcla vertical indico que una muestra representativa del contenido podría estar limitada a la mitad inferior de la bodega. Se uso el contenido de las bodegas de prueba para generar 25 bodegas simuladas y comparar los resultados de tres métodos de muestreo aplicados a estas. Los métodos fueron: (1) muestreo aleatorio (usado también como norma), (2) muestreo extendido, en el cual el proceso de selección fue extendido sobre una porción grande de una bodega, y (3) medición consecutiva de peces durante la descarga de la bodega. EI muestreo repetido con cada método y distintas combinaciones de n y m indico que, puesto que la fuente principal de variación de ta11a ocurría entre las unidades primarias, aumentar n fue la manera mas eficaz de reducir las estimaciones de la varianza de las ta11as de los grupos de edad y el numero total de peces en los desembarcos. El muestreo extendido evito mayormente los efectos de la estratificación por ta11a, y su desempeño fue esencialmente comparable a aquel del muestreo aleatorio. Se recomienda muestrear con este método. El muestreo de peces consecutivos produjo estimaciones mas sesgadas con mayores varianzas. Un análisis de las muestras de frecuencia de ta11a de 1988 indico que, para los grupos de edad que aparecen con mayor frecuencia en la captura, una frecuencia de muestreo minima de una unidad primaria de cada seis para cada estrato de mes-área reduciría los coeficientes de variación (CV) de las estimaciones de ta11a correspondientes a aproximadamente 10% 0 menos. Una estratificación adicional de las muestras por tipo de lance, y no solamente mes-área, redujo aun mas los CV de los grupos de edad escasos, tales como los reclutas, y mejoró potencialmente su precisión. Los CV de las estimaciones del reclutamiento para las cohortes completamente pescadas durante 1981-1984 fueron alrededor de 3-8%. Las estimaciones del reclutamiento y sus varianzas fueron también relativamente insensibles a cambios en las capturas de trimestres individuales y las varianzas, respectivamente, de las cuales fueron derivadas. (PDF contains 70 pages)

Workshop on turtle bycatch mitigation for longline fisheries: experimental design and data analysis, 7-8 November 2007, San Ramón de Alajuela, Costa Rica

Large numbers of fishing vessels operating from ports in Latin America participate in surface longline fisheries in the eastern Pacific Ocean (EPO), and several species of sea turtles inhabit the grounds where these fleets operate. The endangered status of several sea turtle species, and the success of circle hooks (‘treatment’ hooks) in reducing turtle hookings in other ocean areas, as compared to J-hooks and Japanese-style tuna hooks (‘control’ hooks), prompted the initiation of a hook exchange program on the west coast of Latin America, the Eastern Pacific Regional Sea Turtle Program (EPRSTP)1. One of the goals of the EPRSTP is to determine if circle hooks would be effective at reducing turtle bycatch in artisanal fisheries of the EPO without significantly reducing the catch of marketable fish species. Participating fishers were provided with circle hooks at no cost and asked to replace the J/Japanese-style tuna hooks on their longlines with circle hooks in an alternating manner. Data collected by the EPRSTP show differences in longline gear and operational characteristics within and among countries. These aspects of the data, in addition to difficulties encountered with implementation of the alternating-hook design, pose challenges for analysis of these data.

Report on BOBLME Indian mackerel genetic data analysis including final workshop, Phuket, Thailand, 17-19 February 2015

This workshop followed on from two previous workshops held in Colombo, Sri Lanka, 2012 and Kochi, India in 2013. The 14 microsattellite markers had previously been developed for Indian Mackerel (Rastrelliger kanagurta) were used on 31 tissue collections from all eight countries were genotyped in India.

Techniques for spatial analysis and visualization of benthic mapping data: final report

The mapping and geospatial analysis of benthic environments are multidisciplinary tasks that have become more accessible in recent years because of advances in technology and cost reductions in survey systems. The complex relationships that exist among physical, biological, and chemical seafloor components require advanced, integrated analysis techniques to enable scientists and others to visualize patterns and, in so doing, allow inferences to be made about benthic processes. Effective mapping, analysis, and visualization of marine habitats are particularly important because the subtidal seafloor environment is not readily viewed directly by eye. Research in benthic environments relies heavily, therefore, on remote sensing techniques to collect effective data. Because many benthic scientists are not mapping professionals, they may not adequately consider the links between data collection, data analysis, and data visualization. Projects often start with clear goals, but may be hampered by the technical details and skills required for maintaining data quality through the entire process from collection through analysis and presentation. The lack of technical understanding of the entire data handling process can represent a significant impediment to success. While many benthic mapping efforts have detailed their methodology as it relates to the overall scientific goals of a project, only a few published papers and reports focus on the analysis and visualization components (Paton et al. 1997, Weihe et al. 1999, Basu and Saxena 1999, Bruce et al. 1997). In particular, the benthic mapping literature often briefly describes data collection and analysis methods, but fails to provide sufficiently detailed explanation of particular analysis techniques or display methodologies so that others can employ them. In general, such techniques are in large part guided by the data acquisition methods, which can include both aerial and water-based remote sensing methods to map the seafloor without physical disturbance, as well as physical sampling methodologies (e.g., grab or core sampling). The terms benthic mapping and benthic habitat mapping are often used synonymously to describe seafloor mapping conducted for the purpose of benthic habitat identification. There is a subtle yet important difference, however, between general benthic mapping and benthic habitat mapping. The distinction is important because it dictates the sequential analysis and visualization techniques that are employed following data collection. In this paper general seafloor mapping for identification of regional geologic features and morphology is defined as benthic mapping. Benthic habitat mapping incorporates the regional scale geologic information but also includes higher resolution surveys and analysis of biological communities to identify the biological habitats. In addition, this paper adopts the definition of habitats established by Kostylev et al. (2001) as a “spatially defined area where the physical, chemical, and biological environment is distinctly different from the surrounding environment.” (PDF contains 31 pages)

Growth of skipjack, Katsuwonus pelamis, and yellowfin, Thunnus albacares, tunas in the Eastern Pacific Ocean, as estimated from tagging data

ENGLISH: The average linear growth rate of skipjack in the eastern Pacific is less than 1 mm per day except for fish 375 to 424 mm in length at release. The growth rate shows a decrease with increasing length and increasing time at liberty. The growth rate of fish in the length range of about 43 to 57 cm is apparently more rapid in the eastern Pacific than in the western Pacific. Dsing data for the northeastern and southeastern Pacific combined, K and ~ were estimated to be 0.658 (on an annual basis) and 885 mm, respectively, by the ungrouped method and 0.829 and 846 mm, respectively, by the grouped method. Sensitivity analyses have shown however, that the estimates of these parameters are poorly determined by the sum of squares method used to derive them. Estimates of K and ~ for the eastern Pacific tend to be lower and higher, respectively, than those for the western Pacific. The average linear growth rate of yellowfin in the eastern Pacific is a little less than 1 mm per day for fish between about 25 and 100 cm in length at release. The growth appears to be most rapid in Area 2 (Revillagigedo Islands) and slowest in Areas 1 (Baja California), 5 (Central America- Colombia), and 6 (Ecuador-Peru). There is considerable variation in the growth rates of individual fish. The growth does not show a decrease with increasing length or increasing time at liberty so realistic estimates of the parameters of the von Bertalanffy or other similar equations cannot be calculated from these data. If realistic estimates of these parameters are to be secured larger fish must be tagged and released or many more long-term returns from fish to about 100 cm in length at release must be obtained. The growth patterns for the eastern Pacific, central Pacific and eastern Atlantic found by most other investigators differ from one another and from those found in the present study. Some of these differences may be real and others may be due to deficiencies in the data or the methods of analysis. Estimates obtained from tagging data are believed to be realistic provided the tags do not inhibit the growth of the fish. It appears that the growth rates of single- and double-tagged fish are the same; this indicates, though not unequivocally, that the tags do not inhibit the growth. SPANISH: La tasa media de crecimiento lineal del barrilete en el Pacífico oriental es inferior a lmm/día, excepto en el caso de peces de entre 375y 424mm de longitud de liberación. La tasa de crecimiento disminuye a medida que aumenta la longitud y el tiempo en libertad. La tasa de crecimiento de peces de entre unos 43 y 57 cm de longitud parece ser mayor en el Pacífico oriental que en el occidental. A partir de datos del Pacífico nororiental y suroriental combinados, se estimaron K y loo en 0.658 (anual) y 885mm, respectivamente, usando el método no agrupado, y 0.829 y 846mm, respectivamente, usando el método agrupado. Sin embargo, los análisis de sensitividad han demostrado que el método de suma de cuadrados utilizado para derivar las estimaciones de estos parámetros las determina con poca precisión. Las estimaciones de K y loo para el Pacífico oriental suelen ser inferiores y superiores, respectivamente, a los del Pacífico occidental. La tasa media de crecimiento lineal del aleta amarilla en el Pacífico oriental es ligeramente inferior a lmm/día para los peces de entre unos 25y 100cmde longitud de liberación. El crecimiento parece ser más rápido en el Area 2(Islas Revillagigedo),y más lento en las Areas 1(Baja California), 5 (Centroamérica-Colombia), y 6 (Ecuador-Perú). Las tasas de crecimiento de peces individuales varían considerablemente. El crecimiento no muestra una disminuciónconun aumento en la longitud o en el tiempo en libertad, y por consecuencia no se se pueden calcular estimaciones realistas de los parámetros de la ecuación de von Bertalanffy u otras ecuaciones similares a partir de estos datos. Para obtener estimaciones realistas de estos parámetros sería necesario marcar peces mayores u obtener muchas más devoluciones a largo plazo de marcas de peces de unos 100cm de longitud de liberación. Los patrones de crecimiento correspondientes al Pacífico oriental, Pacífico central, y Atlántico oriental descubiertos por la mayoría de los investigadores son diferentes entre síy también de los del presente estudio. Es posibleque algunas de estas diferencias sean verdaderas, mientras que otras se deban a faltas en los datos on en los métodos analíticos utilizados. Se considera que las estimaciones obtenidas a partir de los datos de marcado son realistas, suponiendo siempre que las marcas no impidan el crecimiento de los peces. Parece ser que las tasas de crecimiento de peces con una marca y con dos son idénticas, lo cual indica, aunque sin certeza total, que las marcas no ejercen tal efecto. (PDF contains 76 pages.)

An ecological characterization of the marine resources of Vieques, Puerto Rico. Part I: Historical data synthesis

From the 1940s until 2003, portions of the island of Vieques, a municipality within the Commonwealth of Puerto Rico, were used by the US Navy as a base and training facility, resulting in development and zoning history that differ in comparison to other Caribbean islands. The majority of former Navy lands are now under the jurisdiction of the Department of the Interior’s Fish and Wildlife Service as a National Wildlife Refuge, while a smaller percentage of land was transferred to the Vieques municipality and the Puerto Rico Conservation Trust. An analysis of the distribution and status of the marine resources is timely in light of the recent land transfer, increases in development and tourism, and potential changes in marine zoning around the island. To meet this need, NOAA’s Biogeography Branch, in cooperation with the Office of Response and Restoration and other local and regional partners, conducted Part I of an ecological characterization to integrate historical data and research into a synthesis report. The overall objective of this report is to provide resource managers and residents a comprehensive characterization of the marine resources of Vieques to support research, monitoring, and management. For example, knowledge of the spatial distribution of physical features, habitats, and biological communities is necessary to make an informed decision of the establishment and placement of a marine protected area (MPA). The report is divided into chapters based on the physical environment (e.g., climate, geology, bathymetry), habitat types (e.g., reefs and hardbottom, seagrasses, mangroves) and major faunal groups (e.g. fish, turtles, birds). Each section includes five subsections: an overview, description of the relevant literature, methods of analysis, information on the distribution, status and trends of the particular resource, and a discussion of ecological linkages with other components of the Vieques marine ecosystem and surrounding environment. The physical environment of Vieques is similar to other islands within the Greater Antilles chain, with some distinctions. The warm, tropical climate of Vieques, mediated by the northeasterly trade winds, is characterized by a dry season (December-April) and a rainy season (May-November), the latter of which is characterized by the occasional passage of tropical cyclones. Compared to mainland Puerto Rico, Vieques is characterized by lower elevation, less annual precipitation, and higher average temperatures. The amount of annual precipitation also varies spatially within Vieques, with the western portion of the island receiving higher amounts of rainfall than further east. While the North Equatorial Current dominates the circulation pattern in the Greater Antilles region, small scale current patterns specific to Vieques are not as well characterized. These physical processes are important factors mitigating the distribution and composition of marine benthic habitats around Vieques. In general, the topography of Vieques is characterized by rolling hills. Mt. Pirata, the tallest point at 301 m, is located near the southwest coast. In the absence of island wide sedimentation measurements, information on land cover, slope, precipitation, and soil type were used to estimate relative erosion potential and sediment delivery for each watershed. While slope and precipitation amount are the primary driving factors controlling runoff, land use practices such as urban development, military activity, road construction, and agriculture can increase the delivery of pollution and sediments to coastal waters. Due to the recent land transfer, increased development and tourism is expected, which may result in changes in the input of sediments to the coastal environment.

Herpetofaunal Inventories of the National Parks of South Florida and the Caribbean: Volume III. Big Cypress National Preserve

Amphibian declines and extinctions have been documented around the world, often in protected natural areas. Concern for this trend has prompted the U.S. Geological Survey and the National Park Service to document all species of amphibians that occur within U.S. National Parks and to search for any signs that amphibians may be declining. This study, an inventory of amphibian species in Big Cypress National Preserve, was conducted from 2002 to 2003. The goals of the project were to create a georeferenced inventory of amphibian species, use new analytical techniques to estimate proportion of sites occupied by each species, look for any signs of amphibian decline (missing species, disease, die-offs, and so forth.), and to establish a protocol that could be used for future monitoring efforts. Several sampling methods were used to accomplish these goals. Visual encounter surveys and anuran vocalization surveys were conducted in all habitats throughout the park to estimate the proportion of sites or proportion of area occupied (PAO) by each amphibian species in each habitat. Opportunistic collections, as well as limited drift fence data, were used to augment the visual encounter methods for highly aquatic or cryptic species. A total of 545 visits to 104 sites were conducted for standard sampling alone, and 2,358 individual amphibians and 374 reptiles were encountered. Data analysis was conducted in program PRESENCE to provide PAO estimates for each of the anuran species. All of the amphibian species historically found in Big Cypress National Preserve were detected during this project. At least one individual of each of the four salamander species was captured during sampling. Each of the anuran species in the preserve was adequately sampled using standard herpetological sampling methods, and PAO estimates were produced for each species of anuran by habitat. This information serves as an indicator of habitat associations of the species and relative abundance of sites occupied, but it will also be useful as a comparative baseline for future monitoring efforts. In addition to sampling for amphibians, all encounters with reptiles were documented. The sampling methods used for detecting amphibians are also appropriate for many reptile species. These reptile locations are included in this report, but the number of reptile observations was not sufficient to estimate PAO for reptile species. We encountered 35 of the 46 species of reptiles believed to be present in Big Cypress National Preserve during this study, and evidence exists of the presence of four other reptile species in the Preserve. This study found no evidence of amphibian decline in Big Cypress National Preserve. Although no evidence of decline was observed, several threats to amphibians were identified. Introduced species, especially the Cuban treefrog (Osteopilus septentrionalis), are predators and competitors with several native frog species. The recreational use of off-road vehicles has the potential to affect some amphibian populations, and a study on those potential impacts is currently underway. Also, interference by humans with the natural hydrologic cycle of south Florida has the potential to alter the amphibian community. Continued monitoring of the amphibian species in Big Cypress National Preserve is recommended. The methods used in this study were adequate to produce reliable estimates of the proportion of sites occupied by most anuran species, and are a cost-effective means of determining the status of their populations.

Analysis and adsorption of phosalone on aquatic sediment

Phosalone is a non systematic, wide spectrum organophosphate pesticide which was discovered in 1961 in the laboratories of the Societe des Usines Chimique Rhone-Poulenc in France. It has been approved for commercial use since 1964 in France, in Australia since 1966, in the United Kingdom in 1967 and in many other countries including Japan, Egypt, USSR and the USA. This study provides a full literature review on all aspects of phosalone including its physical, biological and chemical characteristics, and analytical methods of analysis with particular reference to soils/sediments. Furthermore, it aims to develop a method for the determintion of phosalone in aquatic sediments and to determine the adsorption of phosalone onto kaolinite.

Strain comparisons in aquaculture species: a manual

When different strains or breeds of a particular species are available, the best choice is seldom immediately obvious for producers. Scientists are also interested in the relative performance of different strains because it provides a basis for recommendations to producers and it often stimulates the conduct of work aimed at unraveling the underlying biological mechanisms involved in the expression of such differences. Hence, strain or breed comparisons of some sort are frequently conducted. This manual is designed to provide general guidelines for the design of strain comparison trials in aquaculture species. Example analyzes are provided using SAS and SPSS. The manual is intended to serve a wide range of readers from developing countries with limited access to information. The users, however, are expected to have a basic knowledge of quantitative genetics and experience in statistical methods and data analysis as well as familiarity with computer software. The manual mainly focuses on the practical aspects of design and data analysis, and interpretation of results.

Approaches to the evaluation and management of the fish stock in Kainji Lake, Nigeria

Estimates of potential yield for Kainji Lake, and the methods of analysis by earlier workers are discussed. Also summarized is the state of the fishery after impoundment, between 1969 and 1971, based on experimental gillnet catches. Recent sampling of the young of the year along the littoral margin indicates that most of the commercially important species have spawned successful1y in the lake. An intense fishing mortality of juvenile fish, owing to the use of small mesh nets by local fishermen, presents a possible threat to the future establishment of the fish in the lake. The results of gill-net selection studies based on HOLT'S (1957) method are given. The data have been extracted from experimental gill-net catches with graded fleets of nets between 1969 and 1971. Recommendations based on the above studies have been made to ensure a successful establishment of the fish species in the lake and an increase in catch-per~unit effort in subsequent years.