search for: ggraves

I issued the following command to obtain the std dev for each month. psd<-numSummary(Sal, groups=month, statistics=c("sd")) which resulted in > psd sd n NA 1 6.930340 9367 2319 2 7.847003 10827 1008 3 5.962308 12988 404 4

I am trying to create a routine that would take a time series and generate monthly 25%tile and 75%tile limits based on 12 calendar months. I have succeeded to create a do loop to do this, but can't figure out how to initiate the receiving datafile (in this case "limit.list") without sticking an initial set of zeroes in. I'd like to be able to insert a new column

Having given myself carpal tunnel looking for answer to this ... I have a dataset each column of which has 12 rows in it. I created a variable 'z' as follows: z=1:24 Since I have a large number of these plots to make, and they are a bit complex, I want to want to reference the column I want to plot via a variable containing the name of that column. As follows:

My x axis is a series of daily 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

I have a set of 52 weekly values, one is the desired high lake stage for a week, and the other is the desired low lake stage for each week. It looks like this: week High Low 1 16 14.5 2 16 14.5 3 15.95 14.45 4 15.84 14.34 5 15.73 14.23 6 15.61 14.11 7 15.5 14 8 15.38 13.88 9 15.25 13.75 10 15.13 13.63 .... 52 15 13.5 I can plot out two nice curves

I have 2 datafiles 'target' and 'observed' as shown below (I will gladly email these 2 small files to whomever). X25. And X75. Indicate the value of 25th and 75th-percentile of the target ('what should be') and the observed ('what is'). The i.value is simply the month. > target X i.value X25. X75. 1 one.month 1 10.845225 17.87237 2

Newbie here~ I've spent a ton of time on this, but have to admit I am stuck. I've tried various combinations of strptime and the package CHRON, but still can't get there. My simple task is that I want to plot salinity (S1 or S2) on y-axis and time on x-axis, but I want to use a lay-understandable representation of time (not Julian). If I use the following command, I can make a

I have 3 columns: flow, month, and monthname, where month is 1-12, and monthname is name of month. I can't get the plot to replace the 1-12 with monthname using ticks.lab. What am I doing wrong? plot(flow~factor(month),xlab="Month",ylab="Total Flow per Month", ylim=c(0,55000), ticks.lab="monthname") Thanks Gregory A. Graves Lead Scientist REstoration

Among those users of Primer, stress values greater than 0.3 are interpreted as "questionable". Using both isoMDS and metaMDS (vegan package), the stress values returned are much higher using my own data and using examples provided in R Help. For example Rstress = 8.3, and the stressplot r2 = 0.99 indicating (to me) that the ordination is OK. I am guessing that the "stress"

I have a couple of questions about asterisk 1.6: Can you change codecs mid-call upon re-invite? Can you handle the SDP offer-answer in the 200-ACK instead of the usual INVITE-200? Thanks in advance for any insight. Gary -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi, It seems there is random behavior that causes screening mode to be activated when a user calls and the line answered and then forwarded using a dial command such as: EXEC "Dial" "SIP/13365551212 at 8x8|60&SIP/13365541212 at 8x8 |60&SIP/13365531212 at 8x8 |60|dgF(callFlo-in^3^1)M(record^39ff6274-c0f0-453d-aa05-402a7bd6d567^)" -- AGI Script Executing

I am helping a fellow worker get R up and running, and he has run into a peculiar problem I've not encountered in previous install situations. >From the Rconsole menu choice, he can set CRAN mirror to USA CA2, but when he selects "load packages" a very truncated list of packages appears (many packages don't appear at all). If he issues "chooseCRANmirror()" from the

What am I failing to understand here? The script below works fine if the dataset being used is DNase1 <- DNase[ DNase$Run == 1, ] per the example given in help(nlrob). Obviously, I am trying to understand how to use nls and nlrob to fit curves to data using R. #package=DAAG attach(codling) plot(pobs~dose) #next command returns 'step factor reduced below min factor

If you submit these lines, you end up with variable "vare.dis". I want to export vare.dis to csv. Stuck I am. library(vegan,logical.return = TRUE) #return=true verifies package is available library(MASS,logical.return=TRUE) #return=true verifies package is available data(varespec) #varespec is an example data file in the vegan package vare.dis <- vegdist(varespec)

On Mon, Apr 13, 2009 at 12:41 PM, Dan Dube <ddube-at-advisen.com> wrote: > i use tapply and by often, but i always end up banging my head against > the wall with the output. The proposed solution of Dan's problem posted on R-help was: > do.call(rbind,a) When I use this 'solution' I get 'ERROR: second argument must be a list'. So head on wall continues. My

I have a file, each column of which is a separate year, and each row of each column is mean precipitation for that month. Looks like this (except it goes back to 1964). month X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007 X2008 X2009 1 1.600 1.010 4.320 2.110 0.925 3.275 3.460 0.675 1.315 2.920 2 2.960 3.905 3.230 2.380 2.720 1.880 2.430 1.380