Thank you Bert. I will try your suggestion.
You are ever gracious with your knowledge and expertise.
John
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)
>>> Bert Gunter <gunter.ber...@gene.com> 3/28/2012 5:49 PM >>>
Oh ... you may need to specify a single color to avoid a rainbow:
xyplot( Value ~ Time, Group = SubjectID, data = dat, type = "b", col = "black")
-- Bert
On Wed, Mar 28, 2012 at 2:39 PM, John Sorkin
<jsor...@grecc.umaryland.edu> wrote:
> 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 ...{{dropped:30}}
______________________________________________
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.
