similar to: Adding a subset to a glm messes up factors?

Hi, Is there a way to make the filled.contour() function produce a plot similar to the image() function? i.e. not have smooth contours but rather distinct "squares" corresponding to the different cells of the matrix. The reason I ask is because since image doesn't have an option to display colour bars it is a lot less useful. Cheers Muri Soares Department of Maths and Applied

Dear all, I am encountering some issues with my data and need some help. I am trying to run glm analysis with a presence/absence variable as response variable and several explanatory variable (time, location, presence/absence data, abundance data). First I tried to use the glm() function, however I was having 2 warnings concerning glm.fit () : # 1: glm.fit: algorithm did not converge # 2:

Dear Sir/Madam, I am encountering one of those alien computer momements 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]

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

Hi all, I've been having trouble reading in distance matrices (values in lower triangle only) that have been created by other programs. They load with NA values in the upper triangle. Is there an option to read the matrices properly? Thanks, Muri ____________________________________________________________________________________ that gives answers, not web links.

Dear list, I am new with ggplot2 and I have spend quiet some time putting together the following code to create the attached plot. However there's still a few things that I'm having trouble with! I would be grateful if someone can tell me how to fix (1) the colour of my bars into grey scales (2) removing the y-axis (species name) on the right figure to avoid duplication, and (3) fix the

Hola, Vaya desastre al mandar la cabecera de la base de 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

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

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 sampling event) and am having troubld

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

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

Hi everyone, I have a matrix that contains 1000 replicates of a sample of a list of values. I want to group each row (row=replicate) into my defined bin ranges and then calculate the mean and stdev for each of the bin ranges (so I will have 1000 rows but ncol=number of bin ranges). I don't know how to group rows in a matrix matrix according to another vector (which is what I think I need to

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

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

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

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

On 02/11/2014 04:45 PM, Cole Robinson wrote: > On 02/10/2014 06:46 PM, Chris Friesen wrote: >> Hi, >> >> We've run into a problem with libvirt 1.1.2 and are looking for some comments >> on whether this is a bug or design intent. >> >> We're trying to use migrateToURI() but we're using a few things (numatune, >> vcpu mask, etc.) that may need

Dear all, I am trying to use the crr function in the cmprsk package version 2.2 to analyse 198 observations.I have receive the error in solve.default. Can anyone give me some insights into where the problem is? Thanks here is my script : cov=cbind(x1,x2) z<-crr(ftime,fstatus,cov)) and data file: x1 x2 fstatus ftime 0 .02 1 263 0 .03 1 113 0 .03 1 523

Dear all, I am modeling the incidence of recreational anglers along a stretch of coastline, and with a vary large proportion of zeros (>80%) have chosen to use a zero inflated 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: