|
|
| Line 19: |
Line 19: |
| | ==Rationale== | | ==Rationale== |
| | | | |
| − | An rcode for managing the data into a more easily usable format. The produced csv file was then converted to zip and uploaded to [[:File:Monica.zip]]. | + | An rcode for managing the data into a more easily usable format. The produced csv file was then converted to zip and uploaded to [[:File:Monica.zip]]. The code that produced the data file is documented on [[:File:Monica.zip]] |
| − | | |
| − | | |
| − | {{hidden|
| |
| − | <pre>
| |
| − | #File: http://www.tilastotiede.fi/ReplicaX/example_MONICA.r
| |
| − | #Author: Juha Karvanen (juha.karvanen at iki.fi)
| |
| − | #
| |
| − | #Summary: This code demonstrates the creation of data replicas with ReplicaX R functions.
| |
| − | #The R functions needed are available as http://www.tilastotiede.fi/ReplicaX/ReplicaX_functions.r
| |
| − | #For information see http://www.tilastotiede.fi/ReplicaX/ReplicaX_esitys.pdf (in Finnish)
| |
| − | #
| |
| − | #The data used is a published sample of WHO MONICA data.
| |
| − | #Before running the code, you should download the data from
| |
| − | #http://www.thl.fi/publications/monica/monograph_cd/data/form04_3.zip
| |
| − | #
| |
| − | #The data documentation can be found at
| |
| − | #http://www.thl.fi/publications/monica/monograph_cd/formats/form048.htm
| |
| − | #http://www.thl.fi/publications/monica/monograph_cd/formats/survey.htm
| |
| − | #
| |
| − | #For more information on the MONICA project see http://www.thl.fi/monica/
| |
| − | | |
| − | | |
| − | print("ReplicaX demonstration with the MONICA data. Running the code takes a couple minutes.")
| |
| − | | |
| − | #rm(list=ls())
| |
| − | set.seed(20131103) #to reproduce the results
| |
| − | #path <- "/home/jutakarv/Documents/replica_gen/"
| |
| − | #path <- "/home/q/tiede/replica_gen/"
| |
| − | #path <- paste(getwd(),"/",sep="")
| |
| − | #source(paste(path,"ReplicaX_functions.r",sep=""))
| |
| − | | |
| − | path <- "//cesium/jtue$/_Downloads/form04_3/" # This was the temporary file on my own computer when doing the conversion.
| |
| − | if(!file.exists(paste(path,"form04_3.txt",sep="")))
| |
| − | {
| |
| − | stop("Data file form04_3.txt does not exist in the current path.
| |
| − | Please download it from http://www.thl.fi/publications/monica/monograph_cd/data/form04_3.zip
| |
| − | or change the path definition")
| |
| − | }
| |
| − | widths <- c(
| |
| − | 2,1,2,2,6,1,3,8,8,1,1,1,1,2,1,3,1,1,4,1,3,2,1,3,1,3,2,1,1,1,1,1,1,1,1,1,2,1,1,3,3,2,3,3,2,1,2,2,4,2,3,3,8,3,3,8,3,4,2,3,4,
| |
| − | 4,4,2,1,2,1,1,1,1,1,4,4,4,4
| |
| − | )
| |
| − | | |
| − | colnames <- c(
| |
| − | "form","versn","centre","runit","serial","numsur","samunit","dexam","mbirth","agegr","sex","marit","edlevel",
| |
| − | "school","cigs","numcigs","daycigs","evercig","stop","iflyear","maxcigs","cigage","cigarsm","cigar","pipesm","pipe",
| |
| − | "othersm","hibp","drugs","bprecd","hich","chdt","chrx","chrecd","asp","menop","agem","horm","pill","syst1",
| |
| − | "diast1","rz1","syst2","diast2","rz2","cuff","arm","bpcoder","timebp","rtemp","chol","choldl","dchol","hdl",
| |
| − | "hdldl","dhdl","scn","cotin","carbmon","height","weight","waist","hip","whcoder","oversion","eage","eageg",
| |
| − | "cohort1","cohort2","edtert1","edtert2","systm","diastm","chola","hdla"
| |
| − | )
| |
| − | | |
| − | monica <- read.fwf(file=paste(path,"form04_3.txt",sep=""),widths=widths,col.names=colnames)
| |
| − | | |
| − | selectedcolnames <- c("marit","edlevel","school","height","weight","waist","hip")
| |
| − | monica2 <- monica[,selectedcolnames]
| |
| − | vartype <- c("discrete","discrete","discrete","continuous","continuous","continuous","continuous")
| |
| − | | |
| − | #decoding missing values
| |
| − | monica2$height[monica2$height==999] <- NA
| |
| − | monica2$weight[monica2$weight==9999] <- NA
| |
| − | monica2$waist[monica2$waist==9999] <- NA
| |
| − | monica2$hip[monica2$hip==9999] <- NA
| |
| − | monica2$school[monica2$school == 99] <- NA
| |
| − | | |
| − | #unit transformations
| |
| − | monica2$weight <- monica2$weight/10 #transform to kg
| |
| − | monica2$waist <- monica2$waist/10 #transform to cm
| |
| − | monica2$hip <- monica2$hip/10 #transform to cm
| |
| − | | |
| − | write.csv(monica, "//cesium/jtue$/_Documents/Monica.csv")
| |
| − | | |
| − | # This csv file was then zipped and uploaded to [[:File:Monica.zip]]
| |
| − | </pre>
| |
| − | }}
| |
| | | | |
| | ==See also== | | ==See also== |
An rcode for managing the data into a more easily usable format. The produced csv file was then converted to zip and uploaded to File:Monica.zip. The code that produced the data file is documented on File:Monica.zip
