similar to: plot.POSIXct with axes FALSE

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"))

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 try to make hourly average by cut() function, which almost works as *I* expected. What puzled me is that if there is only one item at the end of your data it results in NA. Example will explain what I mean datum<-seq(ISOdate(2004,8,31), ISOdate(2004,9,1), "min") cut(datum[1370:1381],"hour", labels=F) [1] 1 1 1 1 1 1 1 1 1 1 1 NA

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 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

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

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

On 15 Oct 2002 at 18:09, Rieckermann Joerg wrote: > Dear Petr, > > I have been fooling around with this peaks function of yours/Ripley > and don't seem to get the main idea. > > What I intend to do is locate the peaks in a vector and later use thes > values. So the position of the peaks would be of great interest. > > My questions: > * Why does peak retunr a

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 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

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 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 *

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 >

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

Hallo Dan I used your code for ls.objects and it is very usefull. I corrected some typo and it worked well untill now. I observed Error due to class evaluation in lists created by aov fit as the fitted object from aov has two classes (aov and lm) So I made a little change which seems to work. Here it is: # Posted by Dan Putler At 2002-07-15 16:13 # As part of the obveRsive GUI project,

Although the documentation is somewhat sketchy, I() can be used to create objects of class AsIs: > I("a") [1] "a" attr(,"class") [1] "AsIs" "character" > I(4) [1] 4 attr(,"class") [1] "AsIs" "numeric" > I(4 + 0i) [1] 4+0i attr(,"class") [1] "AsIs" "complex" > This

I'd like to make a plot showing frequency of an event. The data is in a data from that includes Year, Month and Day (of month) fields, so I created a Date with ISOdate(Year, Month, Day, tz=''). I can plot frequencies for the year 2002 with > thisyear <- Date[Year==2002] > hist( thisyear, xaxt='n' ) > axis.POSIXct( 1, at=seq(min(thisyear), max(thisyear),

Hello R-experts! I am using R Version 1.7.1 (2003-06-16) on a Debian Linux box and I have discovered an odd result when plotting data involving dates. Please try this minimal example: a = seq(ISOdate(2000,1,1), ISOdate(2001,1,1), "months") b = 1:13 plot(a,b, col="red") What I get is a plot that looks as expected except the x-axis is mostly red. Can anyone reproduce this

I am using R-1.7.0 and have some data which consist of one vector of numbers and a second corresponding vector of dates belonging to the POSIXct class. I would like to plot the numbers against the dates. What is the best way to do this? It almost works to just call `plot.' However if I do this while using the `ylab' parameter I get a warning message: parameter "ylab"