Speaking of open source alternatives to Matlab, Ben Singer, as I
recall, had been working on using the Psycho-Toolbox from octave.
I don't know how far he got with that or what he is up to now but
his current web page (first hit on google for "Ben Singer") is:
http://www.princeton.edu/~bdsinger/
"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 :
http://www.psy.kuleuven.ac.be/lee
?par
Kathryn Wheatley 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 availalable
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
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
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 algorithms'
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 09:
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
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 TukeyHSD
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
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:
http://cran.r-p
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...
Tha
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
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] ma
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 w
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
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 fi
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:
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 perverse,
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
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
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 t
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
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.,
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 =
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
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,
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 - specifically,
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 abo
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 pri
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' (
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:
>
> http://lib.stat.cmu.edu/R/CRAN/src/contrib/Descriptions/H
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
http://lib.stat.cmu.edu/R/CRAN/src/contrib/Des
>> 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
__
[EMAI
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 estim
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 to
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 be
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
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 PROTEC
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 est
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
>
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 in
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
> 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
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 a
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, functi
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, w
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 of
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,
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 pra
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 grave
-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!
|
| Altho
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:
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
b
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<=x<0.5
> 0.5<=x<1
> 1<=x<1.5
> 1.5<=x<2
> and so on...
Hallo
?cut and ?table should be what you a
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
--
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 divid
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<=x<0.5
0.5<=x<1
1<=x<1.5
1.5<=x<2
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
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'
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
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 acce
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,longit
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?
> "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'l
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
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
https://www.stat.math.ethz.
> 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 t
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 src/main/arith
68 matches
Mail list logo