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
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
>
http://www.R-project.org/posting-guide.html<http://www.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]]