Actually, I dput(winter) and is not that big.....so here is the whole thing I
am trying to accomplish.
winter <- structure(list(IDWeek = c(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, 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, 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, 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, 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, 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, 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, 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, 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, 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, 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), Total = c(0L, 0L, 326L, 1735L,
1807L, 2208L, 3883L, 8820L, 6060L, 19326L, 63158L, 100718L, 53015L,
91689L, 152629L, 122708L, 61293L, 15574L, 86538L, 75365L, 303259L,
19691L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 161L, 321L, 1000L, 4425L,
13202L, 19726L, 30518L, 84949L, 157260L, 145691L, 85801L, 62044L,
44439L, 23272L, 22391L, 20159L, 14854L, 35379L, 31142L, 7736L,
13221L, 4894L, 0L, 460L, 1129L, 1757L, 11105L, 27730L, 55769L,
12298L, 32964L, 80110L, 160507L, 74134L, 101736L, 138355L, 107873L,
53806L, 194429L, 7596L, 55497L, 30336L, 18392L, 22291L, 16850L,
16567L, 14799L, 7909L, 93L, 329L, 363L, 2762L, 5573L, 33466L,
96762L, 340947L, 612835L, 347972L, 526506L, 656981L, 344181L,
145465L, 150543L, 935599L, 32789L, 15092L, 7288L, 9889L, 6945L,
5138L, 6299L, 3095L, 108227L, 1191L, 0L, 0L, 132L, 227L, 2157L,
7558L, 16459L, 87220L, 321195L, 209299L, 110807L, 112024L, 128058L,
74373L, 41531L, 15111L, 21890L, 15060L, 22538L, 12443L, 11052L,
3907L, 2640L, 2745L, 3936L, 17593L, 0L, 0L, 10775L, 4166L, 4958L,
16221L, 29401L, 34951L, 33188L, 146044L, 105007L, 185297L, 159682L,
207537L, 140694L, 128275L, 44274L, 27079L, 18928L, 10437L, 29984L,
18395L, 25846L, 4573L, 31995L, 3679L, 707L, 2390L, 8860L, 24430L,
40885L, 72792L, 205521L, 344493L, 662973L, 526409L, 631953L,
850491L, 842678L, 445987L, 558152L, 332032L, 174326L, 80601L,
48696L, 98571L, 103563L, 149469L, 78081L, 182478L, 2158L, 16566L,
0L, 868L, 2044L, 4064L, 6049L, 9399L, 13304L, 45172L, 242155L,
476864L, 712534L, 1058409L, 2115018L, 1510342L, 1138213L, 333192L,
158820L, 94379L, 348882L, 39290L, 29701L, 47258L, 69837L, 7884L,
49338L, 22168L, 0L, 147L, 1231L, 3216L, 5021L, 18462L, 37441L,
64669L, 136590L, 338523L, 559448L, 714837L, 358037L, 180910L,
291943L, 222708L, 163801L, 39109L, 23247L, 15726L, 38701L, 5795L,
12509L, 18721L, 822L, 6674L, 765L, 1753L, 6226L, 9727L, 22033L,
50472L, 67863L, 100909L, 276699L, 239609L, 416465L, 845073L,
1071051L, 571656L, 409021L, 451331L, 167006L, 68069L, 250439L,
104641L, 76062L, 35261L, 55883L, 16064L, 11058L, 8872L, 9346L,
24760L, 41699L, 40032L, 80576L, 152208L, 171888L, 292862L, 522669L,
441518L, 858538L, 906010L, 1464311L, 1334433L, 542204L, 211514L,
64507L, 94830L, 100663L, 116931L, 83139L, 46268L, 14455L, 22811L,
47446L, 54224L), Fry = c(0L, 0L, 326L, 1735L, 1807L, 2208L, 3883L,
8759L, 6060L, 19326L, 63119L, 100524L, 52582L, 88170L, 145564L,
111416L, 38233L, 5248L, 17826L, 11038L, 34008L, 215L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 161L, 321L, 1000L, 4425L, 13055L, 19488L,
30518L, 84818L, 156909L, 144786L, 84207L, 57720L, 31049L, 6858L,
1616L, 719L, 364L, 49L, 0L, 0L, 0L, 0L, 0L, 460L, 1129L, 1757L,
11105L, 27730L, 55769L, 12109L, 32638L, 79547L, 158892L, 72867L,
96249L, 130695L, 88837L, 32531L, 65746L, 1885L, 4188L, 680L,
100L, 109L, 0L, 0L, 0L, 0L, 93L, 329L, 363L, 2762L, 5573L, 33466L,
96694L, 339394L, 611967L, 345807L, 522991L, 645886L, 327353L,
118537L, 109158L, 419571L, 6691L, 456L, 133L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 132L, 227L, 2157L, 7558L, 16459L, 87220L,
320684L, 207921L, 110221L, 110021L, 123153L, 63829L, 23982L,
5090L, 4183L, 374L, 439L, 78L, 67L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
10775L, 4166L, 4958L, 16221L, 29401L, 34851L, 32933L, 145103L,
105007L, 182715L, 153781L, 201622L, 129968L, 87997L, 17235L,
3401L, 1343L, 397L, 316L, 37L, 0L, 0L, 0L, 0L, 707L, 2390L, 8860L,
24430L, 40885L, 72792L, 205521L, 343593L, 658875L, 523360L, 627333L,
847147L, 837189L, 429016L, 523436L, 304785L, 125005L, 28047L,
5141L, 6738L, 2255L, 1065L, 0L, 0L, 0L, 0L, 0L, 868L, 2044L,
4064L, 6049L, 9399L, 13304L, 44961L, 241575L, 474250L, 706084L,
1045237L, 2083214L, 1483109L, 1091328L, 284964L, 96754L, 33605L,
74341L, 1795L, 828L, 36L, 0L, 0L, 0L, 0L, 0L, 147L, 1231L, 3216L,
5021L, 18462L, 37441L, 64508L, 136280L, 337493L, 556858L, 710292L,
353191L, 174072L, 265099L, 166694L, 85538L, 6630L, 2033L, 623L,
196L, 0L, 0L, 0L, 0L, 0L, 765L, 1753L, 6226L, 9727L, 22033L,
50472L, 67863L, 100309L, 273113L, 237617L, 411067L, 838077L,
1059358L, 555461L, 364550L, 387386L, 121865L, 26477L, 25836L,
9169L, 2526L, 368L, 0L, 0L, 0L, 0L, 9346L, 24760L, 41699L, 40032L,
80576L, 152208L, 171844L, 290962L, 519078L, 438136L, 847834L,
891191L, 1445079L, 1298580L, 482887L, 151660L, 27477L, 20607L,
11141L, 2768L, 853L, 32L, 0L, 0L, 0L, 0L), Smolt = c(0L, 0L,
0L, 0L, 0L, 0L, 0L, 62L, 0L, 0L, 38L, 195L, 433L, 3518L, 7067L,
11290L, 23058L, 10327L, 68712L, 64328L, 269248L, 19479L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 147L, 238L, 0L, 131L,
351L, 905L, 1592L, 4324L, 13391L, 16414L, 20774L, 19444L, 14491L,
35330L, 31142L, 7736L, 13221L, 4894L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 188L, 326L, 565L, 1614L, 1268L, 5488L, 7658L, 19034L, 21275L,
128681L, 5710L, 51308L, 29654L, 18292L, 22183L, 16850L, 16567L,
14799L, 7909L, 0L, 0L, 0L, 0L, 0L, 0L, 68L, 1553L, 867L, 2164L,
3514L, 11094L, 16829L, 26927L, 41385L, 516027L, 26097L, 14636L,
7159L, 9889L, 6945L, 5138L, 6299L, 3095L, 108227L, 1191L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 511L, 1378L, 585L, 2004L, 4905L,
10548L, 17549L, 10022L, 17709L, 14686L, 22099L, 12365L, 10985L,
3907L, 2640L, 2745L, 3936L, 17593L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
99L, 255L, 943L, 0L, 2581L, 5902L, 5916L, 10725L, 40276L, 27040L,
23678L, 17586L, 10042L, 29668L, 18359L, 25846L, 4573L, 31995L,
3679L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 900L, 4100L, 3050L, 4619L,
3339L, 5490L, 16971L, 34716L, 27247L, 49321L, 52554L, 43558L,
91836L, 101309L, 148404L, 78081L, 182478L, 2158L, 16566L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 211L, 578L, 2614L, 6453L, 13177L, 31805L,
27233L, 46886L, 48228L, 62066L, 60773L, 274542L, 37494L, 28874L,
47223L, 69837L, 7884L, 49347L, 22192L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 160L, 308L, 1028L, 2583L, 4545L, 4847L, 6841L, 26847L, 56014L,
78264L, 32479L, 21216L, 15104L, 38506L, 5795L, 12509L, 18721L,
822L, 6674L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 601L, 3585L, 1992L,
5398L, 6996L, 11693L, 16197L, 44469L, 63946L, 45142L, 41592L,
224601L, 95474L, 73537L, 34893L, 55883L, 16064L, 11058L, 8872L,
0L, 0L, 0L, 0L, 0L, 0L, 44L, 1900L, 3591L, 3382L, 10704L, 14819L,
19232L, 35853L, 59317L, 59854L, 37030L, 74223L, 89522L, 114163L,
82286L, 46236L, 14455L, 22811L, 47446L, 54224L), FryEq = c(0L,
0L, 326L, 1735L, 1807L, 2208L, 3883L, 8864L, 6060L, 19326L, 63185L,
100854L, 53318L, 94151L, 157576L, 130610L, 77432L, 22805L, 134639L,
120393L, 491733L, 33327L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 161L,
321L, 1000L, 4425L, 13306L, 19894L, 30518L, 85042L, 157506L,
146328L, 86914L, 65073L, 53812L, 34763L, 36931L, 33769L, 24998L,
60110L, 52938L, 13149L, 22476L, 8319L, 0L, 460L, 1129L, 1757L,
11105L, 27730L, 55769L, 12430L, 33192L, 80505L, 161636L, 75022L,
105575L, 143717L, 121197L, 68698L, 284505L, 11592L, 91411L, 51093L,
31198L, 37820L, 28645L, 28162L, 25159L, 13442L, 93L, 329L, 363L,
2762L, 5573L, 33466L, 96810L, 342034L, 613438L, 349489L, 528967L,
664745L, 355963L, 164312L, 179515L, 1296819L, 51058L, 25338L,
12300L, 16808L, 11809L, 8736L, 10708L, 5263L, 163296L, 2024L,
0L, 0L, 132L, 227L, 2157L, 7558L, 16459L, 87220L, 321552L, 210259L,
111217L, 113427L, 131491L, 81757L, 53815L, 22127L, 34285L, 25340L,
38007L, 21098L, 18742L, 6640L, 4486L, 4668L, 6690L, 29906L, 0L,
0L, 10775L, 4166L, 4958L, 16221L, 29401L, 35021L, 33366L, 146703L,
105007L, 187104L, 163812L, 211677L, 148205L, 156467L, 63202L,
43653L, 31239L, 17466L, 50751L, 31248L, 43941L, 7774L, 54393L,
6254L, 707L, 2390L, 8860L, 24430L, 40885L, 72792L, 205521L, 345123L,
665842L, 528545L, 635188L, 852828L, 846522L, 457866L, 582452L,
351103L, 208853L, 117390L, 79187L, 162856L, 174479L, 253350L,
132741L, 310215L, 3671L, 28165L, 0L, 868L, 2044L, 4064L, 6049L,
9399L, 13304L, 45320L, 242559L, 478694L, 717049L, 1067634L, 2137281L,
1529403L, 1171033L, 366950L, 202266L, 136921L, 541064L, 65536L,
49912L, 80312L, 118725L, 13401L, 83935L, 37848L, 0L, 147L, 1231L,
3216L, 5021L, 18462L, 37441L, 64781L, 136805L, 339245L, 561253L,
718017L, 361430L, 185699L, 310734L, 261918L, 218586L, 61845L,
38098L, 26299L, 65656L, 9851L, 21265L, 31828L, 1395L, 11348L,
765L, 1753L, 6226L, 9727L, 22033L, 50472L, 67863L, 101330L, 279209L,
241003L, 420242L, 849967L, 1079237L, 582994L, 440150L, 496096L,
198606L, 97184L, 407660L, 171474L, 127537L, 59687L, 95001L, 27311L,
18797L, 15080L, 9346L, 24760L, 41699L, 40032L, 80576L, 152208L,
171920L, 294192L, 525182L, 443886L, 866027L, 916377L, 1477769L,
1359528L, 583725L, 253410L, 90428L, 146786L, 163327L, 196846L,
140743L, 78635L, 24575L, 39309L, 81900L, 93600L), Year = c(2012L,
2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L,
2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L,
2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2011L, 2011L,
2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L,
2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L,
2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2010L, 2010L, 2010L,
2010L, 2010L, 2010L, 2010L, 2010L, 2010L, 2010L, 2010L, 2010L,
2010L, 2010L, 2010L, 2010L, 2010L, 2010L, 2010L, 2010L, 2010L,
2010L, 2010L, 2010L, 2010L, 2010L, 2009L, 2009L, 2009L, 2009L,
2009L, 2009L, 2009L, 2009L, 2009L, 2009L, 2009L, 2009L, 2009L,
2009L, 2009L, 2009L, 2009L, 2009L, 2009L, 2009L, 2009L, 2009L,
2009L, 2009L, 2009L, 2009L, 2008L, 2008L, 2008L, 2008L, 2008L,
2008L, 2008L, 2008L, 2008L, 2008L, 2008L, 2008L, 2008L, 2008L,
2008L, 2008L, 2008L, 2008L, 2008L, 2008L, 2008L, 2008L, 2008L,
2008L, 2008L, 2008L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L,
2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L,
2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L, 2007L,
2007L, 2007L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L,
2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L,
2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L, 2006L,
2006L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L,
2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L,
2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L, 2005L,
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L,
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L,
2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2004L, 2003L,
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L,
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L,
2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2003L, 2002L, 2002L,
2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L,
2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L, 2002L,
2002L, 2002L, 2002L, 2002L, 2002L, 2002L)), .Names = c("IDWeek",
"Total", "Fry", "Smolt", "FryEq", "Year"), class = "data.frame", row.names =
c(NA,
-286L))
library(ggplot2)
library(plyr)
dim(winter);head(winter)
str(winter)
test <- ggplot(winter,aes(IDWeek,Total)) + geom_line(colour='blue') +
facet_wrap(~Year,as.table=F) +
geom_vline(xintercept=47,colour='red',linetype='dashed') +
geom_hline(yintercept=300000,colour='orange',linetype='dashed',size=1) +
opts(title="Winter Chinook 2002-2012") +
geom_point(data=subset(winter,Year==2012 &
Total==max(Total)),size=7,shape=1,color='green') # Want to be able to subset
using max on the fly
test
# But it doesn't subset at all
# It works if I get the max Total like this:
test <- ggplot(winter,aes(IDWeek,Total)) + geom_line(colour='blue') +
facet_wrap(~Year,as.table=F) +
geom_vline(xintercept=47,colour='red',linetype='dashed') +
geom_hline(yintercept=300000,colour='orange',linetype='dashed',size=1) +
opts(title="Winter Chinook 2002-2012") +
geom_point(data=subset(winter,Year==2012 & Total >
300000),size=7,shape=1,color='green') # it works if get the Total like this
test
# Zooming in to the weeks of interest (weeks 45 to 52) and put a point on
the 2012 max catch which is week 37
windows()
test + coord_cartesian(xlim=c(45,52)) +
scale_x_continuous(breaks=c(45,47,49,51)) +
scale_y_continuous(limits=c(0,350000)) + opts(title="Zooming in to weeks 45
to 52")
Felipe D. Carrillo
Supervisory Fishery Biologist
Department of the Interior
US Fish & Wildlife Service
California, USA
http://www.fws.gov/redbluff/rbdd_jsmp.aspx
>________________________________
>From: arun <smartpink...@yahoo.com>
>To: Felipe Carrillo <mazatlanmex...@yahoo.com>
>Cc: William Dunlap <wdun...@tibco.com>; David Winsemius
><dwinsem...@comcast.net>; R help <r-help@r-project.org>
>Sent: Sunday, December 2, 2012 1:18 PM
>Subject: Re: [R] simple subset question
>
>Hi,
>From the ddply() output, you could get the whole row by:
>
> fish1 <- structure(list(Year = 2002:2012, maxTotal = c(1464311L, 1071051L,
>714837L, 2115018L, 850491L, 207537L, 321195L, 935599L, 194429L,
>157260L, 303259L)), .Names = c("Year", "maxTotal"), row.names = c(NA,
>-11L), class = "data.frame")
>
>
> fish[fish[,2]%in%fish1[,2][fish1[,1]==2012],] #fish (or winter) is your
>original dataset
># IDWeek Total Fry Smolt FryEq Year
>#21 47 303259 34008 269248 491733 2012
>A.K.
>
>
>
>
>
>________________________________
>From: Felipe Carrillo <mazatlanmex...@yahoo.com>
>To: William Dunlap <wdun...@tibco.com>; arun <smartpink...@yahoo.com>
>Cc: R help <r-help@r-project.org>
>Sent: Sunday, December 2, 2012 2:34 PM
>Subject: Re: [R] simple subset question
>
>
>
>Using my whole dataset I get:
>library(plyr)
>ddply(winter,"Year",summarise,maxTotal=max(Total))
>
> fish <- structure(list(Year = 2002:2012, maxTotal = c(1464311L, 1071051L,
>714837L, 2115018L, 850491L, 207537L, 321195L, 935599L, 194429L,
>157260L, 303259L)), .Names = c("Year", "maxTotal"), row.names = c(NA,
>-11L), class = "data.frame")
>
>I only want to extract the max Total for 2012 and want the whole row like this:
> IDWeek Total Fry Smolt FryEq Year
>21 47 303259 34008 269248 491733 2012
>
>My whole dataset is too big to post it so thanks for your help and will try
>to figure out why subset returns an empty row
>
>Felipe D. Carrillo
>Supervisory Fishery Biologist
>Department of the Interior
>US Fish & Wildlife Service
>California, USA
>http://www.fws.gov/redbluff/rbdd_jsmp.aspx
>
>
>
>From: William Dunlap <wdun...@tibco.com>
>>To: Felipe Carrillo <mazatlanmex...@yahoo.com>; arun <smartpink...@yahoo.com>
>>Cc: R help <r-help@r-project.org>
>>Sent: Sunday, December 2, 2012 11:00 AM
>>Subject: RE: [R] simple subset question
>>
>>> I am
>>> still getting an error message
>>> >with :
>>> > x <- subset(fish,Year==2012
>& Total==max(Total));x
>>> >I get:
>>> >[1] IDWeek Total Fry Smolt FryEq Year
>>> ><0 rows> (or 0-length row.names)
>>
>>The above is not an error message. It says that there
>>are no rows satisfying your criteria. Note that Total==max(Total)
>>returns a TRUE for each row in which the Total value
>>equals the maximum Total value over all the years in
>>the data. Are you looking for the maximum value of Total
>>in each year?
>>
>>> tmp <- transform(fish, YearlyMaxTotal = ave(Total, Year, FUN=max))
>>> subset(tmp, Total==YearlyMaxTotal)
>> IDWeek Total Fry Smolt FryEq Year YearlyMaxTotal
>>21 47 303259 34008 269248 491733 2012 303259
>>39 39 157260 156909 351 157506 2011 157260
>>> subset(tmp, Total==YearlyMaxTotal
>& Year==2012)
>> IDWeek Total Fry Smolt FryEq Year YearlyMaxTotal
>>21 47 303259 34008 269248 491733 2012 303259
>>
>>Bill Dunlap
>>Spotfire, TIBCO Software
>>wdunlap tibco.com
>>
>>
>>> -----Original Message-----
>>> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
>>> Behalf
>>> Of Felipe Carrillo
>>> Sent: Sunday, December 02, 2012 10:47 AM
>>> To: arun
>>> Cc: R help
>>> Subject: Re: [R] simple subset question
>>>
>>> Works with the small dataset (2 years) but I get the error message with the
>>> whole
>>> dataset (12 years of data). I am going to have
>>> to check what's
>wrong with it...Thanks
>>>
>>> Felipe D. Carrillo
>>> Supervisory Fishery Biologist
>>> Department of the Interior
>>> US Fish & Wildlife Service
>>> California, USA
>>> http://www.fws.gov/redbluff/rbdd_jsmp.aspx
>>>
>>>
>>> From: arun <smartpink...@yahoo.com>
>>> >To: Felipe Carrillo <mazatlanmex...@yahoo.com>
>>> >Cc: R help <r-help@r-project.org>; R. Michael Weylandt
>>> <michael.weyla...@gmail.com>
>>> >Sent: Sunday, December 2, 2012 10:29 AM
>>> >Subject: Re: [R] simple subset question
>>> >
>>>
>>Hi,
>>> >I am getting this:
>>> >x<-subset(fish,Year==2012 & Total==max(Total))
>>> > x
>>> ># IDWeek Total Fry Smolt FryEq Year
>>> >#21 47 303259 34008 269248 491733 2012
>>> >A.K.
>>> >
>>> >
>>> >
>>> >
>>> >----- Original Message -----
>>> >From: Felipe Carrillo <mazatlanmex...@yahoo.com>
>>> >To: R. Michael Weylandt <michael.weyla...@gmail.com>
>>> >Cc: "r-help@r-project.org" <r-help@r-project.org>
>>> >Sent:
>Sunday, December 2, 2012 1:25 PM
>>> >Subject: Re: [R] simple subset question
>>> >
>>> >Sorry, I was trying it to subset from a bigger dataset called 'winter' and
>>> >forgot to
>>> change the variable names
>>> >when I asked the question. David W suggestion works but the strange part
>>> >is that I am
>>> still getting an error message
>>> >with :
>>> > x <- subset(fish,Year==2012 & Total==max(Total));x
>>> >I get:
>>> >[1] IDWeek Total Fry Smolt FryEq Year
>>> ><0 rows> (or 0-length row.names)
>>> >
>>> >I will start a fresh session to see if that helps...Thank you all
>>> >
>>> >Felipe D. Carrillo
>>> >Supervisory Fishery Biologist
>>> >Department of the Interior
>>> >US Fish & Wildlife Service
>>> >California, USA
>>> >http://www.fws.gov/redbluff/rbdd_jsmp.aspx
>>> >
>>> >
>>> >From: R. Michael Weylandt <michael.weyla...@gmail.com>
>>> >>To: Felipe Carrillo <mazatlanmex...@yahoo.com>
>>> >>Cc: "r-help@r-project.org" <r-help@r-project.org>
>>> >>Sent: Sunday, December 2, 2012 9:42 AM
>>> >>Subject: Re: [R] simple subset question
>>> >>
>>> >>On Sun, Dec 2, 2012 at 5:21 PM, Felipe Carrillo
>>> >><mazatlanmex...@yahoo.com> wrote:
>>> >>> Hi,
>>> >>> Consider the small dataset below, I want to subset by two variables in
>>> >>> one line but it wont work...it works though if I subset separately. I
>>> >>> have
>>> >>> to be missing something obvious that I did not realize before while
>>> >>> using subset..
>>> >>>
>>> >>> fish <- structure(list(IDWeek = c(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, 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), Total = c(0L, 0L, 326L, 1735L,
>>> >>> 1807L, 2208L, 3883L, 8820L, 6060L, 19326L, 63158L, 100718L, 53015L,
>>>
>>>> 91689L, 152629L, 122708L, 61293L, 15574L, 86538L, 75365L, 303259L,
>>> >>> 19691L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 161L, 321L, 1000L, 4425L,
>>> >>> 13202L, 19726L, 30518L, 84949L, 157260L, 145691L, 85801L, 62044L,
>>> >>> 44439L, 23272L, 22391L, 20159L, 14854L, 35379L, 31142L, 7736L,
>>> >>> 13221L, 4894L), Fry = c(0L, 0L, 326L, 1735L, 1807L, 2208L, 3883L,
>>> >>> 8759L, 6060L, 19326L, 63119L, 100524L, 52582L, 88170L, 145564L,
>>> >>> 111416L, 38233L, 5248L, 17826L, 11038L, 34008L, 215L, 0L, 0L,
>>> >>> 0L, 0L, 0L, 0L, 0L, 0L, 161L, 321L, 1000L, 4425L, 13055L, 19488L,
>>> >>> 30518L, 84818L, 156909L, 144786L, 84207L, 57720L, 31049L, 6858L,
>>> >>> 1616L, 719L, 364L, 49L, 0L, 0L, 0L, 0L), Smolt = c(0L, 0L, 0L,
>>> >>> 0L, 0L, 0L, 0L, 62L, 0L, 0L, 38L, 195L, 433L, 3518L, 7067L, 11290L,
>>> >>> 23058L,
>10327L, 68712L, 64328L, 269248L, 19479L, 0L, 0L, 0L,
>>> >>> 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 147L, 238L, 0L, 131L, 351L,
>>> >>> 905L, 1592L, 4324L, 13391L, 16414L, 20774L, 19444L, 14491L, 35330L,
>>> >>> 31142L, 7736L, 13221L, 4894L), FryEq = c(0L, 0L, 326L, 1735L,
>>> >>> 1807L, 2208L, 3883L, 8864L, 6060L, 19326L, 63185L, 100854L, 53318L,
>>> >>> 94151L, 157576L, 130610L, 77432L, 22805L, 134639L, 120393L, 491733L,
>>> >>> 33327L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 161L, 321L, 1000L, 4425L,
>>> >>> 13306L, 19894L, 30518L, 85042L, 157506L, 146328L, 86914L, 65073L,
>>> >>> 53812L, 34763L, 36931L, 33769L, 24998L, 60110L, 52938L, 13149L,
>>> >>> 22476L, 8319L), Year = c(2012L, 2012L, 2012L, 2012L, 2012L, 2012L,
>>> >>> 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L, 2012L,
>>> >>> 2012L, 2012L, 2012L,
>2012L, 2012L, 2012L, 2012L, 2012L, 2012L,
>>> >>> 2012L, 2012L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L,
>>> >>> 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L,
>>> >>> 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L, 2011L,
>>> >>> 2011L)), .Names = c("IDWeek", "Total", "Fry", "Smolt", "FryEq",
>>> >>> "Year"), row.names = c(NA, 52L), class = "data.frame")
>>> >>> fish
>>> >>> # Subset to get the max Total for 2012
>>> >>> x <- subset(winter,Year==2012 & Total==max(Total));b # How come one
>>> >>>line doesn't
>>> work?
>>> >>
>>> >>Works fine for me if I change "winter" to fish here.
>>> >>
>>> >>subset(fish,Year==2012 & Total==max(Total))
>>> >> IDWeek Total Fry Smolt FryEq Year
>>> >>21 47
>303259 34008 269248 491733 2012
>>> >>
>>> >>>
>>> >>> # It works if I subset the year first and then get the Total max from
>>> >>>it
>>> >>> xx <- subset(winter,Year==2012)
>>> >>> xxx <- subset(xx,Total==max(Total));xxx
>>> >>> xxx
>>> >>>
>>> >>> Felipe D. Carrillo
>>> >>> Supervisory Fishery Biologist
>>> >>> Department of the Interior
>>> >>> US Fish & Wildlife Service
>>> >>> California, USA
>>> >>> http://www.fws.gov/redbluff/rbdd_jsmp.aspx
>>> >>>
>>> >>> [[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.
>>> >>>
>>> >>
>>> >>
>>> >>
>>> > [[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.
>>> >
>>> >
>>> >
>>> >
>>> [[alternative HTML version deleted]]
>>
>>
>>
>>
>
>
>
[[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.
