similar to: peaks

Hallo all just a small question I did not find an answer in help pages. Is it possible to use identify() after plotting with plot.POSIX to label points and/or to find out some points? Thanks a lot. Best regards Petr Pikal petr.pikal at precheza.cz p.pik at volny.cz

Dear all I have a character vector of dates something like: timevec<-c("15.5.2003 00:00", "15.5.2003 00:01", "15.5.2003 00:02", "15.5.2003 00:03","15.5.2003 00:04") and I would like to transform it to some more convenient date class. Is there a way how to do it directly without previous reformating to ISO like structure and adding a

Dear all I would like to aggregate a data frame in a three minute interval. I have a time vector: > str(cas.param) `POSIXct', format: chr [1:181] "2003-06-12 09:00:00" "2003-06-12 09:01:00" "2003-06-12 09:02:00" and I want to turn it into a factor by cut. It works if I want to break it in let say hours > str(cut(cas.param,breaks="hour"))

Hello list! I would like to fit a distribution to each of the peaks in a histogram, such as this: http://photos1.blogger.com/blogger/7029/2724/1600/DU145-Bax3-Bcl-xL.png . The peaks are identified using Petr Pikal peaks function ( http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html), but after that I am quite stuck. Any idea as to how I can: Fit a distribution to each peak Integrate the

Dear all Please, is there any function or package for dealing with multimodal distributions? I try to fit multimodal distribution or more precisely to find out mixture of normal distributions which can lead to my actual data. I use optim to find (in that case) two parameters but what I want is to let the function find out arbitrary number of normal distributions underlaying my actual data

Hallo Please I have a question regarding read.table. I wonder if there is any possibility to read a content of 'clipboard' directly (let say as a copy from Excel like programs) or if it is always necessary to save a file in an appropropriate format. I use 1.2.0 W version with W95. Thank you very much. Best regards Petr Pikal Precheza a.s. Nabr.Dr.E.Benese 24 tel: 00420 641 25 2257

Hallo all I went through all files and through some books I have, but still I am not completely sure about the correct aov syntax. Suppose I have a dataset divided to three portions, according to a factor "mereni" which is a measurement day. There is another factor named "frakce" (Each sample is divided into three distinct fractions according to its particle size). samp

Hallo all I have a question regarding POSIX class objects. I try tu use as.POSIXct for creating and plotting some time series. To be able to set labels for time (x) axis I started with plot(..., axes=FALSE) but although it suppressed all other axes, x axis is still plotted together with labels. Is it a bug? see datum.vyber<-seq(ISOdate(2000,1,1), ISOdate(2003,1,1), by="3

Dear all I have a function in which I would like to change some values to NA according to some condition. dropout<-function(y, nahr=FALSE,...) { <some stuff for computing an index> if (nahr) y[index]<<-NA invisible(index) } in case y is a vector all works OK but if it is a part of data frame by calling dropout(df$y) or dropout(df[,number]) no change is done. Please can you

Dear R users I use R (1.3.0) Windows version usually with connection to Excel sheets via RODBC Version: 0.8-2. Everything is OK until I use mixed data types (numeric and character) in one column. In that case I will get a particular column of data frame with NA's. Example vzorek lpas apas NA 51.240 20.26 NA 51.340 20.42 NA 49.770 22.05 RB 12 49.840 19.30 RB 12 49.880

Dear all I would like to ask you about possible bug in using bwplot (from lattice) together with Hmisc library attached. I found it in my actual data, but here is a toy example. It appears only when some levels are missing. library(lattice) library(Hmisc) # preparing data x1<-rnorm(10,5,1) x2<-rnorm(10,5,5) x3<-rnorm(10,1,1) x4<-rnorm(10,1,5) x<-c(x1,x2,x3,x4) x<-c(x,x+5)

Dear all. RE how to estimate parameters of multimodal distribution Thank to prof.Ripley for pointing me to mclust package, although I am not sure I can apply it to my problem. I have another question. I need to change some of my values in data frame to NA. I use something like df[df$v1 < 5, 5:10] <- NA which is OK if there are no NA values in v1. here are some foo attempts >

Dear R users I searched some sources but i did not find an answer.Please give me some hint to following problem. I would like to compute a summary statistic for some vector for different factor levels. I know I can use tapply or aggregate but I do not know if there is a way how to use function with several (two) variable input (like weighted.mean). I wrote a simple a function for factor

I am looking for a function in R that is equivalent to the function "peaks" in Splus. This function gives the local maxima of a vector. I do not find it using the help and I 'm not sure it exits in R. Could someone help me ? Thanks, Herve Cardot -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear all I try to persuade if statement in my function to work as I want (but I am not very successful:( My question is why my function with if statement is evaluated correctly in case of atomic variables and incorrectly in case of vector variables. Here is an example: #function with if statement fff <- function(otac,sklon) { test <- (otac *

I've been asked in private, (and am replying BCC to the asker), >> I saw your post on the R-help archives page about the possibility of >> porting a function from S-Plus called peaks() to R. I am looking for >> some way to locate peaks in a simple x,y data set, and thought that R >> might be the way to go. "of course" it is the way to go, don't get

Hallo all Suppose I have a summary amount of various colours used in different areas colour plastic paint black 15.5 173.8 brown 6.0 523.2 green 2.5 6.4 red 77.1 237.4 yellow 144.6 77.3 It seems that there are some preferable colours in these two areas (yellow in plastic, brown an red in paint). The problem seems to me similar to chi square test but it is intended for testing counts or

In 2004 there was this R-Help posting from Jan 2004: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html R equivalent of Splus peaks() function? The peaks function there has worked well for me on a couple of projects, but some code using "peaks" failed today, which had worked fine in the past. I was looking for a peak in a test case that was a sine curve over one

Hello, Thank you to all those great folks that have helped me in the past (especially Dennis Murphy). I have a new challenge. I often generate time-series data sets that look like the one below, with a variable ("Phase") which has a series of flat-topped peaks (sample data below with 5 "peaks"). I would like to calculate the phase value for each peak. It would be great to

http://r.789695.n4.nabble.com/file/n3312061/x_and_y_values.txt x_and_y_values.txt I have the absorbance values form HPLC Chromatogram. I need to find the peaks in datapoints and area under each peak. I am not sure how how to find the peaks, I tried couple of libraries and peak function but they return me a list of values which when computed from area turns out to be huge nnumbers. In the