>>>>> On Wed, 8 Jan 2020 17:43:29 +0100, Ogbos Okike
<giftedlife2014 at gmail.com> writes:
Ogbos> Dear Enrico,
Ogbos> Thanks for your time.
Ogbos> I have tried to learn how to use dput in R. I have not yet made much
progress.
Ogbos> I have succeeded in using dput to store my data frame. I first
Ogbos> converted my data into a data frame and then used:
Ogbos> dput(dd,file="Ogbos2",control = c("keepNA",
"keepInteger",
Ogbos> "showAttributes")) to output the dput file. dd is my data
frame.
Ogbos> When I opened the file, I didn't like its content as it differs
very
Ogbos> much from my data frame. But I don't know whether that makes
sense to
Ogbos> you. I am attaching the file.
Ogbos> I am thanking you in advance for additional suggestions.
Ogbos> Best wishes
Ogbos> Ogbos
Hello Ogbos
your attempt worked fine: this is the data you sent
dta <- structure(list(dta.year = c(98L, 98L, 98L, 98L, 98L, 98L, 98L,
98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L,
98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L), dta.month = c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), dta.day = c(5L,
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L), dta.hour = c(2L,
3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L,
17L, 18L, 19L, 20L, 21L, 22L, 23L, 0L, 1L, 2L, 3L, 4L, 5L, 6L
), dta.counts = c(6462L, 6450L, 6423L, 6467L, 6480L, 6457L, 6417L,
6467L, 6467L, 6468L, 6500L, 6482L, 6465L, 6465L, 6475L, 6452L,
6440L, 6478L, 6470L, 6422L, 6448L, 6462L, 6485L, 6462L, 6485L,
6470L, 6487L, 6515L, 6488L)), .Names = c("dta.year",
"dta.month",
"dta.day", "dta.hour", "dta.counts"),
row.names = c(NA, -29L),
class = "data.frame")
But then:
head(dta)
## dta.year dta.month dta.day dta.hour dta.counts
## 1 98 1 5 2 6462
## 2 98 1 5 3 6450
## 3 98 1 5 4 6423
## 4 98 1 5 5 6467
## 5 98 1 5 6 6480
## 6 98 1 5 7 6457
The data that you read in has January (1) as month.
So whatever goes wrong, seems to go wrong when you read
the data. Are you quite sure you read the file you
read is the file you have shown?
kind regards
Enrico
Ogbos> On Wed, Jan 8, 2020 at 1:07 PM Enrico Schumann <es at
enricoschumann.net> wrote:
>>
>>
>> Quoting Ogbos Okike <giftedlife2014 at gmail.com>:
>>
>> > Dear Friends,
>> > A sample of my data is:
>> > 98 05 01 02 8541
>> > 98 05 01 03 8548
>> > 98 05 01 04 8512
>> > 98 05 01 05 8541
>> > 98 05 01 06 8509
>> > 98 05 01 07 8472
>> > 98 05 01 08 8454
>> > 98 05 01 09 8461
>> > 98 05 01 10 8462
>> > 98 05 01 11 8475
>> > 98 05 01 12 8433
>> > 98 05 01 13 8479
>> > 98 05 01 14 8417
>> > 98 05 01 15 8463
>> > 98 05 01 16 8473
>> > 98 05 01 17 8450
>> > 98 05 01 18 8433
>> > 98 05 01 19 8437
>> > 98 05 01 20 8437
>> > 98 05 01 21 8438
>> > 98 05 01 22 8421
>> > 98 05 01 23 8420
>> > 98 05 02 00 8371
>> > 98 05 02 01 8338
>> > 98 05 02 02 8251
>> > 98 05 02 03 8204
>> > 98 05 02 04 8183
>> > 98 05 02 05 8231
>> > 98 05 02 06 8242
>> > Columns 1, 2, 3, 4 and 5 stands for year, month, day , hour and
count.
>> >
>> > Using:
>> > Sys.setenv( TZ="GMT" )
>> >
>> >
>> > dta <- read.table("Ohr1may98", col.names =
c("year", "month", "day",
>> > "hour", "counts"))
>> > dta$year <- with( dta, ifelse(year < 50, year + 2000, year
+ 1900))
>> > dta$datetime <- with( dta, as.POSIXct(ISOdatetime(year,
month,day,hour,0,0)))
>> > a = dta$datetime
>> > I converted the datetime and plotted the graph of count vs a.
The plot
>> > was great but I have issues with the date.
>> >
>> > The raw data is for some hours for Ist and second day of may
1998 as
>> > is evident from the sample data. But the result of date stored
in "a"
>> > above shows:
>> >> a
>> > [1] "1998-01-05 02:00:00 GMT" "1998-01-05
03:00:00 GMT"
>> > [3] "1998-01-05 04:00:00 GMT" "1998-01-05
05:00:00 GMT"
>> > [5] "1998-01-05 06:00:00 GMT" "1998-01-05
07:00:00 GMT"
>> > [7] "1998-01-05 08:00:00 GMT" "1998-01-05
09:00:00 GMT"
>> > [9] "1998-01-05 10:00:00 GMT" "1998-01-05
11:00:00 GMT"
>> > [11] "1998-01-05 12:00:00 GMT" "1998-01-05
13:00:00 GMT"
>> > [13] "1998-01-05 14:00:00 GMT" "1998-01-05
15:00:00 GMT"
>> > [15] "1998-01-05 16:00:00 GMT" "1998-01-05
17:00:00 GMT"
>> > [17] "1998-01-05 18:00:00 GMT" "1998-01-05
19:00:00 GMT"
>> > [19] "1998-01-05 20:00:00 GMT" "1998-01-05
21:00:00 GMT"
>> > [21] "1998-01-05 22:00:00 GMT" "1998-01-05
23:00:00 GMT"
>> > [23] "1998-01-06 00:00:00 GMT" "1998-01-06
01:00:00 GMT"
>> > [25] "1998-01-06 02:00:00 GMT" "1998-01-06
03:00:00 GMT"
>> > [27] "1998-01-06 04:00:00 GMT" "1998-01-06
05:00:00 GMT"
>> > [29] "1998-01-06 06:00:00 GMT"
>> > This seems to suggest day 5 and 6 in January 1998 instead of day
1 and
>> > 2 in May of 1998.
>> >
>> > I have spent some time trying to resolve this but I have not
been successful.
>> >
>> > I would be thankful if you could help me to check where I went
astray.
>> >
>> > Thank you.
>> > Best wishes
>> > Ogbos
>> >
>>
>> I cannot reproduce these results. Could you please provide a fully
>> reproducible example, by providing a small example dataset via
'dput(dta)'?
>>
>>
>> --
>> Enrico Schumann
>> Lucerne, Switzerland
>> http://enricoschumann.net
>>
Ogbos> structure(list(dta.year = c(98L, 98L, 98L, 98L, 98L, 98L, 98L,
Ogbos> 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L,
Ogbos> 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L, 98L), dta.month = c(1L,
Ogbos> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
Ogbos> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), dta.day = c(5L,
Ogbos> 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
Ogbos> 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L), dta.hour = c(2L,
Ogbos> 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L,
Ogbos> 17L, 18L, 19L, 20L, 21L, 22L, 23L, 0L, 1L, 2L, 3L, 4L, 5L, 6L
Ogbos> ), dta.counts = c(6462L, 6450L, 6423L, 6467L, 6480L, 6457L, 6417L,
Ogbos> 6467L, 6467L, 6468L, 6500L, 6482L, 6465L, 6465L, 6475L, 6452L,
Ogbos> 6440L, 6478L, 6470L, 6422L, 6448L, 6462L, 6485L, 6462L, 6485L,
Ogbos> 6470L, 6487L, 6515L, 6488L)), .Names = c("dta.year",
"dta.month",
Ogbos> "dta.day", "dta.hour", "dta.counts"),
row.names = c(NA, -29L), class = "data.frame")
--
Enrico Schumann
Lucerne, Switzerland
http://enricoschumann.net