Re: [R] plot with different symbols and colors according to the factor levels

with(lili,plot(y,conc,pch=sample,col=sample,log="y")) gives a better plot Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO

Re: [R] plot with different symbols and colors according to the factor levels

Hanna, lili<-read.table("lili.txt",header=T) # don't forget to label the row number if it's in your data with(lili,plot(y,conc,pch=sample,col=sample)) Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET:

Re: [R] Retrieving data from survey in R Studio

After: survey <- read.csv(file="http://www.maths.usyd.edu.au/u/UG/IM/STAT2012/r/survey.csv";) do: summary(survey) then bone up on "attach" and "with". Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler IN

Re: [R] How to identify runs or clusters of events in time

Mark, I did something similar a couple of year ago by coding non-events as 0, positive events as +1 and negative events as -1 then summing the value through time. In my case the patterns showed up quite clearly and I used other criteria to define the actual periods. Clint Clint Bowman

Re: [R] New installation

I "experiment" with the Fedora distribution at home but am very satisfied with the Scientific Linux distribution here at work--I'm currently using SL7.2. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl..

Re: [R] find high correlated variables in a big matrix

Are you rying to find clusters of variables according to some distance metric? Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600

Re: [R] Make a box-whiskers plot in R with 5 variables, color coded.

Martin, I grew up in the Midwest of the United States--about as native English speaker as you could find. I was taught exactly the same as you have learned. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl

Re: [R] how to work with time of day (independent of date)

t;tz="PST". I hope I've been missing something simple--any suggestions? TIA Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box

Re: [R] Visualization of people's interactions by participation to parties

step. Similarity of vectors could be an indication of interaction. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600

Re: [R] Reading the non delimited file with no particular patterns in the data to R

1 1 0 0 0 0 0 8 8 2015 1 4 0.617127 1 1 0 1 0 1 0 0 0 15 9 9 2015 1 4 0.399207 1 1 0 0 0 0 0 0 0 10 Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Departme

Re: [R] help for lay person assisting R user with disability

utions. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] dplyr - counting a number of specific values in each column - for all columns at once

Thanks, Dimitri. Burt is the real wizard here--I'll bet he can conjure up an elegant solution. For me, just reaching a desired endpoint is enough. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.uta

Re: [R] dplyr - counting a number of specific values in each column - for all columns at once

May want to add headers but the following provides the device number with each set fo sums: for (dev in (unique(md$device))) {cat(colSums(subset(md,md$device==dev)==5,na.rm=T),dev,"\n")} Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler

Re: [R] dplyr - counting a number of specific values in each column - for all columns at once

It would help if I could see beyond my allergy meds. A start could be: colSums(subset(md,md$device==1)==5,na.rm=T) colSums(subset(md,md$device==2)==5,na.rm=T) colSums(subset(md,md$device==3)==5,na.rm=T) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler

Re: [R] dplyr - counting a number of specific values in each column - for all columns at once

Any problem with colSums(md==5, na.rm=T) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407

Re: [R] Calling Function With Arguments In a Script

?source as in source("pairwise-plots-continuous-vars.R") then plotpairs(first,second,third,wise,title) should get you going Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department

Re: [R] cbind question, please

Perhaps: dog <- 1:3 cat <- 2:4 tree <- 5:7 big.char <- cbind(dog,cat,tree) big.char dog cat tree [1,] 1 25 [2,] 2 36 [3,] 3 47 colnames(big.char)<-c("dog","cat","tree") big.char dog cat tree [1,] 1 25

Re: [R] Superimposing 2 curves on the same graph with par(new=TRUE)

Because the range is so different between the two series, I'd suggest using log="y", ylim=c(500,2) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE:

Re: [R] Superimposing 2 curves on the same graph with par(new=TRUE)

Agreed--I neglected to add the secondary y-axis (shouldn't hit send so fast.) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box

Re: [R] Superimposing 2 curves on the same graph with par(new=TRUE)

Try: plot(Date,MORTSBu,lwd=2,lty="dashed",axes=F,xlab="",ylab="") Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE:

Re: [R] Help with Boxplot

Jim, Have you looked at: names: group labels which will be printed under each boxplot. Can be a character vector or an expression (see plotmath). You could use "" where you want a blank. I believe that "at" may work here also. Clint Clint Bowman

Re: [R] Replace the value with 1 and 0

or: with(yourData,table(year,month,yourData[["rain"]]>0)) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600

Re: [R] How to unzip a .gz file

Alexandra, Although you may not have control over the installation of R, 2.15.1 is very old and should be upgraded--the current is 3.1.2 Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology

Re: [R] Help: The Difference Between Workspace, Script, and History

Ian, ls() is your friend. To learn more about ls(), type ?ls Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600

Re: [R] Problem with

I'd also suggest plotting a wind rose for each month (try openair) to understand the statistical test results. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology

Re: [R] lat/log to meter

Eliza, Would transforming to UTM coordinates work? Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX

Re: [R] Split fixed width data in R

?read.fortran Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] Time outside limits

