Re: [R] Help making a simple R code only package

2004-05-03 Thread Prof Brian Ripley
Also, less incidentally, it would have worked in the current version of R, but it seems 1.8.0 was used -- an update of R would be a good idea. On Mon, 3 May 2004, Jason Turner wrote: Stephen Pavlik wrote: Hello, I am new to R and I am trying to figure out how to build and install a

Re: [R] Build problems on Linux SuSE 9.1

2004-05-03 Thread Prof Brian Ripley
It is the result of a bug in the XFree 4.4.0 headers. Please try the patched version r-patched (see the FAQ) where we have a workaround in src/modules/X11/dataentry.c (and quite a few other bug fixes). On Mon, 3 May 2004, Christian Lederer wrote: Hi, did anybody succeed in building R on

[R] multinomial regresion, nls

2004-05-03 Thread Mpiktas
Hi, Does R have any functions implementing such multinomial regression: (S_t^A,S_t^B)~MN(N_t-Y_{t-1},P_t^A,P_t^B) where MN(n,p_1,p_2) is multinomial distribution with parameters n, p_1, p_2. Here P_t^A and P_t^B are nonlinear functions from predictor variables and parameters which need to be

[R] Setting up contrasts

2004-05-03 Thread David J. Netherway
I am using the following model: lm - lm(mydata[[variableName]] ~ Age + Gender + Group, data=mydata) There are 5 groups in Group: nonc (the control), c1,c2,c3 and c4. How do I contrast nonc vs the others? and How do I contrast c1 vs other c's (ie c2,c3,c4 as a subgroup)? I have looked at the

RE: [R] PCA filtering

2004-05-03 Thread Philippe Grosjean
[EMAIL PROTECTED] wrote: Hi, Is there a function under R that does principal component analysis and then use the selected eignevectors to filter the data and transform back to the original coordinates? If there isn't, can someone guide as to how to do this approach. Much appreciated.

RE: [R] Setting up contrasts

