search for: vessel

Dear R users, I am querying an Oracle database using sqlQuery() from the RODBC package. What I would like to do is embed my R Sql query in a for loop. For example suppose the data.frame "vessel" contains two columns: vessel[,1]= vessel id and vessel[,2]=permit year I am using "vessel" as an input in a SQL based function that require inputs of permit year and vessel id as described below: query<-"Select * From vessel Where common.utility(vessel ID,vessel year)&q...

Hello, I want to select same data of a data frame that has information of fish catches: catch, vessel, trip, day, month, year where trip is the fishing operation number done in each day by each vessel, i.e. for each day, any vesel can have several register of catch corresponding to each operation. >From this data frame I want to do: 1. Add the catches of all trip by vessel and day; then 2. Sele...

Hi professors I'm finding the best package in R for manage a big database about vessels. I need to divided the travel of each boat into tracks. I have GPS points. I would join these points in a line on the map. I make that with the function "lines", but I can't divided by vessel or by track, I saw several packages, for tracking animals but they had specifically variab...

...trying to create a plot showing the density distribution of some shipping data. I like the look of violin plots, but my data is not continuous but rather binned and I want to make sure its binned nature (not smooth) is apparent in the final plot. So for example, I have the number of individuals per vessel, but rather than having the actual number of individuals I have data in the format of: 7 values of zero, 11 values between 1-10, 6 values between 10-100, 13 values between 100-1000, etc. To plot this data I generated a new dataset with the first 7 values being 0, representing the 7 values of 0, the...

Hi -I am using R v 2.13.0. I am trying to use the aggregate function to calculate the percent at length for each Trip_id and CommonName. Here is a small subset of the data. Trip_id Vessel CommonName Length Count 1 230 Sunlight Shad,American 19 1 2 230 Sunlight Shad,American 20 1 3 230 Sunlight Shad,American 21 1 4 230 Sunlight Shad,American 23 1 5 230 Sunlight Shad,Americ...

..., 2004 12:47 PM To: r-help at stat.math.ethz.ch Subject: [R] mixed models for analyzing survey data with unequal selection probability Hi, I need the help on this topic because this is out of my statistical trianing as biologist. Here is my brief description of the problem. I have a survey that VESSELs are selected at random with the probability of p(j). Then the tows within the jth VESSEL are sampled at random with probability of p(i|j). I write my model as y = XB + Zb + e where XB is fixed part, Zb is for random effect (VESSEL) and e is within-vessel error. I feel that I should weight the...

...ning messages: 1: glm.fit: algorithm did not converge 2: glm.fit: fitted probabilities numerically 0 or 1 occurred What I am trying to do is model my response variable (CID: correct bird identification) as a function of the predictor variables weather state (WS), sea state (SS), distance from the vessel (DV) and duration of the sighting (DS). I defined both sea state and weather state as factors with three levels (0, 1, or 2). Distance of the vessel values are 100, 80, 60, 40, and 20. Duration of the sighting ranges from 0 to 58 seconds. The output R is giving me is: Deviance Residuals: M...

...nal Medicine University of Virginia, Charlottesville, Virginia ------------------ ######################################################## # Compute the Hosmer-Lemeshow 'goodness-of-fit' test cd.full_model = glm(formula = Collaterals ~ CHF + Age + CABG + relevel (as.factor (num.obst.vessels),"one") + Current.smoker + DM + HTN + ace.inhibitor + MI, family = binomial(link = "logit")) hosmerlem = function(y, yhat, g=10) { cutyhat = cut(yhat, breaks = quantile(yhat, probs=seq(0, 1, 1/g)), include.lowest=TRUE) obs = xtabs(cbind(1 - y, y) ~ cutyhat)...

I am attempting to fit discrete data (daily counts of arrivals of recreational vessels at locks on a river) using the fit.dist package. Some distributions return values of NaN and Inf for certain situations, an example with Inf values is shown below. # of vessels: 1 2 3 4 5 6 7 8 9 10 11 # of days with # of vessels:...

...ed models for analyzing survey data with unequal selection >>probability >> >>Hi, >> >>I need the help on this topic because this is out of my statistical >>trianing as biologist. Here is my brief description of the problem. >> >>I have a survey that VESSELs are selected at random with the probability >>of p(j). Then the tows within the jth VESSEL are sampled at random with >>probability of p(i|j). I write my model as >> >>y = XB + Zb + e >>where XB is fixed part, Zb is for random effect (VESSEL) and e is >>within...