?rle Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] Writing .csv file

?write.csv 2nd line of Usage: example Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534

Re: [R] 8 fast or 4 very fast cores?

I'm in a similar situation and am looking seriously at a pair of E5-2643v3 (6 cores each-hyperthreaded). Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360

Re: [R] Generating a polygon around points

Your question seems to need an answer to, "How do you find a convex hull on a sphere?" Google has many references. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology

Re: [R] data format setting

?sort, ?unique, and subset come to mind. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407

Re: [R] Generate Binary Matrix

A bit kludgey but how about: dimMat <- matrix(0, 1000, 4) for(i in 1:1000){ while(sum(dimMat[i, ] <- sample(c(0,1), 4, replace = TRUE, prob = c(.3, .7)))==0) dimMat[i, ] <- sample(c(0,1), 4, replace = TRUE, prob = c(.3, .7)) } table(rowSums(dimMat)) Clint Bowman

Re: [R] Error: cannot allocate vector of size 5.2 Gb

5.2 won't go into 4 but there may be more problems. 32-bit or 64 bit operating system? RAM is cheap but will your motherboard support more than 4 GB? And don't forget there are other processes that need to run while you are using R. Clint BowmanINTERNET

Re: [R] phlyloclim help

I suspect the space in the filename could be the throwing the parser off--you may need some "" around the filename or fill the space with an underscore. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@mat

Re: [R] shapiro.test

Greg, For some authors the 4th page from the back should be the first page. Not so for you, however. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360

Re: [R] scatter plot to contour plot

Azam, If you also have z ~ f(x,y) you can produce contours. Otherwise, you could look at package hexbin. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE

Re: [R] creating an equivalent of r-help on r.stackexchange.com ? (was: Re: Should there be an R-beginners list?)

Jason, Thanks--I've found an RSS feed from EPA very useful and will check Stackoverflow's. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360

Re: [R] creating an equivalent of r-help on r.stackexchange.com ? (was: Re: Should there be an R-beginners list?)

Liviu, Thanks for the excellent description of the advantages of SE. However, there is a significant fraction of the population that prefers that information be pushed out to them rather than having to pull it to them. The best system is one that accommodates both equally well. Clint

Re: [R] creating an equivalent of r-help on r.stackexchange.com ? (was: Re: Should there be an R-beginners list?)

Don, Thanks for the brilliant summary of my thoughts. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX

Re: [R] Problems downloading and installing nlme under Linux Mint 15

There is the error message: Error : package ‘lattice’ was built before R 3.0.0: please re-install it Perhaps the problem lies there. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology

Re: [R] find max value in each row and return column number and column name

?which.max should start you down the right path Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360

Re: [R] (no subject)

Just guessing, would the following help? list1 <- c("john", "eric", "steve", "john", "eric", "scott", "john") list2 <- c("john", "john", "john", "eric", "eric"

Re: [R] time series has no or less than 2 periods

ind,data.diff,type="b")) abline(h=0) there is also a pretty strong two week signal--is that of any interest? Now you should be able to decide how to proceed. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET:

Re: [R] Renaming variables

or with(score.plot(X1, X2)) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia

Re: [R] Need help to find out the name of my columns and rows in a data file

You are close--think "names", not "titles", as in rownames or colnames (no reason to completely spell out column). Summary already gave you the column names, so type "?rownames" to learn more. Clint BowmanINTERNET: cl...@

Re: [R] XYZ data

perhaps x<-x+0.1*seq(1,length(x)) Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360)

Re: [R] XYZ data

John, That still leaves a string of identical numbers in the vector. Shane, ?jitter perhaps jitter(X,1,0.0001) Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE

Re: [R] How to define desired numbers to a vector based on the present numbers

I suspect the OP may want rep(1:length(unique(x)), rle(x[order(x)])$lengths)[order(order(x))] to allow for variable numbers of unique values. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department

Re: [R] numeric not equal its value

agic will display them as if they are. A little digging into the r-help archives will find further explanation. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE:

Re: [R] Need code for Liddle Exact test?

Second item in Google, searching for "Liddle Exact test McNemar test" was "cran.r-project.org/web/packages/exact2x2/vignettes/exactMcNemar.pdf" Guess you should have looked before asking. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler

Re: [R] NLS results different from Excel -- Tricky fortunes nomination

nalysis System--am I missing something? Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-75

Re: [R] Reading Wyoming radiosonde data files with RadioSonde package

Sure looks as if that second "--..." line is causing another attempt at parsing for varkey and unitkey which for 70.0 and 59.6 just won't make sense to getsonde. It's a pain but I'd experiment by removing that second "..." line. Clint

Re: [R] Changing the order of months within a year

I'm sure the OP has reasons to reorder the months as shown but his data seems to argue that he shouldn't split the wet season (Oct - Apr) but should, instead, tack the 1972 months (Jan - Aug after the 1971 months (Sep - Dec). Clint Clint BowmanINTERNET

Re: [R] R for Windows 32-bit mode versus 64-bit mode

Most Windows computers in my neighborhood have been upgraded to 64 bit. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600

Re: [R] Removing duplicates without a for loop

?duplicated Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] extracting column at regular intervals starting from different columns

1),TRUE, rep(FALSE, 36-start+1))] # [1] 37 74 111 148 185 222 259 296 333 370 407 444 Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815

Re: [R] Tendonitis and R users

I think I've avoided tendonitis by carefully stretching the affected area when I begin to feel discomfort and, as John suggests, "Carefully watch posture and arm/hand actions to reduce strain". Clint BowmanINTERNET: cl...@ecy.wa.gov Air

Re: [R] Remove several numbers from a sequence

yy[!yy%in%xx] Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] Quadrat counting with spatstat

will define your optimum solution. Soap bubbles with micropipets to inflate them may work equally well. Clint Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407

Re: [R] SmoothTrend in OpenAir

I was able to get the plot by converting the dates in the "date" column to dates, e.g., mytable<-read.csv("MCNP-pH.csv") newdate<-with(mytable,strptime(date,"%m/%d/%Y")) mytable$date<-newdate smoothTrend(mytable, pollutant = "pH") C

Re: [R] How to use ggplot to do the binned quantile plots(one type of scatter plot)?

?quantile on the individual bins, make your deciles, then plot the ten series as usual with your x values at the midpoint of the bins. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology

Re: [R] importing NPDB

erate quite a list. Perhaps you can take it from there. Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FA

Re: [R] Gps fraction of degrees function

I see my pre-coffee fingering hit an incorrect key in the final line--"E" is positive, "W" is negative. -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology

Re: [R] Gps fraction of degrees function

#x27;")[[1]][2], "\"")[[1]][1]) nswe<-strsplit(strsplit(strsplit(coords, "?")[[1]][2], "'")[[1]][2], "\"")[[1]][2] deg.frac<-((sec/60+min)/60+deg deg.frac<-ifelse(nswe=="S"|nswe=="E",-deg.frac,deg.frac) Cli

Re: [R] windrose color ramp issue

library(openair) ?windRose -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia

Re: [R] Web based R-help not a list

e more important than number of messages. Reproducible code usually is sufficient to demonstrate the problem. Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE:

Re: [R] How to avoid ifelse statement converting factor to character

most all those posts were right in sourcing the problemits just that nobody actually offered a viable solution the problem is that the new level "C" was not one of the original levels in $social status add "C" as a level and then just do the ol fashioned way and it works just fine do this: d

Re: [R] representing wind date using windrose

I'm also very impressed with openair <http://www.openair-project.org/>, also Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-

Re: [R] round up a number to 10^4

Or does the middle number have two digits switched? 76131.17 would round up to 8 very nicely. -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407

Re: [R] Bug or feature? sum(c(a, b, c)) != (a + b + c)

?all.equal -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] How to use as.Date (or something else) with "31-Jul-2010 23:59:00"

?strptime -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] timeDate with month designated by three letters.

strptime("04-MAY-11 1428",format="%d-%b-%y %H%M") [1] "2011-05-04 14:28:00" -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE:

Re: [R] need to delete by time, not date

You may be happy using: strptime(paste(Date,Time,sep=" "),"%m/%d/%y %H:%M")$hour -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 4

Re: [R] Averaging uneven measurements by time with uneven numbers of measurements

n subtract each succeeding reading from that to obtain the amount of feed dispensed. Now plot that value every 30 minutes. If you are interested in the variation of feed dispensed over a half hour interval, that can be easily obtained by accummulating those half-hour readings.

Re: [R] Averaging uneven measurements by time with uneven numbers of measurements

I'd be tempted to do a robust fit (loess?) to the data with a relatively small span (I'm assuming that there are errors in the measurements and some degree of smoothing is acceptable) then predict the fit at a regular interval (e.g., every 30 minutes). -- Cl