2004-05-03 Thread BXC (Bendix Carstensen)
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of David J. Netherway Sent: Monday, May 03, 2004 10:25 AM To: [EMAIL PROTECTED] Subject: [R] Setting up contrasts I am using the following model: lm - lm(mydata[[variableName]] ~ Age + Gender +

[R] adding a method to the dist function

2004-05-03 Thread Giampiero Salvi
Hi all, I'd like to add the Bhattacharyya method to the dist function. What is the best way to do this? I'm using R 1.9.0 and I was looking for the code that defines the already existing distances, but I didn't manage. As far as I understand, dist is defined in mva that is part of the stats

Re: [R] adding a method to the dist function

2004-05-03 Thread Prof Brian Ripley
On Mon, 3 May 2004, Giampiero Salvi wrote: Hi all, I'd like to add the Bhattacharyya method to the dist function. What is the best way to do this? I'm using R 1.9.0 and I was looking for the code that defines the already existing distances, but I didn't manage. As far as I understand, dist

Re: [R] Least Squares Fit

2004-05-03 Thread Martin Maechler
Jon == Jonathan Baron [EMAIL PROTECTED] on Sun, 2 May 2004 07:10:10 -0400 writes: Jon On 05/02/04 03:26, Fred J. wrote: Hello I need to plot the least squares fit and get the slope of the line that best fit the data. after reading lm and lsfit, since not being

Re: [R] adding a method to the dist function

2004-05-03 Thread Martin Maechler
Giampiero == Giampiero Salvi [EMAIL PROTECTED] on Mon, 3 May 2004 11:50:38 +0200 (CEST) writes: Giampiero I'd like to add the Bhattacharyya method to the Giampiero dist function. What is the best way to do this? Giampiero I'm using R 1.9.0 and I was looking for the code

Re: [R] adding a method to the dist function

2004-05-03 Thread Giampiero Salvi
On Mon, 3 May 2004, Prof Brian Ripley wrote: On Mon, 3 May 2004, Giampiero Salvi wrote: Hi all, I'd like to add the Bhattacharyya method to the dist function. What is the best way to do this? I'm using R 1.9.0 and I was looking for the code that defines the already existing distances,

Re: [R] adding a method to the dist function

2004-05-03 Thread Prof Brian Ripley
On Mon, 3 May 2004, Giampiero Salvi wrote: On Mon, 3 May 2004, Prof Brian Ripley wrote: On Mon, 3 May 2004, Giampiero Salvi wrote: Hi all, I'd like to add the Bhattacharyya method to the dist function. What is the best way to do this? I'm using R 1.9.0 and I was looking for the

Re: [R] adding a method to the dist function

2004-05-03 Thread Giampiero Salvi
On Mon, 3 May 2004, Prof Brian Ripley wrote: dist() compares pairs of rows in the x matrix. How can they have `means and covariances'? -- you have a sample of size one from each of two populations. It seems that (Gaussian) Bhattacharyya is more like mahalanobis(). I had planned to use mean

[R] Multinomial Logistic Models

2004-05-03 Thread Cristian Pattaro
Hi all! Are there any *R* command for fitting multinomial logistic regression models? Thanks! Cristian = Cristian Pattaro = Unit of Epidemiology Medical Statistics Department of Medicines and Public Health

Re: [R] adding a method to the dist function

2004-05-03 Thread Prof Brian Ripley
On Mon, 3 May 2004, Giampiero Salvi wrote: On Mon, 3 May 2004, Prof Brian Ripley wrote: dist() compares pairs of rows in the x matrix. How can they have `means and covariances'? -- you have a sample of size one from each of two populations. It seems that (Gaussian) Bhattacharyya is

Re: [R] adding a method to the dist function

2004-05-03 Thread Martin Maechler
Giampiero == Giampiero Salvi [EMAIL PROTECTED] on Mon, 3 May 2004 14:01:23 +0200 (CEST) writes: ... Giampiero I agree that it would be better to use a new Giampiero function (similar to mahalanobis), as the Giampiero function dist in all the other cases

Re: [R] adding a method to the dist function

2004-05-03 Thread Giampiero Salvi
On Mon, 3 May 2004, Prof Brian Ripley wrote: On Mon, 3 May 2004, Giampiero Salvi wrote: On Mon, 3 May 2004, Prof Brian Ripley wrote: dist() compares pairs of rows in the x matrix. How can they have `means and covariances'? -- you have a sample of size one from each of two

Re: [R] Multinomial Logistic Models

2004-05-03 Thread Prof Brian Ripley
On Mon, 3 May 2004, Cristian Pattaro wrote: Are there any *R* command for fitting multinomial logistic regression models? Yes. help.search(multinomial) will lead you there (if I understand what you mean by the term `multinomial logistic' which strictly a contradiction). -- Brian D. Ripley,

[R] R 1.9.0 on AIX, 64-bit

2004-05-03 Thread Andy Pierce
I'm trying to get R 1.9.0 running on AIX 5.1 with the standard AIX compilers (xlc, xlf) and it is failing 2 of the tests, test-Reg in reg-tests-1.R like this: bash-2.05b$ tail -30 reg-tests-1.Rout.fail [,1] [,2] [1,] 13 [2,] 24 [3,] 13 [4,] 24 stopifnot(typeof(res) == list)

Re: [R] R 1.9.0 on AIX, 64-bit

2004-05-03 Thread Prof Brian Ripley
The second (at least) looks like a broken strptime. Try undefining HAVE_WORKING_STRPTIME in src/include/config.h and recompiling. Unfortunately, it could also be a broken mktime. On Mon, 3 May 2004, Andy Pierce wrote: I'm trying to get R 1.9.0 running on AIX 5.1 with the standard AIX

[R] Lattice: finding out xlim within panel function

2004-05-03 Thread Lutz Prechelt
Dear Lattice bit-meddlers, while within a panel function for xyplot, how can I find out the values of (effectively) xlim and ylim -- no matter whether they have been set explicitly or chosen by Lattice itself? I have just tried for an hour to find out, with no success whatsoever. I looked in

[R] ./configure on SunOS does not run to completion

2004-05-03 Thread wolski
OS: SunOS human 5.9 Generic_112233-12 sun4u sparc SUNW,Sun-Fire-480R I run ./configure --prefix=/non/standard/path ... and it stops at the line checking for xmkmf... /usr/openwin/bin/xmkmf Sincerely Eryk __ [EMAIL PROTECTED] mailing list

Re: [R] Lattice: finding out xlim within panel function

2004-05-03 Thread Deepayan Sarkar
On Monday 03 May 2004 09:36, Lutz Prechelt wrote: Dear Lattice bit-meddlers, while within a panel function for xyplot, how can I find out the values of (effectively) xlim and ylim -- no matter whether they have been set explicitly or chosen by Lattice itself? I have just tried for an hour

Re: [R] numericDeriv

2004-05-03 Thread Ravi Varadhan
Hi All: Along the lines of this thread, I was wondering about the usefulness of putting together a package for numerical differentiation, to perform tasks such as gradient, jacobian, and hessian calculations for exact functions, as well as for noisy functions (via some type of smoothing).

[R] Speed up graphics output?

2004-05-03 Thread dsheuman
Hi all, I've written some code to generate 4 maps per screen and write the output to a jpeg. The output is fairly quick at the start (about 5 jpegs per minute) and then slows down greatly (1-2 jpegs per minute). Is there some way to speed it up? One of my thoughts is to keep the base map

Re: [R] r dev site is down

2004-05-03 Thread Duncan Murdoch
On Sun, 2 May 2004 09:52:16 -0400, zubin [EMAIL PROTECTED] wrote : the R developer site has been down for some time, what gives? http://developer.r-project.org/ It's back online now at a different site. There are a few problems that will be worked out over the next few days: - DNS isn't

[R] RODBC/odbcConnect and Xemacs/ess Problem

2004-05-03 Thread witek
Hi Folks I am trying to use Xemacs/ess 5.2 on windows XP and have a problem with the RODBC package in R1.9.0 library(RODBC) channel-odbcConnect(mysql) Warning messages: 1: [RODBC] ERROR: Could not SQLDriverConnect 2: ODBC connection failed in: odbcDriverConnect(paste(DSN=, dsn, ;UID=, uid,

Re: [R] r dev site is down

2004-05-03 Thread Peter Dalgaard
Duncan Murdoch [EMAIL PROTECTED] writes: On Sun, 2 May 2004 09:52:16 -0400, zubin [EMAIL PROTECTED] wrote : the R developer site has been down for some time, what gives? http://developer.r-project.org/ It's back online now at a different site. There are a few problems that will be

[R] circular correlation

2004-05-03 Thread Joseph Dauer
I have a problem that deals with correlating wind velocity to seed collection data. The problem lies in that I have a wind data set that is contains 2000+ data points and weed collection data on the order of a couple hundred. Both data sets were collected for the same time period, but there is

[R] seek(..., origin='current') doesn't work on gzfile

2004-05-03 Thread Vadim Ogranovich
HI, It seems that seek(..., origin='current') doesn't properly work on gzfile connections, while it does work properly on ordinary files. # open a connection and repeatedly increment current position by 13 chars. Note that the read position does not move con - gzfile(foo.gz); open(con)

Re: [R] ./configure on SunOS does not run to completion

2004-05-03 Thread Don MacQueen
Well, that's interesting, but mine does run to completion. configure:28052: checking for xmkmf configure:28070: found /usr/openwin/bin/xmkmf configure:28082: result: /usr/openwin/bin/xmkmf (SunOS 5.8, but I have that file in that location also on a 5.9 system) I think you need to look at *your*

[R] selecting a coordinate from a graph

2004-05-03 Thread Richard Bonneau
hello everyone, What is the method for selecting a coordinate from a graph with the mouse, if any such function exists? For example, we would like to select outliers in a scatter-plot with the mouse, or add figure legends with the mouse. There was a way to add figure legends with the mouse in

RE: [R] openMosix vs SNOW: redhat kernel causing slowdown?

2004-05-03 Thread Samuelson, Frank*
You're using sockets for connection in snow? or pvm or mpi? There's nothing magical about snow. It just uses the socket connections provided in R, which in turn uses regular BSD sockets. -Original Message- From: Jim Thomas [mailto:[EMAIL PROTECTED] Sent: Thursday, April 29, 2004 6:13

Re: [R] selecting a coordinate from a graph

2004-05-03 Thread Prof Brian Ripley
On Mon, 3 May 2004, Richard Bonneau wrote: What is the method for selecting a coordinate from a graph with the mouse, if any such function exists? For example, we would like to select outliers in a scatter-plot with the mouse, or add figure legends with the mouse. There was a way to add

Re: [R] Speed up graphics output?

2004-05-03 Thread Paul Murrell
Hi Try using par(mfrow=c(2, 2)) instead of split.screen(c(2, 2)) and removing all of the screen(?) commands. I think what is happening is that you are just overwriting previous output every time through the loop so your picture just gets bigger and bigger (more and more overlayed output). With

[R] Legend outside the plotting region

2004-05-03 Thread Miha STAUT
Hello, How would I add a legend to the plot outside the plotting region? I tried different graphical parameters (fig, plt, usr and fig in combination with plt) in par() without success. Thanks in advance, Miha Staut __ [EMAIL PROTECTED] mailing list

Re: [R] Getting parameters from within function

2004-05-03 Thread Jason Turner
[EMAIL PROTECTED] wrote: Hi R-Helpers Is there a way to access the parameters passed to a function all at once? I want something like: misc.params = list(...) except such that it works with formal named parameters. Does match.call() do what you want? foo -

[R] Factor loadings and principal component plots

2004-05-03 Thread ndesnoyers
Hi- Can anyone tell me the command(s) to produce the following plots: -Factor loadings plot for principal components -Plot of principal component scores Also, apart from the prcomp (or princomp) command is there any other way to obtain principal components and if so, how does it/they stack up

Re: [R] Getting parameters from within function (redux)

2004-05-03 Thread Thomas Lumley
On Mon, 3 May 2004 [EMAIL PROTECTED] wrote: Hi R-Helpers Is there a way to access the parameters passed to a function all at once? I want something like: misc.params = list(...) except such that it works with formal named parameters. match.call() or sys.call() -thomas

[R] Getting parameters from within function

2004-05-03 Thread user
Hi R-Helpers Is there a way to access the parameters passed to a function all at once? I want something like: misc.params = list(...) except such that it works with formal named parameters. Context--I am running lots of simulations, want to record how they were specified. Thanks!

RE: [R] Legend outside the plotting region

2004-05-03 Thread John Fox
Dear Miha, Try setting par(xpd=TRUE) or par(xpd=NA). I hope this helps, John -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Miha STAUT Sent: Monday, May 03, 2004 4:15 PM To: [EMAIL PROTECTED] Subject: [R] Legend outside the plotting region

[R] Getting parameters from within function (redux)

2004-05-03 Thread wwsprague
(Sorry about the wack name--don't want to go into your spam buckets...) Hi R-Helpers Is there a way to access the parameters passed to a function all at once? I want something like: misc.params = list(...) except such that it works with formal named parameters. Context--I am running lots of

Re: [R] plotting with R

2004-05-03 Thread Frank E Harrell Jr
On Mon, 3 May 2004 19:14:26 -0400 Liaw, Andy [EMAIL PROTECTED] wrote: From: Karla Meurk (2) If I have three curves on a plot how can I label each curve so the label sits directly above the curve? [A]As to your second question: I don't think there will be a simple answer

[R] I have some problem about save and open

2004-05-03 Thread mathinee
Hello R-help I use R program on linux. and when I would like to quick R, then I type q() after than has question Save workspace image? [y/n/c]: when I ans y , why not to ask about name of file. and when I would like to open last file , how can I do? thanks for your help. and sorry if my English

[R] Epidemiology Tools

2004-05-03 Thread Francis Dimzon
Hi all, Please help on this. We will be teaching epidemiology using opensource software. What are R built-in functions or functions in available packages that are capable of doing these: a) Logistic regression (glm?) b) Conditional logistic regression c) Logistic regression with random

RE: [R] plotting with R

2004-05-03 Thread Liaw, Andy
From: Karla Meurk (2) If I have three curves on a plot how can I label each curve so the label sits directly above the curve? [A]As to your second question: I don't think there will be a simple answer because on has to decide where above is. Could you please be more specific?

Re: [R] circular correlation

2004-05-03 Thread Don MacQueen
I think you're going to have to decide which wind speed values are relevant to which seed collection data points. For example, if you collected a seed on Tuesday, is Monday's wind speed relevant or not? You've given no indication of what the relationships between your wind speed values and

RE: [R] plotting with R

2004-05-03 Thread John Fox
Dear Carla, I often position text on a graph interactively with the mouse via a call to locator(): e.g., text(locator(3), c(one, two, three)). I hope that this helps, John -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Liaw, Andy Sent: Monday,

[R] cointegration

2004-05-03 Thread Erin Hodgess
Dear R People: Is there a function for cointegration in any of the libraries, please? for R Windows, please? thanks, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: [EMAIL PROTECTED]

Re: [R] Lattice: finding out xlim within panel function

2004-05-03 Thread Paul Murrell
Hi Deepayan Sarkar wrote: On Monday 03 May 2004 09:36, Lutz Prechelt wrote: Dear Lattice bit-meddlers, while within a panel function for xyplot, how can I find out the values of (effectively) xlim and ylim -- no matter whether they have been set explicitly or chosen by Lattice itself? I have just