[R] installing knncat package

2004-11-15 Thread Jenny Edmondson
Hi all, I am wondering if anyone has problem installing 'knncat' package. I tried to install using 'Rcmd INSTALL knncat_1.1.1.tar.gz', error message appeared. I attached the message in the end. I allso tried to install from RGui, but it was not listed. Thanks in advance. Regards, jenny

Re: [R] installing knncat package

2004-11-15 Thread Prof Brian Ripley
This is an error in the package. In file src/Makevars it has LAPACK_LIBS=$(R_HOME)/src/modules/lapack/libRlapack.a which is a reference to the R sources. That should not be there: please remove it and this may well work. However, your version of R is rather old, which is why you are not seeing

[R] argument rationalization

2004-11-15 Thread Robin Hankin
Hi I am writing a bunch of functions that take two, three or four arguments. These functions operate on vectors of the same length; but I want the function to behave sensibly if one or more arguments are scalars. + does this for two arguments: +(1:10,3)# interpreted as +(1:10,rep(3,10))

Re: [R] argument rationalization

2004-11-15 Thread Prof Brian Ripley
I think you should implement recycling, ideally at C level. But you could have f - function(x, y, z) { n - max(length(x), length(y), length(z)) .C(something, as.double(rep(x, len=n)), as.double(rep(y, len=n)), as.double(rep(z, len=n)), as.integer(n), ans)$ans } On Mon, 15 Nov 2004,

[R] Problems installing packages on MacOS with R 2.00

2004-11-15 Thread Marco Chiarandini
Dear all, I have a problem installing a package required by Hmisc on MacOS 10.3.5 with R 2.00. g77 -fno-common -g -O2 -c avas.f -o avas.o g77 -fno-common -g -O2 -c rlsmo.f -o rlsmo.o gcc -bundle -flat_namespace -undefined suppress -L/usr/local/lib -o acepack.so ace.o avas.o rlsmo.o

Re: [R] argument rationalization

2004-11-15 Thread Robin Hankin
I think you should implement recycling, ideally at C level. But you could have f - function(x, y, z) { n - max(length(x), length(y), length(z)) .C(something, as.double(rep(x, len=n)), as.double(rep(y, len=n)), as.double(rep(z, len=n)), as.integer(n), ans)$ans } yes!

Re: [R] Problems installing packages on MacOS with R 2.00

2004-11-15 Thread Prof Brian Ripley
It is R 2.0.0! Your problem is that you do not have g77 installed, or at least, not the same version as was used to compile your version of R. (Please do read the posting guide and tell us where you got R from -- I suspect you did not compile it yourself.) I think this should be in the MacOS X

Re: [R] Where has the Debian respository gone?

2004-11-15 Thread Christoph Bier
Dirk Eddelbuettel schrieb am 15.11.2004 02:35 [CRAN Debian respository] It has been turned off by the CRAN masters as the content had slipped further and further behind the Debian content. Current R and CRAN packages are on the Debian archives; you can install these on testing too. To

Re: [R] Where has the Debian respository gone?

2004-11-15 Thread Christoph Bier
Dirk Eddelbuettel schrieb am 15.11.2004 04:12 On Sun, Nov 14, 2004 at 07:35:27PM -0600, Dirk Eddelbuettel wrote: [...] Current R and CRAN packages are on the Debian archives; you can install these on testing too. To the best of my knowledge, there are no backports of current R and Debian CRAN

Re: [R] Problems installing packages on MacOS with R 2.00

2004-11-15 Thread Marco Chiarandini
Dear Prof. Ripley, It is R 2.0.0! Your problem is that you do not have g77 installed, or at least, not the same version as was used to compile your version of R. (Please do read the posting guide and tell us where you got R from -- I suspect you did not compile it yourself.) I took R from

Re: [R] Legend help needed

2004-11-15 Thread Adaikalavan Ramasamy
You have not called legend() in your codes below, so we do not know what your problem is. See other comments below. On Mon, 2004-11-15 at 01:08, Sean David Richards wrote: R : Version 1.9.1 Hi, Am having trouble adding a legend to scatterplot. R code is shown below. I have tried various

Re: [R] Problems installing packages on MacOS with R 2.00

