[R] clock24.plot/radial plot: Fixed

Ogbos Okike giftedlife2014 at gmail.com
Fri Apr 29 21:40:20 CEST 2016


Dear All,
This problem is over. Clock24.plot did the job. Thanks to all those who
assisted me.
Ogbos
On Apr 22, 2016 8:34 PM, "Ogbos Okike" <giftedlife2014 at gmail.com> wrote:

> Dear All,
> One hand. Many thanks!! The code run as soon as I loaded lubridate.
>
> Please can you guide me on how to relate this code to my actual data.
> My actual data is looking like:
> 2005/01/01 00:00   4009
> 2005/01/01 01:00   3969
> 2005/01/01 02:00   3946
> 2005/01/01 03:00   3975
> 2005/01/01 04:00   3960
> 2005/01/01 05:00   3974
> 2005/01/01 06:00   3971
> 2005/01/01 07:00   3970
> 2005/01/01 08:00   3962
> 2005/01/01 09:00   3992
> 2005/01/01 10:00   3955
> 2005/01/01 11:00   3963
> 2005/01/01 12:00   3965
> 2005/01/01 13:00   3947
> 2005/01/01 14:00   3959
> 2005/01/01 15:00   3978
> 2005/01/01 16:00   3967
> and it runs for many years.
>
> Many thanks for time.
> Ogbos
>
> On 4/22/16, David L Carlson <dcarlson at tamu.edu> wrote:
> > Looks like you forgot to load the lubridate package
> >
> > library(lubridate)
> >
> > You are calling functions days(), hours(), minutes(), seconds(), and
> hour()
> > which all come from that package.
> >
> > -------------------------------------
> > David L Carlson
> > Department of Anthropology
> > Texas A&M University
> > College Station, TX 77840-4352
> >
> > -----Original Message-----
> > From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Ogbos
> Okike
> > Sent: Friday, April 22, 2016 2:05 PM
> > To: Ulrik Stervbo
> > Cc: r-help at r-project.org
> > Subject: Re: [R] clock24.plot/radial plot
> >
> > Kind Experts,
> > Many thanks for your guide. I have tried to figure out something that
> > can help me plot my own data using the examples you referred me to. I
> > copied part of the code as:
> >
> > set.seed(44)
> > N=500
> > events <- as.POSIXct("2011-01-01", tz="GMT") +
> >               days(floor(365*runif(N))) +
> >               hours(floor(24*rnorm(N))) +  # using rnorm here
> >               minutes(floor(60*runif(N))) +
> >               seconds(floor(60*runif(N)))
> > hour_of_event <- hour(events)
> > # make a dataframe
> > eventdata <- data.frame(datetime = events, eventhour = hour_of_event)
> > # determine if event is in business hours
> > eventdata$Workday <- eventdata$eventhour %in% seq(9, 17)
> > library(circular)
> > eventdata$eventhour <- circular(hour_of_event%%24, # convert to 24 hrs
> >       units="hours", template="clock24")
> > rose.diag(eventdata$eventhour, bin = 24, col = "lightblue", main =
> > "Events by Hour (sqrt scale)",
> >     prop = 3)
> > I tried to run the above but got an error message: "Error in
> > eval(expr, envir, enclos) : could not find function "days"
> > I was thinking that if I could run this code, I can see what is doing
> > and then start trying to see if I can adapt it to solve my problem.
> >
> > Thank you so much for further assistance.
> > Ogbos
> >
> >
> > On 4/22/16, Ulrik Stervbo <ulrik.stervbo at gmail.com> wrote:
> >> I use ggplot2 for all my plotting needs where you can make plots
> circular
> >> with the coord_polar. Maybe this will help you along:
> >>
> http://rstudio-pubs-static.s3.amazonaws.com/3369_998f8b2d788e4a0384ae565c4280aa47.html
> >>
> >> On Fri, 22 Apr 2016 at 08:31 Ogbos Okike <giftedlife2014 at gmail.com>
> >> wrote:
> >>
> >>> Dear All,
> >>> I am trying to generate a circular/radial plot. The script below has a
> >>> result I am looking for:
> >>> testlen<-rnorm(24)*2+5
> >>>  testpos<-0:23+rnorm(24)/4
> >>>  clock24.plot(testlen,testpos,main="Test Clock24
> >>> (lines)",show.grid=FALSE,
> >>>   line.col="green",lwd=3)
> >>>  if(dev.interactive()) par(ask=TRUE)
> >>>  # now do a 'daylight' plot
> >>>  oldpar<-clock24.plot(testlen[7:19],testpos[7:19],
> >>>   main="Test Clock24 daytime (symbols)",
> >>>   point.col="blue",rp.type="s",lwd=3)
> >>>  # reset everything
> >>>  par(oldpar)
> >>>
> >>> I tried to play with the script to work with my data. I read my data:
> >>> swe<-scan("onedaydata",list(dates="",time="",count=""))
> >>> dates<-swe$dates
> >>> times<-swe$time
> >>> count<-swe$count.
> >>> I tried to replace testlen<-rnorm(24)*2+5 with testlen<-count and
> >>> oldpar<-clock24.plot(testlen[7:19],testpos[7:19], with
> >>> oldpar<-clock24.plot(testlen[0:23],testpos[0:23], but nothing worked.
> >>> The format of my data is 2005/01/01 00:00   4009
> >>> 2005/01/01 01:00   3969
> >>> 2005/01/01 02:00   3946
> >>> 2005/01/01 03:00   3975
> >>> 2005/01/01 04:00   3960
> >>> 2005/01/01 05:00   3974
> >>> 2005/01/01 06:00   3971
> >>> 2005/01/01 07:00   3970
> >>> 2005/01/01 08:00   3962
> >>> 2005/01/01 09:00   3992
> >>> 2005/01/01 10:00   3955
> >>> 2005/01/01 11:00   3963
> >>> 2005/01/01 12:00   3965
> >>> 2005/01/01 13:00   3947
> >>> 2005/01/01 14:00   3959
> >>> 2005/01/01 15:00   3978
> >>> 2005/01/01 16:00   3967
> >>> 2005/01/01 17:00   3978
> >>> 2005/01/01 18:00   3988
> >>> 2005/01/01 19:00   4043
> >>> 2005/01/01 20:00   4026
> >>> 2005/01/01 21:00   3996
> >>> 2005/01/01 22:00   3967
> >>> 2005/01/01 23:00   3969
> >>> 2005/01/02 00:00   3976
> >>> 2005/01/02 01:00   3969
> >>> 2005/01/02 02:00   3955
> >>> 2005/01/02 03:00   3984
> >>> 2005/01/02 04:00   3971
> >>> 2005/01/02 05:00   3960
> >>> 2005/01/02 06:00   3951
> >>> 2005/01/02 07:00   3948
> >>> 2005/01/02 08:00   3954
> >>> 2005/01/02 09:00   3948
> >>> 2005/01/02 10:00   3960
> >>> 2005/01/02 11:00   3964
> >>> 2005/01/02 12:00   3962
> >>> 2005/01/02 13:00   3959
> >>> 2005/01/02 14:00   3950
> >>> 2005/01/02 15:00   3972
> >>> 2005/01/02 16:00   3984
> >>> 2005/01/02 17:00   3983
> >>> 2005/01/02 18:00   3982
> >>> 2005/01/02 19:00   3987
> >>> 2005/01/02 20:00   3989
> >>> 2005/01/02 21:00   3975
> >>> 2005/01/02 22:00   3956
> >>> 2005/01/02 23:00   3975
> >>> 2005/01/03 00:00   3946
> >>> 2005/01/03 01:00   3944
> >>> 2005/01/03 02:00   3915
> >>> 2005/01/03 03:00   3901
> >>> 2005/01/03 04:00   3893
> >>> 2005/01/03 05:00   3854
> >>> 2005/01/03 06:00   3824
> >>> 2005/01/03 07:00   3790
> >>> 2005/01/03 08:00   3770
> >>> 2005/01/03 09:00   3794
> >>> 2005/01/03 10:00   3778
> >>> 2005/01/03 11:00   3803
> >>> 2005/01/03 12:00   3801
> >>> 2005/01/03 13:00   3800
> >>> 2005/01/03 14:00   3783
> >>> 2005/01/03 15:00   3789
> >>> 2005/01/03 16:00   3804
> >>> 2005/01/03 17:00   3781
> >>> 2005/01/03 18:00   3785
> >>> 2005/01/03 19:00   3772
> >>> 2005/01/03 20:00   3777
> >>> 2005/01/03 21:00   3766
> >>> 2005/01/03 22:00   3775
> >>> 2005/01/03 23:00   3779
> >>> 2005/01/04 00:00   3798
> >>> 2005/01/04 01:00   3806
> >>>
> >>> A sample of the plot I want is attached. My data is quite large.
> >>> Thanks for your time.
> >>> Best wishes
> >>> Ogbos
> >>> ______________________________________________
> >>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> >>> https://stat.ethz.ch/mailman/listinfo/r-help
> >>> PLEASE do read the posting guide
> >>> http://www.R-project.org/posting-guide.html
> >>> and provide commented, minimal, self-contained, reproducible code.
> >>
> >
> > ______________________________________________
> > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> > http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> >
>

	[[alternative HTML version deleted]]



More information about the R-help mailing list