Search Archives

You can do that easily with 'try' ?try 'try' is a wrapper to run an expression that might fail and allow the user's code to handle error-recovery. HTH, Andy

See ?barplot and pay attention to the 'height' and 'beside' arguments. They produce very nice stacked plots. Some of the examples using the VA deaths dataset displays stacked plots. HTH, Andy

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20120822/ea3d5624/attachment.pl>

> x <- readLines(...) > tmp <- file() > writeLines(x[substr(x, 83, 86) == "STD"], tmp) > read.fwf(tmp, ...) I wrapped this approach into function and it worked swimmingly. Thanks all. -Andy

Do you just want to out the data vectors together? Try ?scan, ?c, ?data.frame my.data.frame <- data.frame( names = c(OLDa, ALL, OLDc, OLDa, OLDb, NEW......), values = c(1,5,3,1,2,4.....) HTH, -Andy

Try this. Look at ?flsuh.console if you are on windows HTH, Andy #~~~~~~~~~~~~~~~~~~~~~~~~ i<-0 while(i<100){ ## "do a lot of commands" ##"print i" cat(i, "\n") ## if using windows flush.console() i<-i+1 }

What a great package. I imagined the answer before hand, so I guess it's right... R > rParaNorm(1) [1] 0+42i R > This will really help me through a Bayesian analysis I have to do today. Thanks!

Is there is a way to make square bars in xyplot with type="h"? dat <- data.frame(foo = rep(1:10,2), bar = rep(1:10,2)) xyplot(foo~bar, data = dat, type="h",lwd=20) Thanks! Andy

Have you tried layout? layout(matrix(1:10, 5, 2, byrow = TRUE)) for(i in 1:10){ plot(rnorm(100), type = "l", col = "red", main = paste("Your Variable Number", i, sep = " ")) lines(runif(100), col = "blue") } Does this get you...

There is now a package in R called multcomp for general multiple comparisons that does things similar to the Splus library you mentioned. BTW, a search of the help archives for multicomp or "multiple comparisons" brings this up. HTH, Andy

Look at ?Extract myDF <- data.frame(c1 = c(1,1,1,1,2,2,2,2), c2 = c(23,34,45,45,78,65,45,70), c3 = c(12,15,67,87,23,19,90,32)) mean(myDF[myDF$c2...

This is one way: for (j in 1:10) { assign(paste("x", j, sep="."), rnorm(100)) } BTW, there have been many similar posts like this in the past. They are easily found using the search function at http://cran.r...

Hi all, I'm developing a teaching example doing EOF analysis on global temperature data using the clim.pact library. Does anybody know if there is an existing function to do rotated EOF analysis? E.g., a varimax rotation on...

I had a reviewer request a citation for a package that I had neglected to cite. I followed a similar format to that suggested by Thomas Lumley and also referenced a related book by the package's author. The editor...

Is the GLOBE DEM accessible in R via a package? I usually use raster::getData() when I need elev data but I need to sample a SPDF of a few thousand points with latitudes higher than the SRTM covers. Thoughts...

Has anybody implemented a way of calculating multidimensional variance on a multi band image in R? E.g., http://stats.stackexchange.com/questions/49521/how-to-find-variance-between-multidimensional-points Not eager to reinvent the wheel. Regards to you...

You need to look at ?axis. Try this: y <- c(-4,3,-2,1) x <- c("time 1", "time 2", "time 3", "time 4") plot(seq(1,length(x), by = 1), y, type = "b", axes = FALSE, xlab = "Times", ylab = "Stuff...

> Is everybody writing ArcGIS ASCII rasters recently? The GIS community is hopelessly tied to ESRI. So many people have invested their careers in learning Arc that switching to GRASS is an institutional nightmare. Most of what I do now is...

Look at ?all.equal and ?identical as well as searching the archives for those terms. You'll find many an illuminating thread on precision, floating point arithmetic and other wonders. minX <- 4.2 min0 <- 4.1 sigmaG <- 0.1 Diff...

Curious. > I get the following error message: > > Error in cca.default(cnts94, env94): All row and column sums > must be >0 in the community matrix Check the row and column sums if you haven't: colSums(cnts94) rowSums(cnts94) There...