Hello everyone, I was trying to fit an arima model to a daily historical data, but, for some reason, havent been able to. I basically have 212 observations (from 12/1/2012 to 06/30/2013) containing the number of transits for a particular vessel. The following messages are produced by R: dailytrans.fit<-arima(dailytrans$transits, order=c(0,1,2), seasonal=list(order=c(0,1,2), period=365), include.mean=FALSE) Error in arima(dailytrans$transits, order = c(0, 1, 2), seasonal = list(order = c(0, : too few non-missing observations > d...

...2621 hours, 59 minutes ago". I don't know how to convert these data to the time formats used in R (POSIXct). The script is very simple and worked before: library(XML) x <- readHTMLTable('url') where the 'url' is the link to the website with the specification of the vessel. I appreciate any help. Cristina -- Cristina Silva Divis?o de Modela??o e Gest?o de Recursos Pesqueiros Av. Dr. Alfredo Magalh?es Ramalho 1495-165 Lisboa @: csilva at ipma.pt <mailto:csilva at ipma.pt> #: +351 213027096 <phoneto:+351213027096> [[alternative HTML version deleted]]

...odel1,data.frame(year=as.factor(xv),nafdiv=as.factor(rep(" 3N",length(xv))), duration=1000, cfv=as.factor(rep(106166,length(xv))), hooks=rep(1,length(xv))),type="response",se.fit=T) where model1 standardises catch rate and includes nafo area, year, duration of fishing sets, and vessel (cfv) as explanatory variables with an offset of hooks. To predict the years we had data for, we substituted a vector "xv" of years into predict and specified the values to be used for the other explanatory variables. This seemed to work well. When I switched to a glmmPQL model, with ve...

...ship 1 and ship 2? I would think that animals would react similarly to both ships (and the effect would just add up), so I don?t want the model to estimate two covariates that I think are the same. And if I had 5 ships instead of 2, those dfs would really rack up. Note that I can?t just double the vessel values, since their speeds, directions, lengths, etc all differ. Here's a little mock data set for 3 surveys - 2 have a single ship, and 1 has two ships. Note that each survey is only done once, so if there are 2 ships (or more), the number of animals (and all other survey-related info) is jus...

...able to determine the level of recovery that each site has experienced. Recovery is our categorical dependent variable with values of none, low, medium, high which are based upon percent seagrass regrowth into the injury over time. I wish to be able to predict the level of recovery of future vessel grounding sites based upon a number of categorical / continuous predictor variables used here including (but not limited to) such parameters as: sediment grain size, wave exposure, original size (volume) of the injury, injury age, injury location. When I run /rpart/, the data is split into onl...

...t all my mails have been sent to the recycle bin without possibility of recovering 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...

Dear R-helper, I wonder if anyone can help me. I am trying to convert a dataset to a format recognizable by a software onboard a research vessel but I am having problems with some steps. I have a data frame as follows: Conc Lat Lon Depth Point0000 56.25 -5.65 70 Point0001 56.55 -5.35 85 Point0002 56.25 -5.65 65 Point0003 56.37 -5.21 80 Point0004 56.45 -5.23...

...3.0.1) but after loading the package and calling: library(lme4) help(glmer.nb) # I receive the following message: No documentation for ‘glmer.nb’ in specified packages and libraries: you could try ‘??glmer.nb’ # I called the glmer() function successfully, i.e. glmer(CNos ~ Year+Month+(1|Vessel), data=Occ.cr, family=poisson, offset=Occ.cr$logtrapd) # but when I tried this: glmer.nb(CNos ~ Year+Month+(1|Vessel), data=Occ.cr, offset=Occ.cr$logtrapd) # I got the error message: Error: could not find function "glmer.nb" # Can you please provide advice on how I might us...

...ned another. > But to the thoughts train, every user is just one user. > Mapping user to other user is creating a mess you don't want to solve > yourself. > Maybe so. I was merely trying to express a different view. Where Samba is somewhat selfcontained and uses the Linuxhost as a vessel for its purpose; file sharing for Windows. With that in mind, Rowland is right when he sais that the domain adminstrator account becomes an ordinary Unix user on the Linux host. For Samba its good enough. Regards Davor > > -- > With best regards, > Andrey Repin > Sunday, April 19,...

Dear Jeff, I am sorry, I didn't notice that it was not plain text. I hope that it is now in the correct format. I explain the problem again, now with more detais. I am collecting the track positions of our research vessel from www.marinetraffic.com. In the page, the data appear in a table: Timestamp Source Speed (kn) Latitude (?) Longitude (?) Course (?) Show on Map 2017-07-11 10:57 (UTC) Terr-AIS 8.6 37.36228 -9.176811 200 2017-07-11 10:54 (UTC) Terr-AIS 8.4 37.36875 -...