HI Irucka,
Try this:
temp1<-lapply(temp,function(x) x[complete.cases(x),])
temp2<-lapply(temp1,function(tempNew) lapply(tempNew[,-1],function(x)
{x1<-cbind(CYEAR_DECIMAL=tempNew[,1],x)}))
pdf("Irucka2.pdf")
par(mfrow=c(1,2))
lapply(names(temp2),function(i) lapply(temp2[[i]],function(x)
{plot(x[,1],x[,2],main="Fluxmaster versus EGRET/WRTDS \n Seasonal FLux
Sum",sub=i,xlab="Calendar Year Timesteps",ylab="Total Flux
(kg/season)");lines(x[,1],x[,2])}))
dev.off()
A.K.
________________________________
From: Irucka Embry <iruc...@mail2world.com>
To: smartpink...@yahoo.com
Sent: Wednesday, March 6, 2013 5:24 PM
Subject: Re: [R] multiple plots and looping assistance requested (revised codes)
Hi Arun, thanks for the note. Thank you especially for noting the use of the
";" and "{}." I have updated my own code and the possible reproducible example
with those changes.
The attempt at the reproducible code previously was incorrect so I am borrowing
from the data stated in the first comment.
# attempt at reproducible code using data from
http://stackoverflow.com/questions/11548368/making-multiple-plots-in-r-from-one-textfile
> dput(temper)
structure(list(`:Bostoncitydata` = structure(list(Month = 1:3,
Data1 = c(1.5, 12.3, 11.4), Data2 = c(9.1342, 12.31, 3.5),
Data3 = c(8.1231, 1.129, 45.4321)), .Names = c("Month", "Data1",
"Data2", "Data3"), class = "data.frame", row.names = c(NA, -3L
)), `:Chicagocitydata` = structure(list(Month = 1:3, Data1 = c(1.5,
12.3, 0.34), Data2 = c(9.1342, 12.31, 9.0021), Data3 = c(8.1231,
1.129, 24.98)), .Names = c("Month", "Data1", "Data2", "Data3"
), class = "data.frame", row.names = c(NA, -3L))), .Names =
c(":Bostoncitydata",
":Chicagocitydata"))
temper = read.funkyfile("data.txt", "Header", header=TRUE, sep="\t")
tempernow <- lapply(names(temper),function(i) {plot(as.matrix(temp[[i]][1]),
as.matrix(temp[[i]][-1]), main="Fluxmaster versus EGRET/WRTDS \n Seasonal Flux
Sum", sub = i, xlab="Calendar Year Timesteps", ylab="Total Flux (kg/season)");
lines(temp[[i]][1], temp[[i]][-1])})
# below is the revised code that I am trying to implement
> dput(temp)
structure(list(`02143500` = structure(list(Calendar_Year_Timesteps =
c(2000.875,
2001.125, 2001.375, 2001.625, 2001.875, 2002.125, 2002.375, 2002.625,
2002.875, 2003.125, 2003.375, 2003.625, 2003.875, 2004.125, 2004.375,
2004.625, NA), load_00600_W = c(227675.73764, 92777.682029, 84827.680295,
193298.65669, 170799.05034, 103666.8759, 107485.71333, 213765.87505,
472307.65662, 799500.99994, 754868.43185, 454078.02653, 171521.77777,
265827.90007, 120401.25989, 194000.26057, NA), load_00600_F = c(4202.7437226,
16214.840538, 7371.9290209, 3114.1090754, 2464.1114951, 9380.6352081,
3859.2809055, 901.86146915, 22377.413599, 53563.26564, 148264.35049,
103538.36278, 18142.045363, 14672.031667, 18796.93618, 75313.330193,
NA)), .Names = c("Calendar_Year_Timesteps", "load_00600_W", "load_00600_F"
), class = "data.frame", row.names = c(NA, -17L)), `02169000` = structure(list(
Calendar_Year_Timesteps = c(2000.875, 2001.125, 2001.375,
2001.625, 2001.875, 2002.125, 2002.375, 2002.625, 2002.875,
2003.125, 2003.375, 2003.625, 2003.875, 2004.125, 2004.375,
2004.625, NA), load_00600_W = c(61152.4563009734, 3087159.69237187,
83271.1072869206, 78182.9098615795, 10344.3815353876, 2246.4450929669,
989.220241325, 2318.3488271707, 24134.6766402373, 34350.75952138,
39041.6413152116, 28188.925323908, 2281.7849286464, 1756.7094755792,
357.3578073865, 238.1185010405, NA), load_00600_F = c(3940.9979246,
14888.612602, 5551.1377638, 5055.93881, 1689.0656165, 1436.5965696,
849.85774965, 879.76270422, 11452.39104, 62140.01855, 67190.304212,
47562.345621, 3535.1744243, 21923.947402, 2059.0406953, 1685.1384101,
NA)), .Names = c("Calendar_Year_Timesteps", "load_00600_W",
"load_00600_F"), class = "data.frame", row.names = c(NA, -17L
)), `02172300` = structure(list(CYEAR_Decimal = c(2001.875, 2002.125,
2002.375, 2002.625, 2002.875, 2003.125, 2003.375, 2003.625, 2003.875,
2004.125, 2004.375, 2004.625, NA), load_00600_W = c(16150.219744,
18927.53716, 15884.341197, 14102.537116, 11611.582491, 14662.106089,
10576.561683, 7699.1718611, 19837.60763, 22396.338619, 34876.278129,
23822.715775, NA), load_00600_F = c(899.46439731, 1667.869248,
952.12246929, 543.12231908, 1663.8731209, 3123.5771201, 4338.1307758,
4469.2046052, 2642.983822, 2863.529214, 2131.6573312, 2476.4154125,
NA)), .Names = c("CYEAR_Decimal", "load_00600_W", "load_00600_F"
), class = "data.frame", row.names = c(NA, -13L)), `21SC60WQ.CW-206` =
structure(list(
CYEAR_Decimal = c(2001.875, 2002.125, 2002.375, 2002.625),
load_00600_W = c(23693281.5300403, 66710782.7389879, 49036165.5298238,
27086129.9582558), load_00600_F = c(145481.6737, 370915.04829,
282387.90945, 179357.3242)), .Names = c("CYEAR_Decimal",
"load_00600_W", "load_00600_F"), class = "data.frame", row.names = c(NA,
-4L))), .Names = c("02143500", "02169000", "02172300", "21SC60WQ.CW-206"
))
temp = read.funkyfile("sitesmore.csv", "station_id", header=TRUE, sep="\t")
# i have successfully graphed like this for one of the columns, but I need to
do this for both columns in all of the sites (station_ids):
plot(as.matrix(temp[[1]][1]), as.matrix(temp[[1]][2]))
lines(as.matrix(temp[[1]][1]), as.matrix(temp[[1]][2]))
# this is the code that I am using to plot through lapply which is still giving
me problems
resnow <- lapply(names(temp),function(i) {plot(as.matrix(temp[[i]][1]),
as.matrix(temp[[i]][-1]), main="Fluxmaster versus EGRET/WRTDS \n Seasonal Flux
Sum", sub = i, xlab="Calendar Year Timesteps", ylab="Total Flux (kg/season)");
lines(temp[[i]][1], temp[[i]][-1])})
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' and 'y' lengths differ
# I want to have side by side graphs (on the same page) of
Calendar_Year_Timesteps as the x axis versus load_00600_W as the y axis &
Calendar_Year_Timesteps as the x axis versus load_00600_P as the y axis for
each of the data frames. The graphs will be a scatterplot matrix with a line
connecting all of the points.
Thank you.
Irucka
<-----Original Message----->
>From: arun [smartpink...@yahoo.com]
>Sent: 3/6/2013 2:18:52 PM
>To: iruc...@mail2world.com
>Subject: Re: [R] multiple plots and looping assistance requested
>
>Hi,
>Not sure whether I understand it correctly: Your example dataset and code
>gives me scatterplot matrix
>
>
>set.seed(15)
>repro1 <- data.frame(stat1 = rnorm(16, 10), stat2 = rnorm(16, 10), stat3 =
>rnorm(16,10), stat4 =
>rnorm(16, 4))
>set.seed(25)
>repro2 <- data.frame(stat1 = rnorm(16, 10), stat2 = rnorm(16, 10), stat3 =
>rnorm(16,10), stat4 =
>rnorm(16, 4))
>
>set.seed(181)
>repro3 <- data.frame(stat1 = rnorm(16, 10), stat2 = rnorm(16, 10), stat3 =
>rnorm(16,10), stat4 =
>rnorm(16, 4))
>
>set.seed(2052)
>repro4 <- data.frame(stat1 = rnorm(16, 10), stat2 = rnorm(16, 10), stat3 =
>rnorm(16,10), stat4 =
>rnorm(16, 4))
>
>reproList <- list(repro1, repro2, repro3, repro4)
>pdf("Irucka.pdf")
>lapply(reproList,function(x){plot(x);lines(x)}) #here I am getting the
>scatterplot.
>dev.off()
>
>#From your code:
> set.seed(15)
> repro <- data.frame(stat1 = rnorm(16, 10), stat2 = rnorm(16, 10), stat3 =
> rnorm(16,10), stat4 = rnorm(16,
>4))
> reproList <- list(repro, repro, repro, repro)
>pdf("Irucka2.pdf")
>lapply(seq_along(repro),function(i)
>{plot(reproList[[i]]);lines(reproList[[i]])}) #your code didn't had closing
>brackets
> dev.off()
>A.K.
>
>
>
>
>
>
>
>________________________________
>From: Irucka Embry <iruc...@mail2world.com>
>To: smartpink...@yahoo.com
>Sent: Wednesday, March 6, 2013 2:15 PM
>Subject: Re: [R] multiple plots and looping assistance requested
>
>
>Hey Arun, thanks.
>
>I did some back end work in LibreOffice Calc with the previous data set so
>what I had mentioned in the
>previous post does not apply here.
>
>Part of the data set is included in the object temp. I am working on the small
>data set first with the code
>and then I'll use the whole data set.
>
>Thank you.
>
>Irucka
>
>
>> str(temp)
>List of 4
>$ 02143500 :'data.frame': 17 obs. of 3 variables:
>..$ Calendar_Year_Timesteps: num [1:17] 2000.875 2001.125 2001.375 2001.625
>2001.875 ... # the
>decimals
>...$ load_00600_W : num [1:17] 227676 92778 84828 193299 170799 ...
>..$ load_00600_F : num [1:17] 4203 16215 7372 3114 2464 ...
>$ 02169000 :'data.frame': 17 obs. of 3 variables:
>..$ Calendar_Year_Timesteps: num [1:17] 2000.875 2001.125 2001.375 2001.625
>2001.875 ...
>..$ load_00600_W : num [1:17] 61152 3087160 83271 78183 10344 ...
>..$ load_00600_F : num [1:17] 3941 14889 5551 5056 1689 ...
>$ 02172300 :'data.frame': 13 obs. of 3 variables:
>..$ CYEAR_Decimal: num [1:13] 2001.875 2002.125 2002.375 2002.625 2002.875 ...
>..$ load_00600_W : num [1:13] 16150 18928 15884 14103 11612 ...
>..$ load_00600_F : num [1:13] 899 1668 952 543 1664 ...
>$ 21SC60WQ.CW-206:'data.frame': 4 obs. of 3 variables:
>..$ CYEAR_Decimal: num [1:4] 2001.875 2002.125 2002.375 2002.625
>..$ load_00600_W : num [1:4] 23693282 66710783 49036166 27086130
>..$ load_00600_F : num [1:4] 145482 370915 282388 179357
>
>
>temp = read.funkyfile("sitesmore.csv", "station_id", header=TRUE, sep="\t")
>resnow <- lapply(names(temp),function(i) plot(as.matrix(temp[[i]][1]),
>as.matrix(temp[[i]][-1]),
>main="Fluxmaster versus EGRET/WRTDS \n Seasonal Flux Sum", sub = i,
>xlab="Calendar Year Timesteps",
>ylab="Total Flux (kg/season)") lines(temp[[i]][1], temp[[i]][-1]))
>Error in xy.coords(x, y, xlabel, ylabel, log) :
>'x' and 'y' lengths differ
># I want to have side by side graphs (on the same page) of
>Calendar_Year_Timesteps as the x axis versus
>load_00600_W as the y axis & Calendar_Year_Timesteps as the x axis versus
>load_00600_P as the y axis
>for each of the data frames.
>
>
># attempt at a partial reproducible code
>repro <- as.data.frame(stat1 = rnorm(16, 10), stat2 = rnorm(16, 10), stat3 =
>rnorm(16,10), stat4 =
>rnorm(16, 4))
>reproList <- list(repro, repro, repro, repro)
>repronow <- lapply(seq_along(repro) ,function(i) plot(reproList[[i]]
>lines(reproList[[i]])))
>
>
>
>
>
><-----Original Message----->
>>From: arun [smartpink...@yahoo.com]
>>Sent: 3/6/2013 12:53:17 PM
>>To: iruc...@mail2world.com
>>Subject: Re: [R] multiple plots and looping assistance requested
>>
>>
>>
>>HI Irucka,
>>I am fine.
>>Thanks.
>>
>>In the link,I couldn't find the dataset. Or is it the one from the previous
>>email.
>>Arun
>>
>>
>>
>>
>>________________________________
>>From: Irucka Embry <iruc...@mail2world.com>
>>To: smartpink...@yahoo.com
>>Sent: Wednesday, March 6, 2013 1:46 PM
>>Subject: Re: [R] multiple plots and looping assistance requested
>>
>>
>>Hi Arun, how are you?
>>
>>I have revised my previous 2nd question and created this post here:
>>http://r.789695.n4.nabble.com/trouble-with-lapply-plot-labels-and-indexing-in-multiple-plots-
>tp4660512.html
>>
>>Thank-you.
>>
>>Irucka
>>
>><-----Original Message----->
>>>From: arun [smartpink...@yahoo.com]
>>>Sent: 3/5/2013 8:26:10 PM
>>>To: iruc...@mail2world.com
>>>Subject: Re: [R] multiple plots and looping assistance requested
>>>
>>>HI Irucka,
>>>No problem.
>>>
>>>What's your second question?
>>>Arun _______________________________________________________________
Get the Free email that has everyone talking at http://www.mail2world.com
Unlimited Email Storage – POP3 – Calendar – SMS – Translator – Much More!
Irucka2.pdf
Description: Adobe PDF document
______________________________________________ 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.
