Displaying 20 results from an estimated 2000 matches similar to: "Using density() and turnpoint to Identify maxima in data"
2007 May 23
1
Question concerning "pastecs" package
Hi
I just installed the pastecs package and I am wondering: is there an
english (or german) translation of the file pastecs.pdf? If not, is
there an explanation somewhere of the object of type 'turnpoints' as a
result of turnpoints(), especially the "info" field?
Thanks,
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE
2007 Sep 20
2
referencing packages?
Hi
I know how to referenc R in a scientific paper - but is there a
standardised way to reference packages?
Thanks
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH
Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
Biology (UCT)
Plant Conservation Unit
Department of Botany
University of Cape Town
Rondebosch 7701
2007 Oct 01
4
Concatenating one character vector into one string
Hi
I am sure this is simple - but how can I convert one charecter vector
into one string?
example:
x <- c("This ", "is ", "one ", "sentence.")
should become
"This is one entence"
Thanks
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH
Rainer M. Krug, Dipl. Phys.
2008 Jan 27
1
Putting frame around single panels in ggplot 2 and facet_grid()
Hi
I want to highlight two panels in a grid created with facet_grid() by
putting a box around it or usiong another background colour. Is this
possible, and if yes, how?
Thanks
Rainer
--
Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
Biology (UCT)
Plant Conservation Unit
Department of Botany
University of Cape Town
Rondebosch 7701
South Africa
Tel: +27 - (0)21 650 5776 (w)
Fax:
2007 May 28
1
Where to find "nprq"?
Hi
I am trying to install the package "pheno", but it needs the package
"nprq" by Roger Koenker et al. which I can I find this package? It does
not seem to be on CRAN and googling also doesn't give me an URL - is it
still somewhere available?
Thanks,
Rainer
--
NEW EMAIL ADDRESS AND ADDRESS:
Rainer.Krug at uct.ac.za
RKrug at sun.ac.za WILL BE DISCONTINUED END OF
2005 May 16
0
Turnpoints (pastecs): How to specify a limit on the number of tur npoints?
Hello,
I'm trying to get a few turnpoints for a financial time series. There is a
function in pastecs that does that. However, I get a large number of
turnpoints:
library(pastecs)
data(EuStockMarkets)
dax <-EuStockMarkets[,1]
plot(dax)
turnp <-turnpoints(dax)
summary(turnp) #gives 925 peaks/pits!!!
How can specify to get only 30 turnpoints?
Second question: the extract function
2007 Aug 16
6
several plots on several pages
Hi
> version
_
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major 2
minor 5.1
year 2007
month 06
day 27
svn rev 42083
language R
version.string R version 2.5.1 (2007-06-27)
I want to create a pdf withe three graphs on a page and with two pages:
2008 Mar 14
2
Selecting elements in vector
Hi
Consider the following code
> x <- rep(1:13, 13)
> y <- 1:3
I want to select all elements in x which are equal to 1, 2 or 3.
I know that I could use
> sel <- x==y[1] | x==y[2] | x==y[3]
> x[sel]
to obtain the values, but in my analysis, the y-vector is thousands of
elements long.
Is there any way, that I can do that easily?
Thanks
Rainer
--
Rainer M. Krug, Dipl.
2009 Sep 17
1
Turning points in a series
Good morning once more. My problem of yesterday has been addressed. Having
learned a few tricks from that, I wish to ask another question in connection
with that. My data is a cosmic ray data consisting of dates and counts. When
I plot a graph of counts versus dates, the resultant signal shows a number
of maximum and minimum points. These minimum points (turning points) are of
interest to me.
2008 Mar 07
1
How to navigate in layout() created graph?
Hi
I created a complex layout with using layout() and it looks exactly as
I need it. But I don't want to print in the order in which the
subfigure are numbered, but in a different order.
How can I navigate in the layout so that I can specify the subfigure
in which to plot?
At the moment I am using a function which is converting the number to
mfg parameter for par, but it does not seem to
2008 Jul 09
1
version problems of rkward in ubuntu hardy repository
Hi
I tried to install rkward under ubuntu hardy heron, but it tried to
use the one from the cran repository which was newer, but it did not
install. To be able to install rkward, I had to disable the cran
repository, install rkward, lock it's version and enable the
repository again.
I have the feeling, that the dependencies are not correct (it is
looking for libqt4GUI, but I have a
2008 Feb 27
2
Plan to build Package to use GRASS from R
Hi
Sorry for crossposting, but I think this can be of interest for GRASS and R
users.
I am planning to write a package to make the use of GRASS from R easier. The
idea is to wrap the system call to execute the GRASS command into an R
command of the same name.
e.g:
r.to.vect <- function(..., intern=TRUE, ignore.stderr=FALSE)
{
comm <- paste( "r.to.vect ", ...,
2008 Feb 14
1
Analysis with spatstat and Kcross() requires to much memory
Hi
I am running analysis with Kcross from the package spatstat and I am getting
a message that R can not allocate enough memory for a vector of 900MB.
R seems to be running towards the 2GB limit per process.
The dataset is not to big (ca 3000 points) but the mask for the points is
extremely irregular (a buffer around roads which have been sampled) and I
can do trhe analysis if I use a
2008 May 08
0
RMySQL and R 2.7.0 - "Error in field.types$row.names : $ operator is invalid for atomic vectors"
Hi
when executin the following code, I get an "$ operator is invalid for
atomic vectors". I understand the meaning of the error (and have seen
the warnings in earlier R versions). My question is, is there an
updated version of RMySQL which deals with it, or is my on;ly option
to switch to RODBC (which I would not like to do as I am using Linux)?
If I have to use RODBC, how can I easily
2009 Sep 16
2
Turnpoints
Good evening to everybody. I have a data of four columns - three of the
columns represent date while the fourth is counts. Setting 4000 as my
minimum point/value, I wish to find out the number of counts that are less
or equal to 4000 and the dates they occur. I have installed pastecs and have
been trying to use turnpoints function to do this. I have not been much
lucky. I am a still learning. I
2007 Mar 05
4
Identifying last record in individual growth data over different time intervalls
Hi
I have a plist t which contains size measurements of individual plants,
identified by the field "plate". It contains, among other, a field
"year" indicating the year in which the individual was measured and the
"height". The number of measurements range from 1 to 4 measurements in
different years.
My problem is that I would need the LAST measurement. I only
2005 Aug 03
0
possibility of R/Maxima linking
Hello,
I'm aware that recently there has been discussion of R joining forces with
an algebra system in some way. Here is my $0.02.
I suggest linking R with Maxima (http://maxima.sf.net).
Maxima is written in Lisp. Maxima objects are pretty simple.
Maxima works on expressions, and almost every expression
is a Lisp list of the form ((op) arg1 arg2 arg3 ...) which
associates an operator and
2011 Jan 04
1
function masking and gmp questions
Hi,
Here's the problem I ran into: the gmp package has a method for apply()
so it masks the base::apply function. With gmp installed, I tried to
run the function turnpoints() from the pastecs package. It fails
because it calls apply() internally, like this:
apply(mymatrix,1,max,na.rm=TRUE)
,
but the code in the gmp package which sets up the operator overload for
apply() strictly
2010 Feb 11
1
Fwd: Recall function: "evaluation nested too deeply: infinite recursion / options(expressions=)?"
hello all
i searched the archives and couldn't get a solution to the following question.
i have the following function:
F=function(z,v)
{
if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1))
}
and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )
2006 Mar 15
1
Additional arguments in S3 method produces a warning
Hello,
I just notice this:
> x <- c(1:4,0:5, 4, 11)
> library(pastecs)
Loading required package: boot
> tp <- turnpoints(x)
> extract(tp, no.tp = FALSE, peak = TRUE, pit = FALSE)
[1] FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE TRUE FALSE
FALSE
Warning message:
arguments after the first two are ignored in: UseMethod("extract", e, n,
...)
>