I am trying to plot where data points from a give subject are connected by a
line. Each subject is represented by a single row of data. Each subject can have
up to five observations. The first five columns of mydata give the time of
observation, columns 6-10 give the values at each time point. Some subjects have
all data, some are missing values.
The code I wrote to draw the plot is listed below. If you run the code you will
see that (1) The few lines drawn are not drawn on a subject-by-subject basis,
(2) not all lines are drawn.
I hope someone can help me fix my code
Thanks,
John
# Create a matrix to hold our data.
mydata <- matrix(nrow=10,ncol=10)
dimnames(mydata) <- list(NULL,c("time1", "time2", "time3", "time4", "time5",
"clot1", "clot2", "clot3", "clot4", "clot5"))
# load data into the matrix.
# The first five columns represent time of measurement
# Columns 6-10 have the measured values.
mydata[1,]=c( 1, 3, NA,NA, NA, 22, 40, NA, NA, NA)
mydata[2,]=c( 1, NA,NA,NA, NA, 39, NA, NA, NA, NA)
mydata[3,]=c( 1, 6, NA,NA, NA, 47, 43, NA, NA, NA)
mydata[4,]=c( 1, 3, 7, 21, NA, 34, 29, 33, 21, NA)
mydata[5,]=c( NA,NA,NA,NA, NA, NA, NA, NA, NA, NA)
mydata[6,]=c( 1, 3, 7, 19, 35, 42, 60, 86, 56, 39)
mydata[7,]=c( 1, NA,NA,NA, NA, 57, NA, NA, NA, NA)
mydata[8,]=c( 1, 3, 10,17, 25, 71, 55, 57, 62, 47)
mydata[9,]=c( 1, NA,NA,NA, NA, 79, NA, NA, NA, NA)
mydata[10,]=c( 1,3, 9, NA, NA, 60, 64, 68, NA, NA)
# Dispaly data.
mydata
# Plot data.
plot(mydata[,1:5],mydata[,6:10])
# Attempt to connect points for each subject with lines.
lines(mydata[,1:5],mydata[,6:10],type="b")
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/18/GR)
Baltimore, MD 21201-1524
(Phone) 410-605-7119
(Fax) 410-605-7913 (Please call phone number above prior to faxing)
Confidentiality Statement:
This email message, including any attachments, is for th...{{dropped:6}}
______________________________________________
R-help@r-project.org mailing list
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.
