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

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

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] 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] 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] 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] 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] 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] 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] 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] 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] 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)

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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] Best R text editors?

d provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minim

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 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] creating an equivalent of r-help on r.stackexchange.com ? (was: Re: Should there be an R-beginners list?)

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 USPS

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] 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] 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] 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] 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] 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] 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] 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] (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] 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] 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

[R] Dealing with Non-Standard Hours

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

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

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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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 = 1))" ? Clint -- Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler

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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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 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] 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] Comparison between R and MATLAB

The easy way around that is to create an account, "Mathworks", with a common group that all who will use MatLab belong, then su - Mathworks should satisfy the license manager. Clint BowmanINTERNET: [EMAIL PROTECTED] Air Dispersion Modeler INTERNET:

Re: [R] A Maze Generator

I'm familiar with mazes that have a designated starting and ending points but don't see any here (or is this just a maze Clint BowmanINTERNET: [EMAIL PROTECTED] Air Dispersion Modeler INTERNET: [EMAIL PROTECTED] Air Quality Program

Re: [R] things that are difficult/impossible to do in SAS or SPSSbut simple in R

So how does SAS compare with one of the specialty languages such as perl. I've found the combination of perl and R to work quite satisfactorily (as long as I don't confuse the syntax and functions available in each.) Clint Clint BowmanINTERNET: [EMAIL PROT

  1   2   >