similar to: some pairs() questions

Hello, I'm having problems reading a shapefile with read.shape (maptools). I'm absolutely sure my file is there, but I get "no such file". The wd is ok, since read.table for example does find the file. > getwd() [1] "D:/somedirectory/R scripts" > read.table("cities.shp") Error in read.table("cities.shp") : empty beginning of file In addition:

Hi, I'm trying to find a way to determine how many digits a number has. I tried using nchar(paste(number)), but unfortunately paste will reduce 8.00 to "8". Any thoughts? Pieter -- This message was sent on behalf of pieterprovoost at gmail.com at openSubscriber.com http://www.opensubscriber.com/messages/r-help at stat.math.ethz.ch/topic.html

Hello, I'm trying to find out how I can move the y axis label of a horizontal barplot further away from the axis, in order to avoid overlap between the names (which I have rotated using las=1) and the label. I have been fiddling with the graphical parameters but nothing seems to work. Maybe a graphical representation of the graphical parameters would be helpful (cf the CSS box model).

What would be nice is to have something like a "robust" task view... Andy From: Berton Gunter > > Many thanks for this Martin. There now are several packages > with what appear to be overlapping functions (or at least > algorithms). Besides those you mentioned, "robust" and > "roblm" are at least two others. Any recommendations about > how or

Hello, I'm trying to pull off a certain graph using splom, and can't quite get my panel functions right. Basically, the equivalent using pairs would be something like this (using iris data set as an example): panel.corval <- function(x, y, digits=2, prefix="", cex.cor,col,pch) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1))

I have been trying to find a command similiar to bs or ns to add a periodic/or cyclical splie term to a regression model in glm.nb I have been using ns to fit a spline to day of the year (1:365): glm.nb(CNT ~ WKDY + ns(DY,df=6) + HOLIDAY_FLAG + Trend) but I think a periodic spline might be more appropriate. Any suggestions would be appreciated. best, Spencer Jones, M.Stat. NLM Fellow Dept.

> > Leaf Sun wrote: > > The histogram is highly screwed to the right, say, the range > > of the vector is [0, 2], but 95% of the value is squeezed in > > the interval (0.01, 0.2). I guess the histogram is as you wrote. See http://web.maths.unsw.edu.au/~tduong/seminars/intro2kde/ for a short explanation. > -----Original Message----- > From: Berton Gunter

Peter You question is difficult to answer without more information about the distribution of your residuals. Different residual patterns call for different transformations to stabilize the variance. One very common form of heterocedasticity is increasing variance with increasing values of an independent predictor, i.e. the variance of the residuals of y=x increase as x increases. In this case a

I am immediately reminded of something I read which goes "A sufficiently trained statistician can read the vagaries of a Q-Q plot like a sharman can read a chicken's entrails, with a similar recourse to scientific principles. Interpreting Q-Q plots is more a visceral than an intellectual exercise. The uninitiated are often mystified by the process. Experience is the key here."

Hello, I have a customized pairs () fonction as follows that displays correctely my data. ------------------------------------------------------------------------ panel.cor1 <- function (x, y, digits=2, prefix="") { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- cor(x, y,use="pairwise.complete.obs",

Hi, I'm having some problems with a bubble plot (ps package). I don't want tick marks on all four sides (just two), I want to have a smaller font size, and I would like to be able to define bubble sizes shown in the legend (now it shows 0, 0, 0, 9.747 and 4265.757 which is not really convenient. Passing some of the standard plot arguments didn't help (in fact, nothing changed).

Ah... I searched for half an hour for this function... you know, the help function in R could really be a lot better... But wait a minute... looking at this, it appears you have to pass in an expression. What if it is an unknown function, where you only have a handle to the function, but you cannot see it's implementation ? Will this work then ? -----Original Message----- From: Berton Gunter

Hello, I am in need of help in putting histograms on the diagonal of a plot produced with splom(). The plot matrix I am trying to produce is to have standard scatterplots in the upper-left triangle, contour plots in the lower-right triangle, and histograms on the diagonal. I have a function that does the first two, but the histograms on the diagonal has been beyond my ability. Here is my

Thanks Bert for all the help. I got the legend figured out Friday but left early becoz of long weekend so didn't get a chance to reply.. I modified the plot margins a little bit and Here's what I finally had... par(mar=c(c(10, 6, 6, 10) + 0.1)); par(xpd=FALSE); with (dataFrame, stripchart(marbles_buried ~ genotype, method="jitter", vertical=TRUE, col = c('blue',

>On 10/14/2005 1:08 PM, Berton Gunter wrote: >> ?setwd >> e.g. setwd(file.choose()) >> >> BTW, you could have found this on your own via help.search('working >> directory') . Base R has quite good docs -- you should try them first. >> >> -- Bert Gunter > >That won't work in Windows, where the dialogs don't consider directories

Hello list, I am looking to create a figure for my dataset using splom, where there is a splom subplot for each level of a factor within, for example, a 2x2 layout. For each subplot, I wish to put the r-value between each variable pair in the lower panel. The code I have thus far is below, using the iris dataset so that it is reproducible. This almost does what I want, but for some reason

Greetings R experts, I'm having some difficulty recovering lda objects that I've saved within sublists using the save.image() function. I am running a script that exports a variety of different information as a list, included within that list is an lda object. I then take that list and create a list of that with all the different replications I've run. Unfortunately I've been

Hi, The following code produces a plot on X11: splom(~iris[1:4], groups = Species, data = iris, panel = panel.superpose, key = list(title = "Three Varieties of Iris", columns = 3, points = list(pch = super.sym$pch[1:3], col = super.sym$col[1:3]), text = list(c("Setosa", "Versicolor", "Virginica")))) However,

Hi, I'd like to be able to generate a splom plot in R and then use my mouse to click on one of the sub-panels (panel.pairs, specifically) and have R return either the coordinates of that sub-panel, or even better, the names of the corresponding variables plotted in that sub-panel. Here's an example to work with: > library(lattice) > splom(~iris[1:4], groups = Species, data = iris,

Hi, I have a small set of data on which I have tried some normality tests. When I make a histogram of the data the distribution doesn't seem to be normal at all (rather lognormal), but still no matter what test I use (Shapiro, Anderson-Darling,...) it returns a very small p value (which as far as I know means that the distribution is normal). Am I doing something wrong here? Thanks Pieter