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-bounces at r-project.org [mailto:r-help-bounces at 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 <smartpink111 at yahoo.com>
To: Felipe Carrillo <mazatlanmexico at yahoo.com> Cc: R help <r-help at r-project.org>; R. Michael Weylandt
<michael.weylandt at 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 <mazatlanmexico at yahoo.com> To: R. Michael Weylandt <michael.weylandt at gmail.com> Cc: "r-help at r-project.org" <r-help at 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.weylandt at gmail.com>
To: Felipe Carrillo <mazatlanmexico at yahoo.com> Cc: "r-help at r-project.org" <r-help at 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 <mazatlanmexico at 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 at 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 at 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]]