# Re: [R] code reduction (if anyone feels like it)

From: stephen sefick <ssefick_at_gmail.com>
Date: Wed, 16 Jul 2008 06:50:38 -0400

thanks that is a slim version.

Stephen

On Tue, Jul 15, 2008 at 9:03 PM, jim holtman <jholtman_at_gmail.com> wrote:

> Typo in the last one: (resend)
>
> #GPS in Decimal Degrees in the form longitude latitude (raw data)
>
> library(maptools)
>
> # create a list of the coords
> coordList <- list(
> RM215 = matrix(c(-82.1461363, 33.5959109), nrow=1),
>
> SC = matrix(c(-82.025888, 33.606454), nrow=1) ,
>
> RM202 = matrix(c(-81.9906723, 33.5027653), nrow=1),
>
> RM198 = matrix(c(-81.926823, 33.4634678), nrow=1),
>
> HC = matrix(c(-81.920505, 33.46192), nrow=1) ,
>
> RM190 = matrix(c(-81.9317347, 33.3839097), nrow=1),
>
> BC = matrix(c(-81.948189, 33.373043), nrow=1) ,
>
> RM185 = matrix(c(-81.941, 33.3453), nrow=1),
>
> RM179 = matrix(c(-81.890929, 33.317914), nrow=1),
>
> RM148 = matrix(c(-81.7547337, 33.1514072), nrow=1) ,
>
> RM119 = matrix(c(-81.501919, 32.94038), nrow=1),
>
> RM61 = matrix(c(-81.262388, 32.524739), nrow=1))
> #start time for the sequence
> d060101 <- as.POSIXct("2006-01-01", tz="EST")
>
> #sequence in days
> study_seq <- seq(from=d060101, length.out=761, by="days")
>
> # process the list by name so you can use it in the column headings
> times <- lapply(names(coordList), function(.name){
> .sp <- SpatialPoints(coordList[[.name]],
> proj4string=CRS("+proj=longlat +datum=WGS84"))
> .time <- data.frame(
> sunriset(.sp, study_seq, direction="sunrise",
> POSIXct.out=TRUE)\$time,
> sunriset(.sp, study_seq, direction="sunset", POSIXct.out=TRUE)\$time)
> colnames(.time) <- paste(.name, c('sr', 'ss'), sep='')
> .time
> })
> sunriseset <- do.call(cbind, times)
>
>
>
> On Tue, Jul 15, 2008 at 7:20 PM, stephen sefick <ssefick_at_gmail.com> wrote:
> >
> > #GPS in Decimal Degrees in the form longitude latitude (raw data)
> >
> > library(maptools)
> > RM215 <- matrix(c(-82.1461363, 33.5959109), nrow=1)
> >
> > SC <- matrix(c(-82.025888, 33.606454), nrow=1)
> >
> > RM202 <- matrix(c(-81.9906723, 33.5027653), nrow=1)
> >
> > RM198 <- matrix(c(-81.926823, 33.4634678), nrow=1)
> >
> > HC <- matrix(c(-81.920505, 33.46192), nrow=1)
> >
> > RM190 <- matrix(c(-81.9317347, 33.3839097), nrow=1)
> >
> > BC <- matrix(c(-81.948189, 33.373043), nrow=1)
> >
> > RM185 <- matrix(c(-81.941, 33.3453), nrow=1)
> >
> > RM179 <- matrix(c(-81.890929, 33.317914), nrow=1)
> >
> > RM148 <- matrix(c(-81.7547337, 33.1514072), nrow=1)
> >
> > RM119 <- matrix(c(-81.501919, 32.94038), nrow=1)
> >
> > RM61 <- matrix(c(-81.262388, 32.524739), nrow=1)
> >
> >
> > #making the points into a map projection (latlong in WGS84)
> >
> > RM215.sp <- SpatialPoints(RM215, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > #start time for the sequence
> > d060101 <- as.POSIXct("2006-01-01", tz="EST")
> >
> > #sequence in days
> > study_seq <- seq(from=d060101, length.out=761, by="days")
> >
> > #sunrise
> > up.215 <- sunriset(RM215.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > #sunset
> > down.215 <- sunriset(RM215.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > SC.sp <- SpatialPoints(SC, proj4string=CRS("+proj=longlat +datum=WGS84"))
> >
> > up.SC <- sunriset(SC.sp, study_seq, direction="sunrise",
> POSIXct.out=TRUE)
> >
> > down.SC <- sunriset(SC.sp, study_seq, direction="sunset",
> POSIXct.out=TRUE)
> >
> >
> >
> > RM202.sp <- SpatialPoints(RM202, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.202 <- sunriset(RM202.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.202 <- sunriset(RM202.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > RM198.sp <- SpatialPoints(RM198, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.198 <- sunriset(RM198.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.198 <- sunriset(RM198.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > HC.sp <- SpatialPoints(HC, proj4string=CRS("+proj=longlat +datum=WGS84"))
> >
> > up.HC <- sunriset(HC.sp, study_seq, direction="sunrise",
> POSIXct.out=TRUE)
> >
> > down.HC <- sunriset(HC.sp, study_seq, direction="sunset",
> POSIXct.out=TRUE)
> >
> >
> >
> > RM190.sp <- SpatialPoints(RM190, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.190 <- sunriset(RM190.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.190 <- sunriset(RM190.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > BC.sp <- SpatialPoints(BC, proj4string=CRS("+proj=longlat +datum=WGS84"))
> >
> > up.BC <- sunriset(BC.sp, study_seq, direction="sunrise",
> POSIXct.out=TRUE)
> >
> > down.BC <- sunriset(BC.sp, study_seq, direction="sunset",
> POSIXct.out=TRUE)
> >
> >
> >
> > RM185.sp <- SpatialPoints(RM185, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.185 <- sunriset(RM185.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.185 <- sunriset(RM185.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > RM179.sp <- SpatialPoints(RM179, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.179 <- sunriset(RM179.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.179 <- sunriset(RM179.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > RM148.sp <- SpatialPoints(RM148, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.148 <- sunriset(RM148.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.148 <- sunriset(RM148.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > RM119.sp <- SpatialPoints(RM119, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.119 <- sunriset(RM119.sp, study_seq, direction="sunrise",
> > POSIXct.out=TRUE)
> >
> > down.119 <- sunriset(RM119.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> >
> > RM61.sp <- SpatialPoints(RM61, proj4string=CRS("+proj=longlat
> > +datum=WGS84"))
> >
> > up.61 <- sunriset(RM61.sp, study_seq, direction="sunrise",
> POSIXct.out=TRUE)
> >
> > down.61 <- sunriset(RM61.sp, study_seq, direction="sunset",
> > POSIXct.out=TRUE)
> >
> >
> > # put all of the observations side by side sunrise next to sunset then
> put
> > all sites in columns
> >
> > sunriseset <- cbind(as.data.frame(up.215\$time),
> > as.data.frame(down.215\$time), as.data.frame(up.SC\$time),
> > as.data.frame(down.SC\$time),
> >
> > as.data.frame(up.202\$time), as.data.frame(down.202\$time),
> > as.data.frame(up.198\$time), as.data.frame(down.198\$time),
> >
> > as.data.frame(up.HC\$time), as.data.frame(down.HC\$time),
> > as.data.frame(up.190\$time), as.data.frame(down.190\$time),
> > as.data.frame(up.BC\$time),
> >
> > as.data.frame(down.BC\$time), as.data.frame(up.185\$time),
> > as.data.frame(down.185\$time),
> >
> > as.data.frame(up.179\$time), as.data.frame(down.179\$time),
> > as.data.frame(up.148\$time), as.data.frame(down.148\$time),
> >
> > as.data.frame(up.119\$time), as.data.frame(down.119\$time),
> > as.data.frame(up.61\$time), as.data.frame(down.61\$time))
> >
> > #name columns the site and then sunrise (st) and sunset (ss)
> > colnames(sunriseset) <- c("RM215sr", "RM215ss","SCsr", "SCss","RM202sr",
> > "RM202ss" ,"RM198sr", "RM198sr",
> >
> > "HCsr", "HCsr", "RM190sr",
> "RM190ss","BCsr","BCss","RM185sr",
> > "RM185ss","RM179sr", "RM179sss",
> >
> > "RM148sr", "RM148ss","RM119sr", "RM119ss","RM61sr",
> "RM61sr")
> >
> > On Tue, Jul 15, 2008 at 6:10 PM, jim holtman <jholtman_at_gmail.com> wrote:
> >>
> >> I am reasonably sure that it could be more efficient, but it would be
> >> very helpfulto have some comments in the code and a description of
> >> "what is the problem you are trying to solve". It looks like you can
> >> setup some lists and process them more compactly.
> >>
> >> On Tue, Jul 15, 2008 at 3:21 PM, stephen sefick <ssefick_at_gmail.com>
> wrote:
> >> > # I am sure that I could be more efficient than this but how? Thanks
> in
> >> >
> >> > #GPS in Decimal Degrees in the form longitude latitude
> >> > RM215 <- matrix(c(-82.1461363, 33.5959109), nrow=1)
> >> > SC <- matrix(c(-82.025888, 33.606454), nrow=1)
> >> > RM202 <- matrix(c(-81.9906723, 33.5027653), nrow=1)
> >> > RM198 <- matrix(c(-81.926823, 33.4634678), nrow=1)
> >> > HC <- matrix(c(-81.920505, 33.46192), nrow=1)
> >> > RM190 <- matrix(c(-81.9317347, 33.3839097), nrow=1)
> >> > BC <- matrix(c(-81.948189, 33.373043), nrow=1)
> >> > RM185 <- matrix(c(-81.941, 33.3453), nrow=1)
> >> > RM179 <- matrix(c(-81.890929, 33.317914), nrow=1)
> >> > RM148 <- matrix(c(-81.7547337, 33.1514072), nrow=1)
> >> > RM119 <- matrix(c(-81.501919, 32.94038), nrow=1)
> >> > RM61 <- matrix(c(-81.262388, 32.524739), nrow=1)
> >> >
> >> > RM215.sp <- SpatialPoints(RM215, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > d060101 <- as.POSIXct("2006-01-01", tz="EST")
> >> > study_seq <- seq(from=d060101, length.out=761, by="days")
> >> > up.215 <- sunriset(RM215.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.215 <- sunriset(RM215.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > SC.sp <- SpatialPoints(SC, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.SC <- sunriset(SC.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.SC <- sunriset(SC.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM202.sp <- SpatialPoints(RM202, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.202 <- sunriset(RM202.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.202 <- sunriset(RM202.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM198.sp <- SpatialPoints(RM198, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.198 <- sunriset(RM198.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.198 <- sunriset(RM198.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > HC.sp <- SpatialPoints(HC, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.HC <- sunriset(HC.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.HC <- sunriset(HC.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM190.sp <- SpatialPoints(RM190, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.190 <- sunriset(RM190.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.190 <- sunriset(RM190.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > BC.sp <- SpatialPoints(BC, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.BC <- sunriset(BC.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.BC <- sunriset(BC.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM185.sp <- SpatialPoints(RM185, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.185 <- sunriset(RM185.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.185 <- sunriset(RM185.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM179.sp <- SpatialPoints(RM179, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.179 <- sunriset(RM179.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.179 <- sunriset(RM179.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM148.sp <- SpatialPoints(RM148, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.148 <- sunriset(RM148.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.148 <- sunriset(RM148.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM119.sp <- SpatialPoints(RM119, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.119 <- sunriset(RM119.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.119 <- sunriset(RM119.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > RM61.sp <- SpatialPoints(RM61, proj4string=CRS("+proj=longlat
> >> > +datum=WGS84"))
> >> > up.61 <- sunriset(RM61.sp, study_seq, direction="sunrise",
> >> > POSIXct.out=TRUE)
> >> > down.61 <- sunriset(RM61.sp, study_seq, direction="sunset",
> >> > POSIXct.out=TRUE)
> >> >
> >> > sunriseset <- cbind(as.data.frame(up.215\$time),
> >> > as.data.frame(down.215\$time), as.data.frame(up.SC\$time),
> >> > as.data.frame(down.SC\$time),
> >> > as.data.frame(up.202\$time), as.data.frame(down.202\$time),
> >> > as.data.frame(up.198\$time), as.data.frame(down.198\$time),
> >> > as.data.frame(up.HC\$time), as.data.frame(down.HC\$time),
> >> > as.data.frame(up.190\$time), as.data.frame(down.190\$time),
> >> > as.data.frame(up.BC\$time),
> >> > as.data.frame(down.BC\$time), as.data.frame(up.185\$time),
> >> > as.data.frame(down.185\$time),
> >> > as.data.frame(up.179\$time), as.data.frame(down.179\$time),
> >> > as.data.frame(up.148\$time), as.data.frame(down.148\$time),
> >> > as.data.frame(up.119\$time), as.data.frame(down.119\$time),
> >> > as.data.frame(up.61\$time), as.data.frame(down.61\$time))
> >> > colnames(sunriseset) <- c("RM215sr", "RM215ss","SCsr",
> "SCss","RM202sr",
> >> > "RM202ss" ,"RM198sr", "RM198sr",
> >> > "HCsr", "HCsr", "RM190sr",
> >> > "RM190ss","BCsr","BCss","RM185sr",
> >> > "RM185ss","RM179sr", "RM179sss",
> >> > "RM148sr", "RM148ss","RM119sr", "RM119ss","RM61sr",
> >> > "RM61sr")
> >> >
> >> > --
> >> > Let's not spend our time and resources thinking about things that are
> so
> >> > little or so large that all they really do for us is puff us up and
> make
> >> > us
> >> > feel like gods. We are mammals, and have not exhausted the annoying
> >> > little
> >> > problems of being mammals.
> >> >
> >> > -K. Mullis
> >> >
> >> > [[alternative HTML version deleted]]
> >> >
> >> > ______________________________________________
> >> > R-help_at_r-project.org mailing list
> >> > https://stat.ethz.ch/mailman/listinfo/r-help
> >> > http://www.R-project.org/posting-guide.html
> >> > and provide commented, minimal, self-contained, reproducible code.
> >> >
> >>
> >>
> >>
> >> --
> >> Jim Holtman
> >> Cincinnati, OH
> >> +1 513 646 9390
> >>
> >> What is the problem you are trying to solve?
> >
> >
> >
> > --
> > Let's not spend our time and resources thinking about things that are so
> > little or so large that all they really do for us is puff us up and make
> us
> > feel like gods. We are mammals, and have not exhausted the annoying
> little
> > problems of being mammals.
> >
> > -K. Mullis
>
>
>
> --
> Jim Holtman
> Cincinnati, OH
> +1 513 646 9390
>
> What is the problem you are trying to solve?
>

```--
Let's not spend our time and resources thinking about things that are so
little or so large that all they really do for us is puff us up and make us
feel like gods. We are mammals, and have not exhausted the annoying little
problems of being mammals.

-K. Mullis

[[alternative HTML version deleted]]

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help