search for: bathy

...attempted using glm.nb from the MASS package and the negbin function from the AOD package, but both still cause some models to experience errors and warnings. Examples of my 2 largest models, run from the 2 different functions negBin.glm1 <- glm.nb(Count ~ offset(log(Tow.Area)) + Basin + Bathy + Hypoxia + Period + Depth + Basin*Depth + Bathy*Depth + Hypoxia*Depth + Period*Depth + Basin*Period + Bathy*Period + Hypoxia*Period + Hypoxia:Period:Depth + Bathy:Period:Depth + Basin:Period:Depth, control=glm.control(maxit=1000), method="glm.fit", start=coe...

...ed. negbin(formula, random, data, phi.ini = NULL, warnings = FALSE, na.action = na.omit, fixpar = list(), hessian = TRUE, control = list(maxit = 2000), ...) My largest model using glm.nb looks like this: negBin.glm1 <- glm.nb(Count ~ offset(log(Tow.Area)) + Basin + Bathy + Hypoxia + Period + Depth + Basin*Depth + Bathy*Depth + Hypoxia*Depth + Period*Depth + Basin*Period + Bathy*Period + Hypoxia*Period + Hypoxia:Period:Depth + Bathy:Period:Depth + Basin:Period:Depth, control=glm.control(maxit=1000), method="glm.fit",...

...ot; "Finetooth" "Lemon" "Scalloped" "Sandbar" "Spinner" "Abundance" "Diversity" "Depth" "Distance" "Width" "System" "Channel" "Profile" "Bathy" "Slope" "Salinity" X is longitude and Y is Latitude. I typed in the code below and received the reply "Error in eval(expr, envir, enclos) : object 'Atlantic' not found". > library(mgcv) This is mgcv 1.7-6. For overview type 'help("m...

Hi, Is there a way, preferably with R, to read shapefiles and transform them in a format that I could then use with package PBSmapping? I have been able to read such files into R with maptools' read.shape and plot it with plot.Map, but I'd like to bring the data to PBSmapping and plot from there. I also looked at the package shapefile, but it does not seem to do what I want

Hello everyone (and Hadley in particular), I often need to plot data from multiple datasets on the same graph. A common example is when mapping some values: I want to plot the underlying map and then add the points. I currently do it with base graphics, by recording the maximum region in which my map+point will fit, plotting both with these xlim and ylim parameters, adding par (new=T)

...-4 1.012384 -2.00 2.00 3 118 0 5 5 97 429970.1 5572024 1.432096 270.0000 -3 1.004987 -2.00 1.00 -1 -68 0 6 6 98 429960.1 5571904 1.012750 315.0000 0 1.001247 0.00 0.00 0 -21 0 FS_BPI_S Bathy GROUP G1 G2 G3 G4 G5 G6 G7 G8 G9 1 441 -19 8 0 0 0 0 0 0 0 1 0 2 -27 -24 9 0 0 0 0 0 0 0 0 1 3 -27 -24 8 0 0 0 0 0 0 0 1 0 4 -27 -19 8 0 0 0 0 0 0 0 1 0 5 -27 -18 9 0 0 0 0 0 0 0 0 1 6 -27 -18...