R help - Nov 2010 - Extending a plot in a loop (With attachment)

My problem is that I have a data set for every day of measurement in a
seperate file and I want to plot one parameter of the data for all the days
in one graph. I tried to use for loop but only the last data remains in the
program memory, I don`t know how to plot each day`s data continusly after
the others(or how to extending the x axis.) Would you please help me with
it?

This a plot for one day:

radiation.data
<-read.table("C:/updated_CFL_Rad_files/2008/RAD_2008_JD101_0410.dat",
header
= TRUE,sep = ",", quote = " ", dec =
".")
> attach(radiation.data)
The following object(s) are masked from 'radiation.data (position 3)':

    Batt_avg, Batt_st, Day, Hour, Kdown_avg, Kdown_st, LW_in, LW_in_st,
    Minute, Month, PanelT_avg, PanelT_st, PAR_avg, PAR_st, Sec,
    Tcase_avg, Tcase_st, Tdome_avg, Tdome_st, Thermopile_avg,
    Thermopile_st, Tuv_avg, Tuv_st, Uva_avg, Uva_st, Uvb_avg, Uvb_st,
    Year
> names(radiation.data)
 [1] "Year"           "Month"          "Day"      
"Hour"
 [5] "Minute"         "Sec"            "Batt_avg" 
"PanelT_avg"
 [9] "Batt_st"        "PanelT_st"      "Kdown_avg"
"Thermopile_avg"
[13] "Tcase_avg"      "Tdome_avg"      "LW_in"    
"PAR_avg"
[17] "Tuv_avg"        "Uvb_avg"        "Uva_avg"  
"Kdown_st"
[21] "Thermopile_st"  "Tcase_st"       "Tdome_st" 
"LW_in_st"
[25] "PAR_st"         "Tuv_st"         "Uvb_st"   
"Uva_st"

 plot(((PAR_avg*0.216)/Uvb_avg),
main="Par/UVB",xlab="minutes",ylab="Par/UVB")


and this is the algorithm I tried  for plotting all the data in one plot:

x<- matrix( list.files("C:/updated_CFL_Rad_files",full=TRUE)) #
putting all
data sets in a matrix
  for(i in 1:100) {
      if(i < 101) next
     radiation.data <-read.table(x[i], header = TRUE,sep = ",",
quote = " ",
dec = ".")
 attach(radiation.data)
 plot(i*Hour*60+Minute,PAR_avg,main="PAR",xlab="Hour",ylab="Par")
dev.print(device=postscript, "C:/graph5.eps", onefile=FALSE,
horizontal=FALSE)
       }
The plot I see is the last file's plot, I don't know how to keep
previous
data and continue within the same plot.

* I have attached a sample of the data to this email.

-- 
Sincerely

Nasrin  Pak







-- 
Sincerely

Nasrin  Pak

Tena koe Nasrin

Try points() instead of plot() in your second and subsequent calls to plot(). 
points() and lines() adds to the current plot by default.  Of course you may
have difficulties with setting the x and y limits by that's another matter.

HTH ....

Peter Alspach
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of Nasrin Pak
> Sent: Wednesday, 17 November 2010 7:24 a.m.
> Cc: r-help at r-project.org
> Subject: [R] Extending a plot in a loop (With attachment)
> 
> My problem is that I have a data set for every day of measurement in a
> seperate file and I want to plot one parameter of the data for all the
> days in one graph. I tried to use for loop but only the last data
> remains in the program memory, I don`t know how to plot each day`s data
> continusly after the others(or how to extending the x axis.) Would you
> please help me with it?
> 
> This a plot for one day:
> 
> radiation.data
>
<-read.table("C:/updated_CFL_Rad_files/2008/RAD_2008_JD101_0410.dat",
> header = TRUE,sep = ",", quote = " ", dec =
".")
> > attach(radiation.data)
> The following object(s) are masked from 'radiation.data (position
3)':
> 
>     Batt_avg, Batt_st, Day, Hour, Kdown_avg, Kdown_st, LW_in, LW_in_st,
>     Minute, Month, PanelT_avg, PanelT_st, PAR_avg, PAR_st, Sec,
>     Tcase_avg, Tcase_st, Tdome_avg, Tdome_st, Thermopile_avg,
>     Thermopile_st, Tuv_avg, Tuv_st, Uva_avg, Uva_st, Uvb_avg, Uvb_st,
>     Year
> > names(radiation.data)
>  [1] "Year"           "Month"          "Day" 
"Hour"
>  [5] "Minute"         "Sec"           
"Batt_avg"       "PanelT_avg"
>  [9] "Batt_st"        "PanelT_st"     
"Kdown_avg"
> "Thermopile_avg"
> [13] "Tcase_avg"      "Tdome_avg"     
"LW_in"          "PAR_avg"
> [17] "Tuv_avg"        "Uvb_avg"       
"Uva_avg"        "Kdown_st"
> [21] "Thermopile_st"  "Tcase_st"      
"Tdome_st"       "LW_in_st"
> [25] "PAR_st"         "Tuv_st"        
"Uvb_st"         "Uva_st"
> 
>  plot(((PAR_avg*0.216)/Uvb_avg),
> main="Par/UVB",xlab="minutes",ylab="Par/UVB")
> 
> 
> and this is the algorithm I tried  for plotting all the data in one
> plot:
> 
> x<- matrix( list.files("C:/updated_CFL_Rad_files",full=TRUE))
# putting
> all data sets in a matrix
>   for(i in 1:100) {
>       if(i < 101) next
>      radiation.data <-read.table(x[i], header = TRUE,sep =
",", quote > " ", dec = ".")
>  attach(radiation.data)
> 
plot(i*Hour*60+Minute,PAR_avg,main="PAR",xlab="Hour",ylab="Par")
> dev.print(device=postscript, "C:/graph5.eps", onefile=FALSE,
> horizontal=FALSE)
>        }
> The plot I see is the last file's plot, I don't know how to keep
> previous data and continue within the same plot.
> 
> * I have attached a sample of the data to this email.
> 
> --
> Sincerely
> 
> Nasrin  Pak
> 
> 
> 
> 
> 
> 
> 
> --
> Sincerely
> 
> Nasrin  Pak
The contents of this e-mail are confidential and may be subject to legal
privilege.
 If you are not the intended recipient you must not use, disseminate, distribute
or
 reproduce all or any part of this e-mail or attachments.  If you have received
this
 e-mail in error, please notify the sender and delete all material pertaining to
this
 e-mail.  Any opinion or views expressed in this e-mail are those of the
individual
 sender and may not represent those of The New Zealand Institute for Plant and
 Food Research Limited.