[R] Cross build Makefile

2004-06-21 Thread Iago Mosqueira
Hello, I am trying to use Yan and Rossini's Makefile for cross building Windows versions of R packages in Linux with R 1.9.0. When compiling R with the mingw tools I get an error about expm1 being undeclared when first found at src/main/arithmetic.c:1019 If I fiddle a bit with it later on I also

Re: [R] Cross build Makefile

2004-06-21 Thread Prof Brian Ripley
On Mon, 21 Jun 2004, Iago Mosqueira wrote: Hello, I am trying to use Yan and Rossini's Makefile for cross building Windows versions of R packages in Linux with R 1.9.0. When compiling R with the mingw tools I get an error about expm1 being undeclared when first found at

Re: [R] Sweave and echoing R comments

2004-06-21 Thread Friedrich . Leisch
On Sun, 20 Jun 2004 20:34:50 +, Frank E Harrell (FEH) wrote: Charles H. Franklin wrote: Is there any way to echo comments from an R source file into an SWeave- LaTeX document? Example: # Npop is population total # Npoph0..Npoph2 are stratum totals # Npoph is

[R] sunrise, sunset calculation

2004-06-21 Thread Angel Lopez
Are there any functions available to calculate sunrise and sunset times for given latitude,longitude and dates? If not, I'll appreciatte any pointers to C code I could use/port. Thanks, Angel __ [EMAIL PROTECTED] mailing list

Re: [R] sunrise, sunset calculation

2004-06-21 Thread Mike Sumner
At 06:26 PM 6/21/2004, Angel Lopez wrote: Are there any functions available to calculate sunrise and sunset times for given latitude,longitude and dates? If not, I'll appreciatte any pointers to C code I could use/port. Thanks, Angel Hello, I have R code for determining solar elevation for given

Re: [R] sunrise, sunset calculation

2004-06-21 Thread Martin Maechler
Mike == Mike Sumner [EMAIL PROTECTED] on Mon, 21 Jun 2004 19:04:00 +1000 writes: Mike At 06:26 PM 6/21/2004, Angel Lopez wrote: Are there any functions available to calculate sunrise and sunset times for given latitude,longitude and dates? If not, I'll appreciatte any

[R] modifying a glm object

2004-06-21 Thread Adrian Baddeley
I need to modify a glm object directly, by assigning new values to the fitted parameters. (i.e. the canonical parameters stored in object$coefficients and returned by 'coef'). Is there a safe way to do this, e.g. so that 'residuals' and 'predict' give the correct results for the new model?

Re: [R] sunrise, sunset calculation

2004-06-21 Thread Angel Lopez
Martin Maechler wrote: Mike == Mike Sumner [EMAIL PROTECTED] on Mon, 21 Jun 2004 19:04:00 +1000 writes: Mike At 06:26 PM 6/21/2004, Angel Lopez wrote: Are there any functions available to calculate sunrise and sunset times for given latitude,longitude and dates? If

[R] R 1.9.1 is released

2004-06-21 Thread Peter Dalgaard
I've rolled up R-1.9.1.tgz a short while ago. This is a maintenance version mainly to fix a number of minor bugs and issues (the most annoying one seems to have been the change to barplots of tables) and some installation issues. Because of the relocation the CVS archives, there is no direct

RE: [R] [Q] Newbie (continued.. at least I got R running allready :-)

2004-06-21 Thread Petr Pikal
Hi On 18 Jun 2004 at 8:34, Gabor Grothendieck wrote: Check out ?turnpoints in the pastecs library. Try example(turnpoints) str(Nauplii.tp) to see the data structure you get back. It identifies the peaks and pits and other associated information. I am not sure if this

[R] reshape speed improvement?

