I finally see what you and David W mean...I was under the assumption that by subsetting Year==2012 it would also extract the max for that year.....I guess I have to review 'subset' again :) Thank you much
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: David L Carlson <dcarl...@tamu.edu> >To: 'Felipe Carrillo' <mazatlanmex...@yahoo.com>; 'arun' ><smartpink...@yahoo.com> >Cc: 'R help' <r-help@r-project.org> >Sent: Sunday, December 2, 2012 2:54 PM >Subject: RE: [R] simple subset question > >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]] > > > > > [[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.
