Please run your sequence of R statements one at a time so you can tell where
the problem is. Only "run a script" after the code works one line at a time.
There are too many places where things can go wrong otherwise.
Is your file being read in properly? Is the filename correct? was the header
properly recognized?
I see something odd about the input you are providing to the "origin" argument
to the as.Date function.
On June 15, 2024 1:42:45 PM PDT, Jibrin Alhassan <jibrin.alhas...@unn.edu.ng>
wrote:
>Thank you Rui. I ran the following script
>df1 <- read.table("solar_hour", header = TRUE)
>df1$date <- as.Date(paste(df1$year, df1$hour),
> format = "%Y %j",
>origin = "2012-08-01-0")
>df2 <- df1[c("date", "IMF", "SWS", "SSN", "Dst", "f10")]
>head(df1)
>#To display all the rows
> print(df2).
>It gave me this error message
>> source ("script.R")
>Error in `$<-.data.frame`(`*tmp*`, date, value = numeric(0)) :
> replacement has 0 rows, data has 38735
>> print(df2)
>Error: object 'df2' not found
>> My data is an hourly data but desire to have the date as
>year month day hour
>2012 08 01 01
>2012 08 01 02
>2012 08 01 03 etc
>Thanks.
>
>*Jibrin Adejoh Alhassan (Ph.D)*
>Department of Physics and Astronomy,
>University of Nigeria, Nsukka
>
>
>On Sat, Jun 15, 2024 at 8:34 PM Rui Barradas <ruipbarra...@sapo.pt> wrote:
>
>> Às 20:00 de 15/06/2024, Jibrin Alhassan escreveu:
>> > I have solar-geophysical data e.g as blow:
>> > YEAR DOY HR IMF SW SSN Dst f10.7
>> > 2012 214 0 3.4 403. 132 -9 154.6
>> > 2012 214 1 3.7 388. 132 -10 154.6
>> > 2012 214 2 3.7 383. 132 -10 154.6
>> > 2012 214 3 3.7 391. 132 -9 154.6
>> > 2012 214 4 4.2 399. 132 -7 154.6
>> > 2012 214 5 4.1 411. 132 -6 154.6
>> > 2012 214 6 4.0 407. 132 -6 154.6
>> > 2012 214 7 4.2 404. 132 -4 154.6
>> > 2012 214 8 4.3 405. 132 -6 154.6
>> > 2012 214 9 4.4 409. 132 -6 154.6
>> > 2012 214 10 4.4 401. 132 -6 154.6
>> > 2012 214 11 4.5 385. 132 -7 154.6
>> > 2012 214 12 4.7 377. 132 -8 154.6
>> > 2012 214 13 4.7 382. 132 -6 154.6
>> > 2012 214 14 4.3 396. 132 -4 154.6
>> > 2012 214 15 4.1 384. 132 -2 154.6
>> > 2012 214 16 4.0 382. 132 -1 154.6
>> > 2012 214 17 3.9 397. 132 0 154.6
>> > 2012 214 18 3.8 390. 132 1 154.6
>> > 2012 214 19 4.2 400. 132 2 154.6
>> > 2012 214 20 4.6 408. 132 1 154.6
>> > 2012 214 21 4.8 401. 132 -3 154.6
>> > 2012 214 22 4.9 395. 132 -5 154.6
>> > 2012 214 23 5.0 386. 132 -1 154.6
>> > 2012 215 0 5.0 377. 143 -1 138.6
>> > 2012 215 1 4.9 384. 143 -2 138.6
>> > 2012 215 2 4.9 390. 143 -4 138.6
>> > 2012 215 3 4.9 372. 143 -6 138.6
>> > 2012 215 4 5.1 371. 143 -4 138.6
>> > I want to process it to be of the format as shown below
>> > y m d hr imf sws ssn Dst f10.7
>> > 2012-08-01 10 3.4 403. 132 -9 154.6
>> > 2012-08-01 12 3.7 388. 132 -10 154.6
>> > 2012-08-01 15 3.7 383. 132 -10 154.6
>> > 2012-08-01 17 3.7 391. 132 -9 154.6
>> > I want to request an R code to accomplish this task. Thanks for your
>> time.
>> > *Jibrin Adejoh Alhassan (Ph.D)*
>> > Department of Physics and Astronomy,
>> > University of Nigeria, Nsukka
>> >
>> > [[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.
>> Hello,
>>
>> To create a date column, paste the first two columns and coerce to class
>> "Date" with conversion specifications %Y for the 4 digit year and %j for
>> the day of year. See
>>
>> help("strptime")
>>
>>
>>
>> df1 <- read.table(text = "YEAR DOY HR IMF SW SSN Dst f10.7
>> 2012 214 0 3.4 403. 132 -9 154.6
>> 2012 214 1 3.7 388. 132 -10 154.6
>> 2012 214 2 3.7 383. 132 -10 154.6
>> 2012 214 3 3.7 391. 132 -9 154.6
>> 2012 214 4 4.2 399. 132 -7 154.6
>> 2012 214 5 4.1 411. 132 -6 154.6
>> 2012 214 6 4.0 407. 132 -6 154.6
>> 2012 214 7 4.2 404. 132 -4 154.6
>> 2012 214 8 4.3 405. 132 -6 154.6
>> 2012 214 9 4.4 409. 132 -6 154.6
>> 2012 214 10 4.4 401. 132 -6 154.6
>> 2012 214 11 4.5 385. 132 -7 154.6
>> 2012 214 12 4.7 377. 132 -8 154.6
>> 2012 214 13 4.7 382. 132 -6 154.6
>> 2012 214 14 4.3 396. 132 -4 154.6
>> 2012 214 15 4.1 384. 132 -2 154.6
>> 2012 214 16 4.0 382. 132 -1 154.6
>> 2012 214 17 3.9 397. 132 0 154.6
>> 2012 214 18 3.8 390. 132 1 154.6
>> 2012 214 19 4.2 400. 132 2 154.6
>> 2012 214 20 4.6 408. 132 1 154.6
>> 2012 214 21 4.8 401. 132 -3 154.6
>> 2012 214 22 4.9 395. 132 -5 154.6
>> 2012 214 23 5.0 386. 132 -1 154.6
>> 2012 215 0 5.0 377. 143 -1 138.6
>> 2012 215 1 4.9 384. 143 -2 138.6
>> 2012 215 2 4.9 390. 143 -4 138.6
>> 2012 215 3 4.9 372. 143 -6 138.6
>> 2012 215 4 5.1 371. 143 -4 138.6", header = TRUE)
>>
>>
>> with(df1, paste(YEAR, DOY)) |> as.Date(format = "%Y %j")
>> #> [1] "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01"
>> #> [6] "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01"
>> #> [11] "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01"
>> #> [16] "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01"
>> #> [21] "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-01" "2012-08-02"
>> #> [26] "2012-08-02" "2012-08-02" "2012-08-02" "2012-08-02"
>>
>> # now create the column
>> df1$Date <- with(df1, paste(YEAR, DOY)) |> as.Date(format = "%Y %j")
>> # remove the columns no longer needed
>> df1 <- df1[-(1:2)]
>> # relocate the new date column
>> df1 <- df1[c(ncol(df1), 1:(ncol(df1) - 1L))]
>> head(df1)
>> #> Date HR IMF SW SSN Dst f10.7
>> #> 1 2012-08-01 0 3.4 403 132 -9 154.6
>> #> 2 2012-08-01 1 3.7 388 132 -10 154.6
>> #> 3 2012-08-01 2 3.7 383 132 -10 154.6
>> #> 4 2012-08-01 3 3.7 391 132 -9 154.6
>> #> 5 2012-08-01 4 4.2 399 132 -7 154.6
>> #> 6 2012-08-01 5 4.1 411 132 -6 154.6
>>
>>
>> Hope this helps,
>>
>> Rui Barradas
>>
>>
>> --
>> Este e-mail foi analisado pelo software antivírus AVG para verificar a
>> presença de vírus.
>> www.avg.com
>>
>
> [[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.
--
Sent from my phone. Please excuse my brevity.
______________________________________________
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.
