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" <giftedlife2...@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 <dcarl...@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-boun...@r-project.org] On Behalf Of Ogbos > Okike > > Sent: Friday, April 22, 2016 2:05 PM > > To: Ulrik Stervbo > > Cc: r-help@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.ster...@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 <giftedlife2...@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@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@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]] ______________________________________________ R-help@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.
