search for: salin

Hi, I just started with R and I found that there are many options to rearrange the data to do mixed models. I want to use the reshape function. I have 2 between subject variables and one within. I was able to change the data structure but still - the result of the aov functions are calculating everything as a within subject. the table looks like this: SerialNo breed treatment distance_1

Hello all, I am an undergraduate student who is having syntax issues trying to get error bars on my graph. This is the data, which I assigned the name "Saline" to. Time Average SEM 1 -20 0.000000 0.0000000 2 3 30 0.000000 0.0000000 4 45 3.227902 0.7462524 5 60 5.066664 1.1623944 6 80 6.107491 1.5027762 7 110 6.968231 1.3799637 8 140 7.325713 1.2282053 9 2...

Hello, I'm doing the following: library(ncdf) library(fields) library(animation) saline <- open.ncdf("salinity_1990.nc") salt = get.var.ncdf(nc=saline, varid="Salinity") # create an animation of the complete temporal domain in the ncdf file. saveHTML({ for (i in 1:364) { image.plot(salt[, , i]) } }, img.name = "salinity.img&q...

...col = 10, nrow = 10) variables <- c("lightlog", "light", "dstb", "species_count", "masslog", ? ? ? ? ? ? ? "soil_carbon", "soil_organic", "soil_low_flooding", ? ? ? ? ? ? ? "soil_high_flooding", "soil_salinity") rownames(grace) <- colnames(grace) <- variables diag(grace) <- 1 ## Coefficients from the paper. grace.coefficients <- c(0.858, 0.667, -0.251, -0.699, 0.06, 0.012, 0.552, 0.547, ? ? ? ? ? ? ? ? ? ? ? ?0.327, 0.776, -0.404, -0.794, 0.157, 0.120, 0.439, 0.462, ? ? ? ? ? ? ?...

...dates (e.g., 01/01/2000, 01/02/2000, etc.) from 2000 to end of 2008. The default only gives me 4 ticks. I want more. Why doesn't this work? sdate<-as.POSIXct(strptime(date,format="%m/%d/%Y")) plot(ppt~sdate,type="l",ylim=c(0,47),col=1,lwd=1,pch=16,ylab="Salinity, psu",xlab="Year",las=1,main="Duck Key Salinity",xaxt="n") year.text=c("2001","2002","2003","2004", "2005", "2006", "2007","2008", "2009") axis(1, at=c(1:9),labels=year....

...been able to get the break and labels to be correct, however, I can't seem to get the data to match the axis anymore. I must be using the axis.break() in plotrix incorrectly, but I cannot see where my issue is. This is what I have so far. ###### library(plotrix) par(mar=c(6,8,4,4)) ###Data Saline <- structure(list(Time = c(-20L, 0, 30L, 45L, 60L, 80L, 110L,140L,200L, 260L, 320L), Average = c(119.250,118.750,117.500,132.75,151.875,159.75,142.75,160,168,167.125,143),SEM=c(2.211,2.569,2.665,5.435146394,6.208741369,8.363550657,8.51349469,14.30284687,15.93865792,16.76541326,13.796)), .Names...

The University of Miami's Rosenstiel School of Marine and Atmospheric Sciences (RSMAS) is seeking a research associate to work with hydrographic data. The objective of the project is to develop statistical models for estimating salinity and density from observations of temperature together with knowledge of location and time. Because the relationship between salinity and temperature can be highly variable due to meandering fronts, it is important to be able to recognize from which side of a front the data come. Consequently, d...

The University of Miami's Rosenstiel School of Marine and Atmospheric Sciences (RSMAS) is seeking a research associate to work with hydrographic data. The objective of the project is to develop statistical models for estimating salinity and density from observations of temperature together with knowledge of location and time. Because the relationship between salinity and temperature can be highly variable due to meandering fronts, it is important to be able to recognize from which side of a front the data come. Consequently, d...

Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 --------------------------------- Colleagues I'm plotting fish egg densities in temperature-salinity space (i.e. mean temperature ~ mean salinity | fish egg densities) for two survey years. The bivariate plot is overlaid on a temperature ~ salinity plot (not means) showing the distribution of water masses, so I can see the association of the eggs with water masses. It is all working, but I ha...

Hi list. I'm plotting pluviometric (Rain) data as a barplot, and then adding the salinity variable to this plot as lines. Obviously as these Y scales are completely different the salinity appears at the lower part of the graph extremely compacted. I need to plot the line at the exactly same area of the barplot but with its own Y axis (at the right), so the salinity can use the plot a...

Dear R Wizards, After much frustration and days of confusion I have finally broken down and am asking for help, which I don’t like doing, but I just can’t figure this one out on my own. I’ve conducted a laboratory experiment testing the effects of temperature and salinity on whether or not a biological event will occur (Go or NoGo). I’ve coded the factors temperature and salinity as factors for the binomial glm, and I haven’t had any trouble fitting the model and checking assumptions. I am however having trouble with the predict.glm function. I want to create...

Colleagues, I am analyzing data collected during oceanographic cruises. We have conducted many cruises over the last decade. On each cruise we visit ~50 stations. At each station (termed EventNum)we lower an instrument that measures depth, temperature, salinity and oxygen every few seconds as it is lowered through the water. Data from all EventNums on all cruises are stacked, generating a data file that might look like CruiseID Year Season EventNum Region Depth Salinity Temp Oxygen 1 CF0101 2001 Spring 849 Lower 0.056 21.1753 15.5626...

Hi, I have used R for some time, but managed to avoid writing loops. But this time I am afraid there is no way around it. I have a dataframe with time and salinity (see below). I would like to extract the time intervals where salinity changes by less than 0.05. So using the values below this would mean that a subset was made from 2003-07-19 to 2003-07-24, where the change in salinity exceeded 0.05, and then stating again from 2003-07-24. In other words, I...

I have a dataset relating the effects of engineering works on the level of salinity in a river before and after the works. I have modelled this using linear mixed effects models to determine if the significance and level of the response to the works. I am wanting to calculate the se of difference at several points. I am not sure of how to calculate the difference and which es...

...shingle(mn.ts.e.both$eggs2.Pilch.Larv, intervals = int) ## define plot structure out1 <- xyplot(mn.t ~ mn.s | egg.counts, data = mn.ts.e.both, groups = year, xlim = c(35,38), ylim = c(12,24), xlab = "mean salinity", ylab = "mean temperature (deg. C)", main = "2001 and 2002 egg densities in Temperature-Salinity space", aspect = "xy", jitter = T, layout = c(1,6), auto.key=TRUE,...

Hi, We are analizing the relationship between the abundance of groupers in line transects and some variables. We are using the quasipoisson distribution. Do we need to include the length of the transects as an offset if they all have the same length?? Also, can we include in the gam models variables that are measured at different spatial scales? We have done an analysis to see what variables

...sulted a few local statisticians and I’ve not yet received a clear answer. Before I start an experiment I want to conduct an a priori power analysis test for a Generalized Linear Model, family = poisson. The response variables will be counts. I want to look at the effects of 5 temperatures and 3 salinities, and their interaction in 15 different combinations (5 temperatures each testing 3 different salinity levels). I am hoping to learn the power I would achieve if I used a sample size of 18 in each of the combinations, for a total sample size N = 270. Does such a power analysis exist? Can...

Dear Mr. and Mrs. I'm very grateful for these software and this list. My question is: when a use linear multiple regression (lm()) for my data, abundance ichthyoplankton ~ salinity + temperature + month of the year(f1 is a factor: 1 for january, 2 for february, ..., 12 for december), the summary() of results is ... Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.726347 0.285607 2.543 0.011217 * temperature -0.011410...

I have a data assimilation problem that might be amenable to the use of GAMS, but I am not sure how feasible it is to implement. I was told the R mailing list was a great resource. My observations are spatiotemporal salinity in the San Francisco Bay at a number of instruments over a few days. The thing that I want to fit is the initial condition for a salt transport model at the beginning of this time period. The spline basis functions, parameters and curvature penalties would all be purely spatial, though the cross...

...n. Strangely enough I am unable to find a function that could actually pull this off. I am unable to use replications or count the rows using nrow. What I am trying to obtain is the number of observations fitting something along the lines of "data$total[data$spp==1&data$gyro==1&data $salinity==0&data$day==i]" Can anyone help me?