As David W. guessed. The maximum is in year 2005 not 2012 so no row from 2012 matches the maximum.
> subset(winter,Year==2012 & Total==max(Total)) [1] IDWeek Total Fry Smolt FryEq Year <0 rows> (or 0-length row.names) > winter[which(winter$Total==max(winter$Total)),] IDWeek Total Fry Smolt FryEq Year 195 39 2115018 2083214 31805 2137281 2005 Change to > subset(winter, Year==2012 & Total==max(Total[Year==2012])) IDWeek Total Fry Smolt FryEq Year 21 47 303259 34008 269248 491733 2012 ---------------------------------------------- David L Carlson Associate Professor of Anthropology Texas A&M University College Station, TX 77843-4352 > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of Felipe Carrillo > Sent: Sunday, December 02, 2012 4:40 PM > To: arun > Cc: R help > Subject: Re: [R] simple subset question > > 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-bounces@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.