The pipe-delimited text files can be read using the read.table() function in base R. R is powerful and can provide a lot of useful information with very little effort (as demonstrated in the instructions below). Note: These instructions were tested using R 3.1.2 installed on a Linux platform.
If R (and RStudio) not already installed on your local machine, you may install from the following links:
Assuming you have installed/launched R, and have downloaded/unzipped a set of AACT pipe-delimited files...
When getting started, first try reading a smaller file (e.g., ‘id_information.txt’). This may help troubleshoot minor issues. Example R code:
d1 <- read.table(file = "id_information.txt", header = TRUE, sep = "|", na.strings = "", comment.char = "", quote = "\"", fill = FALSE, nrows = 200000) summary(d1) # print a summary of dataframe d1
R should respond with something like:
Notes:
Once you have successfully read a small file, try reading the file ‘designs.txt’. This file is larger, but still contains just simple structured data elements (no free text). Reading this file may help test memory requirements on your system, and is a good file on which to develop strategies for reading larger files. (To try this out, replace id_information.txt with designs.txt in the previous code sample and rerun in R.)
R should respond with something a bit more interesting like:
Example code:
d1 <- read.table(file = "designs.txt", header = TRUE, na.strings = "", stringsAsFactors = FALSE, comment.char = "", quote = "\"", fill = FALSE, nrows = 1000)
Notes:
The following code can be used to estimate the approximate memory needed to read the full dataset
sz <- object.size(d0) # determine approximate memory used to store d0 # multiply by (nrows for full data set / 1000) # to estimate total memory needed to store the full data print( sz * 250000 / 1000, units="Mb") # Estimate that ~50 Mb of memory will be needed. # This is fine since my system has > 4 Gb of RAM.
d1 <- read.table(file = "designs.txt", header = TRUE, sep = "|", na.strings = "", colClasses = c("integer", rep("character",14) ), comment.char = "", quote = "\"", fill = FALSE, nrows = 250000) # a bit larger than needed
Example 1: Read rows 1000 to 1999 in the data file:
d1 <- read.table("designs.txt", skip = 999, # skip first 999 records nrows = 1000, # read the next 1000 records header = FALSE, # turn off header row col.names = names(d0), # use column names determined from d0 na.strings = "", sep = "|", comment.char = "", quote = "\"", fill = FALSE, nrows = 250000) # a bit larger than needed
To identify row number(s) containing string "NCT00000558"...
grep("NCT00000558", readLines("designs.txt")) [1] 222688
The above reveals that NCT00000558 occurs in line 222,688 of the designs.txt file.
# Read just the desired lines with read.table using skip and nrows d1 <- read.table("designs.txt", skip = 214443, # skip first 214443 records nrows = 1, # read 1 record header = FALSE, # turn off header row col.names = names(d0), # use column names determined from d0 na.strings = "", sep = "|", colClasses = c("integer", rep("character",14) ), comment.char = "", quote = "\"", fill = FALSE)
Once you have successfully read the two files mentioned above, you are ready to try reading some other files. The file ‘studies.txt’ contains one record per study (NCT_ID is the unique identifier), and contains many of the study protocol data elements.
The following code reads the first 1000 records from the file ‘studies.txt’.
d0 <- read.table(file = "studies.txt", header = TRUE, sep = "|", na.strings = "", stringsAsFactors = FALSE, comment.char = "", quote = "\"", fill = FALSE, nrows = 1000)
Check memory requirements, assuming ~ 250000 rows: sz <- object.size(d0) print( sz * 250000 / 1000, units="Gb")
Reading the full dataset will require about 0.2 Gb of memory.
# Specify appropriate variable types based on d0 and Data Dictionary Classes = c( "character", rep("Date",3), rep("character",6), "integer", rep("character",6), rep("integer",2), rep("character",11), "Date", rep("character",6) )
# read full dataset of studies d1 <- read.table(file = "studies.txt", header = TRUE, sep = "|", na.strings = "", colClasses = Classes, comment.char = "", nrows = 250000, # a bit larger than needed quote = "\"", fill = FALSE)