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.
>
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.
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. >