2004-06-21 Thread Christian Schulz
Hi, have anybody a hint how i can perform my speed for a big reshape jobs? My long data.frame have dim 1.342.1163and the wide data.frame could occur (if i don't filter the timevar) max. in a 246744 x 1444 data.frame with a lot of NA what could be set to 0. Neverthless i didn't

[R] frequencies

2004-06-21 Thread Silvia Kirkman
Hi I have a set of fish lengths (cm) which I'd like to have divided into bins as specified by myself. I want to classify my bins as: 0=x0.5 0.5=x1 1=x1.5 1.5=x2 and so on... I'd like the frequencies to be given as a vector because I need to use these values for further analysis. Please can

RE: [R] frequencies

2004-06-21 Thread Wayne Jones
See ?hist For example hist(rnorm(100),plot=F,breaks=-3:3)$counts Regards Wayne -Original Message- From: Silvia Kirkman [mailto:[EMAIL PROTECTED] Sent: 21 June 2004 14:01 To: [EMAIL PROTECTED] Subject: [R] frequencies Hi I have a set of fish lengths (cm) which I'd like to have

RE: [R] frequencies

2004-06-21 Thread BXC (Bendix Carstensen)
You want the function cut(), followed by table(). Bendix -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED] www.biostat.ku.dk/~bxc

Re: [R] frequencies

2004-06-21 Thread Petr Pikal
On 21 Jun 2004 at 6:01, Silvia Kirkman wrote: Hi I have a set of fish lengths (cm) which I'd like to have divided into bins as specified by myself. I want to classify my bins as: 0=x0.5 0.5=x1 1=x1.5 1.5=x2 and so on... Hallo ?cut and ?table should be what you are looking for

[R] Visual stimulus presentation using R?

2004-06-21 Thread Christoph Lange
Dear all! Although the Psycho-Toolbox for Matlab is free software, Matlab isn't. I'm planning to do an experiment where it's essentail to travel to the subjects, not let the subjects come to where the Matlab licences are :-( So I need to use a free software for my experiment if I don't want to

Re: [R] frequencies

2004-06-21 Thread Adaikalavan Ramasamy
See the example in help(cut). You will require the option right=FALSE in cut() or you can try hist. x - abs(rnorm(100)) table( cut(x, seq(0, max(ceiling(x)), by=0.5), right=FALSE )) hist(x, breaks=seq(0, max(ceiling(x)), by=0.5), plot=FALSE) On Mon, 2004-06-21 at 14:01, Silvia Kirkman wrote:

Re: [R] Visual stimulus presentation using R?

2004-06-21 Thread Ioannis Dimakos
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Christoph Lange wrote: Have you considered DMDX from the Univ. of Arizona http://www.u.arizona.edu/~kforster/dmdx/overview.htm It is a rather capable program, runs on Windows platforms, and may help you with your situation. Ioannis | Dear all! | |

Re: [R] Visual stimulus presentation using R?

2004-06-21 Thread Spencer Graves
Have you considered open source alternatives to Matlab such as Scilab? I know nothing about their capabilities, but Google produced 2410 hits for Matlab clones, the first of which (http://www.dspguru.com/sw/opendsp/mathclo2.htm) seemed quite useful. hope this helps. spencer graves

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Shin, Daehyok
Glad to know useful functions. How about adding lsv.str function to list only variables bound to values? In my opinion, we are more interested in values than functions in the process of data analysis. In addition, the simple solution of Grothendieck to display only names of objects has its own

RE: [R] binary data

2004-06-21 Thread Shin, Daehyok
I use NetCDF file format using ncdf package for binary data handling. For more complex structure, HDF format will be better, which I do not try in R. Daehyok Shin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ethz.ch]On Behalf Of Moises Hassan Sent: Tuesday,

Re: [R] Visual stimulus presentation using R?

2004-06-21 Thread Christoph Lehmann
Hi Christoph I have never done such stuff with R (and I don't know, e.g. how good the timing would be, in case you are interested in reaction times and so on) but have a look at www.visionegg.org it's a stimulus-presentation framework, written by Andrew Straw entirely in python and has lots

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Prof Brian Ripley
What is a `variable'? R has object,s some of which are function and all are variable or not depending where they are located. On Mon, 21 Jun 2004, Shin, Daehyok wrote: Glad to know useful functions. How about adding lsv.str function to list only variables bound to values? In my opinion, we

Re: [R] A way to list only variables or functions?

2004-06-21 Thread Duncan Murdoch
On Mon, 21 Jun 2004 09:53:35 -0400, Shin, Daehyok [EMAIL PROTECTED] wrote : Glad to know useful functions. How about adding lsv.str function to list only variables bound to values? In my opinion, we are more interested in values than functions in the process of data analysis. In R, functions

Re: [R] Another NEWBIE

2004-06-21 Thread Sven Hartenstein
Hi, F.Kalder [EMAIL PROTECTED] writes: The data entries by ASCII files are strange to me, because I?m so used to work with a (the SPSS) spread sheet (mostly the good old typing in from paper pencil questionnaires), that I don?t know how to handle that yet. Maybe using a SPSS- or at least

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Liaw, Andy
From: Duncan Murdoch On Mon, 21 Jun 2004 09:53:35 -0400, Shin, Daehyok [EMAIL PROTECTED] wrote : Glad to know useful functions. How about adding lsv.str function to list only variables bound to values? In my opinion, we are more interested in values than functions in the process of

Re: [R] A way to list only variables or functions?

2004-06-21 Thread Petr Pikal
On 21 Jun 2004 at 10:39, Duncan Murdoch wrote: On Mon, 21 Jun 2004 09:53:35 -0400, Shin, Daehyok [EMAIL PROTECTED] wrote : Glad to know useful functions. How about adding lsv.str function to list only variables bound to values? In my opinion, we are more interested in values than

RE: [R] problem locfit

2004-06-21 Thread Liaw, Andy
Please do read the posting guide, as the footer suggests, on how to ask questions in a way that induce useful answers. From the little bit that you showed, you stored the result of locfit() to an object named `fitbmt', but then try to look at an object named `fitbmt31'. Is that what you really

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Shin, Daehyok
I am not familar with the unique terminology of R. What I mean with a variable is a sort of a reference which points to a function or a value. Daehyok Shin -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Monday, June 21, 2004 AM 10:26 To: Shin, Daehyok Cc:

[R] Novice question: Smooth interpolation of survival curve

2004-06-21 Thread Mark Nye
Greetings, How might one go about generating a smooth interpolation of a survival curve generated by the survfit function in the survival package? I am able to package my variables by standard methods x-(survfit(Surv(time),data) and plot the step curve, but would like to obtain a smooth

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Shin, Daehyok
Thanks, Petr. The ls.objects function displays a really nice summary about objects. Surely, I will use it more often than ls or ls.str. Daehyok Shin -Original Message- From: Petr Pikal [mailto:[EMAIL PROTECTED] Sent: Monday, June 21, 2004 AM 11:01 To: Duncan Murdoch; [EMAIL

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Prof Brian Ripley
On Mon, 21 Jun 2004, Shin, Daehyok wrote: I am not familar with the unique terminology of R. And that's why we have documentation and hope that people will read it. What I mean with a variable is a sort of a reference which points to a function or a value. In that case, ls() does what you

[R] R 1.9.1 (source) available by rsync

2004-06-21 Thread Martin Maechler
Thanks to Doug Bates (for whom I'm substituting), The source for R-1.9.1 has been available within half an hour of Peter's R-announcement from rsync.R-project.org. You can, for example, use rsync -avC rsync.r-project.org::r-release ./R-1.9.1 to obtain a copy of the sources which [you won't

RE: [R] Visual stimulus presentation using R?

2004-06-21 Thread John Fox
Dear Christoph, Several people have suggested alternative software, and that may well be the way to go, but if you want to stick with R, then it should be possible to use Tcl/Tk via the tcltk package to do what you want. Whether the result would be sufficiently responsive for your needs is hard

[R] Novice question: Smooth interpolation of survival curve

2004-06-21 Thread Mark Nye
Greetings, How might one go about generating a smooth interpolation of a survival curve generated by the survfit function in the survival package? I am able to package my variables by standard methods x-(survfit(Surv(time),data) and plot the step curve, but would like to obtain a smooth

RE: [R] A way to list only variables or functions?

2004-06-21 Thread Gabor Grothendieck
In R, functions often contain useful information about data (in their attached environments). There are several variants of this posted to R-help and R-devel. As an example, search for ls.obj() or ls.object(). Check out ?browseEnv __ [EMAIL

[R] Contrib package downloads?

2004-06-21 Thread drf5n
Where can I download the contrib packages? http://lib.stat.cmu.edu/R/CRAN/src/contrib/PACKAGES.html points to a list of packages available for download. but the packages aren't there: http://lib.stat.cmu.edu/R/CRAN/src/contrib/Descriptions/Hmisc.html

Re: [R] Contrib package downloads?

2004-06-21 Thread drf5n
On Mon, 21 Jun 2004 [EMAIL PROTECTED] wrote: Where can I download the contrib packages? http://lib.stat.cmu.edu/R/CRAN/src/contrib/PACKAGES.html points to a list of packages available for download. but the packages aren't there:

Re: [R] A way to list only variables or functions?

2004-06-21 Thread Barry Rowlingson
Duncan Murdoch wrote: I seem to recall that S-PLUS has such a function, but I forget the name of it. Probably R does too, on CRAN if not in the base packages. objects.summary() I think it was. It always bothered me that the Nth thing you teach Unix users is 'ls' and the N+1th thing is 'ls -l'

[R] question about latex command in Hmisc

2004-06-21 Thread Laura Holt
Hi R People: I got the following error from using the latex command from Hmisc; latex(bbm) 'latex' is not recognized as an internal or external command, operable program or batch file. Warning message: cmd execution failed with error code 1 in: shell(cmd, wait = TRUE, intern = output) I was

RE: [R] question about latex command in Hmisc

2004-06-21 Thread kevin bartz
When you do latex from the command line (Run - cmd), what does DOS say? Maybe latex isn't in your PATH. Kevin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Laura Holt Sent: Monday, June 21, 2004 9:59 AM To: [EMAIL PROTECTED] Subject: [R] question about

[R] visualizing a list of probabilities

2004-06-21 Thread Rajarshi Guha
Hi, I'm using nnet to work on a 2 class classification problem. The result of my code is data.frame of true class, predicted class and associated probability. One way of summarizing the data is by a confusion matrix. However are there any graphical ways I could represent the data -

RE: [R] visualizing a list of probabilities

2004-06-21 Thread Liaw, Andy
Not exactly sure what you want, but you might want to look at the margin() function and the associated plot() method for margin objects in the randokmForest package. (You won't be able to use them directly on an nnet object, but the code should help.) Cheers, Andy From: Rajarshi Guha Hi,

[R] Elementary sapply question

2004-06-21 Thread Ajay Shah
I am discovering sapply! :-) Could you please help me with a very elementary question? Here is what I know. The following two programs generate the same answer. + Loops version| sapply version

Re: [R] question about latex command in Hmisc

2004-06-21 Thread Frank E Harrell Jr
Laura Holt wrote: Hi R People: I got the following error from using the latex command from Hmisc; latex(bbm) 'latex' is not recognized as an internal or external command, operable program or batch file. Warning message: cmd execution failed with error code 1 in: shell(cmd, wait = TRUE, intern =

RE: [R] Elementary sapply question

2004-06-21 Thread Liaw, Andy
At least two ways: 1. Use extra argument in the function being sapply()'ed; e.g., f - function(x, y) x*x + y*y x - 3:5 sapply(x, f, 3) [1] 18 25 34 [See the ... argument in ?sapply.] 2. More generally, if both x and y are vectors (of the same length), then you can use mapply(); e.g., x -

Re: [R] Elementary sapply question

2004-06-21 Thread Ajay Shah
I had asked: My problem is this. Suppose I have: pythagorean - function(x, y) { return(x*x + y*y) } then how do I utilise sapply to replace fixed.x = 3 y.values = c(3,4,5) answers=numeric(3) for (i in 1:3) {

RE: [R] Elementary sapply question

2004-06-21 Thread Brian Desany
Looking in ?mapply, I executed the examples: mapply(rep, 1:4, 4:1) [[1]] [1] 1 1 1 1 [[2]] [1] 2 2 2 [[3]] [1] 3 3 [[4]] [1] 4 mapply(rep, times=1:4, x=4:1) [[1]] [1] 4 [[2]] [1] 3 3 [[3]] [1] 2 2 2 [[4]] [1] 1 1 1 1 I can guess that because these are 2 examples, it is no surprise

Re: [R] Elementary sapply question

2004-06-21 Thread Prof Brian Ripley
You really ought to name ... arguments. sapply(x, f, y=3) makes it clear that f(xx, y=3) is called. But `optional arguments' necessarily come after compulsory ones, which resolves the ambiguity you see. Please note that 1) functions return their values 2) a function body is an expression, so

Re: [R] Elementary sapply question

2004-06-21 Thread Sundar Dorai-Raj
Brian Desany wrote: Looking in ?mapply, I executed the examples: mapply(rep, 1:4, 4:1) [[1]] [1] 1 1 1 1 [[2]] [1] 2 2 2 [[3]] [1] 3 3 [[4]] [1] 4 mapply(rep, times=1:4, x=4:1) [[1]] [1] 4 [[2]] [1] 3 3 [[3]] [1] 2 2 2 [[4]] [1] 1 1 1 1 I can guess that because these are 2 examples, it is no

RE: [R] Elementary sapply question

2004-06-21 Thread Prof Brian Ripley
Named arguments are matched by name. Unnamed arguments are matched by position. S's argument matching rules are quite complex. It normally helps clarity to name arguments in all but the simplest calls. rep(1, 4) and rep(times=1, x=4) are different calls -- the second could be called

Re: [R] Elementary sapply question

2004-06-21 Thread Tony Plate
At Monday 12:57 PM 6/21/2004, Ajay Shah wrote: [...snip...] I am aware of the ... in sapply(). I am unable to understand how sapply will know where to utilise the x[i] values: as the 1st arg or the 2nd arg for f(x, y)? That is, when I say: sapply(x, f, 3) how does sapply know that I mean:

RE: [R] Elementary sapply question

2004-06-21 Thread Liaw, Andy
If you have something like lapply(x, f, ...) What lapply() does (which is the same as what sapply() does, except sapply() tries to `simplify' the result) is roughly: result - vector(mode=list, length=length(x)) for (i in seq(along(x)) { result[[i]] - f(x[i], ...) } I.e., it takes the first

[R] Welch-JM-Test or Brown-Forsythe-Test in R?

2004-06-21 Thread Sven Hartenstein
Hi, I want to test mean differences of 2 groups with heterogenous variances and wonder whether Welch-JM-Test and/or Brown-Forsythe-Test are available in R. (These two are the ones I found in the literature and SPSS provides, I'm open for alternatives.) Any help apreciated! Sven

RE: [R] Welch-JM-Test or Brown-Forsythe-Test in R?

2004-06-21 Thread Wiener, Matthew
Does oneway.test do what you want? Hope this helps, Matt Wiener -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Sven Hartenstein Sent: Monday, June 21, 2004 3:23 PM To: [EMAIL PROTECTED] Subject: [R] Welch-JM-Test or Brown-Forsythe-Test in R? Hi, I

[R] setwd problems

2004-06-21 Thread Joan Serra
Hello, I am using R for Windows and I receive error messages when trying to change my working directory: setwd('C:\BACC_R') Error in setwd(dir) : cannot change working directory I would really appreciate your help, Joan Serra __ [EMAIL PROTECTED]

RE: [R] setwd problems

2004-06-21 Thread Ko-Kang Kevin Wang
Hi, -Original Message- From: [EMAIL PROTECTED] Hello, I am using R for Windows and I receive error messages when trying to change my working directory: setwd('C:\BACC_R') Error in setwd(dir) : cannot change working directory Have you tried: setwd(C:/BACC_R) Or

Re: [R] Welch-JM-Test or Brown-Forsythe-Test in R?

2004-06-21 Thread Sven Hartenstein
Hi, * Wiener, Matthew [EMAIL PROTECTED] writes: Does oneway.test do what you want? Yes, guess so. Thank you! BTW, how could I have found this without asking? I only found lm(), aov() when googling, and in the help() pages of R on these functions there's no reference to oneway.test... Thanks

Re: [R] setwd problems

2004-06-21 Thread Sundar Dorai-Raj
Joan Serra wrote: Hello, I am using R for Windows and I receive error messages when trying to change my working directory: setwd('C:\BACC_R') Error in setwd(dir) : cannot change working directory I would really appreciate your help, Joan Serra Hi Joan, Please review the following:

[R] Bookmarklet for searching R documentation

2004-06-21 Thread Kevin Wright
I recently posted a request for a Mozilla search engine plugin for the R java search applet. Having recieved no response, I pursued this myself and came up with an alternative. I found a bookmarklet that was used to submit searches to Amazon. I modified the code to submit searches to the R

[R] nonparametric two-way structure all-pairwise comparisons

2004-06-21 Thread Marco Chiarandini
Hello, I am conducting a two-way analysis of variance. The ANOVA assumption are not met, hence I need to use non-parametrical methods. In particular I am interested in all-pairwise comparisons between the levels of one of the two factors. If I transform the data in ranks, can I reuse the

Re: [R] sunrise, sunset calculation

2004-06-21 Thread Mike Sumner
It was through Hill et al. (2001) that I eventually found the NOAA website: http://www.soest.hawaii.edu/PFRP/elec.tagdata/tagdata.html Hill, Roger D., and Melinda J. Braun (2001). Geolocation by Light Level - The Next Step: Latitude. In: J.R. Sibert and J. Nielsen (Eds.), Electronic Tagging and

[R] Using xtable with summaries of lm objects

2004-06-21 Thread Peter Dunn
Hi all Suppose I do the following: set.seed(1000) library(xtable) x - runif( 10 ) y - 1 + 2*x + rnorm( length(x) ) test.lm - lm( y ~ x ) summary( test.lm ) xtable ( summary( test.lm ) ) The final xtable output follows: % latex table generated in R 1.8.1 by xtable 1.2-2 package % Tue Jun 22

Re: [R] sunrise, sunset calculation

2004-06-21 Thread Mike Sumner
The Johns Hopkins University/Applied Physics Laboratory also has an IDL library with solar elevation code in a file named sunstuff.pro http://fermi.jhuapl.edu/s1r/idl/s1rlib/local_idl.html At 09:09 PM 6/21/2004, Angel Lopez wrote: but it might be better to say that R code 'uses noaa

Re: [R] Html help does not work in Mac OSX 10.3.4

2004-06-21 Thread Ulises Mora Alvarez
Hello! If I launch R from a console I get: R : Copyright 2004, The R Foundation for Statistical Computing Version 1.9.0 (2004-04-12), ISBN 3-900051-00-3 help.start() Making links in per-session dir ... If /usr/bin/open is already running, it is *not* restarted, and you must switch to its

[R] axis labels for stripcharts

2004-06-21 Thread Kathryn Wheatley
I've produced a stripchart but I would like to change the labelling of the x-axis. However, I can not get the scripts used in the plot() function (to eliminate labels) to work (i.e. xaxt=n or axis=FALSE). Are there other scripts availalable that will work with stripchart? Cheers, Kathryn

Re: [R] axis labels for stripcharts

2004-06-21 Thread Gabor Grothendieck
?par Kathryn Wheatley kew at utas.edu.au writes: : I've produced a stripchart but I would like to change the labelling of the : x-axis. However, I can not get the scripts used in the plot() function (to : eliminate labels) to work (i.e. xaxt=n or axis=FALSE). Are there other : scripts

Re: [R] [Q] Newbie (continued.. at least I got R running allready :-)

2004-06-21 Thread Richard A. O'Keefe
jeroen clarysse [EMAIL PROTECTED] wrote: I'll restate the problem : i got a datafile with 2400 measuerements (every 250msec) of a CO2 measurement device, capturing the breath of a subject. I uploaded such a sample here :