search for: aolinto_r

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

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

Hello, I want to know how do R calculates the number of intervals between tick-marks in the y axis in a plot. I'm making a three y-axes plot and this information would help me a lot. Thanks in advance. Antonio Olinto ---------------------------------------------------------------- Webmail - iBCMG Internet http://www.ibcmg.com.br

Hello, About looping, consider the example: for (i in 1:23194) { dat.stat[i,c(2:8)]<-quantile(dat.bat[BL==block[i],2],prob=c(0,0.025,0.25,0.5,0.75,0.975,1)) print(i) } I'd like to have the value of "i" printed for each loop (step). As I could see the values of "i" are shown on screen only after all the work is done. Thanks in advance for any suggestion. Best

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

...n Dedman Mathematical Sciences Building (Building 27) Australian National University, Canberra ACT 0200. On 16/06/2009, at 2:51 AM, r-help-request@r-project.org wrote: > From: "Gerard M. Keogh" <GMKeogh@justice.ie> > Date: 15 June 2009 8:10:34 PM > To: Antonio Olinto <aolinto_r@bignet.com.br> > Cc: r-help@r-project.org, r-help-bounces@r-project.org > Subject: Re: [R] books on Time series > > > Antonio, > > You basically need the cross-correlation function in R which is very > easy > to use - just look up the examples in ?ccf. > So, if an...

Hi, I have a dataframe with 3 columns: month (1 to 12), length and length class. I'm trying to make a cross table with the counts of lengths per length class and month. I don??t have all classes per month. | 1 2 3 ... ------------- 120 | 0 1 0 150 | 5 6 0 170 | 3 0 7 ... I have already tryed many commands, without success - unfortunately. May someone help me? Thanks in advance!

Hi, I have doubts in using subset command. I have a list of, lets say, 15 species and I want to make a subset with only 2 of them. I??m the command data2 <- subset(data1, species=="sp1"|species=="sp2") Nevertheless, when I ask for the summary (summary(data2)) the others species names still apearing. Also I tried data2<-data1(data1$species %in%

Hi, I have a data set with values (L) varying from 1 to 700, with repeated numbers, and I want to sample them according to the probabilities given by a logistic curve P=1/(1+e(-r*(L-Lc))) where r gives "S" the inclination and Lc the "rotation point". P ranges from 0 to 1. As I could see I cannot use this P as prob in sample function. I'd need the probabilities as given by

Hi all, I used to connect internet via a proxy. Before update packages I wrote in R Sys.putenv("http_proxy"="http://proxy3.redegov.sp.gov.br:80/") Nevertheless the way the connection is done has changed. For example, in the browser the proxy is not indicated and I have to give an username and a password to have access to internet. I read the FAQ and the help for

Hi, I?m trying to make a barplot with the following dataframe, with information on relative frequency per sediment type (ST) for some species: Species ST1 ST2 ST3 SP_A 10 60 30 ... At x-axis are (should be ...) the species names and at y-axis the frequency per sediment, in stacked bars. I tried to use barplot command but with no results. Could anyone help me on this? Thanks in

Hi, I'm using nmds command (library vegan) to analyze some fishing data. I'd like to plot not only points, but also the names of species and stations in a specified position. I used the command text(nmds$points[,1], nmds $points[,2],labels=row.names(nmds $points),pos=3,cex=0.5) But the labels are sometimes overlapped. Is there any way to use identify, or a similar command, to plot the

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

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, I want to make a general routine to draw barplots with numbers plotted above each bar. See the example below. I could not place the numbers on the middle of each bar because I could not calculate the right position of each x-axis tick. axTicks(1) indicated a unitary step, but it does not seem work. I appreciate any help or suggestions. Best regards, Antonio Olinto