inst/data-raw/qwdata_codes_datacreate.R

# source 1
# r - Read multiple tables in from a single text file? - Stack Overflow answered by by0 on Sep 28 2012. See \url{https://stackoverflow.com/questions/7264153/read-multiple-tables-in-from-a-single-text-file}.

# source 2
# loop through all files in a directory, read and save them in an object R - Stack Overflow answered by iembry on Jan 13 15. See \url{https://stackoverflow.com/questions/27911604/loop-through-all-files-in-a-directory-read-and-save-them-in-an-object-r/27912952}.

# source 3
# How to create automatic text file from a list in R? - Stack Overflow answered by Richard Scrivenon Sep 26 14 and edited by Richard Scrivenon Sep 27 14. See \url{https://stackoverflow.com/questions/26068156/how-to-create-automatic-text-file-from-a-list-in-r}.

# Data source \url{https://waterdata.usgs.gov/usa/nwis/qwdata?codes_table26_help}
# Table 26. Parameter codes with fixed values

library("data.table")

# source 1 begins
qwdata_codes <- readLines("./data-raw/qwdata?codes_table26_help.txt")

begin <- grep("^(\\d{5})(\\s)(\\w+)", qwdata_codes)
edge <- vector("integer", length(qwdata_codes))
edge[begin] <- 1
qwdata_codes_names <- qwdata_codes[begin]

edges <- cumsum(edge)

df <- lapply(split(qwdata_codes, edges), function(dat) {

input <- read.delim(textConnection(dat), skip = 1, header = FALSE, fill = TRUE, sep = "$", colClasses = c("character", "numeric", "character"), stringsAsFactors = FALSE)
attr(input, "name") <- dat[1]  # save the name

colnames(input) <- c("Parameter Code", "Fixed Value", "Fixed Text")

input

})

names(df) <- sapply(df, attr, "name")
# source 1 ends

longnames <- names(df)

names(df) <- paste0("pmcode_", substr(names(df), 1, 5))

# source 2 and 3
invisible(lapply(names(df), function(u) {
assign(u, df[[u]])
save(list = u, file = paste0("./data/", u, ".RData"))
}))






parameter <- read.csv("Parameter_Distribution.csv",header=TRUE)
FSTDEC <- parameter$FSTDEC #the FSTDEC parameter from the table
KBOD20 <- parameter$KBOD20
x <- readLines("1.txt")
m <- readLines("1.txt")

for(y in 1:length(KBOD20))
{
  for(i in 1:length(FSTDEC))
{
  regexp <- "KBOD20"
  Lineno <- grep(pattern = regexp, x = x , value = F) #find the place that have the parameter
  newline1 <- paste("    1    3     ",KBOD20[y],"         1     0.101      1.15",sep='')#the newline for the replacement
  m[Lineno+2] <- newline1
  regexp <- "FSTDEC"
  Lineno <- grep(pattern = regexp, x = x , value = F)
  newline1 <- paste("    1   17      ",FSTDEC[i],"         1",sep='')
  x[Lineno+2] <- newline1

filename <- paste("Newfiles\\FSTDEC@",FSTDEC[i],".UCI",sep='')
  writelines(x,m,filename) #generate the file with the new file name
  }
}

original1 <- "00041 Weather"
original2 <- "USGS Weather"
original3 <- "54 rows"
original4 <- "3 variables"
original5 <- "pmcode_00041"

replacement1 <- longnames

dfdim <- sapply(df, dim)
dfdim <- data.frame(sapply(df, dim))
dfrows <- dfdim[1]
dfcols <- dfdim[2]


# https://stackoverflow.com/questions/23474239/writing-multiple-lines-to-a-txt-file-using-r?s=4|0.7866

Try the ie2miscdata package in your browser

Any scripts or data that you put into this service are public.

ie2miscdata documentation built on Aug. 21, 2023, 1:07 a.m.