similar to: data frames

Dear R-list members, Some months ago I wrote a message on the usage of gnls (nlme library) and here I come again. Let me give an example: I have a 10 year length-at-age data set of 10 fishes (see growth.dat at the end of this message) and I want to fit a von Bertalanffy growth model, Li= Linf*(1-exp(-k*(ti-t0))) where Li = length at age i, Linf= asymptotic length, k= curvature parameter, ti=

Thanks for your reply Nassar. What I really want is to overlay (like add=TRUE) the plots 1 and 2. With plots "at the same screen" I meant overlaid and not one beside the other, i.e., only one plot with all the points (1 to 10) in blue (or any other color), the points 5 to 9 in red and a regression line from points 5 to 9. Best regards, Antonio Olinto ----- Original Message -----

Hi, I want to randomly pick 2 fish born the same day but I need those individuals to be from different families. My table includes 1787 fish distributed in 948 families. An example of a subset of fish born in one specific day would look like: >fish fam born spawn 25 46 43 25 46 56 26 46 50 43 46 43 131 46 43 133 46 64 136 46 43 136 46 42 136 46 50 136 46 85 137 46 64 142 46 85 144 46 56

Hi, I'd like to find which values of x will give me a y. In other words, in the example of a gaussian curve, I want to find the values of x that will give me a density, let's say, of 0.02. curve(((1/(sqrt(2*pi)*10))*exp(-((x-50)^2)/(2*10^2))),xlim=c(0,100)) Thanks for any help, Antonio Olinto ------------------------------------------------- WebMail Bignet - O seu provedor do

Dear R-list members, Cynthia M. Jones wrote a paper (Fitting growth curves to retrospective size-at-age data, Fisheries Research 46(2000):123-129; abstract at http://www.elsevier.nl/gej-ng/10/19/44/70/24/37/abstract.html)where the SAS procedure MIXED, Macro NLINMIX (Littell et. al., 1996)was used to estimate the von Bertalanffy growth function parameters assuming that data from the same fish are

Dear R-List members, I want to check if a set of measurements follows better a gamma or a lognormal distribution (see data below). Using shapiro.test I can test for normality (shapiro.test(log (Lt)). To test for gamma (and normal) distribution I would use ks.test but I need to specify its shape and scale. How should I calculate these values in R? I tried > Lt.fit <- glm(Lt ~ 1,

Hello all, Here goes one of my first functions. I want to make a nonparametric multiple sample comparison with unequal sample sizes (see Zar?s Biostatistical Analysis, 3rd. Ed., pg. 201 Example 10.11, pg. 288 Example 11.10). In the real world, I want to compare samples of fish length captured with different fishing gears. After using the Kruskal-Wallis test I want to check the differences

This works fine: fit42<-nls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) But this, identical except using gnls, doesn't converge: fit43<-gnls(Vfs~SSlogis(Months,Asym.Int+Asym.Group*Groupdum,xmid,scal), data=df, start=c(Asym.Int=22,Asym.Group=5,xmid=2,scal=6), na.action=na.omit) Error in gnls(Vfs

Hi, I'm a total beginner in R and this question is probably very simple but I've spent hours reading about it and can't find the answer. I'm trying to reshape a data table from long to wide format. I've tried reshape() and cast() but I get error messages every time and I can't figure why. In my data, I have the length of two fish from each family. My data table (called

Hello R-list members, I’m trying to make a min-avg-max plot (like a boxplot, without box and outliers, i.e., for each class of "x" a vertical line indicating the range (max-min) and an horizontal line or point indicating the mean) with data like: Year Avg Min Max 98 10 7 13 99 12 10 15 01 11 6 14 Which command should be used? Many thanks, Antonio Olinto

Dear list fellows, I want to study time series and use R to analyse time series of fishing data from several species (landings and cpue) investigating the correlation between them and with environmental factors (water temperature, wind, etc.). Searching at Amazon I found three books with examples in R: Time Series Analysis: With Applications in R by Jonathan D. Cryer and Jonathan D. Cryer

Hi to all, So I'm also a new user. I downloaded the program last week and I think it's great. Thanks to those who have developed R. I have a special interest in GLM as a tool to describe fisheries and its variables and I'm just begging to study it. As I could understand there's two types of GLM sun of squares: in "type I" the factors are added in sequence and

Hi, Some time ago I asked about how to use comma as decimal separator in plots and Mr. Paul Murrell wrote: ------- You could try something with axis() and chartr(), like ... par(mfrow=c(2,1)) plot(1:10/11, rep(1, 10), main="Standard X-Axis") plot(1:10/11, rep(1, 10), main="Customised X-Axis", axes=F) axis(1, at=pretty(1:10/11),

Hello, I have a doubt in using the function step (step wise) to select glm models. Usually I apply the gamma distribution to analyze fishery data. To select the terms I use a routine where I first compare single term models to the null model (eg. U~1 vs. U~depth; U~1 vs. U~latitude; etc. ? where U= abundance) and, by means of the result given by a likelihook function applied for each comparison,

Hello all, I have a data frame with year, month, species, fishing gear and catch (Y, M, S, F, C) and I want the sum of C by Y for species "A" and fishing gear "trawl". I tried things like aggregate(C[S=="A" & F=="trawl"], list (Year = Y[S=="A" & F=="trawl"]), fun=sum), but it didn't worked. To overcome this problem I did

As to call files of another softwares(as, Excel, minitab) to be manipulated inside of R-plus? Respectfully Regiane -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20011112/23b21666/attachment.html

Hi, I want to write in x axis label "fitted value of lambda" (lambda in greek letter). xlab=expression(lambda) gives the "lambda", I tryed things like xlab=paste ("fitted value of ", expression(lambda)) but I didn't get the greek letter. Thanks in advance for any hint. Antonio Olinto ------------------------------------------------- WebMail Bignet - O seu

Dear list, I am modelling fish distributions using the glm-function followed by the step-function, and then want to cross-validate the model via the cv.glm-function from the {boot} package. I am working on fish distributions on coral reefs. The code I have works for one fish species. Since I have 227 fishes, I wrote a loop. Now the cv.glm-function comes up with an error message: "Error in

> Background: > OS: Linux Ubuntu Dapper > release: R 2.3.0 > editor: GNU Emacs 21.4.1 > front-end: ESS 5.2.3 > ----------------------------- > Colleagues > I have two spatial datasets (latitude, longitude, fish eggs) and (latitude, longitude, fish larvae) at the same 280 stations (i.e. 280 cases). I want to determine if the 2 datasets are spatially correlated. In other

Hello all, I've come across an online posting http://www.biostat.wustl.edu/archives/html/s-news/2001-10/msg00119.html that described how to get confidence intervals for predicted values from predict.glm. These instructions were meant for S-Plus. Yet, it generally seems to work with R too, but I am encountering some problems. I am explaining my procedure in the following and would be most