clock24.plot/radial plot
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.