similar to: sum(..., na.rm=TRUE) oddity

Hi All, I am having trouble installing the sm package, and only the sm package. My box is RedHat 7.3 on a PIV. Readline 4.3-3 is installed (suited for RH 8, but same error under readline 4.2 for RH 7.3), but it seems to be looking for the readline library in the wrong directory. Will installing readline 4.1 do the trick? What else could be done? Rohan Sadler [root at rsadler R]# R CMD

Hi all, I wish to implement various landscape ecology measures through R, such as: box counting dimension; twist number statistics; contagion and lacuniarity indices; angular second moment; adjacency measures; dominance indices; etc ... Some of the measures can be applied to shape analysis and classification. Is anyone implementing any of these measures? If so I would like to contribute and

Hello all, I'm attempting to conduct spatial analysis of trees within a plot. I want to see if the trees are spatially correlated to soil characteristics, say pH, or moisture content. I think one way to do it is with mpl, however, my soil characteristics were not taken at exactly the same locations as my trees and further, the vectors aren't the same length. I'm getting the

Dear All, This is a question to sound out possibilities. I am with the Faculty of Natural and Agricultural Sciences at the University of Western Australia, representing a few of the more statistically minded in the faculty. Essentially, there have been problems in the past with software support, changing over statistical software, and paying lots of money for it. In R you have an advanced

Dear R-Listers, I have spatial data on an equidistant rectangular grid, similar to topographic data. I know that there are quite a few R-packages or base functions that provide nice iso-contours plot, but I don't want a plot, just the smoothed isocontour line of ONE level (e.g. 10 mm). Data sets are large, so it would be preferable if the availability of regular grid data could be exploited,

Hi all, The suggestions were great. However, after installing readline-devel I had to install ncurses-devel before getting sm to work. Can this be put into the RH 7.3 Readme when downloading R rpms (and other distbns like Mandrake where there may be a similar problem)? Thanks for your help. Rohan Sadler -- Ecosystems Research Group (ERGO) School of Plant Biology (Botany), Faculty of

Definitely 2 camps on this issue; so why not compromise with a drop-down menu for the most frequently used processes, the user responds with the necessary parameters for his choice, and R then writes the source statements on the terminal and executes them. The user follows the familiar gui procedure, is thereby automatically introduced to the command statements involved, and R's

I'm having trouble with tapply, it seems to give me NA for the mean of a set like {2, 3, 4, NA}, can I get it to ignore the NA's and give me a mean like say, three? Am I just doing something stupid? Any help would be apprecitated, Also are there PDF versions of any of the R stats books around that I might be able to get my hands on? Bryn

Hi, We have been running R ver 1.4.1 on a redhat linux (7.1)cluster without problem for quit some time. Now, some user wanted me to install version 1.6.1, and it compiled fine, but it crashes with "segmentation fault" when you run it. I used gcc296 (default) for compiling, Should I change compiler or is there some other thing I could do to avoid this problem. Regards Per-Anders

Hi, Have you seen the tcltk package? I believe it is a better option than the one you outline in your message. Why develop a Windows-only GUI? Maybe in a few years your teaching department will give up Windows for Linux (for the reasons you wrote about R: it is free, advanced, and adaptable). There is an article in R News (vol.1, n?3) by Peter Dalgaard on the tcltk package, and the latter

DeaR comRades: I have a 2D spatial binomial process as shown in the data and code below. I am plotting the number of trials and the number of successes in the spatial binomial experiments and would like to draw the spatial cells were the trials and successes were counted, i.e. a partial grid in the plot only for those cells where there is a number. The cells are 2x2 km cells. The count of Trials

Hello, I am using the following code to create ppp files from csv data and map shape files, but I am getting some errors which I have been unable to fix by searching them online: library(spatstat) library(maps) library(maptools) NYC2<-readShapePoly("nybb.shp") # this is a map of the NYC boroughs without waterways and no census tract divisions (but it does include lines separating

Hi All, After a lengthy battle just to get the package installed, I am not able to actually use Rgraphviz to generate any plots. I tried just using the sample code in the documentation (http://www2.warwick.ac.uk/fac/sci/moac/people/students/peter_cock/r/rgraphviz/) and I get the following: *> > library(Rgraphviz) > test.matrix<-matrix(rep(c(0,1,0,0), 9), ncol=6, nrow=6) >

Dear, I am trying to visualise a time-progressing line (it's supposed to represent spread patterns) using brew package and Google Earth. The idea is to have a function which takes start and end point geographic coordinates, as well as number of intervals to chop the path up, and returns the collection of points segmenting this line. Unfortunately my calculations fail for large distances,

Hi, I noticed the following problem with xy.coords() in R 2.2.1-patched (version info at the foot of this email) and R 2.3.0 unstable (subversion no: r37123): > xy.coords(x = matrix(1:20, ncol = 2)) Error in xy.coords(x = matrix(1:20, ncol = 2)) : argument "y" is missing, with no default > xy.coords(x = matrix(1:20, ncol = 2), y = NULL) $x [1] 1 2 3 4 5 6 7 8 9

SUMMARY: Specific problem: I'm regridding biomass-burning emissions from a global/unprojected inventory to a regional projection (LCC over North America). I need to have boundaries for Canada, Mexico, and US (including US states), but also Caribbean and Atlantic nations (notably the Bahamas). I would also like to add Canadian provinces and Mexican states. How to put these together? General

Dear all,   I have the following written example > coords <- "51°30'48.58\"N" > > as.integer(strsplit(coords, "°")[[1]][1]) [1] 51 > as.integer(strsplit(strsplit(coords, "°")[[1]][2], "'")[[1]][1]) [1] 30 > as.numeric(strsplit(strsplit(strsplit(coords, "°")[[1]][2], "'")[[1]][2],

What is plot.new? and how can i get it to work so i can load other data? > library(splancs) > area = 6*4 > lambda = 1.5 > N = rpois(1,lambda*area) > u = runif(N,-2,4) > v = runif(N,0,4) > plot(u,v,asp=1) > h = chull(u,v) > h = c(h,h[1]) > plot(u[h],v[h],"1",asp=1) Error in plot.xy(xy, type, ...) : invalid plot type '1' >

Hello! My code below creates a data frame and a plot for it. However, I can't figure out how to add a legend that is not ON the plot itself, but outside of it (e.g., to the right of my graph or below it). I tried something: I put a line par(xpd=T, mar=par()$mar+c(0,0,0,4)) right before my plot command), but that screwed up all my gridlines - they covered all graph and do not coincide with

Hello, everybody! I have a code below that creates a data set and then a stacked bar chart based on that data set. No need to look at it - just notice please that my horizontal axis is a date varible (x=my.data$date). I have a question about the last 2 lines of this code: grid(nx=NULL,ny=NULL,col = "lightgray", lty = "dotted",lwd = par("lwd")) axis(1, las = 2) Could