
  Dear R-experts,
   
  I have a dataset of 4 patients and each patient has many records at four different time points. I have done 4 different qqnorm plots on the same graph where each plot represents the records of one patient at each time point. I would like to do the same graph for the remaining patiens, but instead of repeating the below procedure three more times I would like to have a for loop so when I run the loop I would get four graphs of the four patients at the same time, where each graph has 4 different qqnorm plots representing the data at  four time points for each patient. I tried to do the loop but couldn't make it work. 
   
  below is the code used to the graph for one patient. sli-4 - sli_7 contain records at four time points. dat=dataset, y=records, Slide=time point
   
   
  ### Patient 200004
  
  sli_4=dat$y[dat$Slide==4 & dat$Control==0]
  sli_5=dat$y[dat$Slide==5 & dat$Control==0]
  sli_6=dat$y[dat$Slide==6 & dat$Control==0]
  sli_7=dat$y[dat$Slide==7 & dat$Control==0]
   
  ### qq-plot of patient 200004
  q_sli4<-qqnorm(sli_4,plot.it=FALSE)
  q_sli5<-qqnorm(sli_5,plot.it=FALSE)
  q_sli6<-qqnorm(sli_6,plot.it=FALSE)
  q_sli7<-qqnorm(sli_7,plot.it=FALSE)
   
  plot(range(q_sli4,q_sli5,q_sli6,q_sli7),range(q_sli4,q_sli5,q_sli6,q_sli7),type="n",
  xlab = "Theoretical Quantiles",col.main="blue",
  main = "Normal Q-Q Plot of index for patient 200004",ylab = "Sample Quantiles")
  points(q_sli4,col=4,pch=0,cex=1)
  points(q_sli5,col=3,pch=1,cex=1)
  points(q_sli6,col=2,pch=2,cex=1)
  points(q_sli7,col=1,pch=3,cex=1)
  legend("topleft",c("Day 0", "56 days","112 days", "252 days"),col=c(4,3,2,1),
  text.col=c(4,3,2,1),pch=c(0,1,2,3),bg="bisque")
  abline(0,0)
   
   
  Thanks alot for your help,
   
  All the bests,
  Jenny
   
    
---------------------------------
  


 		
---------------------------------


	[[alternative HTML version deleted]]