2004-11-15 Thread Prof Brian Ripley
On Mon, 15 Nov 2004, Marco Chiarandini wrote: Dear Prof. Ripley, It is R 2.0.0! Your problem is that you do not have g77 installed, or at least, not the same version as was used to compile your version of R. (Please do read the posting guide and tell us where you got R from -- I suspect you did

Re: [R] Power curves

2004-11-15 Thread Adaikalavan Ramasamy
What do you mean by power curves ? Is it the power of a study as the effect size varies or power output of a machine with some other parameter ? I usually generate a sequence of numbers (for the x-axis) that spans the range of interest and calculate its output. # Example 1 f - function(x)

Re: [R] Legend help needed

2004-11-15 Thread Adaikalavan Ramasamy
Sorry typo. The last line should read legend(1500, 9000, legend=paste(Data from, sfiles), pch=1:n, col=1:n ) ^^^ On Mon, 2004-11-15 at 11:39, Adaikalavan Ramasamy wrote: You have not called legend() in your codes below, so we do

[R] Somewhat off-topic : Hand and Taylor reference

2004-11-15 Thread Rachel Pearce
I am sorry to trouble the list with this, but I'm hoping that someone will be able to help me track down this reference. I am looking for a good reference on repeated measures, and many of the relevant functions in R refer in their help to: Hand, D. J. and Taylor, C. C. (1987)

[R] how can draw probability density plot?

2004-11-15 Thread LONG Yu
Dear support, I want to draw a probability density plot in R. For example, I provide the mean and variance of a normal distribution, then R can provide me the probability density plot. Now I always generate random numbers of normal distribution and calculate their dnorm(mu, var), finally plot

[R] R: how can draw probability density plot?

2004-11-15 Thread Vito Ricci
I hope this example could help you best vito x-seq(-3.5,3.5,0.1) x [1] -3.5 -3.4 -3.3 -3.2 -3.1 -3.0 -2.9 -2.8 -2.7 -2.6 -2.5 -2.4 -2.3 -2.2 -2.1 [16] -2.0 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1.0 -0.9 -0.8 -0.7 -0.6 [31] -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6

Re: [R] Somewhat off-topic : Hand and Taylor reference

2004-11-15 Thread Peter Dalgaard
Rachel Pearce [EMAIL PROTECTED] writes: Hand, D. J. and Taylor, C. C. (1987) _Multivariate Analysis of Variance and Repeated Measures._ Chapman and Hall. This sounds like a useful text, but I can't find it on Amazon. I can find: Multivariate Analysis of Variance for

[R] R: how can draw probability density plot?

2004-11-15 Thread Vito Ricci
I believe this is better and that you want! x-rnorm(1000,10,1) plot( function(y) dnorm(y, mean(x), sd(x)), from=min(x), to=max(x)) best vito You wrote: Dear support, I want to draw a probability density plot in R. For example, I provide the mean and variance of a normal distribution, then R

Re: [R] Somewhat off-topic : Hand and Taylor reference

2004-11-15 Thread Prof Brian Ripley
The full library info is Multivariate analysis of variance and repeated measures : a practical approach for behavioural scientists Publisher London : Chapman and Hall, 1987. Description xiii, 262 p. : ill. ; 24 cm. Notes Bibliography: p. 208-209. - Includes index. ISBN0412258005

[R] Error when scan()ning dead URL's

2004-11-15 Thread Jonathan Croft
Hi, I often scan web pages directly into R using scan('http://etc...','') however this gives an error if the page/url doesn't exist, or the connection to it is not available. Is it possible to still use scan but have R return something other than an error (which crashes code) when the page or

Re: [R] Error when scan()ning dead URL's

2004-11-15 Thread Peter Dalgaard
Jonathan Croft [EMAIL PROTECTED] writes: I often scan web pages directly into R using scan('http://etc...','') however this gives an error if the page/url doesn't exist, or the connection to it is not available. Is it possible to still use scan but have R return something other than an

[R] help for nls

2004-11-15 Thread BARNERIAS Cyrille
Hello, I am beginning with R and I would like to test a non linear model. But I do not find exactly wath I am looking for in nls packages (or I do not know where to search). I would like to try a model like this : y=b * x exp(n)/(a exp(n) + x exp (n)) Where a = a0 + a1z b= b0 + b1z x and z are

[R] Using R in parallel on a 2 processor machine

2004-11-15 Thread Graham Law
Hi, I have installed R (2.0.0) onto a two processor machine running Windows XP (these two processors have been split into 4 logical processors), with 4Gb of RAM. Rather than multi-threading, I wanted to run 2 instances of R on the machine in parallel, which uses two of the logical processors.

RE: [R] Using R in parallel on a 2 processor machine

2004-11-15 Thread Andy Bunn
How about installing an operating system that knows its way around that much RAM? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Graham Law Sent: Monday, November 15, 2004 9:12 AM To: [EMAIL PROTECTED] Subject: [R] Using R in parallel on a 2 processor

Re: [R] Using R in parallel on a 2 processor machine

2004-11-15 Thread Prof Brian Ripley
Try reading the rw-FAQ! There is no way user processes in Windows XP can access all 4Gb of the address space, but if you tune both R and XP correctly you should be able to get above 2Gb in total. See the rw-FAQ for how to tune R. Nevertheless, as you have already been told, it is a lot easier

Re: [R] help for nls

2004-11-15 Thread Arne Henningsen
Dear Cyrille, type ?nls in the R command line. This will show you how to use the nls function. And please read the posting guide (see bootom of each message). This will help you to get more helpful answers. Best wishes, Arne On Monday 15 November 2004 14:42, BARNERIAS Cyrille wrote: Hello,

[R] R-2.0.1 is released

2004-11-15 Thread Peter Dalgaard
I've rolled up R-2.0.1.tgz a short while ago. This is a maintenance version mainly to fix a number of minor bugs and issues. (Some rather nasty ones were Windows-specific - please check the CHANGES file for those.) You can get it from http://cran.r-project.org/src/base/R-2.0.1.tar.gz or wait

[R] Error whilst building packages

2004-11-15 Thread Dan Bailey
Dear All, I have been working on building a new version of the Wavethresh package for some time now. Having build a working version on Linux, I am getting the following error when checking on Windows: C:\Rpackages\R\rw2000\binRcmd check wavethresh * checking for working latex ... OK * using log

[R] Multivariate Sampling

2004-11-15 Thread stefan . albrecht
Dear all, I am looking for routines which allow multi-variate sampling from non-normal distributions (loglogistic) given correlations among the variables. Unfortunately, I could not find a suitable package for R. Does anybody know one? Many thanks and best regards, Stefan Albrecht

[R] tsdiag() titles

2004-11-15 Thread Andrew Kniss
I am using the ts package to fit ARIMA models, and the tsdiag() function to plot diagnostics. In doing so I'm generating an awful lot of diagnostic plots of different models and different data sets all within the same R session. So my question is, is there an option in tsdiag() similar to

[R] R-2.0.1 is released

2004-11-15 Thread Peter Dalgaard
Whoops! You can get it from http://cran.r-project.org/src/base/R-2/R-2.0.1.tar.gz (Notice the R-2 subdir) I also forgot to sign it on behalf of The R Core Team. -- O__ Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) --

RE: [R] tsdiag() titles

2004-11-15 Thread Andy Bunn
You can use title, but the result is unsatisfying: fit - arima(lh, c(1,0,0)) tsdiag(fit) title(junk) Perhaps mtext with an appropriate par configuration? HTH, Andy -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Andrew Kniss Sent: Monday,

Re: [R] Error whilst building packages

2004-11-15 Thread Prof Brian Ripley
On Mon, 15 Nov 2004, Dan Bailey wrote: Dear All, I have been working on building a new version of the Wavethresh package for some time now. Having build a working version on Linux, I am getting the following error when checking on Windows: C:\Rpackages\R\rw2000\binRcmd check wavethresh *

Re: [R] tsdiag() titles

2004-11-15 Thread Prof Brian Ripley
Did you notice that tsdiag() plots about three plots and gives each a title? You can add a title to the array of plots using title(outer=TRUE), surprisingly enough, but you will need to adjust the outer margins to make room for it. Something like par(oma=c(0,0,2,0)) fit - arima(lh, c(1,0,0))

[R] how to obtain predicted labels for test data using kernelpls

2004-11-15 Thread Huh, Seungho
Dear members, My name is Seungho Huh. I am a statistician who tries to use the Kernel PLS method in a classification problem. I am sending this email to ask you something about the kernelpls function in R (pls.pcr package). I would like to obtain the predicted Y values for test data, using

[R] R and fluxbox: 100% CPU usage

2004-11-15 Thread David Whiting
Hi, In July I reported a problem I was having with R and the fluxbox windows manager (with Linux). The interaction of R and fluxbox causes CPU to go to 100% when trying to create even a simple plot. It is the R process that is at 100%. R with any other window manager was fine, and all other

[R] glim in R?

2004-11-15 Thread Tim F Liao
After some futile searches, I decided to ask the list to see if any of the sages out there would have an answer: I have a function I wrote a few years ago in S, which calls glim numerous times. I'd like to port it to R, but glm works differently from glim, which takes as part of its input an

Re: [R] R and fluxbox: 100% CPU usage

2004-11-15 Thread Marc Schwartz
On Mon, 2004-11-15 at 17:17 +, David Whiting wrote: Hi, In July I reported a problem I was having with R and the fluxbox windows manager (with Linux). The interaction of R and fluxbox causes CPU to go to 100% when trying to create even a simple plot. It is the R process that is at

Re: [R] glim in R?

2004-11-15 Thread Sundar Dorai-Raj
Tim F Liao wrote: After some futile searches, I decided to ask the list to see if any of the sages out there would have an answer: I have a function I wrote a few years ago in S, which calls glim numerous times. I'd like to port it to R, but glm works differently from glim, which takes as

Re: [R] glim in R?

2004-11-15 Thread Rolf Turner
Tim Liao wrote: After some futile searches, I decided to ask the list to see if any of the sages out there would have an answer: I have a function I wrote a few years ago in S, which calls glim numerous times. I'd like to port it to R, but glm works differently from glim, which takes as

RE: [R] R on 64-bit Linux machine

2004-11-15 Thread Vadim Ogranovich
Thanks to everyone for the info. It is very valuable. I am a little bit uneasy about conflicting reports regarding RHEL 3, but I guess at this point I just need to try and see. It's also very soothing to know that there is an official 64-bit build on CRAN. Thanks again for taking time to answer,

Re: [R] Legend help needed

2004-11-15 Thread Sean David Richards
On 15 Nov 2004 at 12:11, Adaikalavan Ramasamy wrote: Have you considered points() or lines() here ? You could simplify to plot(0,1000, type=n, xlim=c(0,1500), ylim=c(1000,9000), xlab=Area above Notch (mm), main=Size Effect Speciments) n - length(sfiles) for (i in 1:n) {

[R] logging messages

2004-11-15 Thread Faber Fedor
Hi all, I'm trying to learn R from an non-statistician's POV. I've got a statistician who uses R, but I'm the schmuck who has to integrate his R functions into an automated process. One of the things I would really like is the ability to log messages to file, specifically using syslog on a Linux

[R] anova, multiple comparison

2004-11-15 Thread Yuandan Zhang
hi, I try to used R to do one-way anova. here is the simple code f1- lm (y ~ block, data=yd) there are 8 levels of factor block, I also want produce multiple pairwise comparisons for the 8 levels of block, inlcuding mean and std err for each of the 8 lelevls. It is tidious to do pair

[R] Barplot difficulties

2004-11-15 Thread Heather J. Branton
Hello. I am an R newbie struggling to learn and use R . I have read many portions of the R Reference Manual, as well as the FAQs. Given that I learn something new each time, I know I might be missing something obvious. But I appeal to your good nature to help me through this initial problem.

[R] lme, two random effects, poisson distribution

2004-11-15 Thread Martina Pavlicova, PhD
Hello, I have a dataset concerning slugs. For each slug, the number of pumps per one time slot was counted. The number of pumps follows Bi(30, p) where p is very small, thus could be approximated by Poisson dist. (# of pumps is very often = 0) The slugs were observed during 12 time slots which

Re: [R] glim in R?

2004-11-15 Thread Tim F Liao
Great suggestions, and it looks like either suggestion should work, although the output from glm may not conform to those from glim, thus some more code there perhaps. Many thanks, Tim Original message Date: Mon, 15 Nov 2004 14:28:50 -0400 (AST) From: Rolf Turner [EMAIL PROTECTED]

Re: [R] Barplot difficulties

2004-11-15 Thread Marc Schwartz
On Mon, 2004-11-15 at 19:03 -0500, Heather J. Branton wrote: Hello. I am an R newbie struggling to learn and use R . I have read many portions of the R Reference Manual, as well as the FAQs. Given that I learn something new each time, I know I might be missing something obvious. But I

RE: [R] how to obtain predicted labels for test data using kerne lpls

2004-11-15 Thread Liaw, Andy
You need to do some extra work if you want to do classification with a regression method. One simple way to do classification with PLS is to code the classes as 0s and 1s (assuming there are only two classes) or -1s and 1s, fit the model, then threshold the prediction; e.g., those with predicted

[R] Re: lme, two random effects, poisson distribution

2004-11-15 Thread Martina Pavlicova, PhD
Hello all, I think that with the help of Mark Irwin, we solved the problem: We fit the model using glmmPQL and instead of using variable, 'state', we model the independent fixed-effect 'state' as I(x10); i.e. it's 0 for resting time slots and 1 for excited times slots. Here is the code:

Re: [R] anova, multiple comparison

2004-11-15 Thread Prof Brian Ripley
On Tue, 16 Nov 2004, Yuandan Zhang wrote: I try to used R to do one-way anova. here is the simple code f1- lm (y ~ block, data=yd) there are 8 levels of factor block, I also want produce multiple pairwise comparisons for the 8 levels of block, inlcuding mean and std err for each of the 8 lelevls.