Re: [R] more on vector vs array

From: Dimitris Rizopoulos <dimitris.rizopoulos_at_med.kuleuven.be>
Date: Tue 09 Aug 2005 - 18:19:37 EST

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


>
> 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@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@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 Received on Tue Aug 09 18:25:52 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:39:45 EST