[R] more on vector vs array
Adaikalavan Ramasamy
ramasamy at cancer.org.uk
Tue Aug 9 12:05:28 CEST 2005
Nice one. But I think you could replace the last line (the one with
do.call) with the simpler
w <- which( dat[ ,2] > dat[ ,3] )
w
[1] 6 11 13 14 16 18 20
dat[ w, ]
station temp.x temp.y
6 b 18 16
11 c 17 15
13 d 16 14
14 d 17 14
16 d 17 14
18 e 16 15
20 e 19 15
Thank you.
Regards, Adai
On Tue, 2005-08-09 at 10:19 +0200, Dimitris Rizopoulos wrote:
> you could use something like this:
>
> dat1 <- data.frame(station = rep(letters[1:5], 4), temp =
> round(rnorm(20, 15, 3)))
> dat2 <- data.frame(station = letters[1:5], temp = round(rnorm(5, 15,
> 4)))
> ################
> dat <- merge(dat1, dat2, by = "station")
> do.call("rbind", lapply(split(dat, dat$station), function(x){
> out <- x[x$temp.x > x$temp.y, ]
> if(nrow(out)) out else rep(NA, length(x))
> }))
>
>
> I hope it helps.
>
> Best,
> Dimitris
>
> ----
> Dimitris Rizopoulos
> Ph.D. Student
> Biostatistical Centre
> School of Public Health
> Catholic University of Leuven
>
> Address: Kapucijnenvoer 35, Leuven, Belgium
> Tel: +32/16/336899
> Fax: +32/16/337015
> Web: http://www.med.kuleuven.be/biostat/
> http://www.student.kuleuven.be/~m0390867/dimitris.htm
>
>
> ----- Original Message -----
> From: "alessandro carletti" <alxmilton at yahoo.it>
> To: "rHELP" <R-help at stat.math.ethz.ch>
> Sent: Tuesday, August 09, 2005 9:58 AM
> Subject: [R] more on vector vs array
>
>
> >
> > Ok, thanks,
> > I'll try with a simplier example:
> >
> > I have a vector with 4 levels
> >
> > dataframe 1
> > station temp
> > aaa 12
> > aaa 13
> > bbb 12
> > bbb 20
> > aaa 23
> > bbb 21
> > ccc 30
> > ccc 18
> > ddd 15
> > aaa 11
> > ddd 15
> > ddd 10
> >
> >
> > and a thresholds vector
> >
> > station thr
> > aaa 20
> > bbb 18
> > ccc 25
> > ddd 10
> >
> >
> > I vant to select from dataframe 1 each value (level by
> > level) > its own threshold value.
> > How to do it automatically? (vector temp and vector
> > thr have different length)
> >
> > Thanks
> >
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide!
> > http://www.R-project.org/posting-guide.html
> >
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>
More information about the R-help
mailing list