Re: [R] Null

with(fcv,shapiro.test(case)) -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534

Re: [R] Compiling a FORTRAN program under Windows 7

You are compiling a subroutine not a program and you compile line should read: gfortran testit.f -c testit.o You then reference that object code testit.o in your final loading stage after compiling other routiens and the main program. -- Clint BowmanINTERNET: cl

Re: [R] filled contour plot with contour lines

Steve, I use filled.contour but have a semicolon, ";", between the two axis calls. Do you need one after "axis(1, seq(45 , 58, by = 10000))" ? Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler

Re: [R] Quick recode of -999 to NA in R

Amen. Ditto for "-999.000", "-999.00" and all of the other ones that various (usually Fortran) programmers have used. Has the most recent Fortran standard come around to understanding NA? -- Clint BowmanINTERNET: cl...@ecy.wa.gov

Re: [R] mean of runoff for several years

Dominique, Let's call your dataframe, df. Then try: with(df, boxplot(runoff ~ month)) which will plot the daily runoff distribution by month. Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.uta

Re: [R] Finding length of unique numbers in a vector

?rle -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7534 Olympia, WA 98504-7600

Re: [R] spacing of color key in filled.contour

ping instead of transforming your data. Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407-7

Re: [R] Remove 100 years from a date object

n paste them together and then as.Date(). Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(

Re: [R] previous business day

May also wish to "or" in a check for holidays. -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600

Re: [R] rounding up (always)

x<-runif(10,1,10) # generate 10 numbers as.integer(log(x,10)+1) -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47

[R] Dealing with Non-Standard Hours

A data set I obtained has the hours running from 01 through 24 rather than the conventional 00 through 23. My favorite, strptime, balks at hour 24. I thought it would be easy to correct but it must be too late on Friday for my brain and caffeine isn't helping. TIA for a hint,

Re: [R] labeling outliers with subject numberss

SubNo[identify(rep(1,8),mydata$score)] -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:(360) 407

Re: [R] Peak Over Threshold values

obably want to ignore the end point at length(walevel). -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX:

Re: [R] Identifying breakpoints/inflection points?

Charlotte, Try: library(msProcess) # you may have to install msProcess year[peaks(birds.pr$fit)] -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407

Re: [R] Identifying breakpoints/inflection points?

Charlotte, Try: birds.lo <- loess(piproute~year,span=.25) # play with span to see your desired pattern birds.pr<-predict(birds.lo, data.frame(year = seq(1967, 2009, 1)), se = FALSE) # plot($year,birds.pr$fit,ylim=c(0,5)) par(new=T) plot(year,birds.pr$fit,pch="+",col=2,ylim=c(

Re: [R] two questions for R beginners

forgotten too much of my fortran and C programming skills to contribute directly to the R Project. Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 4

Re: [R] how to draw image of a data distribution? Thanks

with(p,plot(x,y,pch=20,col=1+as.integer(value))) -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Box 47600FAX

Re: [R] Temperature Prediction Model

that you expect a day-of-week dependence? If so, I'd be more comfortable if you used more than one week to develop it. I fear that you've gotten me quite interested in this analysis, good luck. Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Mod

Re: [R] Temperature Prediction Model

th a sensor. Do we also have to detect the pawing of a "micro-mouse" as well? The collected data also seem to have other parameters which would be valuable--are you limited to just temperature? Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Mo

Re: [R] Best R text editors?

s Curran, of the UK Unix Users Group): "I can touch-type, but I can't touch-mouse" That's a strange disability. It took me several months to learn to touch-type (and years later I'm still not very good at the top-row numbers or the special symbols on them),

Re: [R] Joining Characters in R {issue with paste}

ect.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler INTERNET: cl...@math.utah.edu Department of Ecology VOICE: (360) 407-6815 PO Bo

Re: [R] How do I plot a line followed by two forecast points?

or y values to obtain the desired effect). Clint BowmanINTERNET: cl...@ecy.wa.gov Air Dispersion Modeler INTERNET: cl...@math.utah.edu Air Quality Program VOICE: (360) 407-6815 Department of Ecology FAX:(360) 407

Re: [R] handle graph size in eps

Try adding paper = "special" to your postscript arguments. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Dispersion Modeler INTERNET: cl...@math.utah.edu Air Quality Program VOICE: (360) 407-6815 Department of Ecology

Re: [R] counting strings in a column

Try: a<-c("o","o","a","a","g","o","a","p") table(a) a a g o p 3 1 3 1 Clint BowmanINTERNET: cl...@ecy.wa.gov Air Dispersion Modeler INTERNET: cl...@math.utah.edu A

  1   2   >