similar to: mixed models for analyzing survey data with unequal selec tion probability

Hi, All Thanks to Robert Baskin, Thomas Lumley, and Spencer Graves for the valuable helps. I have learned a lot from this discussion. I put all discussions together without editing, so we can see how things are evolved. Likely, I have a lot of articles to read. As in the discussion, mixed modeling approach is a poosible but may be over-kill in my posted data analyses. I will explore other

Dear All, I have a complex survey data with observations having differing probabilities of selection into the sample. I would like to run a linear mixed effects model and also to use weighting that takes this into consideration. As far as I know, however, the weighting option for the lme command from nlme package (or glmmPQL from MASS, for that matter) is related to heteroscedasticity. Is there a

I am trying to use the LME package to run a multilevel logistic model using the following code: ------------------------------------------------------------------------ ------------------------------------------- Model1 = GLMM(WEAP ~ TSRAT2 , random = ~1 | GROUP , family = binomial, na.action = na.omit ) ------------------------------------------------------------------------

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

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

Hi, I'm 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

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 variables. And I

Hello, When I run the following glm model: modelresult=glm(CID~WS+SS+DV+DS, data=kimu, family=binomial) I get the following warning 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

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

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),

Is there an optimal / minimum sample size for attempting to construct a classification tree using /rpart/? I have 27 seagrass disturbance sites (boat groundings) that have been monitored for a number of years. The monitoring protocol for each site is identical. From the monitoring data, I am able to determine the level of recovery that each site has experienced. Recovery is our

R-help, Sorry for posting again the same question (dated 26-03-2007) but 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

Hi, I am extracting positions data from the marine traffic website. The table has a "Timestamp" column which, in the browser, appears with the format yyyy-mm-dd HH:MM (UTC), e.g. 2017-07-10 14:04 (UTC). When I import the table, the same date "2017-07-10 14:04 (UTC)" appears as "1499696500149969650021 minutes ago", This is the more recent date and time. Older

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

This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X, say, how do the two test specifications test.terms = c("ZB:X","ZC:X") # and test.terms = ~ ZB:X + ZC:X in

Hi all, I have been previously been using scatterplot3d package to create some graphs but unfortunately it does not allow me to rotate the plot on all three axis. The cloud() function in the lattice package does allow me to do so. When I was using scatterplot3d I was using a script (Shown Below) to calculate the mean, quartiles and range limits for all three axis and I was representing that on the

Hi Andrey, 2015-04-19 0:12 GMT+02:00 Andrey Repin <anrdaemon at yandex.ru>: > Greetings, Davor Vusir! > >>> Hi, there are two separate points of view here, map 'Administrator' to the >>> 'root' user, or give 'Administrator' a uidNumber. If you do the first then >>> 'Administrator' can change directory settings on a Unix

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

Not reproducible. [1][2][3] If our answers don't seem to apply to your situation, it will likely be because you did not explain your question clearly. Not plain text. This is a plain text mailing list, and the best-case scenario when you let your email program send HTML is that what you saw is not what we see (worst case is your email is scrambled on our end). Have you read the

Dear David, Thank you very much. You have solved my problem, not of importing in the right format but the date conversion. It works for me, and if I specify the timezone tz='GMT', I get exactly the same time of the webpage. Sincerely yours, Cristina On Tue, 11 Jul 2017 07:15:43 -0700, David Winsemius wrote: >> On Jul 11, 2017, at 6:25 AM, Cristina Silva <csilva at