Hi: Thanks to Dennis and Fernando for your help reordering the levels.
Now I have a different issue:
I am trying to get the?cumulative weekly values using cumsum and it
appears to output the wrong?values. Here's my dataset:
winter <-? structure(list(week = c(26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L,
34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L,
47L, 48L, 49L, 50L, 51L, 52L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L,
9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L,
22L, 23L, 24L, 25L), BY2009 = c(0L, 95L, 436L, 809L, 3668L, 9437L,
42700L, 141135L, 486474L, 1109095L, 1459053L, 1990285L, 2643088L,
2988446L, 3131437L, 3280320L, 4237821L, 4270316L, 4285178L, 4292388L,
4302153L, 4308999L, 4314066L, 4320267L, 4323312L, 4429542L, 4430725L,
4430920L, 4437775L, 4448036L, 4452215L, 4452865L, 4453105L, 4453270L,
4454459L, 4455322L, 4455884L, 4457234L, 4457422L, 4457901L, 4458671L,
4458671L, 4459229L, 4459543L, 4459690L, NA, NA, NA, NA, NA, NA,
NA), BY2008 = c(0L, 0L, 77L, 201L, 360L, 2736L, 12216L, 29530L,
129104L, 452783L, 650994L, 744624L, 873807L, 985627L, 1056887L,
1092128L, 1106148L, 1126926L, 1148620L, 1163636L, 1177036L, 1186223L,
1189830L, 1192634L, 1195051L, 1200342L, 1216558L, 1217456L, 1218014L,
1218723L, 1219695L, 1219756L, 1220128L, 1223214L, 1233322L, 1237617L,
1238499L, 1241092L, 1241128L, 1241361L, 1241604L, 1241604L, 1241674L,
1241946L, 1242254L, 1242388L, 1242428L, 1242428L, 1242428L, 1242428L,
1242428L, 1242428L), BY2007 = c(0L, 0L, 0L, 10775L, 14941L, 19899L,
36120L, 65521L, 100472L, 133660L, 279704L, 384711L, 570008L,
729690L, 937227L, 1077921L, 1206196L, 1250470L, 1277549L, 1296477L,
1306914L, 1336898L, 1355293L, 1381139L, 1385712L, 1417707L, 1421386L,
1422429L, 1432065L, 1433589L, 1434416L, 1441425L, 1441658L, 1442091L,
1442091L, 1443194L, 1443451L, 1443579L, 1443645L, 1443715L, 1444250L,
1444679L, 1444776L, 1444776L, 1444776L, 1444776L, 1444776L, 1444776L,
1444776L, 1444776L, 1444776L, 1444776L), BY2006 = c(0L, 707L,
3097L, 11957L, 36387L, 77272L, 150064L, 355585L, 700078L, 1363051L,
1889460L, 2521413L, 3371904L, 4214582L, 4660569L, 5218721L, 5550753L,
5725079L, 5805680L, 5854376L, 5952947L, 6056510L, 6205979L, 6284060L,
6466538L, 6468696L, 6485262L, 6489289L, 6491944L, 6493056L, 6493623L,
6496218L, 6501194L, 6507530L, 6507824L, 6509582L, 6509873L, 6510076L,
6510526L, 6511624L, 6512412L, 6512607L, 6513139L, 6513139L, 6513139L,
6513139L, 6513139L, 6513139L, 6513139L, 6513306L, 6513306L, 6513306L
), BY2005 = c(0L, 0L, 868L, 2912L, 6976L, 13025L, 22424L, 35728L,
80900L, 323055L, 799919L, 1512453L, 2570862L, 4685880L, 6196222L,
7334435L, 7667627L, 7826447L, 7920826L, 8269708L, 8308998L, 8338699L,
8385957L, 8455794L, 8463678L, 8513016L, 8535184L, 8554581L, 8573978L,
8589962L, 8598650L, 8599850L, 8601473L, 8602764L, 8604120L, 8604827L,
8605702L, 8606577L, 8606799L, 8607682L, 8607682L, 8607682L, 8607811L,
8607811L, 8607811L, 8607811L, 8607811L, 8607811L, 8607811L, 8607811L,
8607811L, 8607811L)), .Names = c("week", "BY2009",
"BY2008",
"BY2007", "BY2006", "BY2005"), class =
"data.frame", row.names = c(NA,
-52L))
apply(winter[,2:6],2,cumsum)
?
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California, USA
It is working correctly, You have NAs in your data, so you get NAs in the cumsum. On Mon, May 10, 2010 at 5:09 PM, Felipe Carrillo <mazatlanmexico@yahoo.com>wrote:> Hi: Thanks to Dennis and Fernando for your help reordering the levels. > Now I have a different issue: > I am trying to get the cumulative weekly values using cumsum and it > appears to output the wrong values. Here's my dataset: > > winter <- structure(list(week = c(26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L, > 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, > 47L, 48L, 49L, 50L, 51L, 52L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, > 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, > 22L, 23L, 24L, 25L), BY2009 = c(0L, 95L, 436L, 809L, 3668L, 9437L, > 42700L, 141135L, 486474L, 1109095L, 1459053L, 1990285L, 2643088L, > 2988446L, 3131437L, 3280320L, 4237821L, 4270316L, 4285178L, 4292388L, > 4302153L, 4308999L, 4314066L, 4320267L, 4323312L, 4429542L, 4430725L, > 4430920L, 4437775L, 4448036L, 4452215L, 4452865L, 4453105L, 4453270L, > 4454459L, 4455322L, 4455884L, 4457234L, 4457422L, 4457901L, 4458671L, > 4458671L, 4459229L, 4459543L, 4459690L, NA, NA, NA, NA, NA, NA, > NA), BY2008 = c(0L, 0L, 77L, 201L, 360L, 2736L, 12216L, 29530L, > 129104L, 452783L, 650994L, 744624L, 873807L, 985627L, 1056887L, > 1092128L, 1106148L, 1126926L, 1148620L, 1163636L, 1177036L, 1186223L, > 1189830L, 1192634L, 1195051L, 1200342L, 1216558L, 1217456L, 1218014L, > 1218723L, 1219695L, 1219756L, 1220128L, 1223214L, 1233322L, 1237617L, > 1238499L, 1241092L, 1241128L, 1241361L, 1241604L, 1241604L, 1241674L, > 1241946L, 1242254L, 1242388L, 1242428L, 1242428L, 1242428L, 1242428L, > 1242428L, 1242428L), BY2007 = c(0L, 0L, 0L, 10775L, 14941L, 19899L, > 36120L, 65521L, 100472L, 133660L, 279704L, 384711L, 570008L, > 729690L, 937227L, 1077921L, 1206196L, 1250470L, 1277549L, 1296477L, > 1306914L, 1336898L, 1355293L, 1381139L, 1385712L, 1417707L, 1421386L, > 1422429L, 1432065L, 1433589L, 1434416L, 1441425L, 1441658L, 1442091L, > 1442091L, 1443194L, 1443451L, 1443579L, 1443645L, 1443715L, 1444250L, > 1444679L, 1444776L, 1444776L, 1444776L, 1444776L, 1444776L, 1444776L, > 1444776L, 1444776L, 1444776L, 1444776L), BY2006 = c(0L, 707L, > 3097L, 11957L, 36387L, 77272L, 150064L, 355585L, 700078L, 1363051L, > 1889460L, 2521413L, 3371904L, 4214582L, 4660569L, 5218721L, 5550753L, > 5725079L, 5805680L, 5854376L, 5952947L, 6056510L, 6205979L, 6284060L, > 6466538L, 6468696L, 6485262L, 6489289L, 6491944L, 6493056L, 6493623L, > 6496218L, 6501194L, 6507530L, 6507824L, 6509582L, 6509873L, 6510076L, > 6510526L, 6511624L, 6512412L, 6512607L, 6513139L, 6513139L, 6513139L, > 6513139L, 6513139L, 6513139L, 6513139L, 6513306L, 6513306L, 6513306L > ), BY2005 = c(0L, 0L, 868L, 2912L, 6976L, 13025L, 22424L, 35728L, > 80900L, 323055L, 799919L, 1512453L, 2570862L, 4685880L, 6196222L, > 7334435L, 7667627L, 7826447L, 7920826L, 8269708L, 8308998L, 8338699L, > 8385957L, 8455794L, 8463678L, 8513016L, 8535184L, 8554581L, 8573978L, > 8589962L, 8598650L, 8599850L, 8601473L, 8602764L, 8604120L, 8604827L, > 8605702L, 8606577L, 8606799L, 8607682L, 8607682L, 8607682L, 8607811L, > 8607811L, 8607811L, 8607811L, 8607811L, 8607811L, 8607811L, 8607811L, > 8607811L, 8607811L)), .Names = c("week", "BY2009", "BY2008", > "BY2007", "BY2006", "BY2005"), class = "data.frame", row.names = c(NA, > -52L)) > apply(winter[,2:6],2,cumsum) > > > Felipe D. Carrillo > Supervisory Fishery Biologist > Department of the Interior > US Fish & Wildlife Service > California, USA > > > > > ______________________________________________ > R-help@r-project.org mailing list > stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > R-project.org/posting-guide.html<r-project.org/posting-guide.html> > and provide commented, minimal, self-contained, reproducible code. >-- Jim Holtman Cincinnati, OH +1 513 646 9390 What is the problem that you are trying to solve? [[alternative HTML version deleted]]