Unpack PDE file as CSV file or ascii screen dump
getPDEcsv(pde = "/home/lees/Site/Santiaguito/pdq.eqs")
getPDEscreen(pde = "/home/lees/Site/Santiaguito/pdq.eqs")
character, file name
list of locations, times and magnitude
Download pde from: http://neic.usgs.gov/neis/epic/epic.html. csv version uses comma separated values. screen versions uses the screen dump and a parser
http://neic.usgs.gov/neis/epic/epic.html
# NOT RUN {
# }
# NOT RUN {
eq1 <- getPDEcsv(pde = "/home/lees/Site/Santiaguito/pdq.eqs")
############################
###### copy/paste from the screen dump at the NEIC web site
##### then:
eq1 <- getPDEscreen(pde = "EQS.txt")
nk <- names(eq1[[1]])
k <- length(nk)
########## get a matrix of values instead
jim <- matrix( unlist(eq1), ncol=k, byrow=TRUE)
colnames(jim) <- nk
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab