I realized that after I hit the send button...I was looking at the wrong
dataset,,,,,Brrrrrr
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California, USA
>
>From: jim holtman <jholt...@gmail.com>
>To: Felipe Carrillo <mazatlanmex...@yahoo.com>
>Cc: r-h...@stat.math.ethz.ch
>Sent: Mon, May 10, 2010 2:33:03 PM
>Subject: Re: [R] cumsum output
>
>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 <mazatlanmex...@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
>>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]]
______________________________________________
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.
