search for: cpue

...ion so I apologize if it's not like it should be and kindly ask you to not hesitate to let me know about it. Thank you for your help Xochitl C. ----------------------------------------------------------------------- Here an extract of my table and the different formula I run : > head (CPUE_table) Year Quarter Subarea Latitude Longitude Presence.S CPUE.S Presence.H CPUE.H Presence.NP CPUE.NP Presence.BW CPUE.BW Presence.C CPUE.C Presence.P CPUE.P Presence.W CPUE.W 1 2000 1 31F1 51.25 1.5 0 0 0 0 0 0 0...

...s one finds every so often in life. See the sequence below: > fish3.fis <-read.csv("emperor2.csv", check.names = TRUE, strip.white = TRUE) > colnames(fish3.fis) [1] "Month" "Year" "FishingArea" "SumOfTotalCatch" "CPUE" [6] "rCPUE" "PA" "Latitude" "Longitude" "Depth" [11] "SST" > hist(CPUE) Error in hist(CPUE) : Object "CPUE" not found So, the system knows CPUE exists, but...

Hello all! I have a data frame with nine variables and 293 cases. (attached goes the csv file). I need to calculate a index using the sum of one variable (N) divided by the length of other variable (Fisherman), but for each day (Date). I tried to use that codes: by (cpue, cpue$Date, function (x) sum (cpue$N)/length(cpue$Fisherman)) tapply (cpue$Data, cpue$Data, function (x) {sum (x)/length(y)}) In both codes, as result, I received the same value for all Date value. And, oddly, this result is the value of calculation for all data frame, as I used just: sum (cpue$...

...of my bars into grey scales (2) removing the y-axis (species name) on the right figure to avoid duplication, and (3) fix the width, such that both figures have the same apparence regardless of the legend. Thank you very much in advance! The code is as follow: (a<-ggplot(Final2, aes(x=spp, y=CPUE, fill=season)) + geom_bar(position=position_dodge()) + geom_errorbar(aes(ymin=CPUE-cpue.se, ymax=CPUE+cpue.se), size=.3, width=.2, position=position_dodge(.9)) + xlab("Species") + scale_x_discrete(breaks=c("100&qu...

...datos. Al enviarlo, se veía bien, pero está visto que por el camino la cosa se torció. Siguiendo el consejo de Javier, he subido un pequeño fragmento al DropBox. He quitado algunas especies para no hacerlo muy largo. Espero que ésta vez no haya problemas. https://www.dropbox.com/s/q7zla50oq7owg8k/CPUE.csv?dl=0 Un saludo y gracias Juan Carlos [[alternative HTML version deleted]]

I have a plotted a stacked panel graph (single x axis and multiple y axis) using the package zoo and would like to add a title for each separate panel. I am using the script: z <- with(mydata,zoo(cbind(mydata$Water.level,mydata$Submerged.plants, mydata$Crayfish.CPUE,mydata$Carp.CPUE),Year)) plot(z,type="b",pch=16,lty=2,xlab="Year",ylab=c("Metres","Realtive density", "CPUE","CPUE"),main="") Any help would be much appreciated. M -- View this message in context: http://old.nabble.com/Titles...

I'm having trouble with some catch per unit effort data (CPUE, fisheries data). Some of the samples were "retained" and some "unretained," and they are entered as 2 separate entries for the same sampling event (Date and time). I want to calculate the total CPUE (so sum the "retained" and "unretained" number for each sam...

Hi R users, I have a dataset with multiple variables and i'm trying to weigh average depths for fish species per year by their abundance (CPUE. I have tried the weighted.mean function but as i have two columns for the x value the lenghts differ with the w (CPUE). How do I solve this problem? So far I have tried this: data<-by(allspecies, list(allspecies$Depth, allspecies$Year), function(allspecies)weighted.mean(allspecies$CPUE)) but...

Hi- I'm wrestling with some of my data apparently not being called into a GLM or an LM. I'm looking at factors affecting fish annual catch rates (ie. CPUE) over 30 years. Two of the factors I'm using are sea surface temperature and sea surface temperature anomaly. A small sample of my data is below: CPUE Year Vessel_ID Base_Port Boat_Lgth Planing SST Anomaly 0.127 1976 2 BOI 40 Planing 19.4 -0.308 0.059 1976 30 BOI 40...

...ing a subset of my data. Whenever I choose a subset, in the summary the factors arent shown (as if the variable was a continuous variable). If I dont use subsets then all the factors are shown. I have copied the output from summary for both cases. Thanks for the help, Muri > model<-glm(log(cpue)~year,family=gaussian) Call: glm(formula = log(cpue) ~ year, family = gaussian) Deviance Residuals: Min 1Q Median 3Q Max -2.0962 -0.5851 -0.1241 0.4805 3.9236 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.8899 0.1844...

Hi, I need to have both italics and standard font on the y axis. This is my script: plot(x,y,pch=16,xlab="metric",ylab="species CPUE") I want species in italics and CPUE in standard text. Any tips would be appreciated. Thanks. -- View this message in context: http://old.nabble.com/Different-fonts-on-the-same-axis-tp26482460p26482460.html Sent from the R help mailing list archive at Nabble.com.

I'm having problems creating an axis label for a plot. y_label = expression(paste(plain('CPUE '),plain('(fish'),plain('x'),plain('h'^{-1}),plain(')'))) I'd like to replace the "x" with an interpunct symbol (dot). Any suggestions how to do this? Thanks in advance. [[alternative HTML version deleted]]

Dear all, I'm trying to create a list of point patterns ppp.object {spatstat} in a loop. My dataset looks like this: > names(OT1);head(OT1);dim(OT1) [1] "EID" "latitude" "longitude" "month" "year" "CPUE" "TSUM" [8] "fTSUM" EID latitude longitude month year CPUE TSUM fTSUM 1 167-1-1996-1135 67.70000 -61.81667 9 1996 0 0 F 2 167-10-1996-1135 67.71667 -59.18333 9 1996 0 0 F 3 167-100-1996-1135 67...

Hi -I'm running a GAM with 7 explanatory variables with a Poisson error structure. All of the variables are continuous so I'm getting error messages in R. cod.fall.full.gam.model<-gam(Kept.CPUE~s(HOUR)+s(LAT_dec)+s(LONG_dec)+s(meantemp_C)+s(meandepth_fa)+s(change_depth)+s(seds), data=cod.fall.version2,family=poisson) In dpois(y, mu, log = TRUE) ... : non-integer x = 5.325517 Thanks Sally -- View this message in context: http://www.nabble.com/GAM-and-Poisson-distribution-tp20508604p2050...

Hi, we have this model glm(CPUE~Year+Month+Length+Year:Month, family=Gamma (link="inverse")) where, CPUE is a numerical variable >0, We use this model for standardization of CPUE and we need the new estimates for each year. How can we obtain it? Thank you

...ed negative binomial (ZINB) distribution. I am using the same variables for both parts of the model, can anyone help me with R code to compute overall marginal effects of each variable? My model is specified as follows: ZINB <- zeroinfl(Tot.Anglers ~ Location + Season + Daytype + Holiday.not + CPUE + ShoreType + Access + Source.pop + WindSpeed + offset(beat_length), dist="negbin", data=anglers) Many thanks, Jeremy [[alternative HTML version deleted]]

...ng and thus I don't know if anyone has answer my query. Here is the original message: I'm applying a gam model (package mgcv) to predict relative abundances of a fish species. The covariates are year, month, vessel and statistical rectangle. The model looks like this: g1 <- gam(log(cpue) ~ s(rekt1) + s(year) + s(mon) + s(reg1), data = dataTest) Once the model is fitted to the data I want to get the mean model estimates by year. I do the following: obsPred <- data.frame(year = dataTest$year, pred = predict(g1, type = "response")) gamFit <- tapply(obsPred$pred,...

Dear list fellows, I want to study time series and use R to analyse time series of fishing data from several species (landings and cpue) investigating the correlation between them and with environmental factors (water temperature, wind, etc.). Searching at Amazon I found three books with examples in R: Time Series Analysis: With Applications in R by Jonathan D. Cryer and Jonathan D. Cryer Time Series Analysis and Its Applicat...

Dear all, I have a serious problem to solve my model. I study over exploitation of fish in the bay of biscay (france). I know only the level of catch and the fishing effort (see data below) by year. My model is composed by the following equations: * the growth function Gt(St) = r*St*(1-St/sbar) with Gt the growth of each period t r intrinsec growth of the stock sbar carriyng capacity of the

...<- delta0 + alpha1 * Month[i] + alpha2 * Lat[i] + beta1 * (SeaW[i] - nSeaW) logit(p[i]) <- gamma0 mu[i, 1] <- 0 mu[i, 2] <- lambda[i] mu.i[i] <- mu[i, index[i]] index[i] ~ dcat(theta[i, 1:2]) theta[i, 1] <- p[i] theta[i, 2] <- 1 - p[i] # mixture YFTCPUE[i] ~ dpois(mu.i[i]) } # recalculate the original intercept term Intercept <- delta0 - beta1 * nSeaW #prior on regression coefficients beta1 ~ dnorm(0,1.0E-6) alpha2 ~ dnorm(0,1.0E-6) alpha1 ~ dnorm(0, 1.0E-6) delta0 ~ dnorm(0, 1.0E-6) gamma0 ~ dnorm(0, 1.0E-6)...